Only rare classical MHC-I alleles are highly expressed in the European house sparrow.

The exceptional polymorphism observed within genes of the major histocompatibility complex (MHC), a core component of the vertebrate immune system, has long fascinated biologists. The highly polymorphic classical MHC class-I (MHC-I) genes are maintained by pathogen-mediated balancing selection (PMBS), as shown by many sites subject to positive selection, while the more monomorphic non-classical MHC-I genes show signatures of purifying selection. In line with PMBS, at any point in time, rare classical MHC alleles are more likely than common classical MHC alleles to confer a selective advantage in host-pathogen interactions. Combining genomic and expression data from the blood of wild house sparrows Passer domesticus, we found that only rare classical MHC-I alleles were highly expressed, while common classical MHC-I alleles were lowly expressed or not expressed. Moreover, highly expressed rare classical MHC-I alleles had more positively selected sites, indicating exposure to stronger PMBS, compared with lowly expressed classical alleles. As predicted, the level of expression was unrelated to allele frequency in the monomorphic non-classical MHC-I alleles. Going beyond previous studies, we offer a fine-scale view of selection on classical MHC-I genes in a wild population by revealing differences in the strength of PMBS according to allele frequency and expression level.

Improved haplotype resolution of highly duplicated MHC genes in a long-read genome assembly using MiSeq amplicons.

Long-read sequencing offers a great improvement in the assembly of complex genomic regions, such as the major histocompatibility complex (MHC) region, which can contain both tandemly duplicated MHC genes (paralogs) and high repeat content. The MHC genes have expanded in passerine birds, resulting in numerous MHC paralogs, with relatively high sequence similarity, making the assembly of the MHC region challenging even with long-read sequencing. In addition, MHC genes show rather high sequence divergence between alleles, making diploid-aware assemblers incorrectly classify haplotypes from the same locus as sequences originating from different genomic regions. Consequently, the number of MHC paralogs can easily be over- or underestimated in long-read assemblies. We therefore set out to verify the MHC diversity in an original and a haplotype-purged long-read assembly of one great reed warbler Acrocephalus arundinaceus individual (the focal individual) by using Illumina MiSeq amplicon sequencing. Single exons, representing MHC class I (MHC-I) and class IIB (MHC-IIB) alleles, were sequenced in the focal individual and mapped to the annotated MHC alleles in the original long-read genome assembly. Eighty-four percent of the annotated MHC-I alleles in the original long-read genome assembly were detected using 55% of the amplicon alleles and likewise, 78% of the annotated MHC-IIB alleles were detected using 61% of the amplicon alleles, indicating an incomplete annotation of MHC genes. In the haploid genome assembly, each MHC-IIB gene should be represented by one allele. The parental origin of the MHC-IIB amplicon alleles in the focal individual was determined by sequencing MHC-IIB in its parents. Two of five larger scaffolds, containing 6-19 MHC-IIB paralogs, had a maternal and paternal origin, respectively, as well as a high nucleotide similarity, which suggests that these scaffolds had been incorrectly assigned as belonging to different loci in the genome rather than as alternate haplotypes of the same locus. Therefore, the number of MHC-IIB paralogs was overestimated in the haploid genome assembly. Based on our findings we propose amplicon sequencing as a suitable complement to long-read sequencing for independent validation of the number of paralogs in general and for haplotype inference in multigene families in particular.

The structure of songbird MHC class I reveals antigen binding that is flexible at the N-terminus and static at the C-terminus.

Long-distance migratory animals such as birds and bats have evolved to withstand selection imposed by pathogens across the globe, and pathogen richness is known to be particularly high in tropical regions. Immune genes, so-called Major Histocompatibility Complex (MHC) genes, are highly duplicated in songbirds compared to other vertebrates, and this high MHC diversity has been hypothesised to result in a unique adaptive immunity. To understand the rationale behind the evolution of the high MHC genetic diversity in songbirds, we determined the structural properties of an MHC class I protein, Acar3, from a long-distance migratory songbird, the great reed warbler Acrocephalus arundinaceus (in short: Acar). The structure of Acar3 was studied in complex with pathogen-derived antigens and shows an overall antigen presentation similar to human MHC class I. However, the peptides bound to Acar3 display an unusual conformation: Whereas the N-terminal ends of the peptides display enhanced flexibility, the conformation of their C-terminal halves is rather static. This uncommon peptide-binding mode in Acar3 is facilitated by a central Arg residue within the peptide-binding groove that fixes the backbone of the peptide at its central position, and potentially permits successful interactions between MHC class I and innate immune receptors. Our study highlights the importance of investigating the immune system of wild animals, such as birds and bats, to uncover unique immune mechanisms which may neither exist in humans nor in model organisms.

Understanding the evolution of immune genes in jawed vertebrates.

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

MHC class II genotype-by-pathogen genotype interaction for infection prevalence in a natural rodent-Borrelia system.

MHC genes are extraordinarily polymorphic in most taxa. Host-pathogen coevolution driven by negative frequency-dependent selection (NFDS) is one of the main hypotheses for the maintenance of such immunogenetic variation. Here, we test a critical but rarely tested assumption of this hypothesis-that MHC alleles affect resistance/susceptibility to a pathogen in a strain-specific way, that is, there is a host genotype-by-pathogen genotype interaction. In a field study of bank voles naturally infected with the tick-transmitted bacterium Borrelia afzelii, we tested for MHC class II (DQB) genotype-by-B. afzelii strain interactions for infection prevalence between 10 DQB alleles and seven strains. One allele (DQB*37) showed an interaction, such that voles carrying DQB*37 had higher prevalence of two strains and lower prevalence of one strain than individuals without the allele. These findings were corroborated by analyses of strain composition of infections, which revealed an effect of DQB*37 in the form of lower ß diversity among infections in voles carrying the allele. Taken together, these results provide rare support at the molecular genetic level for a key assumption of models of antagonistic coevolution through NFDS.

Reciprocal positive effects on parasitemia between coinfecting haemosporidian parasites in house sparrows.

BACKGROUND: Hosts are often simultaneously infected with several parasite species. These co-infections can lead to within-host interactions of parasites, including mutualism and competition, which may affect both virulence and transmission. Birds are frequently co-infected with different haemosporidian parasites, but very little is known about if and how these parasites interact in natural host populations and what consequences there are for the infected hosts. We therefore set out to study Plasmodium and Haemoproteus parasites in house sparrows Passer domesticus with naturally acquired infections using a protocol where the parasitemia (infection intensity) is quantified by qPCR separately for the two parasites. We analysed infection status (presence/absence of the parasite) and parasitemia of parasites in the blood of both adult and juvenile house sparrows repeatedly over the season. RESULTS: Haemoproteus passeris and Plasmodium relictum were the two dominating parasite species, found in 99% of the analyzed Sanger sequences. All birds were infected with both Plasmodium and Haemoproteus parasites during the study period. Seasonality explained infection status for both parasites in the adults: H. passeris was completely absent in the winter while P. relictum was present all year round. Among adults infected with H. passeris there was a positive effect of P. relictum parasitemia on H. passeris parasitemia and likewise among adults infected with P. relictum there was a positive effect of H. passeris parasitemia on P. relictum parasitemia. No such associations on parasitemia were seen in juvenile house sparrows. CONCLUSIONS: The reciprocal positive relationships in parasitemia between P. relictum and H. passeris in adult house sparrows suggests either mutualistic interactions between these frequently occurring parasites or that there is variation in immune responses among house sparrow individuals, hence some individuals suppress the parasitemia of both parasites whereas other individuals suppress neither. Our detailed screening of haemosporidian parasites over the season shows that co-infections are very frequent in both juvenile and adult house sparrows, and since co-infections often have stronger negative effects on host fitness than the single infection, it is imperative to use screening systems with the ability to detect multiple parasites in ecological studies of host-parasite interactions.

MHCtools - an R package for MHC high-throughput sequencing data: Genotyping, haplotype and supertype inference, and downstream genetic analyses in non-model organisms.

The major histocompatibility complex (MHC) plays a central role in the vertebrate adaptive immune system and has been of long-term interest in evolutionary biology. While several protocols have been developed for MHC genotyping, there is a lack of transparent and standardized tools for downstream analysis of MHC data. Here, we present the r package mhctools and demonstrate the use of its functions to (i) assist accurate MHC genotyping from high-throughput amplicon-sequencing data, (ii) infer functional MHC supertypes using bootstrapped clustering analysis, (iii) identify segregating MHC haplotypes from family data, and (iv) analyse functional and genetic distances between MHC sequences. We employed mhctools to analyse MHC class I (MHC-I) amplicon data of 559 great reed warblers (Acrocephalus arundinaceus). We identified 390 MHC-I alleles which clustered into 14 functional supertypes. A phylogenetic analysis and analyses of positive selection suggested that the MHC-I alleles belong to several distinct functional groups. We furthermore identified 107 segregating haplotypes among 116 families, and found substantial variation in diversity with 4-21 MHC-I alleles and 3-13 MHC-I supertypes per haplotype. Finally, we show that the great reed warbler haplotypes harboured combinations of MHC-I supertypes with greater functional divergence than observed in simulated populations of possible haplotypes, a result that is in accordance with the divergent allele advantage hypothesis. Our study demonstrates the power of mhctools to support genotyping and analysis of MHC in non-model species, which we hope will encourage broad implementation among researchers in MHC genetics and evolution.

The genomic architecture of the passerine MHC region: High repeat content and contrasting evolutionary histories of single copy and tandemly duplicated MHC genes.

The major histocompatibility complex (MHC) is of central importance to the immune system, and an optimal MHC diversity is believed to maximize pathogen elimination. Birds show substantial variation in MHC diversity, ranging from few genes in most bird orders to very many genes in passerines. Our understanding of the evolutionary trajectories of the MHC in passerines is hampered by lack of data on genomic organization. Therefore, we assembled and annotated the MHC genomic region of the great reed warbler (Acrocephalus arundinaceus), using long-read sequencing and optical mapping. The MHC region is large (>5.5 Mb), characterized by structural changes compared to hitherto investigated bird orders and shows higher repeat content than the genome average. These features were supported by analyses in three additional passerines. MHC genes in passerines are found in two different chromosomal arrangements, either as single copy MHC genes located among non-MHC genes, or as tandemly duplicated tightly linked MHC genes. Some single copy MHC genes are old and putative orthologues among species. In contrast tandemly duplicated MHC genes are monophyletic within species and have evolved by simultaneous gene duplication of several MHC genes. Structural differences in the MHC genomic region among bird orders seem substantial compared to mammals and have possibly been fuelled by clade-specific immune system adaptations. Our study provides methodological guidance in characterizing complex genomic regions, constitutes a resource for MHC research in birds, and calls for a revision of the general belief that avian MHC has a conserved gene order and small size compared to mammals.

Avian Neo-Sex Chromosomes Reveal Dynamics of Recombination Suppression and W Degeneration.

How the avian sex chromosomes first evolved from autosomes remains elusive as 100 million years (My) of divergence and degeneration obscure their evolutionary history. The Sylvioidea group of songbirds is interesting for understanding avian sex chromosome evolution because a chromosome fusion event ∼24 Ma formed "neo-sex chromosomes" consisting of an added (new) and an ancestral (old) part. Here, we report the complete female genome (ZW) of one Sylvioidea species, the great reed warbler (Acrocephalus arundinaceus). Our long-read assembly shows that the added region has been translocated to both Z and W, and whereas the added-Z has retained its gene order the added-W part has been heavily rearranged. Phylogenetic analyses show that recombination between the homologous added-Z and -W regions continued after the fusion event, and that recombination suppression across this region took several million years to be completed. Moreover, recombination suppression was initiated across multiple positions over the added-Z, which is not consistent with a simple linear progression starting from the fusion point. As expected following recombination suppression, the added-W show signs of degeneration including repeat accumulation and gene loss. Finally, we present evidence for nonrandom maintenance of slowly evolving and dosage-sensitive genes on both ancestral- and added-W, a process causing correlated evolution among orthologous genes across broad taxonomic groups, regardless of sex linkage.

Borrelia Infection in Bank Voles Myodes glareolus Is Associated With Specific DQB Haplotypes Which Affect Allelic Divergence Within Individuals.

The high polymorphism of Major Histocompatibility Complex (MHC) genes is generally considered to be a result of pathogen-mediated balancing selection. Such selection may operate in the form of heterozygote advantage, and/or through specific MHC allele-pathogen interactions. Specific MHC allele-pathogen interactions may promote polymorphism via negative frequency-dependent selection (NFDS), or selection that varies in time and/or space because of variability in the composition of the pathogen community (fluctuating selection; FS). In addition, divergent allele advantage (DAA) may act on top of these forms of balancing selection, explaining the high sequence divergence between MHC alleles. DAA has primarily been thought of as an extension of heterozygote advantage. However, DAA could also work in concert with NFDS though this is yet to be tested explicitly. To evaluate the importance of DAA in pathogen-mediated balancing selection, we surveyed allelic polymorphism of MHC class II DQB genes in wild bank voles (Myodes glareolus) and tested for associations between DQB haplotypes and infection by Borrelia afzelii, a tick-transmitted bacterium causing Lyme disease in humans. We found two significant associations between DQB haplotypes and infection status: one haplotype was associated with lower risk of infection (resistance), while another was associated with higher risk of infection (susceptibility). Interestingly, allelic divergence within individuals was higher for voles with the resistance haplotype compared to other voles. In contrast, allelic divergence was lower for voles with the susceptibility haplotype than other voles. The pattern of higher allelic divergence in individuals with the resistance haplotype is consistent with NFDS favouring divergent alleles in a natural population, hence selection where DAA works in concert with NFDS.

Whole-genome sequencing and genome regions of special interest: Lessons from major histocompatibility complex, sex determination, and plant self-incompatibility.

Whole-genome sequencing of non-model organisms is now widely accessible and has allowed a range of questions in the field of molecular ecology to be investigated with greater power. However, some genomic regions that are of high biological interest remain problematic for assembly and data-handling. Three such regions are the major histocompatibility complex (MHC), sex-determining regions (SDRs) and the plant self-incompatibility locus (S-locus). Using these as examples, we illustrate the challenges of both assembling and resequencing these highly polymorphic regions and how bioinformatic and technological developments are enabling new approaches to their study. Mapping short-read sequences against multiple alternative references improves genotyping comprehensiveness at the S-locus thereby contributing to more accurate assessments of allelic frequencies. Long-read sequencing, producing reads of several tens to hundreds of kilobase pairs in length, facilitates the assembly of such regions as single sequences can span the multiple duplicated gene copies of the MHC region, and sequence through repetitive stretches and translocations in SDRs and S-locus haplotypes. These advances are adding value to short-read genome resequencing approaches by allowing, for example, more accurate haplotype phasing across longer regions. Finally, we assessed further technical improvements, such as nanopore adaptive sequencing and bioinformatic tools using pangenomes, which have the potential to further expand our knowledge of a number of genomic regions that remain challenging to study with classical resequencing approaches.

Trade-offs in expressed major histocompatibility complex diversity seen on a macroevolutionary scale among songbirds.

To survive organisms must defend themselves against pathogens. Classical Major Histocompatibility Complex (MHC) genes play a key role in pathogen defense by encoding molecules involved in pathogen recognition. MHC gene diversity influences the variety of pathogens individuals can recognize and respond to and has consequently been a popular genetic marker for disease resistance in ecology and evolution. However, MHC diversity is predominantly estimated using genomic DNA (gDNA) with little knowledge of expressed diversity. This limits our ability to interpret the adaptive significance of variation in MHC diversity, especially in species with very many MHC genes such as songbirds. Here, we address this issue using phylogenetic comparative analyses of the number of MHC class I alleles (MHC-I diversity) in gDNA and complementary DNA (cDNA), that is, expressed alleles, across 13 songbird species. We propose three theoretical relationships that could be expected between genomic and expressed MHC-I diversity on a macroevolutionary scale and test which of these are best supported. In doing so, we show that significantly fewer MHC-I alleles than the number available are expressed, suggesting that optimal MHC-I diversity could be achieved by modulating gene expression. Understanding the relationship between genomic and expressed MHC diversity is essential for interpreting variation in MHC diversity in an evolutionary context.

Seasonal dynamics of haemosporidian (Apicomplexa, Haemosporida) parasites in house sparrows Passer domesticus at four European sites: comparison between lineages and the importance of screening methods.

Infectious diseases often vary seasonally in a predictable manner, and seasonality may be responsible for geographical differences in prevalence. In temperate regions, vector-borne parasites such as malaria are expected to evolve lower virulence and a time-varying strategy to invest more in transmission when vectors are available. A previous model of seasonal variation of avian malaria described a double peak in prevalence of Plasmodium parasites in multiple hosts resulting from spring relapses and transmission to susceptible individuals in summer. However, this model was rejected by a study describing different patterns of seasonal variation of two Plasmodium spp. at the same site, with the double peak only apparent when these species were combined. Here, we assessed the seasonal variation in prevalence of haemosporidian parasites (Plasmodium, Haemoproteus and Leucocytozoon) in house sparrows (Passer domesticus) sampled across 1 year at four temperate European sites spanning a latitudinal range of 17°. We showed that parasite prevalence and diversity decreased with increasing latitude, but the parasite communities differed between sites, with only one Plasmodium lineage (P_SGS1) occurring at all sites. Moreover, the nested PCR method commonly used to detect and identify haemosporidian parasites strongly underestimated co-infections of Haemoproteus and Plasmodium, significantly biasing the pattern of seasonal variation, so additional molecular methods were used. Finally, we showed that: (i) seasonal variation in prevalence of haemosporidian parasites varied between study sites and parasite lineages/species/genera, describing further cases where the double peak model is not met; (ii) the seasonal dynamics of single lineages (P_SGS1) varied between sites; and (iii) unexpectedly, seasonality was greatest at the southernmost site, a pattern that was mostly driven by lineage H_PADOM05. Limitations of the genotyping methods and consequences of pooling (parasite lineages, sites and years) in studies of haemosporidian parasites are discussed and recommendations proposed, since these actions may obscure the patterns of prevalence and limit ecological inferences.

High MHC gene copy number maintains diversity despite homozygosity in a Critically Endangered single-island endemic bird, but no evidence of MHC-based mate choice.

Small population sizes can, over time, put species at risk due to the loss of genetic variation and the deleterious effects of inbreeding. Losing diversity in the major histocompatibility complex (MHC) could be particularly harmful, given its key role in the immune system. Here, we assess MHC class I (MHC-I) diversity and its effects on mate choice and survival in the Critically Endangered Raso lark Alauda razae, a species restricted to the 7 km2 islet of Raso, Cape Verde, since ~1460, whose population size has dropped as low as 20 pairs. Exhaustively genotyping 122 individuals, we find no effect of MHC-I genotype/diversity on mate choice or survival. However, we demonstrate that MHC-I diversity has been maintained through extreme bottlenecks by retention of a high number of gene copies (at least 14), aided by cosegregation of multiple haplotypes comprising 2-8 linked MHC-I loci. Within-locus homozygosity is high, contributing to low population-wide diversity. Conversely, each individual had comparably many alleles, 6-16 (average 11), and the large and divergent haplotypes occur at high frequency in the population, resulting in high within-individual MHC-I diversity. This functional immune gene diversity will be of critical importance for this highly threatened species' adaptive potential.

Wetter climates select for higher immune gene diversity in resident, but not migratory, songbirds.

Pathogen communities can vary substantially between geographical regions due to different environmental conditions. However, little is known about how host immune systems respond to environmental variation across macro-ecological and evolutionary scales. Here, we select 37 species of songbird that inhabit diverse environments, including African and Palaearctic residents and Afro-Palaearctic migrants, to address how climate and habitat have influenced the evolution of key immune genes, the major histocompatibility complex class I (MHC-I). Resident species living in wetter regions, especially in Africa, had higher MHC-I diversity than species living in drier regions, irrespective of the habitats they occupy. By contrast, no relationship was found between MHC-I diversity and precipitation in migrants. Our results suggest that the immune system of birds has evolved greater pathogen recognition in wetter tropical regions. Furthermore, evolving transcontinental migration appears to have enabled species to escape wet, pathogen-rich areas at key periods of the year, relaxing selection for diversity in immune genes and potentially reducing immune system costs.

Not all birds have a single dominantly expressed MHC-I gene: Transcription suggests that siskins have many highly expressed MHC-I genes.

Passerine birds belong to the most species rich bird order and are found in a wide range of habitats. The extremely polymorphic adaptive immune system of passerines, identified through their major histocompatibility complex class I genes (MHC-I), may explain some of this extreme radiation. Recent work has shown that passerines have higher numbers of MHC-I gene copies than other birds, but little is currently known about expression and function of these gene copies. Non-passerine birds have a single highly expressed MHC-I gene copy, a pattern that seems unlikely in passerines. We used high-throughput sequencing to study MHC-I alleles in siskins (Spinus spinus) and determined gene expression, phylogenetic relationships and sequence divergence. We verified between six and 16 MHC-I alleles per individual and 97% of these were expressed. Strikingly, up to five alleles per individual had high expression. Out of 88 alleles 18 were putatively non-classical with low sequence divergence and expression, and found in a single phylogenetic cluster. The remaining 70 alleles were classical, with high sequence divergence and variable degrees of expression. Our results contradict the suggestion that birds only have a single dominantly expressed MHC-I gene by demonstrating several highly expressed MHC-I gene copies in a passerine.

Genomics of host-pathogen interactions: challenges and opportunities across ecological and spatiotemporal scales.

Evolutionary genomics has recently entered a new era in the study of host-pathogen interactions. A variety of novel genomic techniques has transformed the identification, detection and classification of both hosts and pathogens, allowing a greater resolution that helps decipher their underlying dynamics and provides novel insights into their environmental context. Nevertheless, many challenges to a general understanding of host-pathogen interactions remain, in particular in the synthesis and integration of concepts and findings across a variety of systems and different spatiotemporal and ecological scales. In this perspective we aim to highlight some of the commonalities and complexities across diverse studies of host-pathogen interactions, with a focus on ecological, spatiotemporal variation, and the choice of genomic methods used. We performed a quantitative review of recent literature to investigate links, patterns and potential tradeoffs between the complexity of genomic, ecological and spatiotemporal scales undertaken in individual host-pathogen studies. We found that the majority of studies used whole genome resolution to address their research objectives across a broad range of ecological scales, especially when focusing on the pathogen side of the interaction. Nevertheless, genomic studies conducted in a complex spatiotemporal context are currently rare in the literature. Because processes of host-pathogen interactions can be understood at multiple scales, from molecular-, cellular-, and physiological-scales to the levels of populations and ecosystems, we conclude that a major obstacle for synthesis across diverse host-pathogen systems is that data are collected on widely diverging scales with different degrees of resolution. This disparity not only hampers effective infrastructural organization of the data but also data granularity and accessibility. Comprehensive metadata deposited in association with genomic data in easily accessible databases will allow greater inference across systems in the future, especially when combined with open data standards and practices. The standardization and comparability of such data will facilitate early detection of emerging infectious diseases as well as studies of the impact of anthropogenic stressors, such as climate change, on disease dynamics in humans and wildlife.

Avian MHC Evolution in the Era of Genomics: Phase 1.0.

Birds are a wonderfully diverse and accessible clade with an exceptional range of ecologies and behaviors, making the study of the avian major histocompatibility complex (MHC) of great interest. In the last 20 years, particularly with the advent of high-throughput sequencing, the avian MHC has been explored in great depth in several dimensions: its ability to explain ecological patterns in nature, such as mating preferences; its correlation with parasite resistance; and its structural evolution across the avian tree of life. Here, we review the latest pulse of avian MHC studies spurred by high-throughput sequencing. Despite high-throughput approaches to MHC studies, substantial areas remain in need of improvement with regard to our understanding of MHC structure, diversity, and evolution. Recent studies of the avian MHC have nonetheless revealed intriguing connections between MHC structure and life history traits, and highlight the advantages of long-term ecological studies for understanding the patterns of MHC variation in the wild. Given the exceptional diversity of birds, their accessibility, and the ease of sequencing their genomes, studies of avian MHC promise to improve our understanding of the many dimensions and consequences of MHC variation in nature. However, significant improvements in assembling complete MHC regions with long-read sequencing will be required for truly transformative studies.

Plumage microbiota covaries with the major histocompatibility complex in blue petrels.

To increase fitness, a wide range of vertebrates preferentially mate with partners that are dissimilar at the major histocompatibility complex (MHC) or that have high MHC diversity. Although MHC often can be assessed through olfactory cues, the mechanism by which MHC genes influence odour remains largely unclear. MHC class IIB molecules, which enable recognition and elimination of extracellular bacteria, have been suggested to influence odour indirectly by shaping odour-producing microbiota, i.e. bacterial communities. However, there is little evidence of the predicted covariation between an animal's MHC genotype and its bacterial communities in scent-producing body surfaces. Here, using high-throughput sequencing, we tested the covariation between MHC class IIB genotypes and feather microbiota in the blue petrel (Halobaena caerulea), a seabird with highly developed olfaction that has been suggested to rely on oduor cues during an MHC-based mate choice. First, we show that individuals with similar MHC class IIB profiles also have similar bacterial assemblages in their feathers. Then, we show that individuals with high MHC diversity have less diverse feather microbiota and also a reduced abundance of a bacterium of the genus Arsenophonus, a genus in which some species are symbionts of avian ectoparasites. Our results, showing that feather microbiota covary with MHC, are consistent with the hypothesis that individual MHC genotype may shape the semiochemical-producing microbiota in birds.

Evidence for sexual conflict over major histocompatibility complex diversity in a wild songbird.

Sex differences in parasite load and immune responses are found across a wide range of animals, with females generally having lower parasite loads and stronger immune responses than males. Intrigued by these general patterns, we investigated if there was any sign of sex-specific selection on an essential component of adaptive immunity that is known to affect fitness, the major histocompatibility complex class I (MHC-I) genes, in a 20-year study of great reed warblers. Our analyses on fitness related to MHC-I diversity showed a highly significant interaction between MHC-I diversity and sex, where males with higher, and females with lower, MHC-I diversity were more successful in recruiting offspring. Importantly, mean MHC-I diversity did not differ between males and females, and consequently neither sex reached its MHC-I fitness optimum. Thus, there is an unresolved genetic sexual conflict over MHC-I diversity in great reed warblers. Selection from pathogens is known to maintain MHC diversity, but previous theory ignores that the immune environments are considerably different in males and females. Our results suggest that sexually antagonistic selection is an important, previously neglected, force in the evolution of vertebrate adaptive immunity, and have implications for evolutionary understanding of costs of immune responses and autoimmune diseases.