THE OPPORTUNITY FOR BALANCING SELECTION IN EXPERIMENTAL POPULATIONS OF CAENORHABDITIS ELEGANS

The role of balancing selection in maintaining diversity during the evolution of sexual populations to novel environments is poorly understood. To address this issue, we studied the impact of two mating systems, androdioecy and dioecy, on genotype distributions during the experimental evolution of Caenorhabditis elegans. We analyzed the temporal trajectories of 334 single nucleotide polymorphisms, covering 1/3 of the genome, and found extensive allele frequency changes and little loss of heterozygosities after 100 generations. As modeled with numerical simulations, SNP differentiation was consistent with genetic drift and average fitness effects of 2%, assuming that selection acted independently at each locus. Remarkably, inbreeding by self‐fertilization was of little consequence to SNP differentiation. Modeling selection on deleterious recessive alleles suggests that the initial evolutionary dynamics can be explained by associative overdominance, but not the later stages because much lower heterozygosities would be maintained during experimental evolution. By contrast, models with selection on true overdominant loci can explain the heterozygote excess observed at all periods, particularly when negative epistasis or independent fitness effects were considered. Overall, these findings indicate that selection at single loci, including purging of recessive alleles, underlies most of the genetic differentiation accomplished during the experiment. Nonetheless, they also imply that maintenance of genetic diversity may in large part be due to balancing selection at multiple loci.     

Maintenance of DQB1 polymorphisms in primates.

To understand the evolution of the class II major histocompatibility complex (MHC) DQB1 locus in primates, the second exons of seven DQB1 alleles from five non-human primate species were amplified by polymerase chain reaction. Comparisons of these and other primate sequences show that no between-species diversity is greater than within-species diversity, suggesting maintenance of DQB1 alleles through the history of Old-World primates. There is a preponderance of nonsynonymous nucleotide substitutions at antigen-binding-site codons; this pattern is in marked contrast to what is seen at the closely related, presumably nonfunctional DQB2 gene. The results support the hypothesis that DQB1 polymorphism is maintained by overdominant selection relating to antigen presentation.     

Balancing selection, random genetic drift, and genetic variation at the major histocompatibility complex in two wild populations of guppies (Poecilia reticulata)

Our understanding of the evolution of genes of the major histocompatibility complex (MHC) is rapidly increasing, but there are still enigmatic questions remaining, particularly regarding the maintenance of high levels of MHC polymorphisms in small, isolated populations. Here, we analyze the genetic variation at eight microsatellite loci and sequence variation at exon 2 of the MHC class IIB (DAB) genes in two wild populations of the Trinidadian guppy, Poecilia reticulata. We compare the genetic variation of a small (Ne, 100) and relatively isolated upland population to that of its much larger (Ne approximately 2400) downstream counterpart. As predicted, microsatellite diversity in the upland population is significantly lower and highly differentiated from the population further downstream. Surprisingly, however, these guppy populations are not differentiated by MHC genetic variation and show very similar levels of allelic richness. Computer simulations indicate that the observed level of genetic variation can be maintained with overdominant selection acting at three DAB loci. The selection coefficients differ dramatically between the upland (s > or = 0.2) and lowland (s < or = 0.01) populations. Parasitological analysis on wild-caught fish shows that parasite load is significantly higher on upland than on lowland fish, which suggests that large differences in selection intensity may indeed exist between populations. Based on the infection intensity, a substantial proportion of the upland fish would have suffered direct or indirect fitness consequences as a result of their high parasite loads. Selection by parasites plays a particularly important role in the evolution of guppies in the upland habitat, which has resulted in high levels of MHC diversity being maintained in this population despite considerable genetic drift.     

Diversifying selection underlies the origin of allozyme polymorphism at the phosphoglucose isomerase locus in Tigriopus californicus

The marine copepod Tigriopus californicus lives in intertidal rock pools along the Pacific coast, where it exhibits strong, temporally stable population genetic structure. Previous allozyme surveys have found high frequency private alleles among neighboring subpopulations, indicating that there is limited genetic exchange between populations. Here we evaluate the factors responsible for the diversification and maintenance of alleles at the phosphoglucose isomerase (Pgi) locus by evaluating patterns of nucleotide variation underlying previously identified allozyme polymorphism. Copepods were sampled from eleven sites throughout California and Baja California, revealing deep genetic structure among populations as well as genetic variability within populations. Evidence of recombination is limited to the sample from Pescadero and there is no support for linkage disequilibrium across the Pgi locus. Neutrality tests and codon-based models of substitution suggest the action of natural selection due to elevated non-synonymous substitutions at a small number of sites in Pgi. Two sites are identified as the charge-changing residues underlying allozyme polymorphisms in T. californicus. A reanalysis of allozyme variation at several focal populations, spanning a period of 26 years and over 200 generations, shows that Pgi alleles are maintained without notable frequency changes. Our data suggest that diversifying selection accounted for the origin of Pgi allozymes, while McDonald-Kreitman tests and the temporal stability of private allozyme alleles suggests that balancing selection may be involved in the maintenance of amino acid polymorphisms within populations.     

Evidence for convergent nucleotide evolution and high allelic turnover rates at the complementary sex determiner gene of Western and Asian honeybees

Our understanding of the impact of recombination, mutation, genetic drift, and selection on the evolution of a single gene is still limited. Here we investigate the impact of all these evolutionary forces at the complementary sex determiner (csd) gene that evolves under a balancing mode of selection. Females are heterozygous at the csd gene and males are hemizygous; diploid males are lethal and occur when csd is homozygous. Rare alleles thus have a selective advantage, are seldom lost by the effect of genetic drift, and are maintained over extended periods of time when compared with neutral polymorphisms. Here, we report on the analysis of 17, 19, and 15 csd alleles of Apis cerana, Apis dorsata, and Apis mellifera honeybees, respectively. We observed great heterogeneity of synonymous (piS) and nonsynonymous (piN) polymorphisms across the gene, with a consistent peak in exons 6 and 7. We propose that exons 6 and 7 encode the potential specifying domain (csd-PSD) that has accumulated elevated nucleotide polymorphisms over time by balancing selection. We observed no direct evidence that balancing selection favors the accumulation of nonsynonymous changes at csd-PSD (piN/piS ratios are all <1, ranging from 0.6 to 0.95). We observed an excess of shared nonsynonymous changes, which suggest that strong evolutionary constraints are operating at csd-PSD resulting in the independent accumulation of the same nonsynonymous changes in different alleles across species (convergent evolution). Analysis of csd-PSD genealogy revealed relatively short average coalescence times ( approximately 6 Myr), low average synonymous nucleotide diversity (piS < 0.09), and a lack of trans-specific alleles that substantially contrasts with previously analyzed loci under strong balancing selection. We excluded the possibility of a burst of diversification after population bottlenecking and intragenic recombination as explanatory factors, leaving high turnover rates as the explanation for this observation. By comparing observed allele richness and average coalescence times with a simplified model of csd-coalescence, we found that small long-term population sizes (i.e., N(e) < 10(4)), but not high mutation rates, can explain short maintenance times, implicating a strong historical impact of genetic drift on the molecular evolution of highly social honeybees.     

Positive Darwinian selection promotes charge profile diversity in the antigen-binding cleft of class I major-histocompatibility-complex molecules

Certain major-histocompatibility-complex (MHC) loci are highly polymorphic, and the mechanism of maintenance of this polymorphism remains controversial. Recent studies of the pattern of nucleotide substitution at MHC loci have produced strong evidence that this polymorphism is maintained mainly by positive Darwinian selection that operates on the antigen recognition site (ARS) of the MHC molecule. The ARS of the class I MHC consists of three subregions: (1) the binding cleft, (2) T-cell-receptor-directed residues, and (3) outward-directed residues. Here we report that the rate of nonsynonymous nucleotide substitution is much higher in the binding cleft than in the other ARS subregions. Furthermore, nonsynonymous nucleotide substitutions that result in a change of residue side-chain charge occur significantly more frequently than expected by chance. We conclude that the main target of positive selection on the class I MHC molecules is the binding cleft of the ARS and that this selection acts primarily to promote diversity among alleles with respect to the pattern of residue side-chain charges (charge profile) in the binding cleft. These results provide additional support for the hypothesis that MHC polymorphism is maintained by overdominant selection relating to antigen-binding capacity and thus to disease resistance.     

Overdominance interacts with linkage to determine the rate of adaptation to a new optimum

Overdominance, or a fitness advantage of a heterozygote over both homozygotes, can occur commonly with adaptation to a new optimum phenotype. We model how such overdominant polymorphisms can reduce the evolvability of diploid populations, uncovering a novel form of epistatic constraint on adaptation. The fitness load caused by overdominant polymorphisms can most readily be ameliorated by evolution at tightly linked loci; therefore, traits controlled by multiple loosely linked loci are predicted to be strongly constrained. The degree of constraint is also sensitive to the shape of the relationship between phenotype and fitness, and the constraint caused by overdominance can be strong enough to overcome the effects of clonal interference on the rate of adaptation for a trait. These results point to novel influences on evolvability that are specific to diploids and interact with genetic architecture, and they predict a source of stochastic variability in eukaryotic evolution experiments or cases of rapid evolution in nature.     

The effects of multilocus balancing selection on neutral variability

We studied the effect of multilocus balancing selection on neutral nucleotide variability at linked sites by simulating a model where diallelic polymorphisms are maintained at an arbitrary number of selected loci by means of symmetric overdominance. Different combinations of alleles define different genetic backgrounds that subdivide the population and strongly affect variability. Several multilocus fitness regimes with different degrees of epistasis and gametic disequilibrium are allowed. Analytical results based on a multilocus extension of the structured coalescent predict that the expected linked neutral diversity increases exponentially with the number of selected loci and can become extremely large. Our simulation results show that although variability increases with the number of genetic backgrounds that are maintained in the population, it is reduced by random fluctuations in the frequencies of those backgrounds and does not reach high levels even in very large populations. We also show that previous results on balancing selection in single-locus systems do not extend to the multilocus scenario in a straightforward way. Different patterns of linkage disequilibrium and of the frequency spectrum of neutral mutations are expected under different degrees of epistasis. Interestingly, the power to detect balancing selection using deviations from a neutral distribution of allele frequencies seems to be diminished under the fitness regime that leads to the largest increase of variability over the neutral case. This and other results are discussed in the light of data from the Mhc.     

Contrasting patterns of gene diversity between microsatellites and mitochondrial SNPs in farm and wild Atlantic salmon

Levels of genetic variability at 12 microsatellite loci and 19 single nucleotide polymorphisms in mitochondrial DNA were studied in four farm strains and four wild populations of Atlantic salmon. Within populations, the farm strains showed significantly lower allelic richness and expected heterozygosity than wild populations at the 12 microsatellite loci, but a significantly higher genetic variability with respect to observed number of haplotypes and haplotype diversity in mtDNA. Significant differences in allele- and haplotype-frequencies were observed between farm strains and wild populations, as well as between different farm strains and between different wild populations. The large genetic differentiation at mitochondrial DNA between wild populations (F_ST = 0.24), suggests that the farm strains attained a high mitochondrial genetic variability when created from different wild populations seven generations ago. A large proportion of this variability remains despite an expected lower effective population size for mitochondrial than nuclear DNA. This is best explained by the particular mating schemes in the breeding programmes, with 2–4 females per male. Our observations suggest that for some genetic polymorphisms farm populations might currently hold equal or higher genetic variability than wild populations, but lower overall genetic variability. In the short-term, genetic interactions between escaped farm salmon and wild salmon might increase genetic variability in wild populations, for some, but not most, genetic polymorphisms. In the long term, further losses of genetic variability in farm populations are expected for all genetic polymorphisms, and genetic variability in wild populations will be reduced if escapes of farm salmon continue.     

Adaptation of Drosophila to a novel laboratory environment reveals temporally heterogeneous trajectories of selected alleles

The genomic basis of adaptation to novel environments is a fundamental problem in evolutionary biology that has gained additional importance in the light of the recent global change discussion. Here, we combined laboratory natural selection (experimental evolution) in Drosophila melanogaster with genome‐wide next generation sequencing of DNA pools (Pool‐Seq) to identify alleles that are favourable in a novel laboratory environment and traced their trajectories during the adaptive process. Already after 15 generations, we identified a pronounced genomic response to selection, with almost 5000 single nucleotide polymorphisms (SNP; genome‐wide false discovery rates < 0.005%) deviating from neutral expectation. Importantly, the evolutionary trajectories of the selected alleles were heterogeneous, with the alleles falling into two distinct classes: (i) alleles that continuously rise in frequency; and (ii) alleles that at first increase rapidly but whose frequencies then reach a plateau. Our data thus suggest that the genomic response to selection can involve a large number of selected SNPs that show unexpectedly complex evolutionary trajectories, possibly due to nonadditive effects.     

Genome‐wide signature of local adaptation linked to variable CpG methylation in oak populations

It has long been known that adaptive evolution can occur through genetic mutations in DNA sequence, but it is unclear whether adaptive evolution can occur through analogous epigenetic mechanisms, such as through DNA methylation. If epigenetic variation contributes directly to evolution, species under threat of disease, invasive competition, climate change or other stresses would have greater stores of variation from which to draw. We looked for evidence of natural selection acting on variably methylated DNA sites using population genomic analysis across three climatologically distinct populations of valley oaks. We found patterns of genetic and epigenetic differentiations that indicate local adaptation is operating on large portions of the oak genome. While CHG methyl polymorphisms are not playing a significant role and would make poor targets for natural selection, our findings suggest that CpG methyl polymorphisms as a whole are involved in local adaptation, either directly or through linkage to regions under selection.     

ANTAGONISTIC VERSUS NONANTAGONISTIC MODELS OF BALANCING SELECTION: CHARACTERIZING THE RELATIVE TIMESCALES AND HITCHHIKING EFFECTS OF PARTIAL SELECTIVE SWEEPS

Antagonistically selected alleles‐–those with opposing fitness effects between sexes, environments, or fitness components‐–represent an important component of additive genetic variance in fitness‐related traits, with stably balanced polymorphisms often hypothesized to contribute to observed quantitative genetic variation. Balancing selection hypotheses imply that intermediate‐frequency alleles disproportionately contribute to genetic variance of life‐history traits and fitness. Such alleles may also associate with population genetic footprints of recent selection, including reduced genetic diversity and inflated linkage disequilibrium at linked, neutral sites. Here, we compare the evolutionary dynamics of different balancing selection models, and characterize the evolutionary timescale and hitchhiking effects of partial selective sweeps generated under antagonistic versus nonantagonistic (e.g., overdominant and frequency‐dependent selection) processes. We show that the evolutionary timescales of partial sweeps tend to be much longer, and hitchhiking effects are drastically weaker, under scenarios of antagonistic selection. These results predict an interesting mismatch between molecular population genetic and quantitative genetic patterns of variation. Balanced, antagonistically selected alleles are expected to contribute more to additive genetic variance for fitness than alleles maintained by classic, nonantagonistic mechanisms. Nevertheless, classical mechanisms of balancing selection are much more likely to generate strong population genetic signatures of recent balancing selection.     

The genetic basis and experimental evolution of inbreeding depression in Caenorhabditis elegans

Determining the genetic basis of inbreeding depression is important for understanding the role of selection in the evolution of mixed breeding systems. Here, we investigate how androdioecy (a breeding system characterized by partial selfing and outcrossing) and dioecy (characterized by obligatory outcrossing) influence the experimental evolution of inbreeding depression in Caenorhabditis elegans. We derived inbred lines from ancestral and evolved populations and found that the dioecious lineages underwent more extinction than androdioecious lineages. For both breeding systems, however, there was selection during inbreeding because the diversity patterns of 337 single-nucleotide polymorphisms (SNPs) among surviving inbred lines deviated from neutral expectations. In parallel, we also followed the evolution of embryo to adult viability, which revealed similar starting levels of inbreeding depression in both breeding systems, but also outbreeding depression. Under androdioecy, diversity at a neutral subset of 134 SNPs correlated well with the viability trajectories, showing that the population genetic structure imposed by partial selfing affected the opportunity for different forms of selection. Our findings suggest that the interplay between the disruptions of coevolved sets of loci by outcrossing, the efficient purging of deleterious recessive alleles with selfing and overdominant selection with outcrossing can help explain mixed breeding systems.     

The generation and maintenance of genetic variation by frequency-dependent selection: constructing polymorphisms under the pairwise interaction model

Frequency-dependent selection remains the most commonly invoked heuristic explanation for the maintenance of genetic variation. For polymorphism to exist, new alleles must be both generated and maintained in the population. Here we use a construction approach to model frequency-dependent selection with mutation under the pairwise interaction model. The pairwise interaction model is a general model of frequency-dependent selection at the genotypic level. We find that frequency-dependent selection is able to generate a large number of alleles at a single locus. The construction process generates multiallelic polymorphisms with a wide range of allele-frequency distributions and genotypic fitness relationships. Levels of polymorphism and mean fitness are uncoupled, so constructed polymorphisms remain permanently invasible to new mutants; thus the model never settles down to an equilibrium state. Analysis of constructed fitness sets reveals signatures of heterozygote advantage and positive frequency dependence.     

Frequency-dependent selection and the maintenance of genetic variation: exploring the parameter space of the multiallelic pairwise interaction model

When individuals' fitnesses depend on the genetic composition of the population in which they are found, selection is then frequency dependent. Frequency-dependent selection (FDS) is often invoked as a heuristic explanation for the maintenance of large numbers of alleles at a locus. The pairwise interaction model is a general model of FDS via intraspecific competition at the genotypic level. Here we use a parameter-space approach to investigate the full potential for the maintenance of multiallelic equilibria under the pairwise interaction model. We find that FDS maintains full polymorphism more often than classic constant-selection models and produces more skewed equilibrium allele frequencies. Fitness sets with some degree of rare advantage maintained full polymorphism most often, but a wide variety of nonobvious fitness patterns were also found to have positive potential for polymorphism. An example is put forth suggesting possible explanations for multiallelic polymorphisms maintained despite positive FDS on individual alleles.     

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

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

Sojourn times and substitution rate at overdominant and linked neutral loci

The sojourn times until fixation of an overdominant allele were investigated based on the diffusion equation. Furthermore, the rate of accumulation of mutations, or the substitution rate, was predicted from the mean extinction time of a common overdominant allele. The substitution rate calculated theoretically agreed well with that determined by computer simulation. Overdominant selection enhances the polymorphism at linked loci, while its effect on the sojourn times and the substitution rate at linked loci has not been studied yet. To solve these problems, a model that assumed two linked loci, each with infinite alleles, was examined by computer simulation. A decrease in the recombination rate between two loci markedly changed the distribution of sojourn times of a neutral allele. Although overdominant selection obviously increased the sojourn times and the polymorphism at a linked locus, the rate of nucleotide substitution at the neutral locus was not influenced significantly even if complete linkage was assumed. These results suggest that, in regions containing overdominant genes, linked neutral loci will exhibit elevated levels of polymorphism, but their rate of molecular evolution remains that predicted by neutral theory.     

HLA DNA sequence variation among human populations: molecular signatures of demographic and selective events

Molecular differences between HLA alleles vary up to 57 nucleotides within the peptide binding coding region of human Major Histocompatibility Complex (MHC) genes, but it is still unclear whether this variation results from a stochastic process or from selective constraints related to functional differences among HLA molecules. Although HLA alleles are generally treated as equidistant molecular units in population genetic studies, DNA sequence diversity among populations is also crucial to interpret the observed HLA polymorphism. In this study, we used a large dataset of 2,062 DNA sequences defined for the different HLA alleles to analyze nucleotide diversity of seven HLA genes in 23,500 individuals of about 200 populations spread worldwide. We first analyzed the HLA molecular structure and diversity of these populations in relation to geographic variation and we further investigated possible departures from selective neutrality through Tajima's tests and mismatch distributions. All results were compared to those obtained by classical approaches applied to HLA allele frequencies.Our study shows that the global patterns of HLA nucleotide diversity among populations are significantly correlated to geography, although in some specific cases the molecular information reveals unexpected genetic relationships. At all loci except HLA-DPB1, populations have accumulated a high proportion of very divergent alleles, suggesting an advantage of heterozygotes expressing molecularly distant HLA molecules (asymmetric overdominant selection model). However, both different intensities of selection and unequal levels of gene conversion may explain the heterogeneous mismatch distributions observed among the loci. Also, distinctive patterns of sequence divergence observed at the HLA-DPB1 locus suggest current neutrality but old selective pressures on this gene. We conclude that HLA DNA sequences advantageously complement HLA allele frequencies as a source of data used to explore the genetic history of human populations, and that their analysis allows a more thorough investigation of human MHC molecular evolution.     

QUANTITATIVE GENETIC VARIANCE MAINTAINED BY FLUCTUATING SELECTION WITH OVERLAPPING GENERATIONS: VARIANCE COMPONENTS AND COVARIANCES

The quantitative genetic variance-covariance that can be maintained in a random environment is studied, assuming overlapping generations and Gaussian stabilizing selection with a fluctuating optimum. The phenotype of an individual is assumed to be determined by additive contributions from each locus on paternal and maternal gametes (i.e., no epistasis and no dominance). Recurrent mutation is ignored, but linkage between loci is arbitrary. The genotype distribution in the evolutionarily stable population is generically discrete: only a finite number of polymorphic alleles with distinctly different effects are maintained, even though we allow a continuum of alleles with arbitrary phenotypic contributions to invade. Fluctuating selection maintains nonzero genetic variance in the evolutionarily stable population if the environmental heterogeneity is larger than a certain threshold. Explicit asymptotic expressions for the standing variance-covariance components are derived for the population near the threshold, or for large generational overlap, as a function of environmental variability and genetic parameters (i.e., number of loci, recombination rate, etc.), using the fact that the genotype distribution is discrete. Above the threshold, the population maintains considerable genetic variance in the form of positive linkage disequilibrium and positive gamete covariance (Hardy-Weinberg disequilibrium) as well as allelic variance. The relative proportion of these disequilibrium variances in the total genetic variance increases with the environmental variability.     

Unusually high recombination rate detected in the sex locus region of the honey bee (Apis mellifera)

Sex determination in Hymenoptera is controlled by haplo-diploidy in which unfertilized eggs develop into fertile haploid males. A single sex determination locus with several complementary alleles was proposed for Hymenoptera [so-called complementary sex determination (CSD)]. Heterozygotes at the sex determination locus are normal, fertile females, whereas diploid zygotes that are homozygous develop into sterile males. This results in a strong heterozygote advantage, and the sex locus exhibits extreme polymorphism maintained by overdominant selection. We characterized the sex-determining region by genetic linkage and physical mapping analyses. Detailed linkage and physical mapping studies showed that the recombination rate is <44 kb/cM in the sex-determining region. Comparing genetic map distance along the linkage group III in three crosses revealed a large marker gap in the sex-determining region, suggesting that the recombination rate is high. We suggest that a "hotspot" for recombination has resulted here because of selection for combining favorable genotypes, and perhaps as a result of selection against deleterious mutations. The mapping data, based on long-range restriction mapping, suggest that the Q DNA-marker is within 20,000 bp of the sex locus, which should accelerate molecular analyses.