A sex-linked supergene with large effects on sperm traits has little impact on reproductive traits in female zebra finches.

Despite constituting an essential component of fitness, reproductive success can vary remarkably between individuals and the causes of such variation are not well understood across taxa. In the zebra finch-a model songbird, almost all the variation in sperm morphology and swimming speed is maintained by a large polymorphic inversion (commonly known as a supergene) on the Z chromosome. The relationship between this polymorphism and reproductive success is not fully understood, particularly for females. Here, we explore the effects of female haplotype, and the combination of male and female genotype, on several primary reproductive traits in a captive population of zebra finches. Despite the inversion polymorphism's known effects on sperm traits, we find no evidence that inversion haplotype influences egg production by females or survival of embryos through to hatching. However, our findings do reinforce existing evidence that the inversion polymorphism is maintained by a heterozygote advantage for male fitness. This work provides an important step in understanding the causes of variation in reproductive success in this model species.

Advances in the Evolutionary Understanding of MHC Polymorphism.

Proteins encoded by the classical major histocompatibility complex (MHC) genes incite the vertebrate adaptive immune response by presenting peptide antigens on the cell surface. Here, we review mechanisms explaining landmark features of these genes: extreme polymorphism, excess of nonsynonymous changes in peptide-binding domains, and long gene genealogies. Recent studies provide evidence that these features may arise due to pathogens evolving ways to evade immune response guided by the locally common MHC alleles. However, complexities of selection on MHC genes are simultaneously being revealed that need to be incorporated into existing theory. These include pathogen-driven selection for antigen-binding breadth and expansion of the MHC gene family, associated autoimmunity trade-offs, hitchhiking of deleterious mutations linked to the MHC, geographic subdivision, and adaptive introgression.

High MHC diversity confers no advantage for phenotypic quality and reproductive performance in a wild bird.

Genes of the major histocompatibility complex (MHC) encode antigen-binding molecules and are an integral part of the acquired immune response of vertebrates. In general, high individual MHC diversity is expected to increase fitness by broadening the spectrum of pathogens recognized by the immune system, in accordance with the heterozygote advantage mechanism. On the other hand, the optimality hypothesis assumes that individuals with optimal (intermediate), rather than maximum, diversity of the MHC will achieve the highest fitness because of inherent costs associated with expressing diverse MHC alleles. Here, we tested for associations between individual diversity of the MHC class I and class II genes (binding antigens of intra- and extracellular pathogens respectively) and a range of fitness-related traits (condition, ornament expression and reproduction) in an urban population of the Eurasian coot Fulica atra. Contrary to our expectation, we found that high within-individual allelic diversity of MHC genes (both class I and II) was associated with poorer condition (lower blood haemoglobin concentrations), weaker expression of the putative ornament (smaller frontal shield), later onset of breeding and smaller clutches. An analysis of functional MHC allele clusters (supertypes) provided further support for negative associations of MHC diversity with phenotypic quality and reproductive performance, but most of these relationships could not be explained by the presence of specific maladaptive supertypes. Finally, we found little empirical support for the optimality hypothesis in the Eurasian coot. Our results suggest that the costs of high MHC diversity outweighed any benefits associated with broad MHC repertoire, which could be driven by depauperate pathogen diversity in an urban landscape. To the best of our knowledge, this is one of the first studies providing consistent evidence for negative associations of MHC diversity with a range of fitness-related traits in a natural avian population.

Why Do Sex Chromosomes Stop Recombining?

It is commonly assumed that sex chromosomes evolve recombination suppression because selection favours linkage between sex-determining and sexually antagonistic genes. However, although the role of sexual antagonism during sex chromosome evolution has attained strong support from theory, experimental and observational evidence is rare or equivocal. Here, we highlight alternative, often neglected, hypotheses for recombination suppression on sex chromosomes, which invoke meiotic drive, heterozygote advantage, and genetic drift, respectively. We contrast the hypotheses, the situations when they are likely to be of importance, and outline why it is surprisingly difficult to test them. Lastly, we discuss future research directions (including modelling, population genomics, comparative approaches, and experiments) to disentangle the different hypotheses of sex chromosome evolution.

Models for Eco-Evolutionary Extinction Vortices under Balancing Selection.

AbstractThe smaller a population is, the faster it loses genetic diversity as a result of genetic drift. Loss of genetic diversity can reduce population growth rate, making populations even smaller and more vulnerable to loss of genetic diversity. Ultimately, the population can be driven to extinction by this "eco-evolutionary extinction vortex." While there are already quantitative models for extinction vortices resulting from inbreeding depression and mutation accumulation, to date extinction vortices resulting from loss of genetic diversity at loci under various forms of balancing selection have been mainly described verbally. To understand better when such extinction vortices arise and to develop methods for detecting them, we propose quantitative eco-evolutionary models, both stochastic individual-based simulations and deterministic approximations, linking loss of genetic diversity and population decline. Using mathematical analysis and simulations, we identify parameter combinations that exhibit strong interactions between population size and genetic diversity and match our definition of an eco-evolutionary vortex (i.e., per capita population decline rates and per-locus fixation rates increase with decreasing population size and number of polymorphic loci). We further highlight cues that may be exhibited by such populations but find that classical early-warning signals are of limited use in detecting populations undergoing an eco-evolutionary extinction vortex.

The role of balancing selection in maintaining human RhD blood group polymorphism: A preregistered cross-sectional study.

Maintenance of genetic polymorphism remains one of the big questions of evolutionary biology, which for a long time tended to be explained by balancing selection. This explanation was later criticized, but now is again accepted as an important mechanism in evolution. Human blood group systems seem affected by balancing selection especially strongly. In this preregistered study, we focused on stable coexistence of RhD-positive and RhD-negative subjects in a population. This is an evolutionary conundrum, because carriers of the less frequent negative allele suffer from lower fecundity due to haemolytic disease of the newborn affecting RhD-positive infants born to RhD-negative women. One explanation of persisting stability of RhD polymorphism points to heterozygote advantage. Over the past decade, numerous studies demonstrated that RhD-positive subjects score better than RhD-negative homozygotes in psychomotor tests and physical health-related variables. Still, evidence of better health and performance of heterozygotes is scarce and merely indirect. We compared the physical and mental health of 2,539 subjects whose RhD genotype was estimated based on their and their parents' RhD phenotype. We confirmed that RhD-negative homozygotes fare worse in terms of physical and mental health than subjects with RhD-positive phenotype and that RhD-positive heterozygotes enjoy better health than both homozygotes. For the first time, we demonstrated that RhD-positive homozygotes might suffer from worse health than RhD-negative homozygotes. Our results strongly support the hypothesis that RhD polymorphism is maintained by heterozygote advantage and that balancing selection may have played an important role in human evolution in this context and in general.

Imbalance of Genes Encoding Natural Killer Immunoglobulin-Like Receptors and Human Leukocyte Antigen in Patients With Biliary Cancer.

BACKGROUND & AIMS: Bile duct tumors are rare and have poor prognoses. Natural killer (NK) cells are frequent in human liver and infiltrate these tumors but do not control their progression. Responses of NK cells are regulated by NK immunoglobulin-like receptors (KIRs), which interact with HLA class I ligands. We aimed to characterize the features of the KIR gene loci and their ligands in patients with bile duct cancer (BDC). METHODS: We performed combined multidimensional characterization of genes that encode KIRs and their ligands in blood samples from patients with BDC from Sweden, followed for up to 8 years after diagnosis (n = 148), in 2 geographically matched cohorts of healthy individuals from Northern Europe (n = 204 and n = 900), and in healthy individuals from 6 geographically unrelated populations (n = 2917). We used real-time polymerase chain reaction, RNA sequencing, immunohistochemistry, and flow cytometry to evaluate NK-cell presence, as well as KIR and KIR-ligand expression in bile duct tumors and control tissues. RESULTS: Patients with bile duct tumors had multiple alterations at the KIR gene loci. KIR loci are grouped into genotypes that encode more inhibitory (group A) and more activating (group B) receptors, which can be subdivided into centromeric and telomeric fragments. Patients with BDC had a lower prevalence of KIR2DL3, which was linked to disequilibrium in centromeric A/B and B/B genotypes, compared with control individuals. The associations between KIRs and KIR ligands differed between patients with BDC and control individuals; patients had an altered balance between activating and inhibitory KIRs. KIR-positive NK cells infiltrated biliary tumors that expressed matched KIR ligands. CONCLUSIONS: In a multidimensional analysis of DNA from blood samples of patients with BDC in Europe, we found patients to have multiple alterations at the KIR and HLA gene loci compared with control individuals. These alterations might affect NK-cell tumor surveillance. NK cells from bile duct tumors expressed KIRs and were found in tumors that expressed cognate ligands. This should be considered in development of immune-based therapies for BDC.

Sex and hatching order modulate the association between MHC-II diversity and fitness in early-life stages of a wild seabird.

Genes of the major histocompatibility complex (MHC) play a pivotal role in parasite resistance, and their allelic diversity has been associated with fitness variations in several taxa. However, studies report inconsistencies in the direction of this association, with either positive, quadratic or no association being described. These discrepancies may arise because the fitness costs and benefits of MHC diversity differ among individuals depending on their exposure and immune responses to parasites. Here, we investigated in black-legged kittiwake (Rissa tridactyla) chicks whether associations between MHC class-II diversity and fitness vary with sex and hatching order. MHC-II diversity was positively associated with growth and tick clearance in female chicks, but not in male chicks. Our data also revealed a positive association between MHC-II diversity and survival in second-hatched female chicks (two eggs being the typical clutch size). These findings may result from condition-dependent parasite infections differentially impacting sexes in relation to hatching order. We thus suggest that it may be important to account for individual heterogeneities in traits that potentially exert selective pressures on MHC diversity in order to properly predict MHC-fitness associations.

Genetic and environmental determinants of human TCR repertoire diversity.

T cell discrimination of self and non-self is the foundation of the adaptive immune response, and is orchestrated by the interaction between T cell receptors (TCRs) and their cognate ligands presented by major histocompatibility (MHC) molecules. However, the impact of host immunogenetic variation on the diversity of the TCR repertoire remains unclear. Here, we analyzed a cohort of 666 individuals with TCR repertoire sequencing. We show that TCR repertoire diversity is positively associated with polymorphism at the human leukocyte antigen class I (HLA-I) loci, and diminishes with age and cytomegalovirus (CMV) infection. Moreover, our analysis revealed that HLA-I polymorphism and age independently shape the repertoire in healthy individuals. Our data elucidate key determinants of human TCR repertoire diversity, and suggest a mechanism underlying the evolutionary fitness advantage of HLA-I heterozygosity.

Divergent Allele Advantage at Human MHC Genes: Signatures of Past and Ongoing Selection.

The highly polymorphic genes of the major histocompatibility complex (MHC) play a key role in adaptive immunity. Divergent allele advantage, a mechanism of balancing selection, is proposed to contribute to their exceptional polymorphism. It assumes that MHC genotypes with more divergent alleles allow for broader antigen-presentation to immune effector cells, by that increasing immunocompetence. However, the direct correlation between pairwise sequence divergence and the corresponding repertoire of bound peptides has not been studied systematically across different MHC genes. Here, we investigated this relationship for five key classical human MHC genes (human leukocyte antigen; HLA-A, -B, -C, -DRB1, and -DQB1), using allele-specific computational binding prediction to 118,097 peptides derived from a broad range of human pathogens. For all five human MHC genes, the genetic distance between two alleles of a heterozygous genotype was positively correlated with the total number of peptides bound by these two alleles. In accordance with the major antigen-presentation pathway of MHC class I molecules, HLA-B and HLA-C alleles showed particularly strong correlations for peptides derived from intracellular pathogens. Intriguingly, this bias coincides with distinct protein compositions between intra- and extracellular pathogens, possibly suggesting adaptation of MHC I molecules to present specifically intracellular peptides. Eventually, we observed significant positive correlations between an allele's average divergence and its population frequency. Overall, our results support the divergent allele advantage as a meaningful quantitative mechanism through which pathogen-mediated selection leads to the evolution of MHC diversity.

Genome-wide analysis indicates association between heterozygote advantage and healthy aging in humans.

BACKGROUND: Genetic diversity is known to confer survival advantage in many species across the tree of life. Here, we hypothesize that such pattern applies to humans as well and could be a result of higher fitness in individuals with higher genomic heterozygosity. RESULTS: We use healthy aging as a proxy for better health and fitness, and observe greater heterozygosity in healthy-aged individuals. Specifically, we find that only common genetic variants show significantly higher excess of heterozygosity in the healthy-aged cohort. Lack of difference in heterozygosity for low-frequency variants or disease-associated variants excludes the possibility of compensation for deleterious recessive alleles as a mechanism. In addition, coding SNPs with the highest excess of heterozygosity in the healthy-aged cohort are enriched in genes involved in extracellular matrix and glycoproteins, a group of genes known to be under long-term balancing selection. We also find that individual heterozygosity rate is a significant predictor of electronic health record (EHR)-based estimates of 10-year survival probability in men but not in women, accounting for several factors including age and ethnicity. CONCLUSIONS: Our results demonstrate that the genomic heterozygosity is associated with human healthspan, and that the relationship between higher heterozygosity and healthy aging could be explained by heterozygote advantage. Further characterization of this relationship will have important implications in aging-associated disease risk prediction.

Elevated Heterozygosity in Adults Relative to Juveniles Provides Evidence of Viability Selection on Eagles and Falcons.

Viability selection yields adult populations that are more genetically variable than those of juveniles, producing a positive correlation between heterozygosity and survival. Viability selection could be the result of decreased heterozygosity across many loci in inbred individuals and a subsequent decrease in survivorship resulting from the expression of the deleterious alleles. Alternatively, locus-specific differences in genetic variability between adults and juveniles may be driven by forms of balancing selection, including heterozygote advantage, frequency-dependent selection, or selection across temporal and spatial scales. We use a pooled-sequencing approach to compare genome-wide and locus-specific genetic variability between 74 golden eagle (Aquila chrysaetos), 62 imperial eagle (Aquila heliaca), and 69 prairie falcon (Falco mexicanus) juveniles and adults. Although genome-wide genetic variability is comparable between juvenile and adult golden eagles and prairie falcons, imperial eagle adults are significantly more heterozygous than juveniles. This evidence of viability selection may stem from a relatively smaller imperial eagle effective population size and potentially greater genetic load. We additionally identify ~2000 single-nucleotide polymorphisms across the 3 species with extreme differences in heterozygosity between juveniles and adults. Many of these markers are associated with genes implicated in immune function or olfaction. These loci represent potential targets for studies of how heterozygote advantage, frequency-dependent selection, and selection over spatial and temporal scales influence survivorship in avian species. Overall, our genome-wide data extend previous studies that used allozyme or microsatellite markers and indicate that viability selection may be a more common evolutionary phenomenon than often appreciated.

Multiple Modes of Positive Selection Shaping the Patterns of Incomplete Selective Sweeps over African Populations of Drosophila melanogaster.

It remains a challenge in evolutionary genetics to elucidate how beneficial mutations arise and propagate in a population and how selective pressures on mutant alleles are structured over space and time. By identifying "sweeping haplotypes (SHs)" that putatively carry beneficial alleles and are increasing (or have increased) rapidly in frequency, and surveying the geographic distribution of SH frequencies, we can indirectly infer how selective sweeps unfold in time and thus which modes of positive selection underlie those sweeps. Using population genomic data from African Drosophila melanogaster, we identified SHs from 37 candidate loci under selection. At more than half of loci, we identify single SHs. However, many other loci harbor multiple independent SHs, namely soft selective sweeps, either due to parallel evolution across space or a high beneficial mutation rate. At about a quarter of the loci, intermediate SH frequencies are found across multiple populations, which cannot be explained unless a certain form of frequency-dependent positive selection, such as heterozygote advantage, is invoked given the reasonable range of migration rates between African populations. At one locus, many independent SHs are observed over multiple populations but always together with ancestral haplotypes. This complex pattern is compatible with a large number of mutational targets in a gene and frequency-dependent selection on new variants. We conclude that very diverse modes of positive selection are operating at different sets of loci in D. melanogaster populations.

Transposable Element-Mediated Balancing Selection at Hsp90 Underlies Embryo Developmental Variation.

Understanding the roles of transposable elements (TEs) in the evolution of genome and adaptation is a long-sought goal. Here, we present a new model of TE co-option, in which a TE is harnessed by an essential gene and confers local adaptation through heterozygote advantage. We characterized a human Alu-like TE family, the Lm1 elements, in the genome of the migratory locust Locusta migratoria that harbors 0.7 million copies of the elements. Scanning Lm1 insertions in the natural locust populations revealed the widespread high polymorphism of Lm1. An Lm1 was recruited into the coding region of Heat-shock protein 90 (Hsp90), an important molecular chaperone for diverse signal transduction and developmental pathways. Only heterozygotes of the allele are present in natural populations. Allele frequency increases with decreased latitudes in east coastal China, even increasing up to 76% in southern populations. Regions flanking the Lm1 insertion display clear signatures of a selective sweep linked to Lm1. The Lm1-mediated Hsp90 mutation is consequential for the embryonic development of locust. Heterozygous embryos develop faster than the wild type, particularly when cued by long-day parental photoperiod. The heterozygotes also present a reduced within-population variation in embryonic development, i.e., high developmental synchrony of embryos. The naturally occurring Hsp90 mutation could facilitate multivoltinism and developmental synchronization of the locust in southern tropical region. These results revealed a genetic mechanism behind microevolutionary changes in which balancing selection may have acted to maintain the heterozygote advantage through TE co-option in essential genes.

Long-term balancing selection on chromosomal variants associated with crypsis in a stick insect.

How polymorphisms are maintained within populations over long periods of time remains debated, because genetic drift and various forms of selection are expected to reduce variation. Here, we study the genetic architecture and maintenance of phenotypic morphs that confer crypsis in Timema cristinae stick insects, combining phenotypic information and genotyping-by-sequencing data from 1,360 samples across 21 populations. We find two highly divergent chromosomal variants that span megabases of sequence and are associated with colour polymorphism. We show that these variants exhibit strongly reduced effective recombination, are geographically widespread and probably diverged millions of generations ago. We detect heterokaryotype excess and signs of balancing selection acting on these variants through the species' history. A third chromosomal variant in the same genomic region likely evolved more recently from one of the two colour variants and is associated with dorsal pattern polymorphism. Our results suggest that large-scale genetic variation associated with crypsis has been maintained for long periods of time by potentially complex processes of balancing selection.

Population Differences at MHC Do Not Explain Enhanced Resistance of Song Sparrows to Local Parasites.

Infectious disease represents an emerging threat to natural populations, particularly when hosts are more susceptible to novel parasites (allopatric) than to parasites from the local area (sympatric). This pattern could arise through evolutionary processes (host populations become adapted to their local parasites and genetically differentiated from other populations at immune-related loci) and/or through ecological interactions (host individuals develop resistance to local parasites through previous exposure). The relative importance of these candidate mechanisms remains unclear. In jawed vertebrates, genes of the major histocompatibility complex (MHC) play a fundamental role in immunity and are compelling candidates for spatially varying selection. We recently showed that song sparrows (Melospiza melodia) are more susceptible to allopatric than to sympatric strains of malaria (Plasmodium). In the current study, to determine whether population differences at MHC explain this pattern, we characterized the peptide-binding regions of MHC (classes I and II) of birds that did or did not become infected in the previous experiment. We recovered up to 4 alleles per individual at class I, implying at least 2 loci, and up to 26 alleles per individual at class II, implying at least 13 loci. Individuals with more class I alleles were less likely to become infected by Plasmodium, consistent with parasite-mediated balancing selection. However, we found no evidence for population genetic differentiation at either class of MHC, based on 36 individuals sequenced. Resistance to sympatric parasites previously described for this system likely stems from individuals' prior immune experience, not from population differentiation and locally protective alleles at MHC.

Polymorphisms of the Toll-Like Receptor 2 of Goats (Capra hircus) may be Associated with Somatic Cell Count in Milk.

Toll-like receptor 2 (TLR2) plays an important role in recognition by the innate immune system of Gram-positive bacteria. As Gram-positive bacteria cause mastitis, we examined variations in the region of the TLR2 gene that codes for the extracellular domain. Samples of forty goats from a single dairy herd were collected, half with low SCC (≤200,000 cells/mL), and half with higher SCC. Two synonymous single nucleotide polymorphisms (SNPs) were identified: 840G > A and 1083A > G. One nonsynonymous SNP 739G > A was identified. This coded for valine or isoleucine, which have similar physiochemical properties, and was not in a region coding for a known functional domain. Surprisingly, the least square mean SCC of the heterozygous goats (146,220) was significantly lower than the SCC of homozygous GG goats (537,700; p = 0.004), although these two groups were similar in days in milk (p = 0.984), and there was no significant difference by breed (p = 0.941). Because factors other than mastitis can affect SCC and our sample sizes were limited, additional studies are needed to corroborate an association between TLR2 genotype and SCC or mastitis in goats.

Testing mate choice and overdominance at MH in natural families of Atlantic salmon Salmo salar.

This study aimed to test mate choice and selection during early life stages on major histocompatibility (MH) genotype in natural families of Atlantic salmon Salmo salar spawners and juveniles, using nine microsatellites to reconstruct families, one microsatellite linked to an MH class I gene and one minisatellite linked to an MH class II gene. MH-based mate choice was only detected for the class I locus on the first year, with lower expected heterozygosity in the offspring of actually mated pairs than predicted under random mating. The genotype frequencies of MH-linked loci observed in the juveniles were compared with frequencies expected from Mendelian inheritance of parental alleles to detect selection during early life stages. No selection was detected on the locus linked to class I gene. For the locus linked to class II gene, observed heterozygosity was higher than expected in the first year and lower in the second year, suggesting overdominance and underdominance, respectively. Within family, juveniles' body size was linked to heterozygosity at the same locus, with longer heterozygotes in the first year and longer homozygotes in the second year. Selection therefore seems to differ from one locus to the other and from year to year.

Examining the cause of high inbreeding depression: analysis of whole-genome sequence data in 28 selfed progeny of Eucalyptus grandis.

The genome-wide heterozygosity at 9590 genes, all heterozygous in a single Eucalyptus grandis parent tree, was examined in a group of 28 S1 offspring. Heterozygosity ranged from 52-79%, averaging 65.5%, much higher than the 50% expected under random segregation, supporting the occurrence of strong (47%) selection against homozygosity. The expected pattern of heterozygosity from theoretical calculations and simulations for recessive detrimentals (pseudo-overdominance) and intrinsic heterozygote advantage was examined and compared with that observed. The observed patterns are consistent with at least several detrimental loci with large effects on both parental chromosomes of the 11 pairs. It is likely that 100 or more genes, many with substantial effects on viability, are contributing to this inbreeding depression. Although our genome-wide analysis of nearly 10 000 genes strongly suggested that pseudo-overdominance was responsible for the observed high inbreeding depression, heterozygote advantage could not be excluded. Finding inconvertible evidence of the cause of inbreeding depression still presents a difficult challenge. This study is the first theoretical examination of the genomic effect of inbreeding in a forest tree and provides an approach to analyze these data to determine the extent and cause of inbreeding depression across other plant genomes.

Telomere length reveals cumulative individual and transgenerational inbreeding effects in a passerine bird.

Inbreeding results in more homozygous offspring that should suffer reduced fitness, but it can be difficult to quantify these costs for several reasons. First, inbreeding depression may vary with ecological or physiological stress and only be detectable over long time periods. Second, parental homozygosity may indirectly affect offspring fitness, thus confounding analyses that consider offspring homozygosity alone. Finally, measurement of inbreeding coefficients, survival and reproductive success may often be too crude to detect inbreeding costs in wild populations. Telomere length provides a more precise measure of somatic costs, predicts survival in many species and should reflect differences in somatic condition that result from varying ability to cope with environmental stressors. We studied relative telomere length in a wild population of Seychelles warblers (Acrocephalus sechellensis) to assess the lifelong relationship between individual homozygosity, which reflects genome-wide inbreeding in this species, and telomere length. In juveniles, individual homozygosity was negatively associated with telomere length in poor seasons. In adults, individual homozygosity was consistently negatively related to telomere length, suggesting the accumulation of inbreeding depression during life. Maternal homozygosity also negatively predicted offspring telomere length. Our results show that somatic inbreeding costs are environmentally dependent at certain life stages but may accumulate throughout life.