Patterns of genetic architecture for life-history traits and molecular markers in a subdivided species

Evolutionary Genomics of a Subdivided Species

The ways in which genetic variation is distributed within and among populations is a key determinant of the evolutionary features of a species. However, most comprehensive studies of these features have been restricted to studies of subdivision in settings known to have been driven by local adaptation, leaving our understanding of the natural dispersion of allelic variation less than ideal. Here, we present a geographic population-genomic analysis of 10 populations of the freshwater microcrustacean Daphnia pulex, an emerging model system in evolutionary genomics. These populations exhibit a pattern of moderate isolation-by-distance, with an average migration rate of 0.6 individuals per generation, and average effective population sizes of ∼650,000 individuals. Most populations contain numerous private alleles, and genomic scans highlight the presence of islands of excessively high population subdivision for more common alleles. A large fraction of such islands of population divergence likely reflect historical neutral changes, including rare stochastic migration and hybridization events. The data do point to local adaptive divergence, although the precise nature of the relevant variation is diffuse and cannot be associated with particular loci, despite the very large sample sizes involved in this study. In contrast, an analysis of between-species divergence highlights positive selection operating on a large set of genes with functions nearly nonoverlapping with those involved in local adaptation, in particular ribosome structure, mitochondrial bioenergetics, light reception and response, detoxification, and gene regulation. These results set the stage for using D. pulex as a model for understanding the relationship between molecular and cellular evolution in the context of natural environments.

Does geographic variation in thermal tolerance in Daphnia represent trade-offs or conditional neutrality?

Geographic variation in thermal tolerance in Daphnia seems to represent genetic load at the loci specifically responsible for heat tolerance resulting from conditional neutrality. We see no evidence of trade-offs between fitness-related traits at 25 °C vs. 10 °C or between two algal diets across Daphnia magna clones from a variety of locations representing the opposite ends of the distribution of long-term heat tolerance. Likewise, we found no evidence of within-environment trade-offs between heat tolerance and fitness-related traits in any of the environments. Neither short-term and long-term heat tolerance shows any consistent relationship with lipid fluorescence polarization and lipid peroxidation across clones or environments. Pervasive positive correlations between fitness-related traits indicate differences in genetic load rather than trade-off based local adaptation or thermal specialization. For heat tolerance such differences may be caused by either relaxation of stabilizing selection due to lower exposure to high temperature extremes, i.e., conditional neutrality, or by small effective population size followed by the recent range expansion.

The evolution of reproductive isolation in Daphnia

BACKGROUND

The process by which populations evolve to become new species involves the emergence of various reproductive isolating barriers (RIB). Despite major advancements in understanding this complex process, very little is known about the order in which RIBs evolve or their relative contribution to the total restriction of gene flow during various stages of speciation. This is mainly due to the difficulties of studying reproductive isolation during the early stages of species formation. This study examines ecological and non-ecological RIB within and between Daphnia pulex and Daphnia pulicaria, two recently diverged species that inhabit distinct habitats and exhibit an unusual level of intraspecific genetic subdivision.

RESULTS

We find that while ecological prezygotic barriers are close to completion, none of the non-ecological barriers can restrict gene flow between D. pulex and D. pulicaria completely when acting alone. Surprisingly, we also identified high levels of postzygotic reproductive isolation in 'conspecific' interpopulation crosses of D. pulex.

CONCLUSIONS

While the ecological prezygotic barriers are prevalent during the mature stages of speciation, non-ecological barriers likely dominated the early stages of speciation. This finding indicates the importance of studying the very early stages of speciation and suggests the contribution of postzygotic isolation in initiating the process of speciation.

Genetic control of male production in Daphnia pulex

Daphnia normally reproduce by cyclical parthenogenesis, with offspring sex being determined by environmental cues. However, some females have lost the ability to produce males. Our results demonstrate that this loss of male-producing ability is controlled by a dominant allele at a single locus. We identified the locus by comparing whole-genome sequences of 67 nonmale-producing (NMP) and 100 male-producing (MP) clones from 5 Daphnia pulex populations, revealing 132 NMP-linked SNPs and 59 NMP-linked indels within a single 1.1-Mb nonrecombining region on chromosome I. These markers include 7 nonsynonymous mutations, all of which are located within one unannotated protein-coding gene (gene 8960). Within this single gene, all of the marker-linked NMP haplotypes from different populations form a monophyletic clade, suggesting a single origin of the NMP phenotype, with the NMP haplotype originating by introgression from a sister species, Daphnia pulicaria Methyl farnesoate (MF) is the innate juvenile hormone in daphnids, which induces the production of males and whose inhibition results in female-only production. Gene 8960 is sensitive to treatment by MF in MP clones, but such responsiveness is greatly reduced in NMP clones. Thus, we hypothesize that gene 8960 is located downstream of the MF-signaling pathway in D. pulex, with the NMP phenotype being caused by expression change of gene 8960.

Molecular and quantitative genetic differentiation in Sitobion avenae populations from both sides of the Qinling Mountains

Quantitative trait differences are often assumed to be correlated with molecular variation, but the relationship is not certain, and empirical evidence is still scarce. To address this issue, we sampled six populations of the cereal aphid Sitobion avenae from areas north and south of the Qinling Mountains, and characterized their molecular variation at seven microsatellite loci and quantitative variation at nine life-history traits. Our results demonstrated that southern populations had slightly longer developmental times of nymphs but much higher lifetime fecundity, compared to northern populations. Of the nine tested quantitative characters, eight differed significantly among populations within regions, as well as between northern and southern regions. Genetic differentiation in neutral markers was likely to have been caused by founder events and drift. Increased subdivision for quantitative characters was found in northern populations, but reduced in southern populations. This phenomenon was not found for molecular characters, suggesting the decoupling between molecular and quantitative variation. The pattern of relationships between FST and QST indicated divergent selection and suggested that local adaptation play a role in the differentiation of life-history traits in tested S. avenae populations, particularly in those traits closely related to reproduction. The main role of natural selection over genetic drift was also supported by strong structural differences in G-matrices among S. avenae populations. However, cluster analyses did not result in two groups corresponding to northern and southern regions. Genetic differentiation between northern and southern populations in neutral markers was low, indicating considerable gene flow between them. The relationship between molecular and quantitative variation, as well as its implications for differentiation and evolution of S. avenae populations, was discussed.

An SNP-based second-generation genetic map of Daphnia magna and its application to QTL analysis of phenotypic traits

Background

Although Daphnia is increasingly recognized as a model for ecological genomics and biomedical research, there is, as of yet, no high-resolution genetic map for the genus. Such a map would provide an important tool for mapping phenotypes and assembling the genome. Here we estimate the genome size of Daphnia magna and describe the construction of an SNP array based linkage map. We then test the suitability of the map for life history and behavioural trait mapping. The two parent genotypes used to produce the map derived from D. magna populations with and without fish predation, respectively and are therefore expected to show divergent behaviour and life-histories.

Results

Using flow cytometry we estimated the genome size of D. magna to be about 238 mb. We developed an SNP array tailored to type SNPs in a D. magna F2 panel and used it to construct a D. magna linkage map, which included 1,324 informative markers. The map produced ten linkage groups ranging from 108.9 to 203.6 cM, with an average distance between markers of 1.13 cM and a total map length of 1,483.6 cM (Kosambi corrected). The physical length per cM is estimated to be 160 kb. Mapping infertility genes, life history traits and behavioural traits on this map revealed several significant QTL peaks and showed a complex pattern of underlying genetics, with different traits showing strongly different genetic architectures.

Conclusions

The new linkage map of D. magna constructed here allowed us to characterize genetic differences among parent genotypes from populations with ecological differences. The QTL effect plots are partially consistent with our expectation of local adaptation under contrasting predation regimes. Furthermore, the new genetic map will be an important tool for the Daphnia research community and will contribute to the physical map of the D. magna genome project and the further mapping of phenotypic traits. The clones used to produce the linkage map are maintained in a stock collection and can be used for mapping QTLs of traits that show variance among the F2 clones.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1033) contains supplementary material, which is available to authorized users.

Genomic background and generation time influence deleterious mutation rates in Daphnia

Understanding how genetic variation is generated and how selection shapes mutation rates over evolutionary time requires knowledge of the factors influencing mutation and its effects on quantitative traits. We explore the impact of two factors, genomic background and generation time, on deleterious mutation in Daphnia pulicaria, a cyclically parthenogenic aquatic microcrustacean, using parallel mutation-accumulation experiments. The deleterious mutational properties of life-history characters for individuals from two different populations, and for individuals maintained at two different generation times, were quantified and compared. Mutational properties varied between populations, especially for clutch size, suggesting that genomic background influences mutational properties for some characters. Generation time was found to have a greater effect on mutational properties, with higher per-generation deleterious mutation rates in lines with longer generation times. These results suggest that differences in genetic architecture among populations and species may be explained in part by demographic features that significantly influence generation time and therefore the rate of mutation.

The effect of spontaneous mutations on competitive ability

Understanding the impact of spontaneous mutations on fitness has many theoretical and practical applications in biology. Although mutational effects on individual morphological or life-history characters have been measured in several classic genetic model systems, there are few estimates of the rate of decline due to mutation for complex fitness traits. Here, we estimate the effects of mutation on competitive ability, an important complex fitness trait, in a model system for ecological and evolutionary genomics, Daphnia. Competition assays were performed to compare fitness between mutation-accumulation (MA) lines and control lines from eight different genotypes from two populations of Daphnia pulicaria after 30 and 65 generations of mutation accumulation. Our results show a fitness decline among MA lines relative to controls as expected, but highlight the influence of genomic background on this effect. In addition, in some assays, MA lines outperform controls providing insight into the frequency of beneficial mutations.

Transcontinental phylogeography of the Daphnia pulex species complex

Daphnia pulex is quickly becoming an attractive model species in the field of ecological genomics due to the recent release of its complete genome sequence, a wide variety of new genomic resources, and a rich history of ecological data. Sequences of the mitochondrial NADH dehydrogenase subunit 5 and cytochrome c oxidase subunit 1 genes were used to assess the global phylogeography of this species, and to further elucidate its phylogenetic relationship to other members of the Daphnia pulex species complex. Using both newly acquired and previously published data, we analyzed 398 individuals from collections spanning five continents. Eleven strongly supported lineages were found within the D. pulex complex, and one lineage in particular, panarctic D. pulex, has very little phylogeographical structure and a near worldwide distribution. Mismatch distribution, haplotype network, and population genetic analyses are compatible with a North American origin for this lineage and subsequent spatial expansion in the Late Pleistocene. In addition, our analyses suggest that dispersal between North and South America of this and other species in the D. pulex complex has occurred multiple times, and is predominantly from north to south. Our results provide additional support for the evolutionary relationships of the eleven main mitochondrial lineages of the D. pulex complex. We found that the well-studied panarctic D. pulex is present on every continent except Australia and Antarctica. Despite being geographically very widespread, there is a lack of strong regionalism in the mitochondrial genomes of panarctic D. pulex--a pattern that differs from that of most studied cladocerans. Moreover, our analyses suggest recent expansion of the panarctic D. pulex lineage, with some continents sharing haplotypes. The hypothesis that hybrid asexuality has contributed to the recent and unusual geographic success of the panarctic D. pulex lineage warrants further study.

The effects of reproductive specialization on energy costs and fitness genetic variances in cyclical and obligate parthenogenetic aphids

Organisms with coexisting sexual and asexual populations are ideal models for studying the consequences of either reproductive mode on the quantitative genetic architecture of life-history traits. In the aphid Rhopalosiphum padi, lineages differing in their sex investment coexist but all share a common parthenogenetic phase. Here, we studied multiple genotypes of R. padi specialized either for sexual and asexual reproduction and compared their genetic variation in fitness during the parthenogenetic phase. Specifically, we estimated maintenance costs as standard metabolic rate (SMR), together with fitness (measured as the intrinsic rate of increase and the net reproductive rate). We found that genetic variation (in terms of broad-sense heritability) in fitness was higher in asexual genotypes compared with sexual genotypes. Also, we found that asexual genotypes exhibited several positive genetic correlations indicating that body mass, whole-animal SMR, and apterous individuals production are contributing to fitness. Hence, it appears that in asexual genotypes, energy is fully allocated to maximize the production of parthenogenetic individuals, the simplest possible form of aphid repertoire of life-histories strategies.

The effect of variable frequency of sexual reproduction on the genetic structure of natural populations of a cyclical parthenogen

Cyclical parthenogens are a valuable system in which to empirically test theoretical predictions as to the genetic consequences of sexual reproduction in natural populations, particularly if the frequency of sexual relative to asexual reproduction can be quantified. In this study, we used a series of lake populations of the cyclical parthenogen, Daphnia pulicaria, that vary consistently in their investment in sexual reproduction, to address the questions of whether the ecological variation in investment in sex is detectable at the genetic level, and if so, whether the genetic patterns seen are consistent with theoretical predictions. We show that there is variation in the genetic structure of these populations in a manner consistent with their investment in sexual reproduction. Populations engaging in a high frequency of sex were in Hardy-Weinberg and gametic phase equilibrium, and showed little genotypic differentiation across sampled years. In contrast, populations with a lower frequency of sex deviated widely from equilibrium, had reduced multilocus clonal diversity, and showed significant temporal genotypic deviation.

Speciation with gene flow and the genetics of habitat transitions

Whether speciation can advance to completion in the face of initially high levels of gene flow is a very controversial topic in evolutionary biology. Extensive gene exchange is generally considered to homogenize populations and counteract divergence. Moreover, the role of introgressive hybridization in evolution remains largely unexplored in animals, particularly in freshwater zooplankton in which allopatric speciation is considered to be the norm. Our work investigates the genetic structure of two young ecological species: the pond species, Daphnia pulex and the lake species, Daphnia pulicaria. Phylogenetic and population genetics analyses were conducted on mitochondrial NADH dehydrogenase 5 (ND5) gene, the nuclear Lactate dehydrogenase (Ldh) gene and 21 nuclear microsatellite markers in 416 individuals from habitats with various degrees of permanence. The strong and consistent phylogenetic discordance between nuclear and mitochondrial markers suggests a complex evolutionary history of multiple independent habitat transition events that involved hybridization and introgression between lake and pond Daphnia. On the other hand, the low level of contemporary gene flow between adjacent populations indicates the presence of effective habitat isolating barriers. The Daphnia system provides strong evidence for a divergence-with-gene flow speciation model that involves multiple habitat transition events.

Comparisons between Q(ST) and F(ST) --how wrong have we been?

The comparison between quantitative genetic divergence (Q(ST) ) and neutral genetic divergence (F(ST) ) among populations has become the standard test for historical signatures of selection on quantitative traits. However, when the mutation rate of neutral markers is relatively high in comparison with gene flow, estimates of F(ST) will decrease, resulting in upwardly biased comparisons of Q(ST) vs. F(ST) . Reviewing empirical studies, the difference between Q(ST) and F(ST) is positively related to marker heterozygosity. After refuting alternative explanations for this pattern, we conclude that marker mutation rate indeed has had a biasing effect on published Q(ST) -F(ST) comparisons. Hence, it is no longer clear that populations have commonly diverged in response to divergent selection. We present and discuss potential solutions to this bias. Comparing Q(ST) with recent indices of neutral divergence that statistically correct for marker heterozygosity (Hedrick's G'st and Jost's D) is not advised, because these indices are not theoretically equivalent to Q(ST) . One valid solution is to estimate F(ST) from neutral markers with mutation rates comparable to those of the loci underlying quantitative traits (e.g. SNPs). Q(ST) can also be compared to Φ(ST) (Phi(ST) ) of amova, as long as the genetic distance among allelic variants used to estimate Φ(ST) reflects evolutionary history: in that case, neutral divergence is independent of mutation rate. In contrast to their common usage in comparisons of Q(ST) and F(ST) , microsatellites typically have high mutation rates and do not evolve according to a simple evolutionary model, so are best avoided in Q(ST) -F(ST) comparisons.

Diet restriction and life history trade-offs in short- and long-lived species of Daphnia

The life-extending effects of diet restriction are well documented. One evolutionary model that accounts for this widespread conservation is the resource allocation model, where the selected individuals are those that can delay reproduction during periods of resource limitation. In this study, we use closely related species of a model organism, Daphnia, with widely divergent lifespans to address the relationship between diet restriction and longevity and assess whether the relationships are owing to trade-offs between reproductive and somatic investment. Specifically, we conducted a common garden experiment and constructed reaction norms for lifespan, fecundity, and body size as a function of food concentration. Our study provides evidence that the short-lived species in our study, D. pulex, shows the classically observed relationship of enhanced lifespan in response to reduced diet intake, but does not divert resources to somatic maintenance at the expense of reproduction during chronic diet restriction. In contrast, we find no evidence that the long-lived species in our study, D. pulicaria, gains any life-extending effects through diet restriction. Combined, our results provide evidence that the resource allocation model is not sufficient to explain the evolution of diet-mediated lifespan plasticity.

Does selection or genetic drift explain geographic differentiation of morphological characters in house sparrows Passer domesticus?

Understanding the relative influence of genetic drift and selection is fundamental in evolutionary biology. The theory of neutrality predicts that the genetic differentiation of a quantitative trait (QST) equals the genetic differentiation at neutral molecular markers (FST) if the quantitative trait has not been under selection. Thus, the relative magnitude of observed QST and expected QST under neutral expectations suggests the importance of selection and genetic drift for any observed phenotypic divergence. Because QST is based on additive genetic variance, estimating QST based on phenotypic measurements is problematic due to unknown environmental effects. To account for this, we used a model where the environmental component was allowed to vary when estimating QST. The model was used on data from 14 house sparrow (Passer domesticus) populations in Norway. In accordance with the significant phenotypic inter-population differences our analyses suggested that directional selection may have favoured different optimal phenotypes for some morphological traits across populations. In particular, different body mass and male ornamental phenotypes seemed to have been favoured. The conclusions are, however, dependent on assumptions regarding the proportion of the observed inter-population variation that is due to additive genetic differences, showing the importance of collecting such information in natural populations. By estimating QST, allowing the additive genetic proportion of phenotypic inter-population variation to vary, and by making use of recent statistical methods to compare observed QST with neutral expectations, we can use data that are relatively easy to collect to identify adaptive variation in natural populations.

Adaptation to local ultraviolet radiation conditions among neighbouring Daphnia populations

Understanding the historical processes that generated current patterns of phenotypic diversity in nature is particularly challenging in subdivided populations. Populations often exhibit heritable genetic differences that correlate with environmental variables, but the non-independence among neighbouring populations complicates statistical inference of adaptation. To understand the relative influence of adaptive and non-adaptive processes in generating phenotypes requires joint evaluation of genetic and phenotypic divergence in an integrated and statistically appropriate analysis. We investigated phenotypic divergence, population-genetic structure and potential fitness trade-offs in populations of Daphnia melanica inhabiting neighbouring subalpine ponds of widely differing transparency to ultraviolet radiation (UVR). Using a combination of experimental, population-genetic and statistical techniques, we separated the effects of shared population ancestry and environmental variables in predicting phenotypic divergence among populations. We found that native water transparency significantly predicted divergence in phenotypes among populations even after accounting for significant population structure. This result demonstrates that environmental factors such as UVR can at least partially account for phenotypic divergence. However, a lack of evidence for a hypothesized trade-off between UVR tolerance and growth rates in the absence of UVR prevents us from ruling out the possibility that non-adaptive processes are partially responsible for phenotypic differentiation in this system.

Selection on life-history traits and genetic population divergence in rotifers

A combination of founder effects and local adaptation - the Monopolization hypothesis - has been proposed to reconcile the strong population differentiation of zooplankton dwelling in ponds and lakes and their high dispersal abilities. The role genetic drift plays in genetic differentiation of zooplankton is well documented, but the impact of natural selection has received less attention. Here, we compare differentiation in neutral genetic markers (F(ST)) and in quantitative traits (Q(ST)) in six natural populations of the rotifer Brachionus plicatilis to assess the importance of natural selection in explaining genetic differentiation of life-history traits. Five life-history traits were measured in four temperature x salinity combinations in common-garden experiments. Population differentiation for neutral genetic markers - 11 microsatellite loci - was very high (F(ST) = 0.482). Differentiation in life-history traits was higher in traits related to sexual reproduction than in those related to asexual reproduction. Q(ST) values for diapausing egg production (a trait related to sexual reproduction) were higher than their corresponding F(ST) in some pairs of populations. Our results indicate the importance of divergent natural selection in these populations and suggest local adaptation to the unpredictability of B. plicatilis habitats.

Rapid evolution in response to introduced predators I: rates and patterns of morphological and life-history trait divergence

Background

Introduced species can have profound effects on native species, communities, and ecosystems, and have caused extinctions or declines in native species globally. We examined the evolutionary response of native zooplankton populations to the introduction of non-native salmonids in alpine lakes in the Sierra Nevada of California, USA. We compared morphological and life-history traits in populations of Daphnia with a known history of introduced salmonids and populations that have no history of salmonid introductions.

Results

Our results show that Daphnia populations co-existing with fish have undergone rapid adaptive reductions in body size and in the timing of reproduction. Size-related traits decreased by up to 13 percent in response to introduced fish. Rates of evolutionary change are as high as 4,238 darwins (0.036 haldanes).

Conclusion

Species introductions into aquatic habitats can dramatically alter the selective environment of native species leading to a rapid evolutionary response. Knowledge of the rates and limits of adaptation is an important component of understanding the long-term effects of alterations in the species composition of communities. We discuss the evolutionary consequences of species introductions and compare the rate of evolution observed in the Sierra Nevada Daphnia to published estimates of evolutionary change in ecological timescales.

Phylogenetic clustering and overdispersion in bacterial communities

Very little is known about the structure of microbial communities, despite their abundance and importance to ecosystem processes. Recent work suggests that bacterial biodiversity might exhibit patterns similar to those of plants and animals. However, relative to our knowledge about the diversity of macro-organisms, we know little about patterns of relatedness in free-living bacterial communities, and relatively few studies have quantitatively examined community structure in a phylogenetic framework. Here we apply phylogenetic tools to bacterial diversity data to determine whether bacterial communities are phylogenetically structured. We find that bacterial communities tend to contain lower taxonomic diversity and are more likely to be phylogenetically clustered than expected by chance. Such phylogenetic clustering may indicate the importance of habitat filtering (where a group of closely related species shares a trait, or suite of traits, that allow them to persist in a given habitat) in the assembly of bacterial communities. Microbial communities are especially accessible for phylogenetic analysis and thus have the potential to figure prominently in the integration of evolutionary biology and community ecology.

Recency, range expansion, and unsorted lineages: implications for interpreting neutral genetic variation in the sharp-tailed grouse (Tympanuchus phasianellus)

Both current and historical patterns of variation are relevant to understanding and managing ecological diversity. Recently derived species present a challenge to the reconstruction of historical patterns because neutral molecular data for these taxa are more likely to exhibit effects of recent and ongoing demographic processes. We studied geographical patterns of neutral molecular variation in a species thought to be of relatively recent origin, Tympanuchus phasianellus (sharp-tailed grouse), using mitochondrial control region sequences (CR-I), amplified fragment length polymorphisms (AFLP), and microsatellites. For historical context, we also analysed CR-I in all species of Tympanuchus. Within T. phasianellus, we found evidence for restricted gene flow between eastern and western portions of the species range, generally corresponding with the range boundary of T. p. columbianus and T. p. jamesi. The mismatch distribution and molecular clock estimates from the CR-I data suggested that all Tympanuchus underwent a range expansion prior to sorting of mitotypes among the species, and that sorting may have been delayed as a result of mutation-drift disequilibrium. This study illustrates the challenge of using genetic data to detect historical divergence in groups that are of relatively recent origin, or that have a history dominated by nonequilibrium conditions. We suggest that in such cases, morphological, ecological, and behavioural data may be particularly important adjuncts to molecular data for the recognition of historically or adaptively divergent groups.

Genetic variation and population differentiation in the lichen-forming ascomycete Xanthoria parietina on the island Storfosna, central Norway

Genetic diversity and fine-scale population structure in the lichen-forming ascomycete Xanthoria parietina was investigated using sequence variation in part of the intergenic spacer (IGS) and the complete internal transcribed spacer (ITS) regions of the nuclear ribosomal DNA. Sampling included 213 and 225 individuals, respectively, from seven populations in two different habitats, bark and rock, on the island Storfosna off the central west coast of Norway. Both markers revealed significant variation and a total of 10 IGS and 16 ITS haplotypes were found. There were no signs of significant positive spatial autocorrelation at any spatial size class down to 10% of transect length, nor did we find significant deviations from neutrality or signs of historical population expansion. Analysis of molecular variance (amova) indicated that most of the genetic variance observed was within populations, but when populations were grouped according to habitat, more than a quarter of the variance was explained among groups. Pairwise comparisons of populations (F(ST), exact tests of population differentiation) revealed significant differentiation between populations in different habitats (on bark or rock), but not between populations in the same habitat. Haplotype networks show that internal and presumably old haplotypes are shared between habitats, whereas terminal haplotypes tend to be unique to a habitat, mostly bark. We interpret the observed pattern to mean that there is no evidence of restricted gene flow between populations in the same habitat at the present spatial scale (interpopulation distances one or a few kilometres). On the other hand, differentiation between habitats is considerable, which we attribute to restricted gene flow between habitats (habitat isolation). Evidence suggests that the observed differentiation did not evolve locally. Estimates of divergence time between populations in the respective habitats indicate that an ancestral population started to diverge at least 34,000 years ago but probably much further back in time.

Distinguishing adaptive from nonadaptive genetic differentiation: comparison of QST and FST at two spatial scales

Genetic differentiation in 20 hierarchically sampled populations of wild barley was analyzed with quantitative traits, allozymes and Random Amplified Polymorphic DNAs (RAPDs), and compared for three marker types at two hierarchical levels. Regional subdivision for both molecular markers was much lower than for quantitative traits. For both allozymes and RAPDs, most loci exhibited minor or no regional differentiation, and the relatively high overall estimates of the latter were due to several loci with exceptionally high regional differentiation. The allozyme- and RAPD-specific patterns of differentiation were concordant in general with one another, but not with quantitative trait differentiation. Divergent selection on quantitative traits inferred from very high regional Q(ST) was in full agreement with our previous results obtained from a test of local adaptation and multilevel selection analysis. In contrast, most variation in allozyme and RAPD variation was neutral, although several allozyme loci and RAPD markers were exceptional in their levels of regional differentiation. However, it is not possible to answer the question whether these exceptional loci are directly involved in the response to selection pressure or merely linked to the selected loci. The fact that Q(ST) and F(ST) did not differ at the population scale, that is, within regions, but differed at the regional scale, for which local adaptation has been previously shown, implies that comparison of the level of subdivision in quantitative traits, as compared with molecular markers, is indicative of adaptive population differentiation only when sampling is carried out at the appropriate scale.

Bias and precision in QST estimates: problems and some solutions

Comparison of population differentiation in neutral marker genes and in genes coding quantitative traits by means of F(ST) and Q(ST) indexes has become commonplace practice. While the properties and estimation of F(ST) have been the subject of much interest, little is known about the precision and possible bias in Q(ST) estimates. Using both simulated and real data, we investigated the precision and bias in Q(ST) estimates and various methods of estimating the precision. We found that precision of Q(ST) estimates for typical data sets (i.e., with <20 populations) was poor. Of the methods for estimating the precision, a simulation method, a parametric bootstrap, and the Bayesian approach returned the most precise estimates of the confidence intervals.

Molecular and quantitative genetic divergence among populations of house mice with known evolutionary histories

Evolutionary biologists have long been interested in the processes influencing population differentiation, but separating the effects of neutral and adaptive evolution has been an obstacle for studies of population subdivision. A recently developed method allows tests of whether disruptive (ie, spatially variable) or stabilizing (ie, spatially uniform) selection is influencing phenotypic differentiation among subpopulations. This method, referred to as the F(ST) vs Q(ST) comparison, separates the total additive genetic variance into within- and among-population components and evaluates this level of differentiation against a neutral hypothesis. Thus, levels of neutral molecular (F(ST)) and quantitative genetic (Q(ST)) divergence are compared to evaluate the effects of selection and genetic drift on phenotypic differentiation. Although the utility of such comparisons appears great, its accuracy has not yet been evaluated in populations with known evolutionary histories. In this study, F(ST) vs Q(ST) comparisons were evaluated using laboratory populations of house mice with known evolutionary histories. In this model system, the F(ST) vs Q(ST) comparisons between the selection groups should reveal quantitative trait differentiation consistent with disruptive selection, while the F(ST) vs Q(ST) comparisons among lines within the selection groups should suggest quantitative trait differentiation in agreement with drift. We find that F(ST) vs Q(ST) comparisons generally produce the correct evolutionary inference at each level in the population hierarchy. Additionally, we demonstrate that when strong selection is applied between populations Q(ST) increases relative to Q(ST) among populations diverging by drift. Finally, we show that the statistical properties of Q(ST), a variance component ratio, need further investigation.

Mapping quantitative trait loci in multiple populations of Arabidopsis thaliana identifies natural allelic variation for trichome density

The majority of biological traits are genetically complex. Mapping the quantitative trait loci (QTL) that determine these phenotypes is a powerful means for estimating many parameters of the genetic architecture for a trait and potentially identifying the genes responsible for natural variation. Typically, such experiments are conducted in a single mapping population and, therefore, have only the potential to reveal genomic regions that are polymorphic between the progenitors of the population. What remains unclear is how well the QTL identified in any one mapping experiment characterize the genetics that underlie natural variation in traits. Here we provide QTL mapping data for trichome density from four recombinant inbred mapping populations of Arabidopsis thaliana. By aligning the linkage maps for these four populations onto a common physical map, the results from each experiment were directly compared. Seven of the nine QTL identified are population specific while two were mapped in all four populations. Our results show that many lineage-specific alleles that either increase or decrease trichome density persist in natural populations and that most of this genetic variation is additive. More generally, these findings suggest that the use of multiple populations holds great promise for better understanding the genetic architecture of natural variation.

Ecological and life history characteristics predict population genetic divergence of two salmonids in the same landscape

Ecological and life history characteristics such as population size, dispersal pattern, and mating system mediate the influence of genetic drift and gene flow on population subdivision. Bull trout (Salvelinus confluentus) and mountain whitefish (Prosopium williamsoni) differ markedly in spawning location, population size and mating system. Based on these differences, we predicted that bull trout would have reduced genetic variation within and greater differentiation among populations compared with mountain whitefish. To test this hypothesis, we used microsatellite markers to determine patterns of genetic divergence for each species in the Clark Fork River, Montana, USA. As predicted, bull trout had a much greater proportion of genetic variation partitioned among populations than mountain whitefish. Among all sites, FST was seven times greater for bull trout (FST = 0.304 for bull trout, 0.042 for mountain whitefish. After removing genetically differentiated high mountain lake sites for each species FST, was 10 times greater for bull trout (FST = 0.176 for bull trout; FST = 0.018 for mountain whitefish). The same characteristics that affect dispersal patterns in these species also lead to predictions about the amount and scale of adaptive divergence among populations. We provide a theoretical framework that incorporates variation in ecological and life history factors, neutral divergence, and adaptive divergence to interpret how neutral and adaptive divergence might be correlates of ecological and life history factors.

Latitudinal divergence of common frog (Rana temporaria) life history traits by natural selection: evidence from a comparison of molecular and quantitative genetic data

The relative roles of natural selection and direct environmental induction, as well as of natural selection and genetic drift, in creating clinal latitudinal variation in quantitative traits have seldom been assessed in vertebrates. To address these issues, we compared molecular and quantitative genetic differentiation between six common frog (Rana temporaria) populations along an approximately 1600 km long latitudinal gradient across Scandinavia. The degree of population differentiation (QST approximately 0.81) in three heritable quantitative traits (age and size at metamorphosis, growth rate) exceeded that in eight (neutral) microsatellite loci (FST = 0.24). Isolation by distance was clear for both neutral markers and quantitative traits, but considerably stronger for one of the three quantitative traits than for neutral markers. QST estimates obtained using animals subjected to different rearing conditions (temperature and food treatments) revealed some environmental dependency in patterns of population divergence in quantitative traits, but in general, these effects were weak in comparison to overall patterns. Pairwise comparisons of FST and QST estimates across populations and treatments revealed that the degree of quantitative trait differentiation was not generally predictable from knowledge of that in molecular markers. In fact, both positive and negative correlations were observed depending on conditions where the quantitative genetic variability had been measured. All in all, the results suggest a very high degree of genetic subdivision both in neutral marker genes and genes coding quantitative traits across a relatively recently (< 9000 years) colonized environmental gradient. In particular, they give evidence for natural selection being the primary agent behind the observed latitudinal differentiation in quantitative traits.

Divergent selection maintains adaptive differentiation despite high gene flow between sympatric rainbow smelt ecotypes (Osmerus mordax Mitchill)

In this study, we investigate the relative role of historical factors and evolutionary forces in promoting population differentiation in a new case of sympatric dwarf and normal ecotypes of the rainbow smelt (Osmerus mordax Mitchill) in Lac Saint-Jean (Québec, Canada). Our first objective was to test the hypothesis that the evolution of sympatric smelt ecotypes in Lac Saint-Jean has been contingent upon the secondary contact between two evolutionary lineages in postglacial times. Secondly, the QST method was applied to test the null hypothesis that the extent of phenotypic differences relative to that of neutral marker variation would be similar in comparisons involving populations within and among ecotypes. Thirdly, we applied a quantitative-genetic method as an exploratory assessment as to whether the amount of gene flow observed between populations could affect divergence in adaptive traits under specific conditions. This study revealed a unique situation of dwarf and normal smelt ecotypes that are, respectively, characterized by selmiparous and iteroparous life histories and the occurrence in each of two genetically distinct populations that synchronously use the same spawning habitat in two tributaries. Historical contingency has apparently played little role in the origin of these populations. In contrast, an important role of divergent natural selection in driving their phenotypic divergence was suggested. While divergent selection has apparently been strong enough to maintain phenotypic differentiation in the face of migration, this study suggests that gene flow has been sufficiently important to modulate the extent of adaptive differentiation being achieved between ecotypes, unless the extent of stabilizing selection acting on smelt ecotypes is much more pronounced than usually reported in natural populations.

Locating the sources of an invasive pest, grape phylloxera, using a mitochondrial DNA gene genealogy

Range expansions through human introductions have increased with global commerce and have led to the extinction of native species, alterations in community structure and pest status of the invasive species. Inferring the evolutionary history of invasive species can help to build a firmer footing for management tactics. This study used mitochondrial DNA (mtDNA) sequence comparisons of samples collected from the native and introduced ranges of a pest herbivore of cultivated grapes, grape phylloxera (Daktulosphaira vitifoliae Fitch, Phylloxeridae) to infer the sources and pattern of introductions into worldwide viticulture. Introductions into viticulture from its native North American range first occurred in the mid-19th century. The pattern of spread has suggested a focus of introduction into France, but independent introductions may have occurred elsewhere. The results show that the introduced population represents a limited subsample of the native genetic diversity. The data suggest that most grape phylloxera in viticulture, including all European, have originated in the northeastern USA where the grape species Vitis riparia dominates. There was evidence for independent introductions into South Africa and California. Most California haplotypes were most closely related to native grape phylloxera from the Atlantic Coast on V. vulpina. It is likely that subsequent spread from California into Australia, New Zealand and Peru has occurred.