Non-neutral evolution and reciprocal monophyly of two expressed Mhc class II B genes in Leach’s storm-petrel

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.

Extensive MHC class IIβ diversity across multiple loci in the small-spotted catshark (Scyliorhinus canicula)

The major histocompatibility complex (MHC) is a multigene family responsible for pathogen detection, and initiation of adaptive immune responses. Duplication, natural selection, recombination, and their resulting high functional genetic diversity spread across several duplicated loci are the main hallmarks of the MHC. Although these features were described in several jawed vertebrate lineages, a detailed MHC IIβ characterization at the population level is still lacking for chondrichthyans (chimaeras, rays and sharks), i.e. the most basal lineage to possess an MHC-based adaptive immune system. We used the small-spotted catshark (Scyliorhinus canicula, Carcharhiniformes) as a case-study species to characterize MHC IIβ diversity using complementary molecular tools, including publicly available genome and transcriptome datasets, and a newly developed high-throughput Illumina sequencing protocol. We identified three MHC IIβ loci within the same genomic region, all of which are expressed in different tissues. Genetic screening of the exon 2 in 41 individuals of S. canicula from a single population revealed high levels of sequence diversity, evidence for positive selection, and footprints of recombination. Moreover, the results also suggest the presence of copy number variation in MHC IIβ genes. Thus, the small-spotted catshark exhibits characteristics of functional MHC IIβ genes typically observed in other jawed vertebrates.

Feather chemicals contain information about the major histocompatibility complex in a highly scented seabird

Mate choice informed by the immune genes of the major histocompatibility complex (MHC) may provide fitness benefits including offspring with increased immunocompetence. Olfactory cues are considered the primary mechanism organisms use to evaluate the MHC of potential mates, yet this idea has received limited attention in birds. Motivated by a finding of MHC-dependent mate choice in the Leach's storm-petrel (Oceanodroma leucorhoa), we examined whether the chemical profiles of this highly scented seabird contain information about MHC genes. Whereas previous studies in birds examined non-volatile compounds, we used gas chromatography-mass spectrometry to measure the volatile compounds emitted from feathers that potentially serve as olfactory infochemicals about MHC and coupled this with locus-specific genotyping of MHC IIB genes. We found that feather chemicals reflected individual MHC diversity through interactions with sex and breeding status. Furthermore, similarity in MHC genotype was correlated with similarity in chemical profiles within female-female and male-female dyads. We provide the first evidence that volatile chemicals from bird feathers can encode information about the MHC. Our findings suggest that olfaction likely aids MHC-based mate choice in this species and highlight a role for chemicals in mediating genetic mate choice in birds where this mode of communication has been largely overlooked.

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.

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.

The probability of reciprocal monophyly of gene lineages in three and four species

Reciprocal monophyly, a feature of a genealogy in which multiple groups of descendant lineages each consist of all of the descendants of their respective most recent common ancestors, has been an important concept in studies of species delimitation, phylogeography, population history reconstruction, systematics, and conservation. Computations involving the probability that reciprocal monophyly is observed in a genealogy have played a key role in criteria for defining taxonomic groups and inferring divergence times. The probability of reciprocal monophyly under a coalescent model of population divergence has been studied in detail for groups of gene lineages for pairs of species. Here, we extend this computation to generate corresponding probabilities for sets of gene lineages from three and four species. We study the effects of model parameters on the probability of reciprocal monophyly, finding that it is driven primarily by species tree height, with lesser but still substantial influences of internal branch lengths and sample sizes. We also provide an example application of our results to data from maize and teosinte.

Annual global mean temperature explains reproductive success in a marine vertebrate from 1955 to 2010

The salient feature of anthropogenic climate change over the last century has been the rise in global mean temperature. However, global mean temperature is not used as an explanatory variable in studies of population-level response to climate change, perhaps because the signal-to-noise ratio of this gross measure makes its effect difficult to detect in any but the longest of datasets. Using a population of Leach's storm-petrels breeding in the Bay of Fundy, we tested whether local, regional, or global temperature measures are the best index of reproductive success in the face of climate change in species that travel widely between and within seasons. With a 56-year dataset, we found that annual global mean temperature (AGMT) was the single most important predictor of hatching success, more so than regional sea surface temperatures (breeding season or winter) and local air temperatures at the nesting colony. Storm-petrel reproductive success showed a quadratic response to rising temperatures, in that hatching success increased up to some critical temperature, and then declined when AGMT exceeded that temperature. The year at which AGMT began to consistently exceed that critical temperature was 1988. Importantly, in this population of known-age individuals, the impact of changing climate was greatest on inexperienced breeders: reproductive success of inexperienced birds increased more rapidly as temperatures rose and declined more rapidly after the tipping point than did reproductive success of experienced individuals. The generality of our finding that AGMT is the best predictor of reproductive success in this system may hinge on two things. First, an integrative global measure may be best for species in which individuals move across an enormous spatial range, especially within seasons. Second, the length of our dataset and our capacity to account for individual- and age-based variation in reproductive success increase our ability to detect a noisy signal.

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.

100 million years of multigene family evolution: origin and evolution of the avian MHC class IIB

BACKGROUND

Gene duplication has led to a most remarkable adaptation involved in vertebrates' host-pathogen arms-race, the major histocompatibility complex (MHC). However, MHC duplication history is as yet poorly understood in non-mammalian vertebrates, including birds.

RESULTS

Here, we provide evidence for the evolution of two ancient avian MHC class IIB (MHCIIB) lineages by a duplication event prior to the radiation of all extant birds >100 million years ago, and document the role of concerted evolution in eroding the footprints of the avian MHCIIB duplication history.

CONCLUSIONS

Our results suggest that eroded footprints of gene duplication histories may mimic birth-death evolution and that in the avian MHC the presence of the two lineages may have been masked by elevated rates of concerted evolution in several taxa. Through the presence of a range of intermediate evolutionary stages along the homogenizing process of concerted evolution, the avian MHCIIB provides a remarkable illustration of the erosion of multigene family duplication history.

MHC class I BFIV gene polymorphisms in four Chinese native chicken breeds

The major histocompatibility complex (MHC) includes the most polymorphic genes in vertebrates, and balancing selection has been proposed as a main evolutionary force. Here we present one of the first data sets examining the genetic characteristics of chicken MHC I BFIV molecules in four Chinese native breeds, sourced from different regions in China. In all, 89 BFIV alleles were isolated from 102 individuals sampled, and 13 repeated alleles were observed. No significant correlation was found between genetic differentiation and geographical distance in the phylogenetic tree. BFIV genes exhibited a high level of nucleotide polymorphisms, and most of the polymorphic sites were located in the peptide-binding region (PBR) encoded in exons 2 and 3. A comparison of the three-dimensional structures of PBRs in chicken BFIV and human HLA-A molecules revealed evident structural and functional similarities. The results suggested that MHC I molecules had similar structural features in different species.

Gene duplication and divergence produce divergent MHC genotypes without disassortative mating

Genes of the major histocompatibility complex (MHC) exhibit heterozygote advantage in immune defence, which in turn can select for MHC-disassortative mate choice. However, many species lack this expected pattern of MHC-disassortative mating. A possible explanation lies in evolutionary processes following gene duplication: if two duplicated MHC genes become functionally diverged from each other, offspring will inherit diverse multilocus genotypes even under random mating. We used locus-specific primers for high-throughput sequencing of two expressed MHC Class II B genes in Leach's storm-petrels, Oceanodroma leucorhoa, and found that exon 2 alleles fall into two gene-specific monophyletic clades. We tested for disassortative vs. random mating at these two functionally diverged Class II B genes, using multiple metrics and different subsets of exon 2 sequence data. With good statistical power, we consistently found random assortment of mates at MHC. Despite random mating, birds had MHC genotypes with functionally diverged alleles, averaging 13 amino acid differences in pairwise comparisons of exon 2 alleles within individuals. To test whether this high MHC diversity in individuals is driven by evolutionary divergence of the two duplicated genes, we built a phylogenetic permutation model. The model showed that genotypic diversity was strongly impacted by sequence divergence between the most common allele of each gene, with a smaller additional impact of monophyly of the two genes. Divergence of allele sequences between genes may have reduced the benefits of actively seeking MHC-dissimilar mates, in which case the evolutionary history of duplicated genes is shaping the adaptive landscape of sexual selection.

Next-generation genotyping of hypervariable loci in many individuals of a non-model species: technical and theoretical implications

BACKGROUND

Across species, diversity at the Major Histocompatibility Complex (MHC) is critical to disease resistance and population health; however, use of MHC diversity to quantify the genetic health of populations has been hampered by the extreme variation found in MHC genes. Next generation sequencing (NGS) technology generates sufficient data to genotype even the most diverse species, but workflows for distinguishing artifacts from alleles are still under development. We used NGS to evaluate the MHC diversity of over 300 captive and wild ring-tailed lemurs (Lemur catta: Primates: Mammalia). We modified a published workflow to address errors that arise from deep sequencing individuals and tested for evidence of selection at the most diverse MHC genes.

RESULTS

In addition to evaluating the accuracy of 454 Titanium and Ion Torrent PGM for genotyping large populations at hypervariable genes, we suggested modifications to improve current methods of allele calling. Using these modifications, we genotyped 302 out of 319 individuals, obtaining an average sequencing depth of over 1000 reads per amplicon. We identified 55 MHC-DRB alleles, 51 of which were previously undescribed, and provide the first sequences of five additional MHC genes: DOA, DOB, DPA, DQA, and DRA. The additional five MHC genes had one or two alleles each with little sequence variation; however, the 55 MHC-DRB alleles showed a high dN/dS ratio and trans-species polymorphism, indicating a history of positive selection. Because each individual possessed 1-7 MHC-DRB alleles, we suggest that ring-tailed lemurs have four, putatively functional, MHC-DRB copies.

CONCLUSIONS

In the future, accurate genotyping methods for NGS data will be critical to assessing genetic variation in non-model species. We recommend that future NGS studies increase the proportion of replicated samples, both within and across platforms, particularly for hypervariable genes like the MHC. Quantifying MHC diversity within non-model species is the first step to assessing the relationship of genetic diversity at functional loci to individual fitness and population viability. Owing to MHC-DRB diversity and copy number, ring-tailed lemurs may serve as an ideal model for estimating the interaction between genetic diversity, fitness, and environment, especially regarding endangered species.