Christopher W. Wheat

Rapid transcriptome characterization for a nonmodel organism using 454 pyrosequencing

We present a de novo assembly of a eukaryote transcriptome using 454 pyrosequencing data. The Glanville fritillary butterfly (Melitaea cinxia; Lepidoptera: Nymphalidae) is a prominent species in population biology but had no previous genomic data. Sequencing runs using two normalized complementary DNA collections from a genetically diverse pool of larvae, pupae, and adults yielded 608 053 expressed sequence tags (mean length = 110 nucleotides), which assembled into 48 354 contigs (sets of overlapping DNA segments) and 59 943 singletons. blast comparisons confirmed the accuracy of the sequencing and assembly, and indicated the presence of c. 9000 unique genes, along with > 6000 additional microarray‐confirmed unannotated contigs. Average depth of coverage was 6.5‐fold for the longest 4800 contigs (348–2849 bp in length), sufficient for detecting large numbers of single nucleotide polymorphisms. Oligonucleotide microarray probes designed from the assembled sequences showed highly repeatable hybridization intensity and revealed biological differences among individuals. We conclude that 454 sequencing, when performed to provide sufficient coverage depth, allows de novo transcriptome assembly and a fast, cost‐effective, and reliable method for development of functional genomic tools for nonmodel species. This development narrows the gap between approaches based on model organisms with rich genetic resources vs. species that are most tractable for ecological and evolutionary studies.

The genetic basis of a plant-insect coevolutionary key innovation.

Ehrlich and Raven formally introduced the concept of stepwise coevolution using butterfly and angiosperm interactions in an attempt to account for the impressive biological diversity of these groups. However, many biologists currently envision butterflies evolving 50 to 30 million years (Myr) after the major angiosperm radiation and thus reject coevolutionary origins of butterfly biodiversity. The unresolved central tenet of Ehrlich and Ravens theory is that evolution of plant chemical defenses is followed closely by biochemical adaptation in insect herbivores, and that newly evolved detoxification mechanisms result in adaptive radiation of herbivore lineages. Using one of their original butterfly-host plant systems, the Pieridae, we identify a pierid glucosinolate detoxification mechanism, nitrile-specifier protein (NSP), as a key innovation. Larval NSP activity matches the distribution of glucosinolate in their host plants. Moreover, by using five different temporal estimates, NSP seems to have evolved shortly after the evolution of the host plant group (Brassicales) (≈10 Myr). An adaptive radiation of these glucosinolate-feeding Pierinae followed, resulting in significantly elevated species numbers compared with related clades. Mechanistic understanding in its proper historical context documents more ancient and dynamic plant–insect interactions than previously envisioned. Moreover, these mechanistic insights provide the tools for detailed molecular studies of coevolution from both the plant and insect perspectives.

Rapidly developing functional genomics in ecological model systems via 454 transcriptome sequencing

Next generation sequencing technology affords new opportunities in ecological genetics. This paper addresses how an ecological genetics research program focused on a phenotype of interest can quickly move from no genetic resources to having various functional genomic tools. 454 sequencing and its error rates are discussed, followed by a review of de novo transcriptome assemblies focused on the first successful de novo assembly which happens to be in an ecological model system (the Glanville fritillary butterfly). The potential future developments in 454 sequencing are also covered. Particular attention is paid to the difficulties ecological geneticists are likely to encounter through reviewing relevant studies in both model and non-model systems. Various post-sequencing issues and applications of 454 generated data are presented (e.g. database management, microarray construction, molecular marker and candidate gene development). How to use species with genomic resources to inform study of those without is also discussed. In closing, some of the drawbacks of 454 sequencing are presented along with future prospects of this technology.

Adaptation at specific loci. VII. Natural selection, dispersal and the diversity of molecular–functional variation patterns among butterfly species complexes (Colias: Lepidoptera, Pieridae)

Natural genetic variants at the phosphoglucose isomerase, PGI, gene differ in spatial patterning of their polymorphism among species complexes of Colias butterflies in North America. In both lowland and alpine complexes, molecular–functional properties of the polymorphic genotypes can be used to predict genotype‐specific adult flight performances and resulting large genotypic differences in adult fitness components. In the lowland species complex, there is striking uniformity of PGI polymorph frequencies at a number of sites across the American West; this fits with earlier findings of strong, similar differences in fitness components over this range. In an alpine complex, Colias meadii shows similar uniformity of PGI frequencies within habitat types, either montane steppe or alpine tundra, over several hundred kilometres in the absence of dispersal. At the same time, large shifts (10–20%) in frequency of the most common alleles occur between steppe and tundra populations, whether these are isolated or, as in some cases, are in contact and exchange many dispersing adults each generation. Data on male mating success of common C. meadii PGI genotypes in steppe and tundra show heterozygote advantage in both habitat types, with shifts in relative homozygote disadvantage between habitats which are consistent with observed frequency differences. Nonadaptive explanations for this situation are rejected, and alternative, thermal‐ecology‐based adaptive hypotheses are proposed for later experimental test. These findings show that strong local selection may dominate dispersal as an evolutionary agent, whether or not dispersal is present, and that selection may often be the major force promoting ‘cohesion’ of species over long distances. This case offers new opportunities for integrating studies of molecular structure and function with ecological aspects of natural selection in the wild, both within and among species.

Fitness differences associated with Pgi SNP genotypes in the Glanville fritillary butterfly (Melitaea cinxia)

Allozyme variation at the phosphoglucose isomerase (PGI) locus in the Glanville fritillary butterfly (Melitaea cinxia) is associated with variation in flight metabolic rate, dispersal rate, fecundity and local population growth rate. To map allozyme to DNA variation and to survey putative functional variation in genomic DNA, we cloned the coding sequence of Pgi and identified nonsynonymous variable sites that determine the most common allozyme alleles. We show that these single‐nucleotide polymorphisms (SNPs) exhibit significant excess of heterozygotes in field‐collected population samples as well as in laboratory crosses. This is in contrast to previous results for the same species in which other allozymes and SNPs were in Hardy–Weinberg equilibrium or exhibited an excess of homozygotes. Our results suggest that viability selection favours Pgi heterozygotes. Although this is consistent with direct overdominance at Pgi, we cannot exclude the possibility that heterozygote advantage is caused by the presence of one or more deleterious alleles at linked loci.

Functional genomics of life history variation in a butterfly metapopulation

In fragmented landscapes, small populations frequently go extinct and new ones are established with poorly understood consequences for genetic diversity and evolution of life history traits. Here, we apply functional genomic tools to an ecological model system, the well-studied metapopulation of the Glanville fritillary butterfly. We investigate how dispersal and colonization select upon existing genetic variation affecting life history traits by comparing common-garden reared 2-day adult females from new populations with those from established older populations. New-population females had higher expression of abdomen genes involved in egg provisioning and thorax genes involved in the maintenance of flight muscle proteins. Physiological studies confirmed that new-population butterflies have accelerated egg maturation, apparently regulated by higher juvenile hormone titer and angiotensin converting enzyme mRNA, as well as enhanced flight metabolism. Gene expression varied between allelic forms of two metabolic genes (Pgi and Sdhd), which themselves were associated with differences in flight metabolic rate, population age and population growth rate. These results identify likely molecular mechanisms underpinning life history variation that is maintained by extinction-colonization dynamics in metapopulations.

Evolutionary Origins of a Novel Host Plant Detoxification Gene in Butterflies

Chemical interactions between plants and their insect herbivores provide an excellent opportunity to study the evolution of species interactions on a molecular level. Here, we investigate the molecular evolutionary events that gave rise to a novel detoxifying enzyme (nitrile-specifier protein [NSP]) in the butterfly family Pieridae, previously identified as a coevolutionary key innovation. By generating and sequencing expressed sequence tags, genomic libraries, and screening databases we found NSP to be a member of an insect-specific gene family, which we characterized and named the NSP-like gene family. Members consist of variable tandem repeats, are gut expressed, and are found across Insecta evolving in a dynamic, ongoing birth-death process. In the Lepidoptera, multiple copies of single-domain major allergen genes are present and originate via tandem duplications. Multiple domain genes are found solely within the brassicaceous-feeding Pieridae butterflies, one of them being NSP and another called major allergen (MA). Analyses suggest that NSP and its paralog MA have a unique single-domain evolutionary origin, being formed by intragenic domain duplication followed by tandem whole-gene duplication. Duplicates subsequently experienced a period of relaxed constraint followed by an increase in constraint, perhaps after neofunctionalization. NSP and its ortholog MA are still experiencing high rates of change, reflecting a dynamic evolution consistent with the known role of NSP in plant-insect interactions. Our results provide direct evidence to the hypothesis that gene duplication is one of the driving forces for speciation and adaptation, showing that both within- and whole-gene tandem duplications are a powerful force underlying evolutionary adaptation.

Nucleotide Polymorphism at a Gene (Pgi) under Balancing Selection in a Butterfly Metapopulation

The Glanville fritillary butterfly (Melitaea cinxia, Nymphalidae) has a large, well-studied metapopulation in the Aland Islands in Finland. Previous studies have found that the common allozyme genotypes at the phosphoglucose isomerase (PGI) locus are associated with individual variation in performance and fitness, with phenotypic data suggesting ongoing balancing selection via heterozygote advantage. Here, we analyze nucleotide polymorphism in the coding region of the Pgi gene. Pgi is exceptionally polymorphic, in contrast to three other metabolic genes (Mdh, Idh, and Gapdh) with low levels of polymorphism. Most of the variation is due to two common haplotype clades, which are highly divergent and exhibit extensive linkage disequilibrium. These two clades correspond to the two most common allozyme alleles previously studied. Molecular tests of selection and coalescence simulations indicate that patterns of nucleotide polymorphism depart from neutrality and are consistent with long-term balancing selection. The split between the two main haplotype clades is estimated to predate the last common ancestor of a clade of five extant Melitaea species. Comparative structural analysis of Pgi polymorphism in M. cinxia and the unrelated Colias eurytheme butterfly suggests a similar but not identical target of balancing selection. Our results indicate convergent evolution between these two species at both the phenotypic and molecular levels.

Phosphoglucose isomerase (Pgi) performance and fitness effects among Arthropods and its potential role as an adaptive marker in conservation genetics

The development of conservation genomics will be greatly aided by the use of neutral as well as adaptive molecular markers. Identifying novel adaptive molecular markers that have general application across diverse taxa is challenging, especially in Arthropods where few if any examples of balanced polymorphisms exist that are shared across species. A review of literature on the Pgi gene provides strong evidence for population level fitness consequences of genetic variation in this gene, across very diverse lineages of Arthropods. While these observations demonstrate the potential of using Pgi as an adaptive molecular marker, this gene is fundamentally different from the adaptive markers MHC and SI. Rather than providing insights into individual genetic health, Pgi appears to have a role in conservation genetics by providing insights into gene by environment interactions, local adaptation and evolutionary significant units, and potentially even morphologically cryptic dispersal phenotypes. These findings argue for studying Pgi variation in more species, as it appears central to the goals of conservation genomics.

A mitochondrial-DNA-based phylogeny for some evolutionary-genetic model species of Colias butterflies (Lepidoptera, Pieridae)

We study the phylogenetic relationships among some North American Colias (sulfur) butterflies, using mitochondrial gene sequences (ribosomal RNA, cytochrome oxidase I+II) totaling about 20% of the mitochondrial genome. We find that (1) the lowland species complex shows a branching order different from earlier views; (2) several montane and northern taxa may be more distinct than in earlier views; (3) one morphologically conservative Holarctic assemblage, C. hecla, is differentiated at the molecular-genetic level into at least three taxa which occupy distinct positions in the phylogeny and are sisters to diverse other taxa. These conclusions, constituting phylogenetic hypotheses, are supported by parsimony, maximum-likelihood, and Bayesian reconstruction algorithms. They are tested formally, by interior branch tests and paired-site tests, against alternative hypotheses derived from conventional species and subspecies naming combinations. In all cases our hypotheses are supported by these tests and the conventional alternatives are rejected. The barcoding subset of cytochrome oxidase I sequence identifies only some of the taxa supported by our full data set. Comparison of genetic divergence values among Colias taxa with those among related Pierid butterflies suggests that species radiations within Colias are comparatively younger. This emerging Colias phylogeny facilitates comparisons of genetic polymorphism and other adaptive mechanisms among taxa, thereby connecting micro- and macro-evolutionary processes.