Natural Selection and Neutral Evolution Jointly Drive Population Divergence between Alpine and Lowland Ecotypes of the Allopolyploid Plant Anemone multifida (Ranunculaceae)

Population differentiation can be driven in large part by natural selection, but selectively neutral evolution can play a prominent role in shaping patters of population divergence. The decomposition of the evolutionary history of populations into the relative effects of natural selection and selectively neutral evolution enables an understanding of the causes of population divergence and adaptation. In this study, we examined heterogeneous genomic divergence between alpine and lowland ecotypes of the allopolyploid plant, Anemone multifida. Using peak height and dominant AFLP data, we quantified population differentiation at non-outlier (neutral) and outlier loci to determine the potential contribution of natural selection and selectively neutral evolution to population divergence. We found 13 candidate loci, corresponding to 2.7% of loci, with signatures of divergent natural selection between alpine and lowland populations and between alpine populations (Fst  = 0.074–0.445 at outlier loci), but neutral population differentiation was also evident between alpine populations (F_ST  = 0.041–0.095 at neutral loci). By examining population structure at both neutral and outlier loci, we determined that the combined effects of selection and neutral evolution are associated with the divergence of alpine populations, which may be linked to extreme abiotic conditions and isolation between alpine sites. The presence of outlier levels of genetic variation in structured populations underscores the importance of separately analyzing neutral and outlier loci to infer the relative role of divergent natural selection and neutral evolution in population divergence.     

Genetic diversity at neutral and adaptive loci determines individual fitness in a long-lived territorial bird

There is compelling evidence about the manifest effects of inbreeding depression on individual fitness and populations' risk of extinction. The majority of studies addressing inbreeding depression on wild populations are generally based on indirect measures of inbreeding using neutral markers. However, the study of functional loci, such as genes of the major histocompatibility complex (MHC), is highly recommended. MHC genes constitute an essential component of the immune system of individuals, which is directly related to individual fitness and survival. In this study, we analyse heterozygosity fitness correlations of neutral and adaptive genetic variation (22 microsatellite loci and two loci of the MHC class II, respectively) with the age of recruitment and breeding success of a decimated and geographically isolated population of a long-lived territorial vulture. Our results indicate a negative correlation between neutral genetic diversity and age of recruitment, suggesting that inbreeding may be delaying reproduction. We also found a positive correlation between functional (MHC) genetic diversity and breeding success, together with a specific positive effect of the most frequent pair of cosegregating MHC alleles in the population. Globally, our findings demonstrate that genetic depauperation in small populations has a negative impact on the individual fitness, thus increasing the populations' extinction risk.     

A program to compare genetic differentiation statistics across loci using resampling of individuals and loci

Comparisons of genetic differentiation across populations based on different loci can provide insight into the evolutionary patterns acting on various regions of genomes. Here, we develop a program to statistically compare population genetic differentiation statistics (F(ST) or G'(ST) ) calculated from different loci. The program employs a routine that resamples either or both of individuals and loci and calculates a bootstrap confidence interval in the statistics. Resampling individuals is important when fewer than 25 individuals are sampled per population and when confidence intervals are required for individual loci. Resampling loci provides confidence intervals for sets of loci, such as a set presumed to be neutral, but can be anticonservative if fewer than 20 loci are analysed. We demonstrate the program using previously published data on the genetic differentiation at a major histocompatibility complex locus and at microsatellite loci across 10 populations of the guppy (Poecilia reticulata).     

The effect of linkage and population size on inbreeding depression due to mutational load.

Using a stochastic model of a finite population in which there is mutation to partially recessive detrimental alleles at many loci, we study the effects of population size and linkage between the loci on the population mean fitness and inbreeding depression values. Although linkage between the selected loci decreases the amount of inbreeding depression, neither population size nor recombination rate have strong effects on these quantities, unless extremely small values are assumed. We also investigate how partial linkage between the loci that determine fitness affects the invasion of populations by alleles at a modifier locus that controls the selfing rate. In most of the cases studied, the direction of selection on modifiers was consistent with that found in our previous deterministic calculations. However, there was some evidence that linkage between the modifier locus and the selected loci makes outcrossing less likely to evolve; more losses of alleles promoting outcrossing occurred in runs with linkage than in runs with free recombination. We also studied the fate of neutral alleles introduced into populations carrying detrimental mutations. The times to loss of neutral alleles introduced at low frequency were shorter than those predicted for alleles in the absence of selected loci, taking into account the reduction of the effective population size due to inbreeding. Previous studies have been confined to outbreeding populations, and to alleles at frequencies close to one-half, and have found an effect in the opposite direction. It therefore appears that associations between neutral and selected loci may produce effects that differ according to the initial frequencies of the neutral alleles.     

Microsatellite analysis of genetic diversity of the Vietnamese sika deer (Cervus nippon pseudaxis)

The Vietnamese sika deer (Cervus nippon pseudaxis) is an endangered subspecies of economic and traditional value in Vietnam. Most living individuals are held in traditional farms in central Vietnam, others being found in zoos around the world. Here we study the neutral genetic diversity and population structure of this subspecies using nine microsatellite loci in order to evaluate the consequences of the limited number of individuals from which this population was initiated and of the breeding practices (i.e., possible inbreeding). Two hundred individuals were sampled from several villages. Our data show both evidence for limited local inbreeding and isolation by distance with a mean F(ST) value of 0.02 between villages. This suggests that exchange of animals occurs at a local scale, at a rate such that highly inbred mating is avoided. However, the genetic diversity, with an expected heterozygosity (H(e)) of 0.60 and mean number of alleles (k) of 5.7, was not significantly larger than that estimated from zoo populations of much smaller census size (17 animals sampled; H(e) = 0.65, k = 4.11). Our results also suggest that the Vietnamese population might have experienced a slight bottleneck. However, this population is sufficiently variable to constitute a source of individuals for reintroduction in the wild in Vietnam.     

Using population genetics and demographic reconstruction to predict outcomes of genetic rescue for an endangered songbird

Genetic rescue can be a successful way to restore species genetic diversity, but it can also lead to outbreeding depression (decreases in hybrid fitness) if attempted in incompatible populations. Thus, population genetic profiles and demographic history are needed to evaluate the feasibility of translocation. We used population genetic analyses and Approximate Bayesian Computation (ABC) to assess genetic rescue as an option for two populations of the yellow-shouldered blackbird (Agelaius xanthomus), an endangered Puerto Rico endemic. The candidate recipient population, a managed population in Pitahaya (southwestern Puerto Rico), had been characterized previously for its mating system and population genetics. Here, we used nine microsatellite loci to measure the genetic diversity of a candidate source population, a subspecies (A. x. monensis) on Mona Island, 66 km west of Puerto Rico. We compared genetic diversity and inferred historical and contemporary gene flow between the two populations. We found clear population structure and no migration between populations, as well as evidence that the Mona population descended from the Pitahaya population approximately 95 generations ago. Despite the historical gene flow, the degree of contemporary genetic and environmental divergence means the Mona population may not be suitable for immediate use as a source population. We recommend (a) further investigating whether the observed population structure is due to adaptive or neutral forces, (b) testing the Mona population for behavioral plasticity in different environments, and (c) evaluating other source populations in addition to the Mona population for genetic rescue.     

Patterns of differentiation in a colour polymorphism and in neutral markers reveal rapid genetic changes in natural damselfly populations

The existence and mode of selection operating on heritable adaptive traits can be inferred by comparing population differentiation in neutral genetic variation between populations (often using F(ST) values) with the corresponding estimates for adaptive traits. Such comparisons indicate if selection acts in a diversifying way between populations, in which case differentiation in selected traits is expected to exceed differentiation in neutral markers [F(ST )(selected) > F(ST )(neutral)], or if negative frequency-dependent selection maintains genetic polymorphisms and pulls populations towards a common stable equilibrium [F(ST) (selected) < F(ST) (neutral)]. Here, we compared F(ST) values for putatively neutral data (obtained using amplified fragment length polymorphism) with estimates of differentiation in morph frequencies in the colour-polymorphic damselfly Ischnura elegans. We found that in the first year (2000), population differentiation in morph frequencies was significantly greater than differentiation in neutral loci, while in 2002 (only 2 years and 2 generations later), population differentiation in morph frequencies had decreased to a level significantly lower than differentiation in neutral loci. Genetic drift as an explanation for population differentiation in morph frequencies could thus be rejected in both years. These results indicate that the type and/or strength of selection on morph frequencies in this system can change substantially between years. We suggest that an approach to a common equilibrium morph frequency across all populations, driven by negative frequency-dependent selection, is the cause of these temporal changes. We conclude that inferences about selection obtained by comparing F(ST) values from neutral and adaptive genetic variation are most useful when spatial and temporal data are available from several populations and time points and when such information is combined with other ecological sources of data.     

Interpretation of variation across marker loci as evidence of selection

Population structure and history have similar effects on the genetic diversity at all neutral loci. However, some marker loci may also have been strongly influenced by natural selection. Selection shapes genetic diversity in a locus-specific manner. If we could identify those loci that have responded to selection during the divergence of populations, then we may obtain better estimates of the parameters of population history by excluding these loci. Previous attempts were made to identify outlier loci from the distribution of sample statistics under neutral models of population structure and history. Unfortunately these methods depend on assumptions about population structure and history that usually cannot be verified. In this article, we define new population-specific parameters of population divergence and construct sample statistics that are estimators of these parameters. We then use the joint distribution of these estimators to identify outlier loci that may be subject to selection. We found that outlier loci are easier to recognize when this joint distribution is conditioned on the total number of allelic states represented in the pooled sample at each locus. This is so because the conditional distribution is less sensitive to the values of nuisance parameters.     

Beyond the beneficial effects of translocations as an effective tool for the genetic restoration of isolated populations

Translocations are becoming increasingly popular as appropriate management strategies for the genetic restoration of endangered species and populations. Although a few studies have shown that the introduction of novel alleles has reversed the detrimental effects of inbreeding over the short-term (i.e., genetic rescue), it is not clear how effective such translocations are for both maintaining neutral variation that may be adaptive in the future (i.e., genetic restoration) and increasing population viability over the long-term. In addition, scientists have expressed concerns regarding the potential genetic swamping of locally adapted populations, which may eliminate significant components of genetic diversity through the replacement of the target population by the source individuals used for translocations. Here we show that bird translocations into a wild population of greater prairie-chickens (Tympanuchus cupido pinnatus) in southeastern Illinois were effective in both removing detrimental variation associated with inbreeding depression as well as restoring neutral genetic variation to historical levels. Furthermore, we found that although translocations resulted in immediate increases in fitness, the demographic recovery and long-term viability of the population appears to be limited by the availability of suitable habitat. Our results demonstrate that although translocations can be effective management tools for the genetic restoration of wild populations on the verge of extinction, their long-term viability may not be guaranteed unless the initial conditions that led to most species declines (e.g., habitat loss) are reversed.     

Genetic rescue in interconnected populations of small and large size of the self-incompatible Ranunculus reptans

Small populations of our study species Ranunculus reptans have reduced fitness because of inbreeding, genetic load, and reduced mate availability; that is, they suffer from a three-fold genetic Allee effect. Here, we investigate how the effect of interpopulation outbreeding on offspring fitness depends on population size. We performed within- and between-population crosses with plants originating from 15 populations, and measured offspring performance in a common environment. Interpopulation outbreeding led to an increase in population means of clonal performance, which was defined as the number of rooted offspring rosettes produced per maternal ovule. This fitness gain mainly occurred at the life stage of seed set. It was especially pronounced for populations with a long-term history of small size inferred from their low genetic diversity, estimated from eight allozyme loci. We conclude that in a self-incompatible plant such as R. reptans, interpopulation outbreeding can lead to an immediate genetic rescue effect due to increased cross-compatibility and heterosis, and that this rescue effect is increased as population size decreases.     

TEMPORAL VARIATION AT THE MHC CLASS IIB IN WILD POPULATIONS OF THE GUPPY (POECILIA RETICULATA)

Understanding genetic diversity in natural populations is a fundamental objective of evolutionary biology. The immune genes of the major histocompatibility complex (MHC) are excellent candidates to study such diversity because they are highly polymorphic in populations. Although balancing selection may be responsible for maintaining diversity at these functionally important loci, temporal variation in selection pressure has rarely been examined. We examine temporal variation in MHC class IIB diversity in nine guppy (Poecilia reticulata) populations over two years. We found that five of the populations changed significantly more at the MHC than at neutral (microsatellite) loci as measured by F_ST, which suggests that the change at the MHC was due to selection and not neutral processes. Additionally, pairwise population differentiation measures at the MHC were higher in 2007 than in 2006, with the signature of selection changing from homogenizing to diversifying selection or neutral evolution. Interestingly, within the populations the magnitude of the change at the MHC between years was related to the change in the proportion of individuals infected by a common parasite, indicating a link between genetic structure and the parasite. Our data thereby implicate temporal variation in selective pressure as an important mechanism maintaining diversity at the MHC in wild populations.     

Inbreeding depression and low between-population heterosis in recently diverged experimental populations of a selfing species

In fragmented populations, genetic drift and selection reduce genetic diversity, which in turn results in a loss of fitness or in a loss of evolvability. Genetic rescue, that is, controlled input of diversity from distant populations, may restore evolutionary potential, whereas outbreeding depression might counteract the positive effect of this strategy. We carried out self-pollination and crosses within and between populations in an experimental subdivided population of a selfing species, Triticum aestivum L., to estimate the magnitude of these two phenomena. Surprisingly, for a self-fertilizing species, we found significant inbreeding depression within each population for four of the six traits studied, indicating that mildly deleterious mutations were still segregating in these populations. The progeny of within- and between-population crosses was very similar, indicating low between-population heterosis and little outbreeding depression. We conclude that relatively large population effective sizes prevented fixation of a high genetic load and that local adaptation was limited in these recently diverged populations. The kinship coefficient estimated between the parents using 20 neutral markers was a poor predictor of the progeny phenotypic values, indicating that there was a weak link between neutral diversity and genes controlling fitness-related traits. These results show that when assessing the viability of natural populations and the need for genetic rescue, the use of neutral markers should be complemented with information about the presence of local adaptation in the subdivided population.     

The probability distribution under a population divergence model of the number of genetic founding lineages of a population or species

The composition of genetic variation in a population or species is shaped by the number of events that led to the founding of the group. We consider a neutral coalescent model of two populations, where a derived population is founded as an offshoot of an ancestral population. For a given locus, using both recursive and nonrecursive approaches, we compute the probability distribution of the number of genetic founding lineages that have given rise to the derived population. This number of genetic founding lineages is defined as the number of ancestral individuals that contributed at the locus to the present-day derived population, and is formulated in terms of interspecific coalescence events. The effects of sample size and divergence time on the probability distribution of the number of founding lineages are studied in detail. For 99.99% of the loci in the derived population to each have one founding lineage, the two populations must be separated for 9.9N generations. However, only approximately 0.87N generations must pass since divergence for 99.99% of the loci to have <6 founding lineages. Our results are useful as a prior expectation on the number of founding lineages in scenarios that involve the evolution of one population from the splitting of an ancestral group, such as in the colonization of islands, the formation of polyploid species, and the domestication of crops and livestock from wild ancestors.     

Genetic analysis of microsatellite polymorphism in the Elliot's Pheasant (Syrmaticus ellioti) in China

In this study, we reported the population genetic analyses in the Elliot's Pheasant(Syrnaticus ellioti) using seven polymorphism microsatellite loci based on 105 individuals from 4 geographical populations. Departures from Hardy-Weinberg equilibrium were found in four geographical populations. The average number of alleles was 8.86, with a total of 62 alleles across 7 loci; observed heterozygosity (HO) was generally low and the average number was 0.504. For the seven microsatellite loci, the polymorphism information content ranged from 0.549 to 0.860, with an average number 0.712. Population bottlenecks of the four geographical populations were tested by infinite allele mutation model, step-wise mutation model and two-phase mutation model, which found that each population had experienced bottleneck effect during the recent period. Fst analysis across all geographical populations indicated that the genetic differentiaton between the Guizhou geographical population and the Hunan geographical population was highly significant (P<0.001), a finding supported by the far genetic relationship showed by the neighbor-joining tree of four geographical populations based on Nei's unbiased genetic distances. Using hierarchical analysis of molecular variance (Guizhou geographical population relative to all others pooled), we found a low level of the genetic variation among geographical populations and that between groups. However, differences among populations relative to the total sample explained most of the genetic variance (92.84%), which was significant.     

Population genetic analysis of 15 autosomal STR loci in the Russian population of northeastern Inner-Mongolia, China

In the present study, we investigated the genetic polymorphisms of 15 autosomal STR loci in the Russian population of northeastern Inner-Mongolia, China as well as its genetic relationships with other populations. DNA typing for 15 autosomal STR loci was performed on 148 randomly selected healthy individuals from the Russian population living in Eerguna, northeastern Inner-Mongolia. Allelic frequencies of these loci were calculated by direct counting. The genotype data of this Russian population was moreover compared to other populations using neighbor-joining method, as such constructing a phylogenic tree. A total of 143 alleles were found in the Russian population with corresponding allele frequencies in the range from 0.0034 to 0.5372. Among all the 15 loci, D18S51 had the highest polymorphism (PIC = 0.8632), whereas TPOX had the lowest (PIC = 0.5179). In the phylogenic tree, this Russian population has a close relationship with the populations of South Siberia and northeastern Asia. This study may increase our understanding of the genetic background of the Russian population in Eerguna, China as such providing useful information for anthropological research, forensic sciences as well as disease-association studies.     

Heterozygosity and survival in blue tits (Cyanistes caeruleus): contrasting effects of presumably functional and neutral loci

The relationship between genetic diversity and fitness has important implications in evolutionary and conservation biology. This relationship has been widely investigated at the individual level in studies of heterozygosity-fitness correlations (HFC). General effects caused by inbreeding and/or local effects at single loci have been used as explanations of HFC, but the debate about the causes of HFC in open, natural populations is still ongoing. Study designs that control for variation in the inbreeding level of the individuals, and knowledge on the function and location of the markers used to measure heterozygosity, are fundamental to understand the causes of HFC. Here we investigated correlations between individual heterozygosity and estimates of survival at different life-history stages in an open population of blue tits (Cyanistes caeruleus). For survival at the embryo, nestling and fledgling stage, we used a full-sibling approach, i.e. we controlled for the level of inbreeding. We genotyped 1496 individuals with 79 microsatellites mapped across 25 chromosomes in the zebra finch (Taeniopygia guttata) that were classified either as potentially functional (58 loci) or as neutral (21 loci). We found different effects of standardized multilocus heterozygosity (SH): SH(functional) had a negative effect on the probability of hatching and local recruitment of females, whereas SH(neutral) had a positive effect on adult survival. The negative effects of functional loci are better explained by local effects, whereas the positive effects of neutral markers could reflect inbreeding effects in the population. Our results highlight the importance of considering the characteristics of the markers used in HFC studies and confirm the mixed effects of heterozygosity in different contexts (e.g. sex and life-history stage).     

The relationship of nucleotide polymorphism, recombination rate and selection in wild tomato species

We analyzed the effects of mating system and recombination rate on single nucleotide polymorphisms using 14 single-copy nuclear loci from single populations of five species of wild tomatoes (Solanum section Lycopersicon). The taxa investigated comprise two self-compatible (SC) and three self-incompatible (SI) species. The observed reduction in nucleotide diversity in the SC populations compared to the SI populations is much stronger than expected under the neutral effects of the mating system on effective population size. Importantly, outgroup sequences available for 11 of the 14 loci yield strong positive correlations between silent nucleotide diversity and silent divergence, indicative of marked among-locus differences in mutation rates and/or selective constraints. Furthermore, using a physical estimate of local recombination rates, we find that silent nucleotide diversity (but not divergence) is positively correlated with recombination rate in two of the SI species. However, this correlation is not nearly as strong as in other well-characterized species (in particular, Drosophila). We propose that nucleotide diversity in Lycopersicon is dominated mainly by differences in neutral mutation rates and/or selective constraints among loci, demographic processes (such as population subdivision), and background selection. In addition, we hypothesize that the soil seed bank plays an important role in the maintenance of the large genetic diversity in the SI species (in particular L. peruvianum).     

Reduced MHC class II diversity in island compared to mainland populations of the black-footed rock-wallaby (<Emphasis Type="Italic">Petrogale lateralis lateralis</Emphasis>)

Many animal populations that are endangered in mainland areas exist in stable island populations, which have the potential to act as an “ark” in case of mainland population declines. Previous studies have found neutral genetic variation in such species to be up to an order of magnitude lower in island compared to mainland populations. If low genetic variation is prevalent across fitness-related loci, this would reduce the effectiveness of island populations as a source of individuals to supplement declining mainland populations or re-establish extinct mainland populations. One such species, the black-footed rock-wallaby (Petrogale lateralis lateralis), exists within fragmented mainland populations and small island populations off Western Australia. We examined sequence variation in this species within a fitness-related locus under positive selection, the MHC class II DAB β1 locus. The mainland populations displayed greater levels of allelic diversity (4–7 alleles) than the island population, despite being small and isolated, and contained at least two DAB gene copies. The island population displayed low allelic diversity (2 alleles) and fewer alleles per individual in comparison to mainland populations, and probably possesses only one DAB gene copy. The patterns of DAB diversity suggested that the island population has a markedly lower level of genetic variation than the mainland populations, in concordance with results from microsatellites (genotyped in a previous study), but preserved unique alleles which were not found in mainland populations. Where possible, conservation actions should pool individuals from multiple populations, not only island populations, for translocation programs, and focus on preventing further declines in mainland populations.     

Can balancing selection on MHC loci counteract genetic drift in small fragmented populations of black grouse?

The ability of natural populations to adapt to new environmental conditions is crucial for their survival and partly determined by the standing genetic variation in each population. Populations with higher genetic diversity are more likely to contain individuals that are better adapted to new circumstances than populations with lower genetic diversity. Here, we use both neutral and major histocompatibility complex (MHC) markers to test whether small and highly fragmented populations hold lower genetic diversity than large ones. We use black grouse as it is distributed across Europe and found in populations with varying degrees of isolation and size. We sampled 11 different populations; five continuous, three isolated, and three small and isolated. We tested patterns of genetic variation in these populations using three different types of genetic markers: nine microsatellites and 21 single nucleotide polymorphisms (SNPs) which both were found to be neutral, and two functional MHC genes that are presumably under selection. The small isolated populations displayed significantly lower neutral genetic diversity compared to continuous populations. A similar trend, but not as pronounced, was found for genotypes at MHC class II loci. Populations were less divergent at MHC genes compared to neutral markers. Measures of genetic diversity and population genetic structure were positively correlated among microsatellites and SNPs, but none of them were correlated to MHC when comparing all populations. Our results suggest that balancing selection at MHC loci does not counteract the power of genetic drift when populations get small and fragmented.     

A strong genetic footprint of the re‐introduction history of Alpine ibex (Capra ibex ibex)

A population’s neutral genetic variation is a composite of its size, degree of isolation and demographic history. Bottlenecks and founder events increase genetic drift, leading to the loss of genetic variation and increased genetic differentiation among populations. Gene flow has the opposite effects. Thus, gene flow can override the genetic patterns caused by founder events. Using 37 microsatellite loci, we investigated the effects of serial bottlenecks on genetic variation and differentiation among 42 Alpine ibex populations in Switzerland with known re‐introduction histories. We detected a strong footprint of re‐introduction events on contemporary genetic structure, with re‐introduction history explaining a substantial part of the genetic differentiation among populations. As a result of the translocation of a considerable number of individuals from the sole formerly surviving population in northern Italy, most of the genetic variation of the ancestral population is now present in the combined re‐introduced Swiss populations. However, re‐introductions split up the genetic variation among populations, such that each contemporary Swiss population showed lower genetic variation than the ancestral population. As expected, serial bottlenecks had different effects on the expected heterozygosity (He) and standardized number of alleles (sNa). While loss of sNa was higher in the first bottlenecks than in subsequent ones, He declined to a similar degree with each bottleneck. Thus, genetic drift was detected with each bottleneck, even when no loss of sNa was observed. Overall, more than a hundred years after the beginning of this successful re‐introduction programme, re‐introduction history was the main determinant of today’s genetic structure.