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

Genetic variation at innate and adaptive immune genes - contrasting patterns of differentiation and local adaptation in a wild gull

Immunogenetic variation in natural vertebrate populations is expected to respond to spatial and temporal fluctuations in pathogen assemblages. While spatial heterogeneity in pathogen-driven selection enhances local immunogenetic adaptations and population divergence, different immune genes may yield contrasting responses to the environment. Here, we investigated population differentiation at the key pathogen recognition genes of the innate and adaptive immune system in a colonial bird species, the black-headed gull Chroicocephalus ridibundus. We assessed genetic variation at three toll-like receptor (TLR) genes (innate immunity) and the major histocompatibility complex (MHC) class I and II genes (adaptive immunity) in gulls from seven colonies scattered across Poland. As expected, we found much greater polymorphism at the MHC than TLRs. Population differentiation at the MHC class II, but not MHC-I, was significantly stronger than at neutral microsatellite loci, suggesting local adaptation. This could reflect spatial variation in the composition of extracellular parasite communities (e.g., helminths), possibly driven by sharp differences in habitat structure between colonies. Despite contrasting patterns of population differentiation, both MHC classes showed similar regimes of diversifying selection. Some significant population differentiation was also observed at TLRs, suggesting that innate immune receptors may respond to fine-scale spatial variation in pathogen pressure, although this pattern could have been enhanced by drift. Our results suggested that local adaptation at the pathogen recognition immune genes can be maintained at relatively small or moderate spatial scales in species with high dispersal potential and they highlighted the complexity of immunogenetic responses of animals to heterogeneous environments.

High functional allelic diversity and copy number in both MHC classes in the common buzzard

BACKGROUND

The major histocompatibility complex (MHC), which encodes molecules that recognize various pathogens and parasites and initiates the adaptive immune response in vertebrates, is renowned for its exceptional polymorphism and is a model of adaptive gene evolution. In birds, the number of MHC genes and sequence diversity varies greatly among taxa, believed due to evolutionary history and differential selection pressures. Earlier characterization studies and recent comparative studies suggest that non-passerine species have relatively few MHC gene copies compared to passerines. Additionally, comparative studies that have looked at partial MHC sequences have speculated that non-passerines have opposite patterns of selection on MHC class I (MHC-I) and class II (MHC-II) loci than passerines: namely, greater sequence diversity and signals of selection on MHC-II than MHC-I. However, new sequencing technology is revealing much greater MHC variation than previously expected while also facilitating full sequence variant detection directly from genomic data. Our study aims to take advantage of high-throughput sequencing methods to fully characterize both classes and domains of MHC of a non-passerine bird of prey, the common buzzard (Buteo buteo), to test predictions of MHC variation and differential selection on MHC classes.

RESULTS

Using genetic, genomic, and transcriptomic high-throughput sequencing data, we established common buzzards have at least three loci that produce functional alleles at both MHC classes. In total, we characterize 91 alleles from 113 common buzzard chicks for MHC-I exon 3 and 41 alleles from 125 chicks for MHC-IIB exon 2. Among these alleles, we found greater sequence polymorphism and stronger diversifying selection at MHC-IIB exon 2 than MHC-I exon 3, suggesting differential selection pressures on MHC classes. However, upon further investigation of the entire peptide-binding groove by including genomic data from MHC-I exon 2 and MHC-IIA exon 2, this turned out to be false. MHC-I exon 2 was as polymorphic as MHC-IIB exon 2 and MHC-IIA exon 2 was essentially invariant. Thus, comparisons between MHC-I and MHC-II that included both domains of the peptide-binding groove showed no differences in polymorphism nor diversifying selection between the classes. Nevertheless, selection analysis indicates balancing selection has been acting on common buzzard MHC and phylogenetic inference revealed that trans-species polymorphism is present between common buzzards and species separated for over 33 million years for class I and class II.

CONCLUSIONS

We characterize and confirm the functionality of unexpectedly high copy number and allelic diversity in both MHC classes of a bird of prey. While balancing selection is acting on both classes, there is no evidence of differential selection pressure on MHC classes in common buzzards and this result may hold more generally once more data for understudied MHC exons becomes available.

Both selection and drift drive the spatial pattern of adaptive genetic variation in a wild mammal

The major histocompatibility complex (MHC) has been intensively studied for the relative effects of different evolutionary forces in recent decades. Pathogen-mediated balancing selection is generally thought to explain the high polymorphism observed in MHC genes, but it is still unclear to what extent MHC diversity is shaped by selection relative to neutral drift. In this study, we genotyped MHC class II DRB genes and 15 neutral microsatellite loci across 26 geographic populations of European badgers (Meles meles) covering most of their geographic range. By comparing variation of microsatellite and diversity of MHC at different levels, we demonstrate that both balancing selection and drift have shaped the evolution of MHC genes. When only MHC allelic identity was investigated, the spatial pattern of MHC variation was similar to that of microsatellites. By contrast, when functional aspects of the MHC diversity (e.g., immunological supertypes) were considered, balancing selection appears to decrease genetic structuring across populations. Our comprehensive sampling and analytical approach enable us to conclude that the likely mechanisms of selection are heterozygote advantage and/or rare-allele advantage. This study is a clear demonstration of how both balancing selection and genetic drift simultaneously affect the evolution of MHC genes in a widely distributed wild mammal.

Functional immune diversity in reindeer reveals a high Arctic population at risk

Climate changes the geographic range of both species as well as pathogens, causing a potential increase in the vulnerability of populations or species with limited genetic diversity. With advances in high throughput sequencing (HTS) technologies, we can now define functional expressed genetic diversity of wild species at a larger scale and identify populations at risk. Previous studies have used genomic DNA to define major histocompatibility complex (MHC) class II diversity in reindeer. Varying numbers of expressed genes found in many ungulates strongly argues for using cDNA in MHC typing strategies to ensure that diversity estimates relate to functional genes. We have used available reindeer genomes to identify candidate genes and established an HTS approach to define expressed MHC class I and class II diversity. To capture a broad diversity we included samples from wild reindeer from Southern Norway, semi-domesticated reindeer from Northern Norway and reindeer from the high Artic archipelago Svalbard. Our data show a medium MHC diversity in semi-domesticated and wild Norwegian mainland reindeer, and low MHC diversity reindeer in Svalbard reindeer. The low immune diversity in Svalbard reindeer provides a potential risk if the pathogenic pressure changes in response to altered environmental conditions due to climate change, or increased human-related activity.

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

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

Using de novo genome assembly and high-throughput sequencing to characterize the MHC region in a non-model bird, the Eurasian coot

Genes of the Major Histocompatibility Complex (MHC) form a key component of vertebrate adaptive immunity, as they code for molecules which bind antigens of intra- and extracellular pathogens (MHC class I and II, respectively) and present them to T cell receptors. In general, MHC genes are hyper-polymorphic and high MHC diversity is often maintained within natural populations (via balancing selection) and within individuals (via gene duplications). Because of its complex architecture with tandems of duplicated genes, characterization of MHC region in non-model vertebrate species still poses a major challenge. Here, we combined de novo genome assembly and high-throughput sequencing to characterize MHC polymorphism in a rallid bird species, the Eurasian coot Fulica atra. An analysis of genome assembly indicated high duplication rate at MHC-I, which was also supported by targeted sequencing of peptide-binding exons (at least five MHC-I loci genotyped). We found high allelic richness at both MHC-I and MHC-II, although signature of diversifying selection and recombination (gene conversion) was much stronger at MHC-II. Our results indicate that Eurasian coot retains extraordinary polymorphism at both MHC classes (when compared to other non-passerine bird species), although they may be subject to different evolutionary mechanism.

Bioaccumulation of trace elements affects chick body condition and gut microbiome in greater flamingos

Elevated concentrations of trace elements represent a major concern to wetland ecosystems, since river estuaries are geochemical endpoints that accumulate pollution. Although the negative impact of environmental exposure of highly toxic elements such as Pb and Hg has received substantial attention, we still lack a comprehensive understanding of the effects that these and other common trace elements have on natural populations. We used greater flamingos as a study system within three sites that represent a gradient of pollution. Controlling for environmental sediment exposure, we assessed if signatures of bioaccumulation in feathers for ten trace elements (As, Cd, Cr, Cu, Hg, Ni, Pb, Se, Sn and Zn) are associated with two known proxies of health: body condition and the gut bacterial microbiome. We found evidence of an adverse effect of Se, Hg, and Pb bioaccumulation on body condition. Furthermore, bioaccumulation of the elements As, Cu, Se, Pb and Zn influenced different aspects of the gut microbiome. Bioaccumulation of Se led to a shift in the microbiome composition, largely driven by an enrichment of Bacteroides plebeius, which is linked to the breakdown of sulphated polysaccharides of algae. Bacteroides plebeius was negatively associated with chick body condition, suggesting an adverse effect of a microalgae diet rich in Se. Pb bioaccumulation was linked with a decrease in microbial diversity (adjusted-R² = 10.4%) and an increase in heterogeneity of the microbial community (adjusted-R² = 10.5%), an indication of impaired gut homeostasis. As, Cu and Zn had more nuanced effects on gut microbiome heterogeneity according to breeding site and bioaccumulation concentration. Our results therefore suggest that in addition to well-studied elements, bioaccumulation of poorly studied elements also adversely affect health of natural populations.

A novel workflow to improve genotyping of multigene families in wildlife species: An experimental set-up with a known model system

Genotyping complex multigene families in novel systems is particularly challenging. Target primers frequently amplify simultaneously multiple loci leading to high PCR and sequencing artefacts such as chimeras and allele amplification bias. Most genotyping pipelines have been validated in nonmodel systems whereby the real genotype is unknown and the generation of artefacts may be highly repeatable. Further hindering accurate genotyping, the relationship between artefacts and genotype complexity (i.e. number of alleles per genotype) within a PCR remains poorly described. Here, we investigated the latter by experimentally combining multiple known major histocompatibility complex (MHC) haplotypes of a model organism (chicken, Gallus gallus, 43 artificial genotypes with 2-13 alleles per amplicon). In addition to well-defined 'optimal' primers, we simulated a nonmodel species situation by designing 'cross-species' primers based on sequence data from closely related Galliform species. We applied a novel open-source genotyping pipeline (ACACIA; https://gitlab.com/psc_santos/ACACIA), and compared its performance with another, previously published pipeline (AmpliSAS). Allele calling accuracy was higher when using ACACIA (98.5% versus 97% and 77.8% versus 75% for the 'optimal' and 'cross-species' data sets, respectively). Systematic allele dropout of three alleles owing to primer mismatch in the 'cross-species' data set explained high allele calling repeatability (100% when using ACACIA) despite low accuracy, demonstrating that repeatability can be misleading when evaluating genotyping workflows. Genotype complexity was positively associated with nonchimeric artefacts, chimeric artefacts (nonlinearly by levelling when amplifying more than 4-6 alleles) and allele amplification bias. Our study exemplifies and demonstrates pitfalls researchers should avoid to reliably genotype complex multigene families.

Host-parasite interactions in non-native invasive species are dependent on the levels of standing genetic variation at the immune locus

Background

Parasites may mediate the success of biological invasions through their effect on host fitness and thus, on host population growth and stability. However, a release from the pressure of parasites is strongly related to the genetic differentiation of the host. In invasive host populations, the number of available genetic variants, allowing them to ‘fight’ the infection, are likely to be influenced by founder events and genetic drift. The level standing genetic variation of invasive populations may be crucial in successfully adapting to new environments and resisting diseases. We studied invasive populations of raccoon that experienced a random reduction in genetic diversity during the establishment and evaluated the relationship between host immune genetic diversity and intestinal parasites infection.

Results

We distinguished two different genetic clusters that are characterized by different sets of functionally relevant MHC-DRB alleles. Both clusters were characterized by considerably different allele-parasite associations and different levels of parasite infection. The specific resistance MHC-DRB alleles explained the lower prevalence of Digenea parasites. An increased infection intensity was related to the presence of two MHC-DRB alleles. One of these alleles significantly decreased in frequency over time, causing a decrease of Digenea abundance in raccoons in consecutive years.

Conclusions

Our findings suggest that intestinal parasites can exert selective pressure on an invasive host with lowered levels of immune genetic diversity and contribute to promoting local adaptation over time. The random genetic drift that created the two different genetic clusters in the invasive raccoon range imposed completely different MHC-parasite associations, strongly associated with the infection status of populations. Our findings underline the role of standing genetic variation in shaping host-parasite relationships and provide empirical support that functional genetic variation may be, at least partly, responsible for differences in the success of invasive populations.

Can extreme MHC class I diversity be a feature of a wide geographic range? The example of Seba's short-tailed bat (Carollia perspicillata)

The major histocompatibility complex (MHC) is one of the most diverse genetic regions under pathogen-driven selection because of its central role in antigen binding and immunity. The highest MHC variability, both in terms of the number of individual alleles and gene copies, has so far been found in passerine birds; this is probably attributable to passerine adaptation to both a wide geographic range and a diverse array of habitats. If extraordinary high MHC variation and duplication rates are adaptive features under selection during the evolution of ecologically and taxonomically diverse species, then similarly diverse MHC architectures should be found in bats. Bats are an extremely species-rich mammalian group that is globally widely distributed. Many bat species roost in multitudinous groups and have high contact rates with pathogens, conspecifics, and allospecifics. We have characterized the MHC class I diversity in 116 Panamanian Seba's short-tailed bats (Carollia perspicillata), a widely distributed, generalist, neotropical species. We have detected a remarkable individual and population-level diversity of MHC class I genes, with between seven and 22 alleles and a unique genotype in each individual. This diversity is comparable with that reported in passerine birds and, in both taxonomic groups, further variability has evolved through length polymorphisms. Our findings support the hypothesis that, for species with a geographically broader range, high MHC class I variability is particularly adaptive. Investigation of the details of the underlying adaptive processes and the role of the high MHC diversity in pathogen resistance are important next steps for a better understanding of the role of bats in viral evolution and as carriers of several deadly zoonotic viruses.

Offspring Microbiomes Differ Across Breeding Sites in a Panmictic Species

High dispersal rates are known to homogenize host’s population genetic structure in panmictic species and to disrupt host local adaptation to the environment. Long-distance dispersal might also spread micro-organisms across large geographical areas. However, so far, to which extent selection mechanisms that shape host’s population genetics are mirrored in the population structure of the enteric microbiome remains unclear. High dispersal rates and horizontal parental transfer may homogenize bacterial communities between breeding sites (homogeneous hypothesis). Alternatively, strong selection from the local environment may differentiate bacterial communities between breeding sites (heterogeneous hypothesis). Furthermore, selection from age-specific environmental or physiological factors may differentiate the microbiome between juveniles and adults. Here, we analyzed the cloacal bacterial 16S rRNA gene of fledgling greater flamingos, Phoenicopterus roseus, across nine western Mediterranean breeding sites and four breeding seasons (n = 731) and adult birds (n = 27) from a single site. We found that fledgling cloacal microbiome, as measured by alpha diversity, beta diversity, the relative abundance of assigned sequence variants (ASVs) belonging to a phylum and genus composition within phylum, varied significantly between sampling sites and across time within site despite high adult dispersal rates. The spatio-temporal effects were stronger on individual ASV absence/presence than on ASV abundance (i.e., than on core microbiome composition). Spatial effects had a stronger effect than temporal effects, particularly on ASV abundance. Our study supports the heterogeneous hypothesis whereby local environmental conditions select and differentiate bacterial communities, thus countering the homogenizing effects of high-dispersing host species. In addition, differences in core microbiome between adult vs. fledgling samples suggests that differences in age-specific environmental and/or physiological factors result in differential selection pressure of core enteric microbiome between age classes, even within the same environment. In particular, the genus Corynebacterium, associated with both seasonal fat uptake and migration in previous studies, was much more abundant in high-dispersing fledglings than in more resident adults. To conclude, selection mechanisms that shape the host’s genetic structure cannot be extended to the genetic structure of the enteric microbiome, which has important implications regarding our understanding of both host local adaptation mechanisms and enteric microbiome population genetics.

Positive selection and climatic effects on MHC class II gene diversity in hares (Lepus capensis) from a steep ecological gradient

In natural populations, allelic diversity of the major histocompatibility complex (MHC) is commonly interpreted as resulting from positive selection in varying spatiotemporal pathogenic landscapes. Composite pathogenic landscape data are, however, rarely available. We studied the spatial distribution of allelic diversity at two MHC class II loci (DQA, DQB) in hares, Lepus capensis, along a steep ecological gradient in North Africa and tested the role of climatic parameters for the spatial distribution of DQA and DQB proteins. Climatic parameters were considered to reflect to some extent pathogenic landscape variation. We investigated historical and contemporary forces that have shaped the variability at both genes, and tested for differential selective pressure across the ecological gradient by comparing allelic variation at MHC and neutral loci. We found positive selection on both MHC loci and significantly decreasing diversity from North to South Tunisia. Our multinomial linear models revealed significant effects of geographical positions that were correlated with mean annual temperature and precipitation on the occurrence of protein variants, but no effects of co-occurring DQA or DQB proteins, respectively. Diversifying selection, recombination, adaptation to local pathogenic landscapes (supposedly reflected by climate parameters) and neutral demographic processes have shaped the observed MHC diversity and differentiation patterns.

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.