Contrasting heterozygosity-fitness correlations between populations of a self-compatible shrub in a fragmented landscape

The mechanisms underlying heterozygosity-fitness correlations (HFCs) are subject of intense debates, especially about how important population features such as size or degree of isolation influence HFCs. Here, we report variation in HFCs between Large and Small populations of a self-compatible shrub (Myrtus communis) occurring within an extremely fragmented landscape. In each of the five study populations, we obtained data on both heterozygosity and fitness for 9-12 maternal families (i.e. offspring from the same mother plant). Whereas heterozygosity explained most of the variance (60-86?%) in growth rate of seedling families within Large populations, this relationship was absent within Small populations. Our results suggest that inbreeding may explain the observed HFCs within Large populations, and that different genetic processes (such as genetic drift and/or selection) could have overridden HFCs within Small populations. While it is difficult to draw general conclusions from five populations, we think our results open new research perspectives on how different genetic processes underlie variation in HFCs under different population contexts. Our study also points to a need for further attention on the complex relationships between heterozygosity in self-compatible plants and their progeny in relation to mating system variation. Finally, our results provide interesting new insights into how population genetic diversity is maintained or lost in a highly fragmented landscape.     

Towards genetic markers in animal populations as biomonitors for human-induced environmental change

Genetic markers provide potentially sensitive indicators of changes in environmental conditions because the genetic constitution of populations is normally altered well before populations become extinct. Genetic indicators in populations include overall genetic diversity, genetic changes in traits measured at the phenotypic level, and evolution at specific loci under selection. While overall genetic diversity has rarely been successfully related to environmental conditions, genetically based changes in traits have now been linked to the presence of toxins and both local and global temperature shifts. Candidate loci for monitoring stressors are emerging from information on how specific genes influence traits, and from screens of random loci across environmental gradients. Drosophila research suggests that chromosomal regions under recent intense selection can be identified from patterns of molecular variation and a high frequency of transposable element insertions. Allele frequency changes at candidate loci have been linked to pesticides, pollutants and climate change. Nevertheless, there are challenges in interpreting allele frequencies in populations, particularly when a large number of loci control a trait and when interactions between alleles influence trait expression. To meet these challenges, population samples should be collected for longitudinal studies, and experimental programmes should be undertaken to link variation at candidate genes to ecological processes.     

DOES LINKAGE DISEQUILIBRIUM GENERATE HETEROZYGOSITY‐FITNESS CORRELATIONS IN GREAT REED WARBLERS?

Heterozygosity-fitness correlations (HFCs) at noncoding genetic markers are commonly assumed to reflect fitness effects of heterozygosity at genomewide distributed genes in partially inbred populations. However, in populations with much linkage disequilibrium (LD), HFCs may arise also as a consequence of selection on fitness loci in the local chromosomal vicinity of the markers. Recent data suggest that relatively high levels of LD may prevail in many ecological situations. Consequently, LD may be an important factor, together with partial inbreeding, in causing HFCs in natural populations. In the present study, we evaluate whether LD can generate HFCs in a small and newly founded population of great reed warblers (Acrocephalus arundinaceus). For this purpose dyads of full siblings of which only one individual survived to adult age (i.e., returned to breed at the study area) were scored at 19 microsatellite loci, and at a gene region of hypothesized importance for survival, the major histocompatibility complex (MHC). By examining siblings, we controlled for variation in the inbreeding coefficient and thus excluded genome-wide fitness effects in our analyses. We found that recruited individuals had significantly higher multilocus heterozygosity (MLH), and mean d2 (a microsatellite-specific variable), than their nonrecruited siblings. There was a tendency for the survivors to have a more diverse MHC than the nonsurvivors. Single-locus analyses showed that the strength of the genotype-survival association was especially pronounced at four microsatellite loci. By using genotype data from the entire breeding population, we detected significant LD between five of 162 pairs of microsatellite loci after accounting for multiple tests. Our present finding of a significant within-family multilocus heterozygosity-survival association in a nonequilibrium population supports the view that LD generates HFCs in natural populations.     

Local heterozygosity-fitness correlations with global positive effects on fitness in threespine stickleback

The complex interactions between genetic diversity and evolution have important implications in many biological areas including conservation, speciation, and mate choice. A common way to study these interactions is to look at heterozygosity-fitness correlations (HFCs). Until recently, HFCs based on noncoding markers were believed to result primarily from global inbreeding effects. However, accumulating theoretical and empirical evidence shows that HFCs may often result from genes being linked to the markers used (local effect). Moreover, local effect HFCs could differ from global inbreeding effects in their direction and occurrence. Consequently, the investigation of the structure and consequences of local HFCs is emerging as a new important goal in evolutionary biology. In this study of a wild threespine stickleback (Gasterosteus aculeatus) population, we first tested the presence of significant positive or negative local effects of heterozygosity at 30 microsatellites loci on five fitness components: survival, mating success, territoriality, length, and body condition. Then, we evaluated the direction and shape of total impact of local HFCs, and estimated the magnitude of the impacts on fitness using regression coefficients and selection differentials. We found that multilocus heterozygosity was not a reliable estimator of individual inbreeding coefficient, which supported the relevance of single-locus based analyses. Highly significant and temporally stable local HFCs were observed. These were mainly positive, but negative effects of heterozygosity were also found. Strong and opposite effects of heterozygosity are probably present in many populations, but may be blurred in HFC analyses looking for global effects only. In this population, both negative and positive HFCs are apparently driving mate preference by females, which is likely to contribute to the maintenance of both additive and nonadditive genetic variance.     

Unraveling the Effects of Selection and Demography on Immune Gene Variation in Free-Ranging Plains Zebra (Equus quagga) Populations

Demography, migration and natural selection are predominant processes affecting the distribution of genetic variation among natural populations. Many studies use neutral genetic markers to make inferences about population history. However, the investigation of functional coding loci, which directly reflect fitness, is critical to our understanding of species'' ecology and evolution. Immune genes, such as those of the Major Histocompatibility Complex (MHC), play an important role in pathogen recognition and provide a potent model system for studying selection. We contrasted diversity patterns of neutral data with MHC loci, ELA-DRA and -DQA, in two southern African plains zebra (Equus quagga) populations: Etosha National Park, Namibia, and Kruger National Park, South Africa. Results from neutrality tests, along with observations of elevated diversity and low differentiation across populations, supported previous genus-level evidence for balancing selection at these loci. Despite being low, MHC divergence across populations was significant and may be attributed to drift effects typical of geographically separated populations experiencing little to no gene flow, or alternatively to shifting allele frequency distributions driven by spatially variable and fluctuating pathogen communities. At the DRA, zebra exhibited geographic differentiation concordant with microsatellites and reduced levels of diversity in Etosha due to highly skewed allele frequencies that could not be explained by demography, suggestive of spatially heterogeneous selection and local adaptation. This study highlights the complexity in which selection affects immune gene diversity and warrants the need for further research on the ecological mechanisms shaping patterns of adaptive variation among natural populations.     

Genetic architecture and balancing selection: the life and death of differentiated variants

Balancing selection describes any form of natural selection, which results in the persistence of multiple variants of a trait at intermediate frequencies within populations. By offering up a snapshot of multiple co‐occurring functional variants and their interactions, systems under balancing selection can reveal the evolutionary mechanisms favouring the emergence and persistence of adaptive variation in natural populations. We here focus on the mechanisms by which several functional variants for a given trait can arise, a process typically requiring multiple epistatic mutations. We highlight how balancing selection can favour specific features in the genetic architecture and review the evolutionary and molecular mechanisms shaping this architecture. First, balancing selection affects the number of loci underlying differentiated traits and their respective effects. Control by one or few loci favours the persistence of differentiated functional variants by limiting intergenic recombination, or its impact, and may sometimes lead to the evolution of supergenes. Chromosomal rearrangements, particularly inversions, preventing adaptive combinations from being dissociated are increasingly being noted as features of such systems. Similarly, due to the frequency of heterozygotes maintained by balancing selection, dominance may be a key property of adaptive variants. High heterozygosity and limited recombination also influence associated genetic load, as linked recessive deleterious mutations may be sheltered. The capture of deleterious elements in a locus under balancing selection may reinforce polymorphism by further promoting heterozygotes. Finally, according to recent genomewide scans, balanced polymorphism might be more pervasive than generally thought. We stress the need for both functional and ecological studies to characterize the evolutionary mechanisms operating in these systems.     

Fluctuating selection and the determinants of genetic variation

Recent studies of cosmopolitan Drosophila populations have found hundreds to thousands of genetic loci with seasonally fluctuating allele frequencies, bringing temporally fluctuating selection to the forefront of the historical debate surrounding the maintenance of genetic variation in natural populations. Numerous mechanisms have been explored in this longstanding area of research, but these exciting empirical findings have prompted several recent theoretical and experimental studies that seek to better understand the drivers, dynamics, and genome-wide influence of fluctuating selection. In this review, we evaluate the latest evidence for multilocus fluctuating selection in Drosophila and other taxa, highlighting the role of potential genetic and ecological mechanisms in maintaining these loci and their impacts on neutral genetic variation.     

Rapid genome‐wide evolution in Brassica rapa populations following drought revealed by sequencing of ancestral and descendant gene pools

There is increasing evidence that evolution can occur rapidly in response to selection. Recent advances in sequencing suggest the possibility of documenting genetic changes as they occur in populations, thus uncovering the genetic basis of evolution, particularly if samples are available from both before and after selection. Here, we had a unique opportunity to directly assess genetic changes in natural populations following an evolutionary response to a fluctuation in climate. We analysed genome‐wide differences between ancestors and descendants of natural populations of Brassica rapa plants from two locations that rapidly evolved changes in multiple phenotypic traits, including flowering time, following a multiyear late‐season drought in California. These ancestor‐descendant comparisons revealed evolutionary shifts in allele frequencies in many genes. Some genes showing evolutionary shifts have functions related to drought stress and flowering time, consistent with an adaptive response to selection. Loci differentiated between ancestors and descendants (F_ST outliers) were generally different from those showing signatures of selection based on site frequency spectrum analysis (Tajima's D), indicating that the loci that evolved in response to the recent drought and those under historical selection were generally distinct. Very few genes showed similar evolutionary responses between two geographically distinct populations, suggesting independent genetic trajectories of evolution yielding parallel phenotypic changes. The results show that selection can result in rapid genome‐wide evolutionary shifts in allele frequencies in natural populations, and highlight the usefulness of combining resurrection experiments in natural populations with genomics for studying the genetic basis of adaptive evolution.     

Controlling for P‐value inflation in allele frequency change in experimental evolution and artificial selection experiments

Experimental evolution studies can be used to explore genomic response to artificial and natural selection. In such studies, loci that display larger allele frequency change than expected by genetic drift alone are assumed to be directly or indirectly associated with traits under selection. However, such studies report surprisingly many loci under selection, suggesting that current tests for allele frequency change may be subject to P‐value inflation and hence be anticonservative. One factor known from genomewide association (GWA) studies to cause P‐value inflation is population stratification, such as relatedness among individuals. Here, we suggest that by treating presence of an individual in a population after selection as a binary response variable, existing GWA methods can be used to account for relatedness when estimating allele frequency change. We show that accounting for relatedness like this effectively reduces false‐positives in tests for allele frequency change in simulated data with varying levels of population structure. However, once relatedness has been accounted for, the power to detect causal loci under selection is low. Finally, we demonstrate the presence of P‐value inflation in allele frequency change in empirical data spanning multiple generations from an artificial selection experiment on tarsus length in two free‐living populations of house sparrow and correct for this using genomic control. Our results indicate that since allele frequencies in large parts of the genome may change when selection acts on a heritable trait, such selection is likely to have considerable and immediate consequences for the eco‐evolutionary dynamics of the affected populations.     

Experimental evidence for ecological selection on genome variation in the wild

Understanding natural selection's effect on genetic variation is a major goal in biology, but the genome‐scale consequences of contemporary selection are not well known. In a release and recapture field experiment we transplanted stick insects to native and novel host plants and directly measured allele frequency changes within a generation at 186 576 genetic loci. We observed substantial, genome‐wide allele frequency changes during the experiment, most of which could be attributed to random mortality (genetic drift). However, we also documented that selection affected multiple genetic loci distributed across the genome, particularly in transplants to the novel host. Host‐associated selection affecting the genome acted on both a known colour‐pattern trait as well as other (unmeasured) phenotypes. We also found evidence that selection associated with elevation affected genome variation, although our experiment was not designed to test this. Our results illustrate how genomic data can identify previously underappreciated ecological sources and phenotypic targets of selection.     

Genetic differentiation between alpine and lowland populations of a butterfly is related to PGI enzyme genotype

Understanding the ecological process of population differentiation and identifying the molecular changes that contribute to adaptation is a central aspect of evolutionary biology. In this study we analyzed the geographic variation in allozyme allele frequencies (based on 15 enzyme systems representing 18 loci) across 18 populations of the butterfly Lycaena tityrus from different altitudes. Population genetic analyses showed that within population genetic diversity (e.g. mean number of alleles per loci: 1.8; expected heterozygosity: 12%) was within the typical value range for the Lepidoptera. The populations of L. tityrus investigated showed a remarkable genetic differentiation (F_ST: 0.065), being clearly divided into an alpine (high-altitude) and a non-alpine (low-altitude) cluster. This differentiation was almost entirely caused by variation at a single enzyme locus, PGI. Although the involvement of historical events cannot be ruled out, several lines of evidence suggest that the specific pattern of allozyme (and in this case particularly PGI) variation found is caused by thermal selection. For example, genetic variation was highly locus-specific rather than relatively uniform, as should be expected for effects of natural selection. Further, the PGI 2–2 genotype dominating in alpine (in contrast to low-altitude) populations is known to exhibit increased cold stress resistance and relatively long development times typical of alpine populations. Thus, PGI (or possibly a closely linked gene) is an obvious target for thermal selection in L. tityrus . This study exemplifies how allozyme analyses can be used to detect candidate loci likely to be under selection.     

Allele frequency divergence reveals ubiquitous influence of positive selection in Drosophila

Resolving the role of natural selection is a basic objective of evolutionary biology. It is generally difficult to detect the influence of selection because ubiquitous non-selective stochastic change in allele frequencies (genetic drift) degrades evidence of selection. As a result, selection scans typically only identify genomic regions that have undergone episodes of intense selection. Yet it seems likely such episodes are the exception; the norm is more likely to involve subtle, concurrent selective changes at a large number of loci. We develop a new theoretical approach that uncovers a previously undocumented genome-wide signature of selection in the collective divergence of allele frequencies over time. Applying our approach to temporally resolved allele frequency measurements from laboratory and wild Drosophila populations, we quantify the selective contribution to allele frequency divergence and find that selection has substantial effects on much of the genome. We further quantify the magnitude of the total selection coefficient (a measure of the combined effects of direct and linked selection) at a typical polymorphic locus, and find this to be large (of order 1%) even though most mutations are not directly under selection. We find that selective allele frequency divergence is substantially elevated at intermediate allele frequencies, which we argue is most parsimoniously explained by positive—not negative—selection. Thus, in these populations most mutations are far from evolving neutrally in the short term (tens of generations), including mutations with neutral fitness effects, and the result cannot be explained simply as an ongoing purging of deleterious mutations.     

The effects of genetic drift during range expansion on geographical patterns of variation: a computer simulation of the colonization of Australia by Bufo marinus

A computer simulation is performed of allele frequency changes resulting from genetic drift at eleven loci during the probable course of colonization of Australia by populations of the Giant Toad, Bufo marinus. The history of twelve populations for which allele frequency data are available is modelled. Account is taken of the likely pattern of relationship among the populations, the effective size of the populations (as indicated by the observed variance of allele frequencies) and the probable isolation of the populations from each other following their separation. In all simulations, allele frequencies at some loci show a significant association with latitude while others do not. In six of ten simulations, there is a significant association between degree of variation at a locus and the presence of a latitudinal cline of allele frequencies. There are also indications of this kind of association in simulations of a uni-directional range expansion. These results demonstrate that such associations, which were also observed in the data from the actual populations, can result from genetic drift during a range expansion, and therefore cannot be taken as evidence of the action of natural selection.     

Evolution of quantitative traits in the wild: mind the ecology

Recent advances in the quantitative genetics of traits in wild animal populations have created new interest in whether natural selection, and genetic response to it, can be detected within long-term ecological studies. However, such studies have re-emphasized the fact that ecological heterogeneity can confound our ability to infer selection on genetic variation and detect a population''s response to selection by conventional quantitative genetics approaches. Here, I highlight three manifestations of this issue: counter gradient variation, environmentally induced covariance between traits and the correlated effects of a fluctuating environment. These effects are symptomatic of the oversimplifications and strong assumptions of the breeder''s equation when it is applied to natural populations. In addition, methods to assay genetic change in quantitative traits have overestimated the precision with which change can be measured. In the future, a more conservative approach to inferring quantitative genetic response to selection, or genomic approaches allowing the estimation of selection intensity and responses to selection at known quantitative trait loci, will provide a more precise view of evolution in ecological time.     

Temporal changes in allele frequencies in two reciprocally selected maize populations

The effects of breeding on allele frequency changes at 82 restriction fragment length polymorphism (RFLP) loci were examined in two maize (Zea mays L.) populations undergoing reciprocal recurrent selection, Iowa Stiff Stalk Synthetic and Iowa Corn Borer Synthetic #1. After 12 cycles of selection, approximately 30% of the alleles were extinct and 10% near fixation in each population. A test of selective neutrality identified several loci in each population whose allele frequency changes cannot be explained by genetic drift; interpopulation mean expected heterozygosity increased for that subset of 28 loci but not for the remaining 54 loci. Mean expected heterozygosity within the two subpopulations decreased 39%, while the between-population component of genetic variation increased from 0.5% to 33.4% of the total. Effective population size is a key parameter for discerning allele frequency changes due to genetic drift versus those resulting from selection and genetic hitchhiking. Empirical estimates of effective population size for each population were within the range predicted by the breeding method. Received: 10 June 1998 / Accepted: 29 April 1999     

Heterozygosity at a single locus explains a large proportion of variation in two fitness‐related traits in great tits: a general or a local effect?

In natural populations, mating between relatives can have important fitness consequences due to the negative effects of reduced heterozygosity. Parental level of inbreeding or heterozygosity has been also found to influence the performance of offspring, via direct and indirect parental effects that are independent of the progeny own level of genetic diversity. In this study, we first analysed the effects of parental heterozygosity and relatedness (i.e. an estimate of offspring genetic diversity) on four traits related to offspring viability in great tits (Parus major) using 15 microsatellite markers. Second, we tested whether significant heterozygosity–fitness correlations (HFCs) were due to ‘local’ (i.e. linkage to genes influencing fitness) and/or ‘general’ (genome‐wide heterozygosity) effects. We found a significant negative relationship between parental genetic relatedness and hatching success, and maternal heterozygosity was positively associated with offspring body size. The characteristics of the studied populations (recent admixture, polygynous matings) together with the fact that we found evidence for identity disequilibrium across our set of neutral markers suggest that HFCs may have resulted from genome‐wide inbreeding depression. However, one locus (Ase18) had disproportionately large effects on the observed HFCs: heterozygosity at this locus had significant positive effects on hatching success and offspring size. It suggests that this marker may lie near to a functional locus under selection (i.e. a local effect) or, alternatively, heterozygosity at this locus might be correlated to heterozygosity across the genome due to the extensive ID found in our populations (i.e. a general effect). Collectively, our results lend support to both the general and local effect hypotheses and reinforce the view that HFCs lie on a continuum from inbreeding depression to those strictly due to linkage between marker loci and genes under selection.     

Demographic Genetics of Brown Trout (Salmo Trutta) and Estimation of Effective Population Size from Temporal Change of Allele Frequencies

We studied temporal allele frequency shifts over 15 years and estimated the genetically effective size of four natural populations of brown trout (Salmo trutta L.) on the basis of the variation at 14 polymorphic allozyme loci. The allele frequency differences between consecutive cohorts were significant in all four populations. There were no indications of natural selection, and we conclude that random genetic drift is the most likely cause of temporal allele frequency shifts at the loci examined. Effective population sizes were estimated from observed allele frequency shifts among cohorts, taking into consideration the demographic characteristics of each population. The estimated effective sizes of the four populations range from 52 to 480 individuals, and we conclude that the effective size of natural brown trout populations may differ considerably among lakes that are similar in size and other apparent characteristics. In spite of their different effective sizes all four populations have similar levels of genetic variation (average heterozygosity) indicating that excessive loss of genetic variability has been retarded, most likely because of gene flow among neighboring populations.     

Heterozygosity is unrelated to adult fitness measures in a large, noninbred population of great tits (Parus major)

The extent to which heterozygosity-fitness correlations (HFCs) are expected in wild populations is an important and unresolved question in evolutionary biology, because it relates to our understanding of the genetic architecture of fitness. Here, we report a study of HFCs in a wild, noninbred population of great tits (Parus major), based on a sample comprising 281 individuals typed at 26 markers, resulting in a data set comprising over 5600 genotypes. We regressed pedigree-derived f-score and multilocus genetic diversity against eight life-history traits known to be associated with fitness in this population, including lifetime reproductive success (LRS), as well as several morphological traits under weak selection. We found no evidence for either multilocus or single-locus HFCs for any morphological or fitness trait, and further found no evidence that effect sizes were stronger for those life-history traits more closely associated with reproductive fitness. This result may, in part, be explained by the fact that we found no evidence that our set of 26 markers had any power to infer genome-wide heterozygosity in this population and that marker-derived heterozygosity was uncorrelated with pedigree-derived f-score. Overall, these results emphasize the fact that the often-reported strong HFCs detected in small, inbred populations do not reflect a general phenomenon of increasing individual reproductive fitness with increasing heterozygosity.     

Detecting selection along environmental gradients: analysis of eight methods and their effectiveness for outbreeding and selfing populations

Thanks to genome‐scale diversity data, present‐day studies can provide a detailed view of how natural and cultivated species adapt to their environment and particularly to environmental gradients. However, due to their sensitivity, up‐to‐date studies might be more sensitive to undocumented demographic effects such as the pattern of migration and the reproduction regime. In this study, we provide guidelines for the use of popular or recently developed statistical methods to detect footprints of selection. We simulated 100 populations along a selective gradient and explored different migration models, sampling schemes and rates of self‐fertilization. We investigated the power and robustness of eight methods to detect loci potentially under selection: three designed to detect genotype–environment correlations and five designed to detect adaptive differentiation (based on F_ST or similar measures). We show that genotype–environment correlation methods have substantially more power to detect selection than differentiation‐based methods but that they generally suffer from high rates of false positives. This effect is exacerbated whenever allele frequencies are correlated, either between populations or within populations. Our results suggest that, when the underlying genetic structure of the data is unknown, a number of robust methods are preferable. Moreover, in the simulated scenario we used, sampling many populations led to better results than sampling many individuals per population. Finally, care should be taken when using methods to identify genotype–environment correlations without correcting for allele frequency autocorrelation because of the risk of spurious signals due to allele frequency correlations between populations.     

Cline coupling and uncoupling in a stickleback hybrid zone

Strong ecological selection on a genetic locus can maintain allele frequency differences between populations in different environments, even in the face of hybridization. When alleles at divergent loci come into tight linkage disequilibrium, selection acts on them as a unit and can significantly reduce gene flow. For populations interbreeding across a hybrid zone, linkage disequilibria between loci can force clines to share the same slopes and centers. However, strong ecological selection on a locus can also pull its cline away from the others, reducing linkage disequilibrium and weakening the barrier to gene flow. We looked for this “cline uncoupling” effect in a hybrid zone between stream resident and anadromous sticklebacks at two genes known to be under divergent natural selection (Eda and ATP1a1) and five morphological traits that repeatedly evolve in freshwater stickleback. These clines were all steep and located together at the top of the estuary, such that we found no evidence for cline uncoupling. However, we did not observe the stepped shape normally associated with steep concordant clines. It thus remains possible that these clines cluster together because their individual selection regimes are identical, but this would be very surprising given their diverse roles in osmoregulation, body armor, and swimming performance.