Balancing selection maintains intraspecific diversity in a deep-sea fish.

Genes involved in immune defense against pathogens provide some of the most well-known examples of both directional and balancing selection. Antimicrobial peptides (AMPs) are innate immune effector genes, playing a key role in pathogen clearance in many species, including Drosophila. Conflicting lines of evidence have suggested that AMPs may be under directional, balancing, or purifying selection. Here, we use both a linear model and control-gene-based approach to show that balancing selection is an important force shaping AMP diversity in Drosophila. In Drosophila melanogaster, this is most clearly observed in ancestral African populations. Furthermore, the signature of balancing selection is even more striking once background selection has been accounted for. Balancing selection also acts on AMPs in Drosophila mauritiana, an isolated island endemic separated from D. melanogaster by about 4 Myr of evolution. This suggests that balancing selection may be broadly acting to maintain adaptive diversity in Drosophila AMPs, as has been found in other taxa.

Balancing selection has been shown to be common in plants for several different types of traits, such as self-incompatibility and heterostyly. Generally, for these traits balancing selection is generated by interactions among individuals or between individuals and other species (e.g., pathogens or pollinators). However, there are phenotypic polymorphisms in plants that do not obviously involve types of interactions that generate balancing selection. Little is known about the extent to which balancing selection also acts to preserve these polymorphisms. Here we ask whether balancing selection preserves an anther-color polymorphism in Erythronium umbilicatum (Liliaceae). We identified a major gene underlying this polymorphism. We then attempted to detect signatures of balancing selection on that gene by developing a new coalescence test for balancing selection. We found that variation in anther color is in large part caused by variation in a paralog of EuMYB3, an anthocyanin-regulating R2R3-MYB transcription factor. However, we found little evidence for balancing selection having acted historically on EuMYB3. Our results thus suggest that plant polymorphisms, especially those not involved in interactions that are likely to generate negative frequency-dependent selection, may reflect a transient state in which one morph will eventually be fixed by either genetic drift or directional selection. Our results also suggest that regulation of the anthocyanin pathway is more evolutionarily labile than is generally believed.

Understanding the processes that maintain genetic variation in natural populations is a major focus of evolutionary biology. Theoretical investigations have shown that more than a dozen processes can explain the maintenance of genetic variation (see Hartl 1980 for a summary), whereas experiments with laboratory populations indicate that at least some of these processes can operate under some conditions (e.g., Dobzhansky 1948; Dobzhansky and Pavlovsky 1953; Ehrman 1967; Powell and Wistrand 1978; Jones and Probert 1980; Cavender and Clegg 1981). Nevertheless, detailed studies of the selective forces acting to maintain genetic variation in natural populations are rare. The study reported here is part of an ongoing examination of the selective forces acting to maintain genetic variation at a locus influencing floral pigment intensity in the common morning glory, Ipomoea purpurea. This locus (designated W) is polymorphic in natural populations throughout the southeastern United States. Plants homozygous for the w allele have white flowers with pigmented rays (i.e., whites). Heterozygous plants have lightly pigmented flowers with dark rays (i.e., lights), and homozygotes for the W allele have darkly pigmented flowers with even darker rays (i.e., darks) (Ennos and Clegg 1983; Epperson and Clegg 1988). In natural populations, frequencies of the white allele generally range from 0 to 0.4, with a mean of about 0.1 (Epperson and Clegg 1986). Previous work has shown that when white (ww) plants are in the minority, they are undervisited by bees (Brown and Clegg 1984; Epperson and Clegg 1986; Rausher et al. 1993). This undervisitation results in higher selfing rates for white plants, without an associated reduction in the contribution of whites to the outcross pollen pool (i.e., no pollen discounting, Rausher et al. 1993; Iwao 1995). Moreover, inbreeding depression is minimal or absent in North American populations that have been examined (Pear 1983; S. M. Chang, pers. comm. 1996). These observations indicate that the white allele should enjoy a transmission advantage, and thus increase in frequency, when rare (Fisher 1941) and may explain the protection of the white allele in natural populations. It is unclear, however, what forces operate to protect the dark allele. In the absence of any other selective forces acting on the W locus, a combination of genetic drift and active protection of the white allele is expected to result in the fixation of the white allele. The high frequency of the dark allele (0.6-1) in most natural populations therefore suggests the action of other selective forces favoring the dark allele.

Editorial and retrospective 2010

Sexual antagonism (SA) occurs when an allele that is beneficial to one sex, is detrimental to the other. This conflict can result in balancing, directional, or disruptive selection acting on SA alleles. A body of theory predicts the conditions under which sexually antagonistic mutants will invade and be maintained in stable polymorphism under balancing selection. There remains, however, considerable debate over the distribution of SA genetic variation across autosomes and sex chromosomes, with contradictory evidence coming from data and theory. In this article, we investigate how the interplay between selection and genetic drift will affect the genomic distribution of sexually antagonistic alleles. The effective population sizes can differ between the autosomes and the sex chromosomes due to a number of ecological factors and, consequently, the distribution of SA genetic variation in genomes. In general, we predict the interplay of SA selection and genetic drift should lead to the accumulation of SA alleles on the X in male heterogametic (XY) species and, on the autosomes in female heterogametic (ZW) species, especially when sexual competition is strong among males.

The role of balancing selection is a long-standing evolutionary puzzle. Balancing selection is a crucial evolutionary process that maintains genetic variation (polymorphism) over extended periods of time; however, detecting it poses a significant challenge. Building upon the Polymorphism-aware phylogenetic Models (PoMos) framework rooted in the Moran model, we introduce a PoMoBalance model. This novel approach is designed to disentangle the interplay of mutation, genetic drift, and directional selection (GC-biased gene conversion), along with the previously unexplored balancing selection pressures on ultra-long timescales comparable with species divergence times by analyzing multi-individual genomic and phylogenetic divergence data. Implemented in the open-source RevBayes Bayesian framework, PoMoBalance offers a versatile tool for inferring phylogenetic trees as well as quantifying various selective pressures. The novel aspect of our approach in studying balancing selection lies in polymorphism-aware phylogenetic models' ability to account for ancestral polymorphisms and incorporate parameters that measure frequency-dependent selection, allowing us to determine the strength of the effect and exact frequencies under selection. We implemented validation tests and assessed the model on the data simulated with SLiM and a custom Moran model simulator. Real sequence analysis of Drosophila populations reveals insights into the evolutionary dynamics of regions subject to frequency-dependent balancing selection, particularly in the context of sex-limited color dimorphism in Drosophila erecta.

Genetic linkage and molecular evolution

Strong balancing selection on the major histocompatibility complex (MHC) can lead to different patterns in gene frequencies and neutral genomic variation within species. We investigated diversity and geographic structure of MHC genes DQA and DQB, as well as their inferred functional haplotypes, from 2 regional populations (East and West Coast) of the endangered Hector's dolphin (Cephalorhynchus hectori hectori) and the critically endangered Māui dolphin (Cephalorhynchus hectori maui) (West Coast, North Island), and contrasted these results with patterns from neutral microsatellites. The Māui had the lowest number of alleles for DQA (2) and DQB (3), consistent with strong genetic drift acting on this remnant population. However, the 2 retained DQA alleles are among the most divergent combinations of all 4 alleles found across the Hector's metapopulation, potentially reflecting the retention of divergent alleles due to balancing selection. The high frequency of the divergent DQB*04 allele also gave this population the highest nucleotide diversity for DQB. Strong differentiation was evident for DQA, DQB, and DQA-DQB haplotypes between the regional populations of Hector's dolphins (FST > 0.213) and both subspecies (FST > 0.311). Differentiation was generally greater than observed at neutral microsatellite loci, suggesting the influence of selection between geographically proximate East and West Coast populations. This might be the result of spatial differences in directional selection on those opposite coastlines. In addition, measures of the ratio of nonsynonymous to synonymous substitutions (dN/dS) were consistent with balancing selection over evolutionary time. Together, these results suggest a complex interplay of balancing selection, directional selection, local fidelity, and genetic drift.

Natural selection is expected to leave an imprint on the neutral polymorphisms at the adjacent genomic regions of a selected gene. While directional selection tends to reduce within-population genetic diversity and increase among-population differentiation, the reverse is expected under balancing selection. To identify targets of natural selection in the three-spined stickleback (Gasterosteus aculeatus) genome, 103 microsatellite and two indel markers including expressed sequence tags (EST) and quantitative trait loci (QTL)-associated loci, were genotyped in four freshwater and three marine populations. The results indicated that a high proportion of loci (14.7%) might be affected by balancing selection and a lower proportion (2.8%) by directional selection. The strongest signatures of directional selection were detected in a microsatellite locus and two indel markers located in the intronic regions of the Eda-gene coding for the number of lateral plates. Yet, other microsatellite loci previously found to be informative in QTL-mapping studies revealed no signatures of selection. Two novel microsatellite loci (Stn12 and Stn90) located in chromosomes I and VIII, respectively, showed signals of directional selection and might be linked to genomic regions containing gene(s) important for adaptive divergence. Although the coverage of the total genomic content was relatively low, the predominance of balancing selection signals is in agreement with the contention that balancing, rather than directional selection is the predominant mode of selection in the wild.

JOHNSON1 has suggested that the relationship observed between the actual and effective numbers of alleles in certain natural populations is not consistent with the selectively neutral allele model of Kimura and Crow2. We show here, however, that the data from a large number of natural populations are in fact compatible with this model.

Tandemly arranged paralogous genes lbe and lbl are members of the Drosophila NK homeobox family. We analyzed population samples of Drosophila melanogaster from Africa, Europe, North and South America, and single strains of D. sechellia, D. simulans, and D. yakuba within two linked regions encompassing partial sequences of lbe and lbl. The evolution of lbe and lbl is highly constrained due to their important regulatory functions. Despite this, a variety of forces have shaped the patterns of variation in lb genes: recombination, intragenic gene conversion and natural selection strongly influence background variation created by linkage disequilibrium and dimorphic haplotype structure. The two genes exhibited similar levels of nucleotide diversity and positive selection was detected in the noncoding regions of both genes. However, synonymous variability was significantly higher for lbe: no nonsynonymous changes were observed in this gene. We argue that balancing selection impacts some synonymous sites of the lbe gene. Stability of mRNA secondary structure was significantly different between the lbe (but not lbl) haplotype groups and may represent a driving force of balancing selection in epistatically interacting synonymous sites. Balancing selection on synonymous sites may be the first, or one of a few such observations, in Drosophila. In contrast, recurrent positive selection on lbl at the protein level influenced evolution at three codon sites. Transcription factor binding-site profiles were different for lbe and lbl, suggesting that their developmental functions are not redundant. Combined with our previous results on nucleotide variation in esterase and other homeobox genes, these results suggest that interplay of balancing and directional selection may be a general feature of molecular evolution in Drosophila and other eukaryote genomes.

Natural selection not only affects the actual loci under selection but also leaves "footprints" in patterns of genetic variation in linked genetic regions. This offers exciting opportunities for inferring selection and for understanding the processes shaping levels of genetic variation in natural populations. Here, we develop analytical approximations based on coalescent theory to characterize the genetic footprint of a complex, but potentially common type of natural selection: balancing selection with seasonally fluctuating allele frequencies. As we show analytically and confirm with stochastic simulations, seasonal allele frequency fluctuations can have important (and partly unexpected) consequences for the genetic footprint of balancing selection. Fluctuating balancing selection generally leads to an increase in genetic diversity close to the selected site, the effect of balancing selection, but reduces diversity further away from the selected site, which is a consequence of the allele-frequency fluctuations effectively producing recurrent bottlenecks of allelic backgrounds. This medium- and long-range reduction usually outweighs the short-range increase when averaging diversity levels across the entire chromosome. Strong fluctuating balancing selection even induces a loss of genetic variation in unlinked regions, e.g. on different chromosomes. If many loci in the genome are simultaneously under fluctuating balancing selection this can lead to substantial genome-wide reductions in genetic diversity, even when allele-frequency fluctuations are small and local footprints are difficult to detect. Thus, together with genetic drift, selective sweeps and background selection, fluctuating selection could be a major force shaping levels of genetic diversity in natural populations.

Antagonistic selection--where alleles at a locus have opposing effects on male and female fitness ("sexual antagonism") or between components of fitness ("antagonistic pleiotropy")--might play an important role in maintaining population genetic variation and in driving phylogenetic and genomic patterns of sexual dimorphism and life-history evolution. While prior theory has thoroughly characterized the conditions necessary for antagonistic balancing selection to operate, we currently know little about the evolutionary interactions between antagonistic selection, recurrent mutation, and genetic drift, which should collectively shape empirical patterns of genetic variation. To fill this void, we developed and analyzed a series of population genetic models that simultaneously incorporate these processes. Our models identify two general properties of antagonistically selected loci. First, antagonistic selection inflates heterozygosity and fitness variance across a broad parameter range--a result that applies to alleles maintained by balancing selection and by recurrent mutation. Second, effective population size and genetic drift profoundly affect the statistical frequency distributions of antagonistically selected alleles. The "efficacy" of antagonistic selection (i.e., its tendency to dominate over genetic drift) is extremely weak relative to classical models, such as directional selection and overdominance. Alleles meeting traditional criteria for strong selection (N(e)s >> 1, where N(e) is the effective population size, and s is a selection coefficient for a given sex or fitness component) may nevertheless evolve as if neutral. The effects of mutation and demography may generate population differences in overall levels of antagonistic fitness variation, as well as molecular population genetic signatures of balancing selection.

How phenotypic diversity originates and persists within populations are classic puzzles in evolutionary biology. While balanced polymorphisms segregate within many species, it remains rare for both the genetic basis and the selective forces to be known, leading to an incomplete understanding of many classes of traits under balancing selection. Here, we uncover the genetic architecture of a balanced sexual mimicry polymorphism and identify behavioral mechanisms that may be involved in its maintenance in the swordtail fish Xiphophorus birchmanni. We find that ∼40% of X.birchmanni males develop a "false gravid spot," a melanic pigmentation pattern that mimics the "pregnancy spot" associated with sexual maturity in female live-bearing fish. Using genome-wide association mapping, we detect a single intergenic region associated with variation in the false gravid spot phenotype, which is upstream of kitlga, a melanophore patterning gene. By performing long-read sequencing within and across populations, we identify complex structural rearrangements between alternate alleles at this locus. The false gravid spot haplotype drives increased allele-specific expression of kitlga, which provides a mechanistic explanation for the increased melanophore abundance that causes the spot. By studying social interactions in the laboratory and in nature, we find that males with the false gravid spot experience less aggression; however, they also receive increased attention from other males and are disdained by females. These behavioral interactions may contribute to the maintenance of this phenotypic polymorphism in natural populations. We speculate that structural variants affecting gene regulation may be an underappreciated driver of balanced polymorphisms across diverse species.

Animals display incredibly diverse colour patterns, a testament to evolution's endless innovation in shaping life. In many species, the interplay between males and females in the pursuit of mates has driven the evolution of a myriad of colour forms, from the flashy peacock tail feathers to the tiniest colour markings in damselflies. In others, colour provides crypsis by allowing to blend into the background and to escape the eyes of predators. While the obvious benefits of this dazzling diversity for reproduction and survival seem straightforward, its maintenance is not. Theory predicts that genetic drift and various forms of selection reduce variation over time, making the persistence of colour variants over generations a puzzle. In this issue of Molecular Ecology, Lindtke etal. () study the cryptic colour morphs of Timema cristinae walking sticks to shed light on the genetic architecture and mechanisms that allow colour polymorphism maintenance over long timescales. By combining genome-wide data with phenotyping information from natural populations, they were able to map the green and melanistic colour to one genomic region with highly reduced effective recombination rate between two main chromosomal variants, consistent with an inversion polymorphism. These two main chromosomal variants showed geographically widespread heterozygote excess, and genomic signatures consistent with long-term balancing selection. A younger chromosomal variant was detected for the third morph, the green-striped colour morphs, in the same genomic regions as the melanistic and the green-unstriped morphs. Together, these results suggest that the genetic architecture of cryptic T.cristinae morphs is caused by nonrecombining genomic blocks that have been maintained over extended time periods by balancing selection making this study one of the few available empirical examples documenting that balancing selection of various forms may play an important role in maintaining adaptive genetic variation in nature.