Genotyping strategy matters when analyzing hypervariable major histocompatibility complex-Experience from a passerine bird

Associations among MHC genes, latitude, and avian malaria infections in the rufous-collared sparrow (Zonotrichia capensis)

The major histocompatibility complex (MHC) is a genetic region in jawed vertebrates that contains key genes involved in the immune response. Associations between the MHC and avian malaria infections in wild birds have been observed and mainly explored in the Northern Hemisphere, while a general lack of information remains in the Southern Hemisphere. Here, we investigated the associations between the MHC genes and infections with Plasmodium and Haemoproteus blood parasites along a latitudinal gradient in South America. We sampled 93 rufous-collared sparrows (Zonotrichia capensis) individuals from four countries, Colombia, Ecuador, Peru, and Chile, and estimated MHC-I and MHC-II allele diversity. We detected between 1-4 (MHC-I) and 1-6 (MHC-II) amino acidic alleles per individual, with signs of positive selection. We obtained generalized additive mixed models to explore the associations between MHC-I and MHC-II diversity and latitude. We also explored the relationship between infection status and latitude/biome. We found a non-linear association between the MHC-II amino acidic allele diversity and latitude. Individuals from north Chile presented a lower MHC genetic diversity than those from other locations. We also found an association between deserts and xeric shrublands and a lower prevalence of Haemoproteus parasites. Our results support a lower MHC genetic in arid or semi-arid habitats in the region with the lower prevalence of Haemoproteus parasites.

A pan-cetacean MHC amplicon sequencing panel developed and evaluated in combination with genome assemblies

The major histocompatibility complex (MHC) is a highly polymorphic gene family that is crucial in immunity, and its diversity can be effectively used as a fitness marker for populations. Despite this, MHC remains poorly characterised in non-model species (e.g., cetaceans: whales, dolphins and porpoises) as high gene copy number variation, especially in the fast-evolving class I region, makes analyses of genomic sequences difficult. To date, only small sections of class I and IIa genes have been used to assess functional diversity in cetacean populations. Here, we undertook a systematic characterisation of the MHC class I and IIa regions in available cetacean genomes. We extracted full-length gene sequences to design pan-cetacean primers that amplified the complete exon 2 from MHC class I and IIa genes in one combined sequencing panel. We validated this panel in 19 cetacean species and described 354 alleles for both classes. Furthermore, we identified likely assembly artefacts for many MHC class I assemblies based on the presence of class I genes in the amplicon data compared to missing genes from genomes. Finally, we investigated MHC diversity using the panel in 25 humpback and 30 southern right whales, including four paternity trios for humpback whales. This revealed copy-number variable class I haplotypes in humpback whales, which is likely a common phenomenon across cetaceans. These MHC alleles will form the basis for a cetacean branch of the Immuno-Polymorphism Database (IPD-MHC), a curated resource intended to aid in the systematic compilation of MHC alleles across several species, to support conservation initiatives.

Template-specific optimization of NGS genotyping pipelines reveals allele-specific variation in MHC gene expression

Using high-throughput sequencing for precise genotyping of multi-locus gene families, such as the major histocompatibility complex (MHC), remains challenging, due to the complexity of the data and difficulties in distinguishing genuine from erroneous variants. Several dedicated genotyping pipelines for data from high-throughput sequencing, such as next-generation sequencing (NGS), have been developed to tackle the ensuing risk of artificially inflated diversity. Here, we thoroughly assess three such multi-locus genotyping pipelines for NGS data, the DOC method, AmpliSAS and ACACIA, using MHC class IIβ data sets of three-spined stickleback gDNA, cDNA and "artificial" plasmid samples with known allelic diversity. We show that genotyping of gDNA and plasmid samples at optimal pipeline parameters was highly accurate and reproducible across methods. However, for cDNA data, the gDNA-optimal parameter configuration yielded decreased overall genotyping precision and consistency between pipelines. Further adjustments of key clustering parameters were required tο account for higher error rates and larger variation in sequencing depth per allele, highlighting the importance of template-specific pipeline optimization for reliable genotyping of multi-locus gene families. Through accurate paired gDNA-cDNA typing and MHC-II haplotype inference, we show that MHC-II allele-specific expression levels correlate negatively with allele number across haplotypes. Lastly, sibship-assisted cDNA-typing of MHC-I revealed novel variants linked in haplotype blocks, and a higher-than-previously-reported individual MHC-I allelic diversity. In conclusion, we provide novel genotyping protocols for the three-spined stickleback MHC-I and -II genes, and evaluate the performance of popular NGS-genotyping pipelines. We also show that fine-tuned genotyping of paired gDNA-cDNA samples facilitates amplification bias-corrected MHC allele expression analysis.

Do female bluethroats without extra-pair offspring have more MHC-compatible social mates?

Abstract

Genes of the major histocompatibility complex (MHC) are crucial for adaptive immunity in jawed vertebrates, and theory predicts that there should be mate choice for optimizing MHC constitution in the offspring. In a previous study, we demonstrated a non-random female choice of extra-pair males in the bluethroat (Luscinia svecica), yielding offspring that was closer to an intermediate MHC class II (MHCII) allele count than their within-pair halfsiblings. The present study tests whether social pairs with only within-pair young (WPY) in their brood, in the same study population, had a combined MHC-constitution closer to a presumed intermediate optimum, than social pairs with extra-pair young (EPY), with a corresponding pattern in their offspring. As expected, we found that WPY from pure WPY-broods were more MHC-optimal than WPY from mixed broods, but only in broods of young (second year) males. Correspondingly, there was a tendency for social pairs with only WPY in their brood to be more MHC-compatible than social pairs with EPY in their brood, when the male was young. Older bluethroat males have considerably larger testes than young males, and their higher sperm competitiveness could help them secure paternity in their own brood, also when they are not MHC-compatible. In other words, in the sexual conflict over paternity, females may be more likely to realise their preference for a MHC-compatible mate when paired to a young male. As a possible fitness indicator, immune responsiveness to an injected antigen (PHA) was elevated for offspring closer to “the golden mean” in MHCII allele count.

Significance statement

This study contributes to our understanding of MHC-based mate choice in extra-pair mating systems, by showing that female bluethroats (Luscinia svecica) with an MHCII-compatible social mate tend to have no extra-pair young in their brood, but only when the social male is young. This elucidates a possible sexual conflict, in which older social males are able to override female preferences and prevent other males from gaining paternity in their brood through higher sperm production. Studying systems in which extra-pair paternity occurs offers an insight into the genetic benefits of mate choice, as extra-pair males, in contrast to social males, generally contribute only sperm. Further, the strict and thorough genotyping scheme applied in this study enabled us to demonstrate a preference for “the golden mean” in MHC-diversity in a species with one of the highest MHC class II-diversity known to date.

Divergent allele advantage in the MHC and amphibian emerging infectious disease

Genetic variation in the major histocompatibility complex (MHC) may be associated with resistance to the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd). The pathogen originated in Asia, then spread worldwide, causing amphibian population declines and species extinctions. We compared the expressed MHC IIβ1 alleles of a Bd-resistant toad species, Bufo gargarizans, from South Korea with those of a Bd-susceptible Australasian frog species, Litoria caerulea. We found at least six expressed MHC IIβ1 loci in each of the two species. Amino acid diversity encoded by these MHC alleles was similar between species, but the genetic divergence of those alleles known for broader pathogen-derived peptide binding was greater in the Bd-resistant species. In addition, we found a potentially rare allele in one resistant individual from the Bd-susceptible species. Deep next-generation sequencing recovered approximately triple the genetic resolution accessible from traditional cloning-based genotyping. Targeting more than one MHC IIβ1 expressed locus enables us to better understand how host MHC may adapt to emerging infectious diseases.

Diversifying selection and climatic effects on major histocompatibility complex class II gene diversity in the greater horseshoe bat

Heterogeneous pathogenic stress can shape major histocompatibility complex (MHC) diversity by influencing the functional plasticity of the immune response. Therefore, MHC diversity could reflect environmental stress, demonstrating its importance in uncovering the mechanisms of adaptive genetic variation. In this study, we combined neutral microsatellite loci, an immune-related MHC II-DRB locus, and climatic factors to unravel the mechanisms affecting the diversity and genetic differentiation of MHC genes in the greater horseshoe bat (Rhinolophus ferrumequinum), a species with a wide geographical distribution that has three distinct genetic lineages in China. First, increased genetic differentiation at the MHC locus among populations compared using microsatellites indicated diversifying selection. Second, the genetic differentiation of MHC and microsatellites were significantly correlated, suggesting that demographic processes exist. However, MHC genetic differentiation was significantly correlated with geographical distance among populations, even after controlling for the neutral markers, suggesting a major effect of selection. Third, although the MHC genetic differentiation was larger than that for microsatellites, there was no significant difference in the genetic differentiation between the two markers among genetic lineages, indicating the effect of balancing selection. Fourth, combined with climatic factors, MHC diversity and supertypes showed significant correlations with temperature and precipitation, but not with the phylogeographic structure of R. ferrumequinum, suggesting an effect of local adaptation driven by climate on MHC diversity. Moreover, the number of MHC supertypes varied between populations and lineages, suggesting regional characteristics and support for local adaptation. Taken together, the results of our study provide insights into the adaptive evolutionary driving forces at different geographic scales in R. ferrumequinum. In addition, climate factors may have played a vital role in driving adaptive evolution in this species.

Performance Comparison of Different Approaches in Genotyping MHC-DRB: The Contrast between Single-Locus and Multi-Locus Species

Major histocompatibility complex (MHC) genes are widely recognised as valuable markers for wildlife genetic studies given their extreme polymorphism and functional importance in fitness-related traits. Newly developed genotyping methods, which rely on the use of next-generation sequencing (NGS), are gradually replacing traditional cloning and Sanger sequencing methods in MHC genotyping studies. Allele calling in NGS methods remains challenging due to extreme polymorphism and locus multiplication in the MHC coupled with allele amplification bias and the generation of artificial sequences. In this study, we compared the performance of molecular cloning with Illumina and Ion Torrent NGS sequencing in MHC-DRB genotyping of single-locus species (roe deer) and species with multiple DRB loci (red deer) in an attempt to adopt a reliable and straightforward method that does not require complex bioinformatic analyses. Our results show that all methods work similarly well in roe deer, but we demonstrate non-consistency in results across methods in red deer. With Illumina sequencing, we detected a maximum number of alleles in 10 red deer individuals (42), while other methods were somewhat less accurate as they scored 69–81% of alleles detected with Illumina sequencing.

Adaptive evolution of major histocompatibility complex class I immune genes and disease associations in coastal juvenile sea turtles

Characterizing polymorphism at the major histocompatibility complex (MHC) genes is key to understanding the vertebrate immune response to disease. Despite being globally afflicted by the infectious tumour disease fibropapillomatosis (FP), immunogenetic variation in sea turtles is minimally explored. We sequenced the α 1 peptide-binding region of MHC class I genes (162 bp) from 268 juvenile green (Chelonia mydas) and 88 loggerhead (Caretta caretta) sea turtles in Florida, USA. We recovered extensive variation (116 alleles) and trans-species polymorphism. Supertyping analysis uncovered three functional MHC supertypes corresponding to the three well-supported clades in the phylogeny. We found significant evidence of positive selection at seven amino acid sites in the class I exon. Random forest modelling and risk ratio analysis of Ch. mydas alleles uncovered one allele weakly associated with smooth FP tumour texture, which may be associated with disease outcome. Our study represents the first characterization of MHC class I diversity in Ch. mydas and the largest sample of sea turtles used to date in any study of adaptive genetic variation, revealing tremendous genetic variation and high adaptive potential to viral pathogen threats. The novel associations we identified between MHC diversity and FP outcomes in sea turtles further highlight the importance of evaluating genetic predictors of disease, including MHC and other functional markers.

Improved high-throughput MHC typing for non-model species using long-read sequencing

The major histocompatibility complex (MHC) plays a critical role in the vertebrate immune system. Accurate MHC typing is critical to understanding not only host fitness and disease susceptibility, but also the mechanisms underlying host‐pathogen co‐evolution. However, due to the high degree of gene duplication and diversification of MHC genes, it is often technically challenging to accurately characterise MHC genetic diversity in non‐model species. Here we conducted a systematic review to identify common issues associated with current widely used MHC typing approaches. Then to overcome these challenges, we developed a long‐read based MHC typing method along with a new analysis pipeline. Our approach enables the sequencing of fully phased MHC alleles spanning all key functional domains and the separation of highly similar alleles as well as the removal of technical artefacts such as PCR heteroduplexes and chimeras. Using this approach, we performed population‐scale MHC typing in the Tasmanian devil (Sarcophilus harrisii), revealing previously undiscovered MHC functional diversity in this endangered species. Our new method provides a better solution for addressing research questions that require high MHC typing accuracy. Since the method is not limited by species or the number of genes analysed, it will be applicable for studying not only the MHC but also other complex gene families.

Avian blood parasite richness decreases with major histocompatibility complex class I loci number

Major histocompatibility complex (MHC) genes are among the most polymorphic in the vertebrate genome. The high allele diversity is believed to be maintained primarily by sexual and pathogen-mediated balancing selection. The number of MHC loci also varies greatly across vertebrates, most notably across birds. MHC proteins play key roles in presenting antigens on the cell surface for recognition by T cells, with class I proteins specifically targeting intracellular pathogens. Here, we explore the hypothesis that MHC class I diversity (measured as loci number) coevolves with haemosporidian parasite burden of the host. Using data on 54 bird species, we demonstrate that high-MHC class I diversity is associated with significantly lower richness of Plasmodium, Haemoproteus as well as overall haemosporidian parasite lineages, the former thus indicating more efficient protection against intracellular pathogens. Nonetheless, the latter associations were only detected when MHC diversity was assessed using cloning and not 454 pyrosequencing-based studies, nor across all genotyping methods combined. Our results indicate that high-MHC class I diversity might play a key role in providing qualitative resistance against diverse haemosporidian parasites in birds, but further clarification is needed for the origin of contrasting results when using different genotyping methods for MHC loci quantification.

MHC Genotyping by SSCP and Amplicon-Based NGS Approach in Chamois

Genes of the major histocompatibility complex (MHC) code for cell surface proteins essential for adaptive immunity. They show the most outstanding genetic diversity in vertebrates, which has been connected with various fitness traits and thus with the long-term persistence of populations. In this study, polymorphism of the MHC class II DRB locus was investigated in chamois with Single-Strand Conformation Polymorphism (SSCP)/Sanger genotyping and Ion Torrent S5 next-generation sequencing (NGS). From eight identified DRB variants in 28 individuals, five had already been described, and three were new, undescribed alleles. With conventional SSCP/Sanger sequencing, we were able to detect seven alleles, all of which were also detected with NGS. We found inconsistencies in the individual genotypes between the two methods, which were mainly caused by allelic dropout in the SSCP/Sanger method. Six out of 28 individuals were falsely classified as homozygous with SSCP/Sanger analysis. Overall, 25% of the individuals were identified as genotyping discrepancies between the two methods. Our results show that NGS technologies are better performing in sequencing highly variable regions such as the MHC, and they also have a higher detection capacity, thus allowing a more accurate description of the genetic composition, which is crucial for evolutionary and population genetic studies.

Lack of statistical power as a major limitation in understanding MHC-mediated immunocompetence in wild vertebrate populations

Disentangling the sources of variation in developing an effective immune response against pathogens is of major interest to immunoecology and evolutionary biology. To date, the link between immunocompetence and genetic variation at the major histocompatibility complex (MHC) has received little attention in wild animals, despite the key role of MHC genes in activating the adaptive immune system. Although several studies point to a link between MHC and immunocompetence, negative findings have also been reported. Such disparate findings suggest that limited statistical power might be affecting studies on this topic, owing to insufficient sample sizes and/or a generally small effect of MHC on the immunocompetence of wild vertebrates. To clarify this issue, we investigated the link between MHC variation and seven immunocompetence proxies in a large sample of barn owls and estimated the effect sizes and statistical power of this and published studies on this topic. We found that MHC poorly explained variation in immunocompetence of barn owls, with small-to-moderate associations between MHC and immunocompetence in owls (effect size: .1 ≥ r ≤ .3) similar to other vertebrates studied to date. Such small-to-moderate effects were largely associated with insufficient power, which was only sufficient (>0.8) to detect moderate-to-large effect sizes (r ≥ .3). Thus, studies linking MHC variation with immunocompetence in wild populations are underpowered to detect MHC effects, which are likely to be of generally small magnitude. Larger sample sizes (>200) will be required to achieve sufficient power in future studies aiming to robustly test for a link between MHC variation and immunocompetence.

Endless forms of sexual selection

In recent years, the field of sexual selection has exploded, with advances in theoretical and empirical research complementing each other in exciting ways. This perspective piece is the product of a "stock-taking" workshop on sexual selection and sexual conflict. Our aim is to identify and deliberate on outstanding questions and to stimulate discussion rather than provide a comprehensive overview of the entire field. These questions are organized into four thematic sections we deem essential to the field. First we focus on the evolution of mate choice and mating systems. Variation in mate quality can generate both competition and choice in the opposite sex, with implications for the evolution of mating systems. Limitations on mate choice may dictate the importance of direct vs. indirect benefits in mating decisions and consequently, mating systems, especially with regard to polyandry. Second, we focus on how sender and receiver mechanisms shape signal design. Mediation of honest signal content likely depends on integration of temporally variable social and physiological costs that are challenging to measure. We view the neuroethology of sensory and cognitive receiver biases as the main key to signal form and the 'aesthetic sense' proposed by Darwin. Since a receiver bias is sufficient to both initiate and drive ornament or armament exaggeration, without a genetically correlated or even coevolving receiver, this may be the appropriate 'null model' of sexual selection. Thirdly, we focus on the genetic architecture of sexually selected traits. Despite advances in modern molecular techniques, the number and identity of genes underlying performance, display and secondary sexual traits remains largely unknown. In-depth investigations into the genetic basis of sexual dimorphism in the context of long-term field studies will reveal constraints and trajectories of sexually selected trait evolution. Finally, we focus on sexual selection and conflict as drivers of speciation. Population divergence and speciation are often influenced by an interplay between sexual and natural selection. The extent to which sexual selection promotes or counteracts population divergence may vary depending on the genetic architecture of traits as well as the covariance between mating competition and local adaptation. Additionally, post-copulatory processes, such as selection against heterospecific sperm, may influence the importance of sexual selection in speciation. We propose that efforts to resolve these four themes can catalyze conceptual progress in the field of sexual selection, and we offer potential avenues of research to advance this progress.