Neutral Genetic Differentiation Research Articles

Adaptation and colonization history affect the evolution of clines in two introduced species

Phenotypic and genetic clines have long been synonymous with adaptive evolution. However, other processes (for example, migration, range expansion, invasion) may generate clines in traits or loci across geographical and environmental gradients. It is therefore important to distinguish between clines that represent adaptive evolution and those that result from selectively neutral demographic or genetic processes. We tested for the differentiation of phenotypic traits along environmental gradients using two species in the genus Silene, whilst statistically controlling for colonization history and founder effects. We sampled seed families from across the native and introduced ranges, genotyped individuals and estimated phenotypic differentiation in replicated common gardens. The results suggest that post-glacial expansion of S. vulgaris and S. latifolia involved both neutral and adaptive genetic differentiation (clines) of life history traits along major axes of environmental variation in Europe and North America. Phenotypic clines generally persisted when tested against the neutral expectation, although some clines disappeared (and one cline emerged) when the effects of genetic ancestry were statistically removed. Colonization history, estimated using genetic markers, is a useful null model for tests of adaptive trait divergence, especially during range expansion and invasion when selection and gene flow may not have reached equilibrium.

The relative role of drift and selection in shaping the human skull

Human populations across the world vary greatly in cranial morphology. It is highly debated to what extent this variability has accumulated through neutral processes (genetic drift) or through natural selection driven by climate. By taking advantage of recent work showing that geographic distance along landmasses is an excellent proxy for neutral genetic differentiation, we quantify the relative role of drift versus selection in an exceptionally large dataset of human skulls. We show that neutral processes have been much more important than climate in shaping the human cranium. We further demonstrate that a large proportion of the signal for natural selection comes from populations from extremely cold regions. More generally, we show that, if drift is not explicitly accounted for, the effect of natural selection can be greatly overestimated.

Comparison of quantitative and molecular genetic variation of native vs. invasive populations of purple loosestrife (Lythrum salicaria L., Lythraceae)

Study of adaptive evolutionary changes in populations of invasive species can be advanced through the joint application of quantitative and population genetic methods. Using purple loosestrife as a model system, we investigated the relative roles of natural selection, genetic drift and gene flow in the invasive process by contrasting phenotypical and neutral genetic differentiation among native European and invasive North American populations (Q(ST) - F(ST) analysis). Our results indicate that invasive and native populations harbour comparable levels of amplified fragment length polymorphism variation, a pattern consistent with multiple independent introductions from a diverse European gene pool. However, it was observed that the genetic variation reduced during subsequent invasion, perhaps by founder effects and genetic drift. Comparison of genetically based quantitative trait differentiation (Q(ST)) with its expectation under neutrality (F(ST)) revealed no evidence of disruptive selection (Q(ST) > F(ST)) or stabilizing selection (Q(ST) < F(ST)). One exception was found for only one trait (the number of stems) showing significant sign of stabilizing selection across all populations. This suggests that there are difficulties in distinguishing the effects of nonadaptive population processes and natural selection. Multiple introductions of purple loosestrife may have created a genetic mixture from diverse source populations and increased population genetic diversity, but its link to the adaptive differentiation of invasive North American populations needs further research.

Are QST–FST comparisons for natural populations meaningful?

Comparisons between putatively neutral genetic differentiation amongst populations, F(ST), and quantitative genetic variation, Q(ST), are increasingly being used to test for natural selection. However, we find that approximately half of the comparisons that use only data from wild populations confound phenotypic and genetic variation. We urge the use of a clear distinction between narrow-sense Q(ST), which can be meaningfully compared with F(ST), and phenotypic divergence measured between populations, P(ST), which is inadequate for comparisons in the wild. We also point out that an unbiased estimate of Q(ST) can be found using the so-called 'animal model' of quantitative genetics.

Development of variable microsatellite loci and range-wide characterization of nuclear genetic diversity in the important dryland shrub antelope bitterbrush (Purshia tridentata)

Antelope bitterbrush (Purshia tridentata Pursh DC; Rosaceae) is an arid-land shrub that occupies an important ecological niche in various fire-dominated communities across much of the western United States. Because of its importance as a browse for large mammals and a food source for granivores, P. tridentata is frequently planted by Federal agencies in arid-land revegetation. We are currently analyzing the range-wide genetic diversity of this species as part of a larger effort to develop seed movement guidelines. In this study, we describe the development of eight novel nuclear microsatellite loci and characterize the amount and apportionment of range-wide nuclear genetic diversity. The eight microsatellite loci exhibited a high level of polymorphism (13–33 alleles per locus) and in general, observed levels of heterozygosity did not deviate form Hardy–Weinberg expectations. An initial screen of 196 individuals from 12 widely distributed populations revealed a moderate amount of genetic differentiation (FST=0.09, p<0.001). Moreover, these loci successfully produced PCR products in cross-species amplifications with two closely related Purshia species. These results demonstrate the utility of these markers and provide useful measures of neutral genetic diversity and population differentiation.

A Genome-Scan Method to Identify Selected Loci Appropriate for Both Dominant and Codominant Markers: A Bayesian Perspective

Identifying loci under natural selection from genomic surveys is of great interest in different research areas. Commonly used methods to separate neutral effects from adaptive effects are based on locus-specific population differentiation coefficients to identify outliers. Here we extend such an approach to estimate directly the probability that each locus is subject to selection using a Bayesian method. We also extend it to allow the use of dominant markers like AFLPs. It has been shown that this model is robust to complex demographic scenarios for neutral genetic differentiation. Here we show that the inclusion of isolated populations that underwent a strong bottleneck can lead to a high rate of false positives. Nevertheless, we demonstrate that it is possible to avoid them by carefully choosing the populations that should be included in the analysis. We analyze two previously published data sets: a human data set of codominant markers and a Littorina saxatilis data set of dominant markers. We also perform a detailed sensitivity study to compare the power of the method using amplified fragment length polymorphism (AFLP), SNP, and microsatellite markers. The method has been implemented in a new software available at our website (http://www-leca.ujf-grenoble.fr/logiciels.htm).

Speciation with gene flow in the large white-headed gulls: does selection counterbalance introgression?

We investigated the role of selection in generating and maintaining species distinctness in spite of ongoing gene flow, using two zones of secondary contact between large gull species in Europe (Larus argentatus and Larus cachinnans) and North America (Larus glaucescens and Larus occidentalis). We used the pattern of neutral genetic differentiation at nine microsatellite loci (F(ST)) as an indicator of expected changes under neutral processes and compared it with phenotypic differentiation (P(ST)) for a large number of traits (size, plumage melanism and coloration of bare parts). Even assuming very low heritability, interspecific divergence between L. glaucescens and L. occidentalis in plumage melanism and orbital ring colour clearly exceeded neutral differentiation. Similarly, melanism of the central primaries was highly divergent between L. argentatus and L. cachinnans. Such divergence is unlikely to have arisen randomly and is therefore attributed to spatially varying selection. Variation in plumage melanism in both transects agrees with Gloger's rule, which suggests that latitude (and associated sun and humidity gradients) could be the selective pressure shaping differentiation in plumage melanism. We suggest that strong species differentiation in orbital ring colour results from sexual selection. We conclude that these large gull species, along with other recently diverged species that hybridize after coming into secondary contact, may differ only in restricted regions of the genome that are undergoing strong disruptive selection because of their phenotypic effects.

HETEROGENEOUS GENOMIC DIFFERENTIATION BETWEEN WALKING-STICK ECOTYPES: "ISOLATION BY ADAPTATION" AND MULTIPLE ROLES FOR DIVERGENT SELECTION

Genetic differentiation can be highly variable across the genome. For example, loci under divergent selection and those tightly linked to them may exhibit elevated differentiation compared to neutral regions. These represent "outlier loci" whose differentiation exceeds neutral expectations. Adaptive divergence can also increase genome-wide differentiation by promoting general barriers to neutral gene flow, thereby facilitating genomic divergence via genetic drift. This latter process can yield a positive correlation between adaptive phenotypic divergence and neutral genetic differentiation (described here as "isolation-by-adaptation"). Here, we examine both these processes by combining an AFLP genome scan of two host plant ecotypes of Timema cristinae walking-sticks with existing data on adaptive phenotypic divergence and ecological speciation in these insects. We found that about 8% of loci are outliers in multiple population comparisons. Replicated comparisons between population-pairs using the same versus different host species revealed that 1-2% of loci are subject to host-related selection specifically. Locus-specific analyses revealed that up to 10% of putatively neutral (nonoutlier) AFLP loci exhibit significant isolation-by-adaptation. Our results suggest that selection may affect differentiation directly, via linkage, or by facilitating genetic drift. They thus illustrate the varied and sometimes nonintuitive contributions of selection to heterogeneous genomic differentiation.

Genetic isolation of fragmented populations is exacerbated by drift and selection

Reduced genetic variation at marker loci in small populations has been well documented, whereas the relationship between quantitative genetic variation and population size has attracted little empirical investigation. Here we demonstrate that both neutral and quantitative genetic variation are reduced in small populations of a fragmented plant metapopulation, and that both drift and selective change are enhanced in small populations. Measures of neutral genetic differentiation (F(ST)) and quantitative genetic differentiation (Q(ST)) in two traits were higher among small demes, and Q(ST) between small populations exceeded that expected from drift alone. This suggests that fragmented populations experience both enhanced genetic drift and divergent selection on phenotypic traits, and that drift affects variation in both neutral markers and quantitative traits. These results highlight the need to integrate natural selection into conservation genetic theory, and suggests that small populations may represent reservoirs of genetic variation adaptive within a wide range of environments.

Adaptive vs. neutral genetic diversity: implications for landscape genetics

Genetic diversity is important for the maintenance of the viability and the evolutionary or adaptive potential of populations and species. However, there are two principal types of genetic diversity: adaptive and neutral – a fact widely neglected by non-specialists. We introduce these two types of genetic diversity and critically point to their potential uses and misuses in population or landscape genetic studies. First, most molecular-genetic laboratory techniques analyse neutral genetic variation. This means that the gene variants detected do not have any direct effect on fitness. This type of genetic variation is thus selectively neutral and tells us nothing about the adaptive or evolutionary potential of a population or a species. Nevertheless, neutral genetic markers have great potential for investigating processes such as gene flow, migration or dispersal. Hence, they allow us to empirically test the functional relevance of spatial indices such as connectivity used in landscape ecology. Second, adaptive genetic variation, i.e. genetic variation under natural selection, is analysed in quantitative genetic experiments under controlled and uniform environmental conditions. Unfortunately, the genetic variation (i.e. heritability) and population differentiation at quantitative, adaptive traits is not directly linked with neutral genetic diversity or differentiation. Thus, neutral genetic data cannot serve as a surrogate of adaptive genetic data. In summary, neutral genetic diversity is well suited for the study of processes within landscapes such as gene flow, while the evolutionary or adaptive potential of populations or species has to be assessed in quantitative genetic experiments. Landscape ecologists have to mind these differences between neutral and adaptive genetic variation when interpreting the results of landscape genetic studies.

Distinct population structure in a phenotypically homogeneous rock‐dwelling cichlid fish from Lake Tanganyika

Several lineages of cichlid fishes in the East African Great Lakes display stunning levels of morphological diversification. The rapid evolution of rock-dwelling polygynous mouthbrooders in Lake Malawi, for example, was in part ascribed to their allopatric distribution on disjunct stretches of rocky coast, where even short habitat discontinuities reduce gene flow effectively. However, as seen in other cichlids, ecological barriers do not always prevent gene flow, whereas genetic structure can develop along continuous habitat, and morphological diversification does not necessarily accompany genetic differentiation. The present study investigates the population structure of Variabilichromis moorii, a monogamous substrate-brooding lamprologine of rocky coasts in Lake Tanganyika, which occurs over about 1000 km of shoreline almost without phenotypic variation. Phylogeographic analyses of mitochondrial DNA sequences indicated that dispersal is infrequent and generally occurs between adjacent locations only. Exceptions to this pattern are closely related haplotypes from certain locations on opposite lakeshores, a phenomenon which has been observed in other species and is thought to reflect lake crossing along an underwater ridge in times of low water level. Genetic population differentiation, estimated from mitochondrial DNA and microsatellite data in six adjacent populations, was equally high across localities separated by sandy shores and along uninterrupted stretches of rocky shore. Our results suggest that ecological barriers are not required to induce philopatric behavior in Variabilichromis, and that morphological stasis persists in the face of high levels of neutral genetic differentiation.

CONTRASTING PATTERNS OF FLORAL AND MOLECULAR VARIATION ACROSS A CLINE IN MIMULUS AURANTIACUS

Steep clines in ecologically important traits may be caused by divergent natural selection. However, processes that do not necessarily invoke ongoing selection, such as secondary contact or restricted gene flow, can also cause patterns of phenotypic differentiation over short spatial scales. Distinguishing among all possible scenarios is difficult, but an attainable goal is to establish whether scenarios that imply selection need to be invoked. We compared the extent of morphological and genetic differentiation between geographically structured red and yellow floral races of Mimulus aurantiacus (bush monkeyflower; Phrymaceae). Flower color was assessed in a common garden as well as in the field to determine whether variation was genetic and to quantify the extent of geographical differentiation. Population genetic differentiation at marker loci was measured for both chloroplast and nuclear genomes, and the degree of population structure within and among the floral races was evaluated. Flower color shows both a strong genetic basis and a sharp geographic transition, with pure red-flowered populations in western San Diego County and pure yellow-flowered populations to the east. In the zone of contact, both pure and intermediate phenotypes occur. Patterns of genetic differentiation at marker loci are far less pronounced, as little of the variation is partitioned according to the differences in flower color. Phenotypic differentiation (Q(ST)) between populations with different flower colors is much greater than neutral genetic differentiation (F(ST)). When comparisons are made between populations of the same flower color, the opposite trend is evident. Limited neutral genetic structure between the floral races, combined with sharp differentiation in flower color, is consistent with the hypothesis that current or recent natural selection maintains the cline in flower color.

Maternal genetic effects on adaptive divergence between anadromous and resident brook charr during early life history

The importance of directional selection relative to neutral evolution may be determined by comparing quantitative genetic variation in phenotype (Q(ST)) to variation at neutral molecular markers (F(ST)). Quantitative divergence between salmonid life history types is often considerable, but ontogenetic changes in the significance of major sources of genetic variance during post-hatch development suggest that selective differentiation varies by developmental stage. In this study, we tested the hypothesis that maternal genetic differentiation between anadromous and resident brook charr (Salvelinus fontinalis Mitchill) populations for early quantitative traits (embryonic size/growth, survival, egg number and developmental time) would be greater than neutral genetic differentiation, but that the maternal genetic basis for differentiation would be higher for pre-resorption traits than post-resorption traits. Quantitative genetic divergence between anadromous (seawater migratory) and resident Laval River (Québec) brook charr based on maternal genetic variance was high (Q(ST) > 0.4) for embryonic length, yolk sac volume, embryonic growth rate and time to first response to feeding relative to neutral genetic differentiation [F(ST) = 0.153 (0.071-0.214)], with anadromous females having positive genetic coefficients for all of the above characters. However, Q(ST) was essentially zero for all traits post-resorption of the yolk sac. Our results indicate that the observed divergence between resident and anadromous brook charr has been driven by directional selection, and may therefore be adaptive. Moreover, they provide among the first evidence that the relative importance of selective differentiation may be highly context-specific, and varies by genetic contributions to phenotype by parental sex at specific points in offspring ontogeny. This in turn suggests that interpretations of Q(ST)-F(ST) comparisons may be improved by considering the structure of quantitative genetic architecture by age category and the sex of the parent used in estimation.

CONTRASTING PATTERNS OF FLORAL AND MOLECULAR VARIATION ACROSS A CLINE IN MIMULUS AURANTIACUS

Steep clines in ecologically important traits may be caused by divergent natural selection. However, processes that do not necessarily invoke ongoing selection, such as secondary contact or restricted gene flow, can also cause patterns of phenotypic differentiation over short spatial scales. Distinguishing among all possible scenarios is difficult, but an attainable goal is to establish whether scenarios that imply selection need to be invoked. We compared the extent of morphological and genetic differentiation between geographically structured red and yellow floral races of Mimulus aurantiacus (bush monkeyflower; Phrymaceae). Flower color was assessed in a common garden as well as in the field to determine whether variation was genetic and to quantify the extent of geographical differentiation. Population genetic differentiation at marker loci was measured for both chloroplast and nuclear genomes, and the degree of population structure within and among the floral races was evaluated. Flower color shows both a strong genetic basis and a sharp geographic transition, with pure red-flowered populations in western San Diego County and pure yellow-flowered populations to the east. In the zone of contact, both pure and intermediate phenotypes occur. Patterns of genetic differentiation at marker loci are far less pronounced, as little of the variation is partitioned according to the differences in flower color. Phenotypic differentiation (QST) between populations with different flower colors is much greater than neutral genetic differentiation (FST). When comparisons are made between populations of the same flower color, the opposite trend is evident. Limited neutral genetic structure between the floral races, combined with sharp differentiation in flower color, is consistent with the hypothesis that current or recent natural selection maintains the cline in flower color.

CONTRASTING PATTERNS OF QUANTITATIVE AND NEUTRAL GENETIC VARIATION IN LOCALLY ADAPTED POPULATIONS OF THE NATTERJACK TOAD, BUFO CALAMITA

The relative importance of natural selection and genetic drift in determining patterns of phenotypic diversity observed in nature is still unclear. The natterjack toad (Bufo calamita) is one of a few amphibian species capable of breeding in saline ponds, even though water salinity represents a considerable stress for them. Results from two common-garden experiments showed a pattern of geographic variation in embryonic salinity tolerance among populations from either fresh or brackish environments, consistent with the hypothesis of local adaptation. Full-sib analysis showed increased variation in survival among sibships within population for all populations as osmotic stress was increased (broad-sense heritability increased as salinity raised). Nevertheless, toads native to the brackish water environment had the highest overall survival under brackish conditions. Levels of population genetic differentiation for salinity tolerance were higher than those of neutral genetic differentiation, the latter obtained through the analysis of eight microsatellite loci. Microsatellite markers also revealed little population differentiation, lack of an isolation-by-distance pattern, and moderate gene flow connecting the populations. Therefore, environmental stress tolerance appears to have evolved in absence of geographic isolation, and consequently we reject the null hypothesis of neutral differentiation.

CONTRASTING PATTERNS OF QUANTITATIVE AND NEUTRAL GENETIC VARIATION IN LOCALLY ADAPTED POPULATIONS OF THE NATTERJACK TOAD, BUFO CALAMITA

The relative importance of natural selection and genetic drift in determining patterns of phenotypic diversity observed in nature is still unclear. The natterjack toad (Bufo calamita) is one of a few amphibian species capable of breeding in saline ponds, even though water salinity represents a considerable stress for them. Results from two common-garden experiments showed a pattern of geographic variation in embryonic salinity tolerance among populations from either fresh or brackish environments, consistent with the hypothesis of local adaptation. Full-sib analysis showed increased variation in survival among sibships within population for all populations as osmotic stress was increased (broad-sense heritability increased as salinity raised). Nevertheless, toads native to the brackish water environment had the highest overall survival under brackish conditions. Levels of population genetic differentiation for salinity tolerance were higher than those of neutral genetic differentiation, the latter obtained through the analysis of eight microsatellite loci. Microsatellite markers also revealed little population differentiation, lack of an isolation-by-distance pattern, and moderate gene flow connecting the populations. Therefore, environmental stress tolerance appears to have evolved in absence of geographic isolation, and consequently we reject the null hypothesis of neutral differentiation.

Neutral genetic differentiation in an island model with cyclical parthenogenesis

Abstract Unusually high levels of genetic differentiation are often observed between populations of Daphnia (Crustacea: Cladocera) and other cladocerans. Selection and departure from migration–mutation–drift equilibrium have been invoked to explain this fact. However, the null model of neutral genes at equilibrium has not been explicitly stated. Using a simple island model that takes into account the main characteristics of the cyclical parthenogenetic life-cycle, we show that a high level of differentiation can be obtained for neutral genes at equilibrium if the migration rate during the asexual phase is low and the effective size over the season is reduced, or if mating within clones is common. Recurrent population bottlenecks have a large effect on differentiation and the time to reach equilibrium. Such fluctuating clonal dynamics could be brought about by limitations in the number of clones that hatch and establish themselves each season, or by random associations between neutral and selected genes, as previously suggested in the literature.

Species distribution in marine phytoplankton : Braarud T., 1962. J. Oceanogr. Soc., Jap., 20th Ann. Vol.: 628–649

Understanding the extent to which neutral processes and adaptive divergence shape the spatial structure of natural populations is a major goal in evolutionary biology and is especially important for the identification of significant levels of biodiversity. Our results identified replicated habitat-specific (adaptive) phenotypic divergence in the brown macroalga Fucus vesiculosus that is independent of population (neutral) genetic structure. F. vesiculosus inhabits contiguous and contrasting marine to estuarine intertidal habitats. Combining analyses of genetic and phenotypic traits of populations living under differential selective regimes (estuaries and open coast), we investigated levels of neutral genetic differentiation and adaptive physiological responses to emersion stress. In southwest England (SW UK) and northern Iberia (N. Iberia), populations living in estuaries and marine coastal habitats were genetically characterized at six microsatellite loci. In N. Iberia, two clades with limited admixture were recovered, each including one open coast site and the adjacent estuarine location. In contrast, SW UK samples clustered according to habitat and formed three distinct groups of genotypes; one including the two open coast locations and the other two representing each of the estuarine sites. Temperature loggers revealed distinct emersion regimes that characterized each habitat type independently of the region, while water and air temperature profiles showed site-specific trends. Despite acclimation under usual conditions, trait means of emersion stress resilience showed a strong phenotypic divergence between habitats, consistent with environmental clines in exposure time observed in the different habitats. We demonstrate that neutral genetic clusters do not reflect locally adapted population units. Our results identified replicated habitat-specific (adaptive) phenotypic divergence that is independent of population (neutral) genetic structure in F. vesiculosus. The significance of such findings extends beyond the theoretical evolutionary and ecological interest of discovering parallel adaptive responses to the broader implications for conservation of intraspecific biodiversity.