Gene duplication and divergence produce divergent MHC genotypes without disassortative mating

Diversity of major histocompatibility complex of II B gene and mate choice in a monogamous and long-lived seabird, the Little Auk (Alle alle)

The major histocompatibility complex (MHC) plays a key role in the adaptive immune system of vertebrates, and is known to influence mate choice in many species. In birds, the MHC has been extensively examined but mainly in galliforms and passerines while other taxa that represent specific ecological and evolutionary life-histories, like seabirds, are underexamined. Here, we characterized diversity of MHC Class II B exon 2 in a colonial pelagic seabird, the Little Auk (or Dovekie Alle alle). We further examined whether MHC variation could be maintained through balancing selection and disassortative mating. We found high polymorphism at the genotyped MHC fragment, characterizing 99 distinct alleles across 140 individuals from three populations. The alleles frequencies exhibited a similar skewed distribution in both sexes, with the four most commonly occurring alleles representing approximately 35% of allelic variation. The results of a Bayesian site-by-site selection analysis suggest evidence of balancing selection and no direct evidence for MHC-dependent disassortative mating preferences in the Little Auk. The latter result might be attributed to the high overall polymorphism of the examined fragment, which itself may be maintained by the large population size of the species.

Female and male Leach's Storm Petrels (Hydrobates leucorhous) pursue different foraging strategies during the incubation period

Reproduction in procellariiform birds is characterized by a single egg clutch, slow development, a long breeding season and obligate biparental care. Female Leach's Storm Petrels Hydrobates leucorhous, nearly monomorphic members of this order, produce eggs that are between 20 and 25% of adult body weight. We tested whether female foraging behaviour differs from male foraging behaviour during the ~ 44-day incubation period across seven breeding colonies in the Northwest Atlantic. Over six breeding seasons, we used a combination of Global Positioning System and Global Location Sensor devices to measure characteristics of individual foraging trips during the incubation period. Females travelled significantly greater distances and went farther from the breeding colony than did males on individual foraging trips. For both sexes, the longer the foraging trip, the greater the distance. Independent of trip duration, females travelled farther, and spent a greater proportion of their foraging trips prospecting widely as defined by behavioural categories derived from a Hidden Markov Model. For both sexes, trip duration decreased with date. Sex differences in these foraging metrics were apparently not a consequence of morphological differences or spatial segregation. Our data are consistent with the idea that female foraging strategies differed from male foraging strategies during incubation in ways that would be expected if females were still compensating for egg formation.

Intra-specific copy number variation of MHC class II genes in the Siamese fighting fish

Duplicates of genes for major histocompatibility complex (MHC) molecules can be subjected to selection independently and vary markedly in their evolutionary rates, sequence polymorphism, and functional roles. Therefore, without a thorough understanding of their copy number variation (CNV) in the genome, the MHC-dependent fitness consequences within a species could be misinterpreted. Studying the intra-specific CNV of this highly polymorphic gene, however, has long been hindered by the difficulties in assigning alleles to loci and the lack of high-quality genomic data. Here, using the high-quality genome of the Siamese fighting fish (Betta splendens), a model for mate choice studies, and the whole-genome sequencing (WGS) data of 17 Betta species, we achieved locus-specific amplification of their three classical MHC class II genes - DAB1, DAB2, and DAB3. By performing quantitative PCR and depth-of-coverage analysis using the WGS data, we revealed intra-specific CNV at the DAB3 locus. We identified individuals that had two allelic copies (i.e., heterozygous or homozygous) or one allele (i.e., hemizygous) and individuals without this gene. The CNV was due to the deletion of a 20-kb-long genomic region harboring both the DAA3 and DAB3 genes. We further showed that the three DAB genes were under different modes of selection, which also applies to their corresponding DAA genes that share similar pattern of polymorphism. Our study demonstrates a combined approach to study CNV within a species, which is crucial for the understanding of multigene family evolution and the fitness consequences of CNV.

Characterization, Selection, and Trans-Species Polymorphism in the MHC Class II of Heermann’s Gull (Charadriiformes)

The major histocompatibility complex (MHC) enables vertebrates to cope with pathogens and maintain healthy populations, thus making it a unique set of loci for addressing ecology and evolutionary biology questions. The aim of our study was to examine the variability of Heermann’s Gull MHC class II (MHCIIB) and compare these loci with other Charadriiformes. Fifty-nine MHCIIB haplotypes were recovered from sixty-eight Heermann’s Gulls by cloning, of them, twelve were identified as putative true alleles, forty-five as unique alleles, and two as pseudogenes. Intra and interspecific relationships indicated at least two loci in Heermann’s Gull MHCIIB and trans-species polymorphism among Charadriiformes (coinciding with the documented evidence of two ancient avian MHCIIB lineages, except in the Charadriidae family). Additionally, sites under diversifying selection revealed a better match with peptide-binding sites inferred in birds than those described in humans. Despite the negative anthropogenic activity reported on Isla Rasa, Heermann’s Gull showed MHCIIB variability consistent with population expansion, possibly due to a sudden growth following conservation efforts. Duplication must play an essential role in shaping Charadriiformes MHCIIB variability, buffering selective pressures through balancing selection. These findings suggest that MHC copy number and protected islands can contribute to seabird conservation.

Many loci make light work: High individual diversity despite low population diversity and random mating at class I MHC in a Critically Endangered island songbird

Multi-cellular organisms are under constant attack from parasites, making immune defence a critical aspect of fitness. In vertebrate animals, genes of the major histocompatibility complex (MHC) determine the breadth of pathogens to which individuals can respond. Having many MHC alleles can confer better protection against infectious disease, and balancing selection at MHC is widespread. Indeed, MHC loci are famously variable, with some populations harbouring thousands of alleles (Biedrzycka et al., 2018; Robinson, Soormally, Hayhurst, & Marsh, 2016). MHC has also long fascinated behavioural ecologists because mate choice-for example, preferring MHC-dissimilar partners-may amplify the effects of natural selection (Penn & Potts, 1999). But despite keen interest in the evolutionary ecology of MHC, extensive duplication (Minias, Pikus, Whittingham, & Dunn, 2019) has made these genes challenging to study. In a From the Cover article in this issue of Molecular Ecology, Stervander, Dierickx, Thorley, Brooke, and Westerdahl (2020) characterizes class I MHC in a Critically Endangered songbird, relating genotype to mate choice and survivorship. By inferring copy number and patterns of allelic co-segregation, the authors pave the way to elucidating the genomic architecture of MHC in this bottlenecked population. These insights help reconcile apparently counterintuitive findings: no effect of MHC genotype on mate choice or survival, and high MHC diversity within individuals despite low diversity at the population level. The latter finding is cause for optimism regarding conservation prospects. Moreover, these results suggest that ancient duplication events can have longstanding effects on the adaptive landscapes of natural and sexual selection.

Divergence between genes but limited allelic polymorphism in two MHC class II A genes in Leach's storm-petrels Oceanodroma leucorhoa

The major histocompatibility complex (MHC) is critical to host-pathogen interactions. Class II MHC is a heterodimer, with α and β subunits encoded by different genes. The peptide-binding groove is formed by the first domain of both subunits (α₁ and β₁), but studies of class II variation or natural selection focus primarily on the β subunit and II B genes. We explored MHC II A in Leach's storm-petrel, a seabird with two expressed, polymorphic II B genes. We found two II A genes, Ocle-DAA and Ocle-DBA, in contrast to the single II A gene in chicken and duck. In exon 2 which encodes the α₁ domain, the storm-petrel II A genes differed strongly from each other but showed little within-gene polymorphism in 30 individuals: just one Ocle-DAA allele, and three Ocle-DBA alleles differing from each other by single non-synonymous substitutions. In a comparable sample, the two II B genes had nine markedly diverged alleles each. Differences between the α₁ domains of Ocle-DAA and Ocle-DBA showed signatures of positive selection, but mainly at non-peptide-binding site (PBS) positions. In contrast, positive selection within and between the II B genes corresponded to putative PBS codons. Phylogenetic analysis of the conserved α₂ domain did not reveal deep or well-supported lineages of II A genes in birds, in contrast to the pronounced differentiation of DQA, DPA, and DRA isotypes in mammals. This uncertain homology complicates efforts to compare levels of functional variation and modes of evolution of II A genes across taxa.

Evolution of major histocompatibility complex gene copy number

MHC genes, which code for proteins responsible for presenting pathogen-derived antigens to the host immune system, show remarkable copy-number variation both between and within species. However, the evolutionary forces driving this variation are poorly understood. Here, we use computer simulations to investigate whether evolution of the number of MHC variants in the genome can be shaped by the number of pathogen species the host population encounters (pathogen richness). Our model assumed that while increasing a range of pathogens recognised, expressing additional MHC variants also incurs costs such as an increased risk of autoimmunity. We found that pathogen richness selected for high MHC copy number only when the costs were low. Furthermore, the shape of the association was modified by the rate of pathogen evolution, with faster pathogen mutation rates selecting for increased host MHC copy number, but only when pathogen richness was low to moderate. Thus, taking into account factors other than pathogen richness may help explain wide variation between vertebrate species in the number of MHC genes. Within population, variation in the number of unique MHC variants carried by individuals (INV) was observed under most parameter combinations, except at low pathogen richness. This variance gave rise to positive correlations between INV and host immunocompetence (proportion of pathogens recognised). However, within-population variation in host immunocompetence declined with pathogen richness. Thus, counterintuitively, pathogens can contribute more to genetic variance for host fitness in species exposed to fewer pathogen species, with consequences to predictions from “Hamilton-Zuk” theory of sexual selection.

Highly polymorphic genes of the Major Histocompatibility Complex (MHC) code for proteins responsible for presenting antigens to lymphocytes, thus initiating adaptive immune response. The polymorphism is driven by coevolution with parasites which are selected to evade recognition by MHC proteins. Expressing many MHC molecules could ensure that an individual could present antigens of most pathogen species encountered, but this comes at a cost, such as enhanced negative selection on lymphocytes leading to holes in T-cell receptor repertoire. Our simulations showed that evolution of the number of MHC genes in the genome is driven by a complex interaction between three factors we explored: pathogen richness, the intrinsic cost of expressing additional MHC variants, and pathogen mutation rate. In contrast to verbal arguments, our results indicate that pathogen richness does not always selects for MHC gene family expansion. Taking into account factors other than pathogen richness, in particular costs of expressing additional MHC variants which are still poorly understood, may help explain striking interspecific variation in the number of MHC genes. Counterintuitively, our results also demonstrated that opportunity for selection on immunocompetence should decrease with MHC gene family expansion.

Ecology can inform genetics: Disassortative mating contributes to MHC polymorphism in Leach's storm-petrels (Oceanodroma leucorhoa)

Studies of MHC-based mate choice in wild populations often test hypotheses on species exhibiting female choice and male-male competition, which reflects the general prevalence of females as the choosy sex in natural systems. Here, we examined mutual mate-choice patterns in a small burrow-nesting seabird, the Leach's storm-petrel (Oceanodroma leucorhoa), using the major histocompatibility complex (MHC). The life history and ecology of this species are extreme: both partners work together to fledge a single chick during the breeding season, a task that requires regularly travelling hundreds of kilometres to and from foraging grounds over a 6- to 8-week provisioning period. Using a 5-year data set unprecedented for this species (n = 1078 adults and 925 chicks), we found a positive relationship between variation in the likelihood of female reproductive success and heterozygosity at Ocle-DAB2, a MHC class IIB locus. Contrary to previous reports rejecting disassortative mating as a mechanism for maintaining genetic polymorphism in this species, here we show that males make significant disassortative mate-choice decisions. Variability in female reproductive success suggests that the most common homozygous females (Ocle-DAB2*01/Ocle-DAB2*01) may be physiologically disadvantaged and, therefore, less preferred as lifelong partners for choosy males. The results from this study support the role of mate choice in maintaining high levels of MHC variability in a wild seabird species and highlight the need to incorporate a broader ecological framework and sufficient sample sizes into studies of MHC-based mating patterns in wild populations in general.

Lack of evidence for selection favouring MHC haplotypes that combine high functional diversity

High rates of gene duplication and the highest levels of functional allelic diversity in vertebrate genomes are the main hallmarks of the major histocompatibility complex (MHC), a multigene family with a primordial role in pathogen recognition. The usual tight linkage among MHC gene duplicates may provide an opportunity for the evolution of haplotypes that associate functionally divergent alleles and thus grant the transmission of optimal levels of diversity to coming generations. Even though such associations may be a crucial component of disease resistance, this hypothesis has been given little attention in wild populations. Here, we leveraged pedigree data from a barn owl (Tyto alba) population to characterize MHC haplotype structure across two MHC class I (MHC-I) and two MHC class IIB (MHC-IIB) duplicates, in order to test the hypothesis that haplotypes' genetic diversity is higher than expected from randomly associated alleles. After showing that MHC loci are tightly linked within classes, we found limited evidence for shifts towards MHC haplotypes combining high diversity. Neither amino acid nor functional within-haplotype diversity were significantly higher than in random sets of haplotypes, regardless of MHC class. Our results therefore provide no evidence for selection towards high-diversity MHC haplotypes in barn owls. Rather, high rates of concerted evolution may constrain the evolution of high-diversity haplotypes at MHC-I, while, in contrast, for MHC-IIB, fixed differences among loci may provide barn owls with already optimized functional diversity. This suggests that at the MHC-I and MHC-IIB respectively, different evolutionary dynamics may govern the evolution of within-haplotype diversity.

Reduced population size does not affect the mating strategy of a vulnerable and endemic seabird

Bottleneck episodes may occur in small and isolated animal populations, which may result in decreased genetic diversity and increased inbreeding, but also in mating strategy adjustment. This was evaluated in the vulnerable and socially monogamous Monteiro's Storm-petrel Hydrobates monteiroi, a seabird endemic to the Azores archipelago which has suffered a dramatic population decline since the XVth century. To do this, we conducted a genetic study (18 microsatellite markers) in the population from Praia islet, which has been monitored over 16 years. We found no evidence that a genetic bottleneck was associated with this demographic decline. Monteiro's Storm-petrels paired randomly with respect to genetic relatedness and body measurements. Pair fecundity was unrelated to genetic relatedness between partners. We detected only two cases of extra-pair parentage associated with an extra-pair copulation (out of 71 offspring). Unsuccessful pairs were most likely to divorce the next year, but genetic relatedness between pair mates and pair breeding experience did not influence divorce. Divorce enabled individuals to improve their reproductive performances after re-mating only when the new partner was experienced. Re-pairing with an experienced partner occurred more frequently when divorcees changed nest than when they retained their nest. This study shows that even in strongly reduced populations, genetic diversity can be maintained, inbreeding does not necessarily occur, and random pairing is not risky in terms of pair lifetime reproductive success. Given, however, that we found no clear phenotypic mate choice criteria, the part played by non-morphological traits should be assessed more accurately in order to better understand seabird mating strategies.

Morphological and genetic factors shape the microbiome of a seabird species (Oceanodroma leucorhoa) more than environmental and social factors

BACKGROUND

The microbiome provides multiple benefits to animal hosts that can profoundly impact health and behavior. Microbiomes are well-characterized in humans and other animals in controlled settings, yet assessments of wild bird microbial communities remain vastly understudied. This is particularly true for pelagic seabirds with unique life histories that differ from terrestrial bird species. This study was designed to examine how morphological, genetic, environmental, and social factors affect the microbiome of a burrow-nesting seabird species, Leach's storm petrel (Oceanodroma leucorhoa). These seabirds are highly olfactory and may rely on microbiome-mediated odor cues during mate selection. Composition and structure of bacterial communities associated with the uropygial gland and brood patch were assessed using 16S rRNA amplicon-based Illumina Mi-Seq analysis and compared to burrow-associated bacterial communities. This is the first study to examine microbial diversity associated with multiple body sites on a seabird species.

RESULTS

Results indicate that sex and skin site contribute most to bacterial community variation in Leach's storm petrels and that major histocompatibility complex (MHC) genotype may impact the composition of bacterial assemblages in males. In contrast to terrestrial birds and other animals, environmental and social interactions do not significantly influence storm petrel-associated bacterial assemblages. Thus, individual morphological and genetic influences outweighed environmental and social factors on microbiome composition.

CONCLUSIONS

Contrary to observations of terrestrial birds, microbiomes of Leach's storm petrels vary most by the sex of the bird and by the body site sampled, rather than environmental surroundings or social behavior.

Very high MHC Class IIB diversity without spatial differentiation in the mediterranean population of greater Flamingos

Background

Selective pressure from pathogens is thought to shape the allelic diversity of major histocompatibility complex (MHC) genes in vertebrates. In particular, both local adaptation to pathogens and gene flow are thought to explain a large part of the intraspecific variation observed in MHC allelic diversity. To date, however, evidence that adaptation to locally prevalent pathogens maintains MHC variation is limited to species with limited dispersal and, hence, reduced gene flow. On the one hand high gene flow can disrupt local adaptation in species with high dispersal rates, on the other hand such species are much more likely to experience spatial variation in pathogen pressure, suggesting that there may be intense pathogen mediated selection pressure operating across breeding sites in panmictic species. Such pathogen mediated selection pressure operating across breeding sites should therefore be sufficient to maintain high MHC diversity in high dispersing species in the absence of local adaptation mechanisms. We used the Greater Flamingo, Phoenicopterus roseus, a long-lived colonial bird showing a homogeneous genetic structure of neutral markers at the scale of the Mediterranean region, to test the prediction that higher MHC allelic diversity with no population structure should occur in large panmictic populations of long-distance dispersing birds than in other resident species.

Results

We assessed the level of allelic diversity at the MHC Class IIB exon 2 from 116 individuals born in four different breeding colonies of Greater Flamingo in the Mediterranean region. We found one of the highest allelic diversity (109 alleles, 2 loci) of any non-passerine avian species investigated so far relative to the number of individuals and loci genotyped. There was no evidence of population structure between the four major Mediterranean breeding colonies.

Conclusion

Our results suggest that local adaptation at MHC Class IIB in Greater Flamingos is constrained by high gene flow and high MHC diversity appears to be maintained by population wide pathogen-mediated selection rather than local pathogen-mediated selection. Further understanding of how pathogens vary across space and time will be crucial to further elucidate the mechanisms maintaining MHC diversity in species with large panmictic populations and high dispersal rates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0905-3) contains supplementary material, which is available to authorized users.