Population genetics and independently replicated evolution of predator-associated burst speed ecophenotypy in mosquitofish.

Many species show replicated ecophenotypy due to recurring patterns of natural selection. Based on the presence or absence of pursuit predators, at least 17 species of fish repeatedly differentiated in body shape in a manner that increases burst swimming speed and the likelihood of predator escape. The predator-associated burst speed (PABS) ecophenotype is characterized by a small head and trunk and enlarged caudal region. Mechanisms promoting replicated phenotype-environment association include selection (without evolution), a single instance of adaptive evolution followed by biased habitat occupation, repeated instances of local adaptation, or adaptive phenotypic plasticity. Common garden rearing of mosquitofish, Gambusia affinis, demonstrated a likely heritable basis for PABS phenotypy, but it is unknown whether populations are otherwise genetically distinct or whether replicated ecophenotypy represents a single or replicated instances of adaptation. To genetically characterize the populations and test hypotheses of single or multiple adaptations, we characterized variation in 12 polymorphic DNA microsatellites in the previously studied G. affinis populations. Populations were genetically distinct by multilocus analysis, exhibited high allelic diversity, and were heterozygote deficient, which effects were attributed to G. affinis's shoaling nature and habitat patchiness. Genetic and phenotypic distances among populations were correlated for non-PABS but not PABS morphology. Multilocus analysis demonstrated ecophenotype polyphyly and scattered multivariate genetic structure which support only the replicated-adaptation model. As all of the diverse tests performed demonstrated lack of congruence between patterns of molecular genetic and PABS differentiation, it is likely that divergent natural selection drove multiple instances of adaptive evolution.

Adapterama III: Quadruple-indexed, double/triple-enzyme RADseq libraries (2RAD/3RAD).

Molecular ecologists frequently use genome reduction strategies that rely upon restriction enzyme digestion of genomic DNA to sample consistent portions of the genome from many individuals (e.g., RADseq, GBS). However, researchers often find the existing methods expensive to initiate and/or difficult to implement consistently, especially because it is difficult to multiplex sufficient numbers of samples to fill entire sequencing lanes. Here, we introduce a low-cost and highly robust approach for the construction of dual-digest RADseq libraries that build on adapters and primers designed in Adapterama I. Major features of our method include: (1) minimizing the number of processing steps; (2) focusing on a single strand of sample DNA for library construction, allowing the use of a non-phosphorylated adapter on one end; (3) ligating adapters in the presence of active restriction enzymes, thereby reducing chimeras; (4) including an optional third restriction enzyme to cut apart adapter-dimers formed by the phosphorylated adapter, thus increasing the efficiency of adapter ligation to sample DNA, which is particularly effective when only low quantity/quality DNA samples are available; (5) interchangeable adapter designs; (6) incorporating variable-length internal indexes within the adapters to increase the scope of sample indexing, facilitate pooling, and increase sequence diversity; (7) maintaining compatibility with universal dual-indexed primers and thus, Illumina sequencing reagents and libraries; and, (8) easy modification for the identification of PCR duplicates. We present eight adapter designs that work with 72 restriction enzyme combinations. We demonstrate the efficiency of our approach by comparing it with existing methods, and we validate its utility through the discovery of many variable loci in a variety of non-model organisms. Our 2RAD/3RAD method is easy to perform, has low startup costs, has increased utility with low-concentration input DNA, and produces libraries that can be highly-multiplexed and pooled with other Illumina libraries.

Adapterama II: universal amplicon sequencing on Illumina platforms (TaggiMatrix).

Next-generation sequencing (NGS) of amplicons is used in a wide variety of contexts. In many cases, NGS amplicon sequencing remains overly expensive and inflexible, with library preparation strategies relying upon the fusion of locus-specific primers to full-length adapter sequences with a single identifying sequence or ligating adapters onto PCR products. In Adapterama I, we presented universal stubs and primers to produce thousands of unique index combinations and a modifiable system for incorporating them into Illumina libraries. Here, we describe multiple ways to use the Adapterama system and other approaches for amplicon sequencing on Illumina instruments. In the variant we use most frequently for large-scale projects, we fuse partial adapter sequences (TruSeq or Nextera) onto the 5' end of locus-specific PCR primers with variable-length tag sequences between the adapter and locus-specific sequences. These fusion primers can be used combinatorially to amplify samples within a 96-well plate (8 forward primers + 12 reverse primers yield 8 × 12 = 96 combinations), and the resulting amplicons can be pooled. The initial PCR products then serve as template for a second round of PCR with dual-indexed iTru or iNext primers (also used combinatorially) to make full-length libraries. The resulting quadruple-indexed amplicons have diversity at most base positions and can be pooled with any standard Illumina library for sequencing. The number of sequencing reads from the amplicon pools can be adjusted, facilitating deep sequencing when required or reducing sequencing costs per sample to an economically trivial amount when deep coverage is not needed. We demonstrate the utility and versatility of our approaches with results from six projects using different implementations of our protocols. Thus, we show that these methods facilitate amplicon library construction for Illumina instruments at reduced cost with increased flexibility. A simple web page to design fusion primers compatible with iTru primers is available at: http://baddna.uga.edu/tools-taggi.html. A fast and easy to use program to demultiplex amplicon pools with internal indexes is available at: https://github.com/lefeverde/Mr_Demuxy.

A High-Quality Reference Genome for the Invasive Mosquitofish Gambusia affinis Using a Chicago Library.

The western mosquitofish, Gambusia affinis, is a freshwater poecilid fish native to the southeastern United States but with a global distribution due to widespread human introduction. Gambusia affinis has been used as a model species for a broad range of evolutionary and ecological studies. We sequenced the genome of a male G. affinis to facilitate genetic studies in diverse fields including invasion biology and comparative genetics. We generated Illumina short read data from paired-end libraries and in vitro proximity-ligation libraries. We obtained 54.9× coverage, N50 contig length of 17.6 kb, and N50 scaffold length of 6.65 Mb. Compared to two other species in the Poeciliidae family, G. affinis has slightly fewer genes that have shorter total, exon, and intron length on average. Using a set of universal single-copy orthologs in fish genomes, we found 95.5% of these genes were complete in the G. affinis assembly. The number of transposable elements in the G. affinis assembly is similar to those of closely related species. The high-quality genome sequence and annotations we report will be valuable resources for scientists to map the genetic architecture of traits of interest in this species.

Insertions/Deletions-Associated Nucleotide Polymorphism in Arabidopsis thaliana.

Although high levels of within-species variation are commonly observed, a general mechanism for the origin of such variation is still lacking. Insertions and deletions (indels) are a widespread feature of genomes and we hypothesize that there might be an association between indels and patterns of nucleotide polymorphism. Here, we investigate flanking sequences around 18 indels (>100 bp) among a large number of accessions of the plant, Arabidopsis thaliana. We found two distinct haplotypes, i.e., a nucleotide dimorphism, present around each of these indels and dimorphic haplotypes always corresponded to the indel-present/-absent patterns. In addition, the peaks of nucleotide diversity between the two divergent alleles were closely associated with these indels. Thus, there exists a close association between indels and dimorphisms. Further analysis suggests that indel-associated substitutions could be an important component of genetic variation shaping nucleotide polymorphism in Arabidopsis. Finally, we suggest a mechanism by which indels might generate these highly divergent haplotypes. This study provides evidence that nucleotide dimorphisms, which are frequently regarded as evidence of frequency-dependent selection, could be explained simply by structural variation in the genome.

High degree of clonal reproduction and lack of large-scale geographic patterning mark the introduced range of the invasive vine, kudzu (Pueraria montana var. lobata), in North America.

PREMISE OF THE STUDY: Pueraria montana var. lobata, or kudzu, is an invasive species whose invasion in North America is not genetically well characterized. The clonality of kudzu introduces challenges to population genetic analyses that can bias the assessment of spatial patterns of genotypes. Assessing patterns of genetic diversity while considering clonality is necessary to understand the invasion and spread of kudzu in its invasive range. METHODS: We screened 1747 individuals from 87 populations across the invasive range with 15 microsatellite markers and a 789 bp chloroplast region. We performed detailed clonal analyses and tested levels of genetic diversity, population structure, and phylogeographic relationships. KEY RESULTS: Kudzu exhibited a clonal rate of 80%, and was more heterozygous than other long-lived perennials. We detected only 353 distinct clonal lineages, with over 60% sharing a maternal haplotype. Populations were established with few genotypes, many consisting of only a single clone. We found no isolation by distance. Despite high genetic diversity, we found little geographic patterning. CONCLUSIONS: Kudzu is highly clonal with few genetically distinct lineages and haplotypes existing in the introduced range. Our data are consistent with a large single introduction, or a few at most. Introduced lineages are geographically randomly distributed but isolated, suggesting that genotypes rarely expand into already established populations. No route of expansion was detectable from an original introduction. The invasion of kudzu does not seem to have been dominated by a single genotype, thus standing genetic variation and phenotypic plasticity are more likely mechanisms explaining kudzu's invasion success.

RADcap: sequence capture of dual-digest RADseq libraries with identifiable duplicates and reduced missing data.

Molecular ecologists seek to genotype hundreds to thousands of loci from hundreds to thousands of individuals at minimal cost per sample. Current methods, such as restriction-site-associated DNA sequencing (RADseq) and sequence capture, are constrained by costs associated with inefficient use of sequencing data and sample preparation. Here, we introduce RADcap, an approach that combines the major benefits of RADseq (low cost with specific start positions) with those of sequence capture (repeatable sequencing of specific loci) to significantly increase efficiency and reduce costs relative to current approaches. RADcap uses a new version of dual-digest RADseq (3RAD) to identify candidate SNP loci for capture bait design and subsequently uses custom sequence capture baits to consistently enrich candidate SNP loci across many individuals. We combined this approach with a new library preparation method for identifying and removing PCR duplicates from 3RAD libraries, which allows researchers to process RADseq data using traditional pipelines, and we tested the RADcap method by genotyping sets of 96-384 Wisteria plants. Our results demonstrate that our RADcap method: (i) methodologically reduces (to <5%) and allows computational removal of PCR duplicate reads from data, (ii) achieves 80-90% reads on target in 11 of 12 enrichments, (iii) returns consistent coverage (≥4×) across >90% of individuals at up to 99.8% of the targeted loci, (iv) produces consistently high occupancy matrices of genotypes across hundreds of individuals and (v) costs significantly less than current approaches.

Eleven microsatellites in an emerging invader, Phytolacca americana (Phytolaccaceae), from its native and introduced ranges.

PREMISE OF THE STUDY: To facilitate population genetic analyses, microsatellite markers were developed for pokeweed (Phytolacca americana), a large, weedy, perennial herb native to eastern North America that is emerging as a significant invasive species in China. METHODS AND RESULTS: We mined 1,100,538 Illumina MiSeq reads from genomic DNA for microsatellites and identified 58 primer pairs. We screened these primers for polymorphism in two native and two invasive populations. We identified 11 loci that amplified consistently. The number of alleles per locus ranged from two to six, and observed heterozygosity ranged from 0.00 to 1.00. All loci were largely monomorphic within populations but different among populations. The primers were of very limited use in the congener P. acinosa. CONCLUSIONS: These loci will provide a valuable resource to study the population genetics and invasion history of P. americana.

Genetic factors associated with mating system cause a partial reproductive barrier between two parapatric species of leavenworthia (brassicaceae).

Reproductive barriers play a major role in the origin and maintenance of biodiversity by restricting gene flow between species. Although both pre- and postzygotic barriers often isolate species, prezygotic barriers are thought to contribute more to reproductive isolation. We investigated possible reproductive barriers between Leavenworthia alabamica and L. crassa, parapatric species with high morphological and ecological similarity and the ability to hybridize. Using greenhouse and field experiments, we tested for habitat isolation and genetic incompatibilities. From controlled crosses, we identified unilateral incompatibility (a partial prezygotic barrier associated with the self-incompatibility system), but no evidence of other genetic incompatibilities. We found a small reduction in pollen viability of F(1) hybrids and early germination of F(1), F(2), and BC hybrids relative to L. alabamica and L. crassa in a common garden experiment, but the effect on fitness was not tested. Field studies of hybrid pollen viability and germination are needed to determine if they contribute to reproductive isolation. In a reciprocal transplant, we found no evidence of habitat isolation or reduced hybrid survival (from seedling to adult stage) or reproduction. These data suggest unilateral incompatibility partially reproductively isolates L. alabamica and L. crassa, but no other reproductive barriers could be detected.

Constraints on the evolution of tolerance to herbicide in the common morning glory: resistance and tolerance are mutually exclusive.

Evolutionary biologists explain the maintenance of intermediate levels of defense in plant populations as being due to trade-offs, or negative genetic covariances among ecologically important traits. Attempts at detecting trade-offs as constraints on the evolution of defense have not always been successful, leading some to conclude that such trade-offs rarely explain current levels of defense in the population. Using the agricultural pest Ipomoea purpurea, we measured correlations between traits involved in defense to glyphosate, the active ingredient in Roundup, a widely used herbicide. We found significant allocation costs of tolerance, as well as trade-offs between resistance and two measures of tolerance to glyphosate. Selection on resistance and tolerance exhibited differing patterns: tolerance to leaf damage was under negative directional selection, whereas resistance was under positive directional selection. The joint pattern of selection on resistance and tolerance to leaf damage indicated the presence of alternate peaks in the fitness landscape such that a combination of either high tolerance and low resistance, or high resistance and low tolerance was favored. The widespread use of this herbicide suggests that it is likely an important selective agent on weed populations. Understanding the evolutionary dynamics of herbicide defense traits is thus of increasing importance in the context of human-mediated evolution.

QTL-based evidence for the role of epistasis in evolution.

The extent to which epistasis contributes to adaptation and speciation has been a controversial topic in evolutionary genetics. One experimental approach to study epistasis is based on quantitative trait locus (QTL) mapping using molecular markers. Comparisons can be made among all possible pair-wise combinations of the markers, irrespective of whether an additive QTL is associated with a marker; several software packages have been developed that facilitate this. We review several examples of using this approach to identify epistatic QTLs for traits of evolutionary or ecological interest. While there is variability in the results, the number of epistatic QTL interactions is often greater than or equal to the number of additive QTLs. The magnitude of epistatic effects can be larger than the additive effects. Thus, epistatic interactions seem to be an important part of natural genetic variation. Future studies of epistatic QTLs could lead to descriptions of the genetic networks underlying variation for fitness-related traits.

Epistasis for fitness-related quantitative traits in Arabidopsis thaliana grown in the field and in the greenhouse.

The extent to which epistasis contributes to adaptation, population differentiation, and speciation is a long-standing and important problem in evolutionary genetics. Using recombinant inbred (RI) lines of Arabidopsis thaliana grown under natural field conditions, we have examined the genetic architecture of fitness-correlated traits with respect to epistasis; we identified both single-locus additive and two-locus epistatic QTL for natural variation in fruit number, germination, and seed length and width. For fruit number, we found seven significant epistatic interactions, but only two additive QTL. For seed germination, length, and width, there were from two to four additive QTL and from five to eight epistatic interactions. The epistatic interactions were both positive and negative. In each case, the magnitude of the epistatic effects was roughly double that of the effects of the additive QTL, varying from -41% to +29% for fruit number and from -5% to +4% for seed germination, length, and width. A number of the QTL that we describe participate in more than one epistatic interaction, and some loci identified as additive also may participate in an epistatic interaction; the genetic architecture for fitness traits may be a network of additive and epistatic effects. We compared the map positions of the additive and epistatic QTL for germination, seed width, and seed length from plants grown in both the field and the greenhouse. While the total number of significant additive and epistatic QTL was similar under the two growth conditions, the map locations were largely different. We found a small number of significant epistatic QTL x environment effects when we tested directly for them. Our results support the idea that epistatic interactions are an important part of natural genetic variation and reinforce the need for caution in comparing results from greenhouse-grown and field-grown plants.

Ontogenetics of QTL: the genetic architecture of trichome density over time in Arabidopsis thaliana.

Although much is known about the molecular genetic basis of trichome development in Arabidopsis thaliana, less is known about the underlying genetic basis of continuous variation in a trait known to be of adaptive importance: trichome density. The density of leaf trichomes is known to be a major determinant of herbivore damage in natural populations of A. thaliana and herbivores are a significant selective force on genetic variation for trichome density. A number of developmental changes occur during ontogeny in A. thaliana, including changes in trichome density. I used multiple interval mapping (MIM) analysis to identify QTL responsible for trichome density on both juvenile leaves and adult leaves in replicate, independent trials and asked whether those QTL changed with ontogeny. In both juvenile and adult leaves, I detected a single major QTL on chromosome 2 that explained much of the genetic variance. Although additional QTL were detected, there were no consistent differences in the genetic architecture of trichome density measured on juvenile and adult leaves. The finding of a single QTL of major effect for a trait of known adaptive importance suggests that genes of major effect may play an important role in adaptation.

Hybridization as a source of evolutionary novelty: leaf shape in a Hawaiian composite.

Hybridization is increasingly recognized as a significant creative force in evolution. Interbreeding among species can lead to the creation of novel genotypes and morphologies that lead to adaptation. On the Hawaiian island of O'ahu, populations of two species of plants in the endemic genus Lipochaeta grow at similar elevations in the northern Wai'anae Mountains. These two species represent extremes of the phenotypic distribution of leaf shape: the leaves of Lipochaeta tenuifolia individuals are compound and highly dissected while leaves of L. tenuis are simple. Based primarily on leaf shape morphology, a putative hybrid population of Lipochaeta located at Pu'u Kawiwi was identified. Individuals in this population exhibit a range of leaf shapes intermediate in varying degrees between the leaf shapes of the putative parental species. We analyzed individuals from pure populations of L. tenuifolia, L. tenuis and the putative hybrids using 133 AFLP markers. Genetic analysis of these neutral markers provided support for the hybrid origin of this population. The correlation between genetic background and leaf morphology in the hybrids suggested that the genome of the parental species with simple leaves might have significantly contributed to the evolution of a novel, compound leaf morphology.

Can ecology help genomics: the genome as ecosystem?

Ecologists study the rules that govern processes influencing the distribution and abundance of organisms, particularly with respect to the interactions of organisms with their biotic and abiotic environments. Over the past decades, using a combination of sophisticated mathematical models and rigorous experiments, ecologists have made considerable progress in understanding the complex web of interactions that constitute an ecosystem. The field of genomics runs on a path parallel to ecology. Like ecology, genomicists seek to understand how each gene in the genome interacts with every other gene and how each gene interacts with multiple, environmental factors. Gene networks connect genes as complex as the 'webs' that connect the species in an ecosystem. In fact, genes exist in an ecosystem we call the genome. The genome as ecosystem is more than a metaphor--it serves as the conceptual foundation for an interdisciplinary approach to the study of complex systems characteristic of both genomics and ecology. Through the infusion of genomics into ecology and ecology into genomics both fields will gain fresh insight into the outstanding major questions of their disciplines.

Fitness costs and benefits of novel herbicide tolerance in a noxious weed.

Glyphosate, the active ingredient in the herbicide RoundUp, has increased dramatically in use over the past decade and constitutes a potent anthropogenic source of selection. In the southeastern United States, weedy morning glories have begun to develop tolerance to glyphosate, representing a unique opportunity to examine the evolutionary genetics of a novel trait. We found genetic variation for tolerance, indicating the potential for the population to respond to selection by glyphosate. However, the following significant evolutionary constraint exists: in the absence of glyphosate, tolerant genotypes produced fewer seeds than susceptible genotypes. The combination of strong positive directional selection in the presence of glyphosate and strong negative directional selection in its absence may indicate that the selective landscape of land use could drive the evolutionary trajectory of glyphosate tolerance. Understanding these evolutionary forces is imperative for devising comprehensive management strategies to help slow the rate of the evolution of tolerance.

Natural selection for polymorphism in the disease resistance gene Rps2 of Arabidopsis thaliana.

Pathogen resistance is an ecologically important phenotype increasingly well understood at the molecular genetic level. In this article, we examine levels of avrRpt2-dependent resistance and Rps2 locus DNA sequence variability in a worldwide sample of 27 accessions of Arabidopsis thaliana. The rooted parsimony tree of Rps2 sequences drawn from a diverse set of ecotypes includes a deep bifurcation separating major resistance and susceptibility clades of alleles. We find evidence for selection maintaining these alleles and identify the N-terminal part of the leucine-rich repeat region as a probable target of selection. Additional protein variants are found within the two major clades and correlate well with measurable differences among ecotypes in resistance to the avirulence gene avrRpt2 of the pathogen Pseudomonas syringae. Long-lived polymorphisms have been observed for other resistance genes of A. thaliana; the Rps2 data suggest that the long-term maintenance of phenotypic variation in resistance genes may be a general phenomenon and are consistent with diversifying selection acting in concert with selection to maintain variation.

Testing for environmentally induced bias in phenotypic estimates of natural selection: theory and practice.

Measuring natural selection has been a fundamental goal of evolutionary biology for more than a century, and techniques developed in the last 20 yr have provided relatively simple means for biologists to do so. Many of these techniques, however, share a common limitation: when applied to phenotypic data, environmentally induced covariances between traits and fitness can lead to biased estimates of selection and misleading predictions about evolutionary change. Utilizing estimates of breeding values instead of phenotypic data with these methods can eliminate environmentally induced bias, although this approach is more difficult to implement. Despite this potential limitation to phenotypic methods and the availability of a potential solution, little empirical evidence exists on the extent of environmentally induced bias in phenotypic estimates of selection. In this article, we present a method for detecting bias in phenotypic estimates of selection and demonstrate its use with three independent data sets. Nearly 25% of the phenotypic selection gradients estimated from our data are biased by environmental covariances. We find that bias caused by environmental covariances appears mainly to affect quantitative estimates of the strength of selection based on phenotypic data and that the magnitude of these biases is large. As our estimates of selection are based on data from spatially replicated field experiments, we suggest that our findings on the prevalence of bias caused by environmental covariances are likely to be conservative.