Balancing and directional selection at exon-2 of the MHC DQB1 locus among populations of odontocete cetaceans

Similar patterns of genetic diversity and linkage disequilibrium in Western chimpanzees (Pan troglodytes verus) and humans indicate highly conserved mechanisms of MHC molecular evolution

BACKGROUND

Many species are threatened with extinction as their population sizes decrease with changing environments or face novel pathogenic threats. A reduction of genetic diversity at major histocompatibility complex (MHC) genes may have dramatic effects on populations' survival, as these genes play a key role in adaptive immunity. This might be the case for chimpanzees, the MHC genes of which reveal signatures of an ancient selective sweep likely due to a viral epidemic that reduced their population size a few million years ago. To better assess how this past event affected MHC variation in chimpanzees compared to humans, we analysed several indexes of genetic diversity and linkage disequilibrium across seven MHC genes on four cohorts of chimpanzees and we compared them to those estimated at orthologous HLA genes in a large set of human populations.

RESULTS

Interestingly, the analyses uncovered similar patterns of both molecular diversity and linkage disequilibrium across the seven MHC genes in chimpanzees and humans. Indeed, in both species the greatest allelic richness and heterozygosity were found at loci A, B, C and DRB1, the greatest nucleotide diversity at loci DRB1, DQA1 and DQB1, and both significant global linkage disequilibrium and the greatest proportions of haplotypes in linkage disequilibrium were observed at pairs DQA1 ~ DQB1, DQA1 ~ DRB1, DQB1 ~ DRB1 and B ~ C. Our results also showed that, despite some differences among loci, the levels of genetic diversity and linkage disequilibrium observed in contemporary chimpanzees were globally similar to those estimated in small isolated human populations, in contrast to significant differences compared to large populations.

CONCLUSIONS

We conclude, first, that highly conserved mechanisms shaped the diversity of orthologous MHC genes in chimpanzees and humans. Furthermore, our findings support the hypothesis that an ancient demographic decline affecting the chimpanzee populations - like that ascribed to a viral epidemic - exerted a substantial effect on the molecular diversity of their MHC genes, albeit not more pronounced than that experienced by HLA genes in human populations that underwent rapid genetic drift during humans' peopling history. We thus propose a model where chimpanzees' MHC genes regenerated molecular variation through recombination/gene conversion and/or balancing selection after the selective sweep.

Sex-specific impact of inbreeding on pathogen load in the striped dolphin

The impact of inbreeding on fitness has been widely studied and provides consequential inference about adaptive potential and the impact on survival for reduced and fragmented natural populations. Correlations between heterozygosity and fitness are common in the literature, but they rarely inform about the likely mechanisms. Here, we investigate a pathology with a clear impact on health in striped dolphin hosts (a nematode infection that compromises lung function). Dolphins varied with respect to their parasite burden of this highly pathogenic lung nematode (Skrjabinalius guevarai). Genetic diversity revealed by high-resolution restriction-associated DNA (43 018 RADseq single nucleotide polymorphisms) analyses showed a clear association between heterozygosity and pathogen load, but only for female dolphins, for which the more heterozygous individuals had lower Sk. guevarai burden. One locus identified by RADseq was a strong outlier in association with parasite load (heterozygous in all uninfected females, homozygous for 94% of infected females), found in an intron of the citron rho-interacting serine/threonine kinase locus (associated with milk production in mammals). Allelic variation at the Class II major histocompatability complex DQB locus was also assessed and found to be associated with both regional variation and with pathogen load. Both sex specificity and the identification of associating functional loci provide insight into the mechanisms by which more inbred individuals may be more susceptible to the infection of this parasite. This provides important insight towards our understanding of the impact of inbreeding in natural populations, relevant to both evolutionary and practical conservation considerations.

Is MHC diversity a better marker for conservation than neutral genetic diversity? A case study of two contrasting dolphin populations

Genetic diversity is essential for populations to adapt to changing environments. Measures of genetic diversity are often based on selectively neutral markers, such as microsatellites. Genetic diversity to guide conservation management, however, is better reflected by adaptive markers, including genes of the major histocompatibility complex (MHC). Our aim was to assess MHC and neutral genetic diversity in two contrasting bottlenose dolphin (Tursiops aduncus) populations in Western Australia-one apparently viable population with high reproductive output (Shark Bay) and one with lower reproductive output that was forecast to decline (Bunbury). We assessed genetic variation in the two populations by sequencing the MHC class II DQB, which encompasses the functionally important peptide binding regions (PBR). Neutral genetic diversity was assessed by genotyping twenty-three microsatellite loci. We confirmed that MHC is an adaptive marker in both populations. Overall, the Shark Bay population exhibited greater MHC diversity than the Bunbury population-for example, it displayed greater MHC nucleotide diversity. In contrast, the difference in microsatellite diversity between the two populations was comparatively low. Our findings are consistent with the hypothesis that viable populations typically display greater genetic diversity than less viable populations. The results also suggest that MHC variation is more closely associated with population viability than neutral genetic variation. Although the inferences from our findings are limited, because we only compared two populations, our results add to a growing number of studies that highlight the usefulness of MHC as a potentially suitable genetic marker for animal conservation. The Shark Bay population, which carries greater adaptive genetic diversity than the Bunbury population, is thus likely more robust to natural or human-induced changes to the coastal ecosystem it inhabits.

Positive selection in coding regions and motif duplication in regulatory regions of bottlenose dolphin MHC class II genes

The vertebrate immune response is mediated through highly adaptive, quickly evolving cell surface receptors, the major histocompatibility complex (MHC). MHC molecules bind and present a diverse array of pathogenic molecules and trigger a cascade of defenses. Use of MHC variation as a marker for population health has also evolved quickly following advances in sequencing methods. We applied a combination of traditional and next generation sequencing methodology to characterize coding (peptide binding region) and regulatory (proximal promoter) sequence variation in MHC Class II DQA and DQB genes between estuarine and coastal populations of the bottlenose dolphin, Tursiops truncatus, an apex predator whose health status is indicative of anthropogenic impacts on the ecosystem. The coding regions had 10 alleles each at DQA and DQB; the promoters had 6 and 7 alleles at DQA and DQB, respectively with variation within key regulatory motifs. Positive selection was observed for the coding regions of both genes while both coding and promoter regions exhibited geographic differences in allele composition that likely indicates diversifying selection across habitats. Most notable was the discovery of a complete duplication of a 14-bp T-box motif in the DQA promoter. Four class II promoter regions (DQA, DQB, DRA, DRB) were characterized in species from four cetacean families (Delphinidae, Monodontidae, Lipotidae, and Physeteridae) and revealed substantial promoter structural diversity across this order. Peptide binding regions may not be the only source of adaptive potential within cetacean MHC for responding to pathogenic threats. These findings are the first analysis of cetacean MHC regulatory motifs, which may divulge unique immunogenetic strategies among cetaceans and reveal how MHC transcriptional control continues to evolve. The combined MHC regulatory and coding data provide new genetic context for distinct vulnerability profiles between coastal and estuarine populations, which are key concerns for health and risk management.

Bridging immunogenetics and immunoproteomics: Model positional scanning library analysis for Major Histocompatibility Complex class II DQ in Tursiops truncatus

The Major Histocompatibility Complex (MHC) is a critical element in mounting an effective immune response in vertebrates against invading pathogens. Studies of MHC in wildlife populations have typically focused on assessing diversity within the peptide binding regions (PBR) of the MHC class II (MHC II) family, especially the DQ receptor genes. Such metrics of diversity, however, are of limited use to health risk assessment since functional analyses (where changes in the PBR are correlated to recognition/pathologies of known pathogen proteins), are difficult to conduct in wildlife species. Here we describe a means to predict the binding preferences of MHC proteins: We have developed a model positional scanning library analysis (MPSLA) by harnessing the power of mixture based combinatorial libraries to probe the peptide landscapes of distinct MHC II DQ proteins. The algorithm provided by NNAlign was employed to predict the binding affinities of sets of peptides generated for DQ proteins. These binding affinities were then used to retroactively construct a model Positional Scanning Library screen. To test the utility of the approach, a model screen was compared to physical combinatorial screens for human MHC II DP. Model library screens were generated for DQ proteins derived from sequence data from bottlenose dolphins from the Indian River Lagoon (IRL) and the Atlantic coast of Florida, and compared to screens of DQ proteins from Genbank for dolphin and three other cetaceans. To explore the peptide binding landscape for DQ proteins from the IRL, combinations of the amino acids identified as active were compiled into peptide sequence lists that were used to mine databases for representation in known proteins. The frequency of which peptide sequences predicted to bind the MHC protein are found in proteins from pathogens associated with marine mammals was found to be significant (p values <0.0001). Through this analysis, genetic variation in MHC (classes I and II) can now be associated with the binding repertoires of the expressed MHC proteins and subsequently used to identify target pathogens. This approach may be eventually applied to evaluate individual population and species risk for outbreaks of emerging diseases.

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.

Adaptive tolerance to a pathogenic fungus drives major histocompatibility complex evolution in natural amphibian populations

Amphibians have been affected globally by the disease chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), and we are just now beginning to understand how immunogenetic variability contributes to disease susceptibility. Lineages of an expressed major histocompatibility complex (MHC) class II locus involved in acquired immunity are associated with chytridiomycosis susceptibility in controlled laboratory challenge assays. Here, we extend these findings to natural populations that vary both in exposure and response to Bd. We find that MHC alleles and supertypes associated with Bd survival in the field show a molecular signal of positive selection, while those associated with susceptibility do not, supporting the hypothesis that heritable Bd tolerance is rapidly evolving. We compare MHC supertypes to neutral loci to demonstrate where selection versus demography is shaping MHC variability. One population with Bd tolerance in nature shows a significant signal of directional selection for the same allele (allele Q) that was significantly associated with survival in an earlier laboratory study. Our findings indicate that selective pressure for Bd survival drives rapid immunogenetic adaptation in some natural populations, despite differences in environment and demography. Our field-based analysis of immunogenetic variation confirms that natural amphibian populations have the evolutionary potential to adapt to chytridiomycosis.

Low Major Histocompatibility Complex Class II Variation in the Endangered Indo-Pacific Humpback Dolphin (Sousa chinensis): Inferences About the Role of Balancing Selection

It has been widely reported that the major histocompatibility complex (MHC) is under balancing selection due to its immune function across terrestrial and aquatic mammals. The comprehensive studies at MHC and other neutral loci could give us a synthetic evaluation about the major force determining genetic diversity of species. Previously, a low level of genetic diversity has been reported among the Indo-Pacific humpback dolphin (Sousa chinensis) in the Pearl River Estuary (PRE) using both mitochondrial marker and microsatellite loci. Here, the expression and sequence polymorphism of 2 MHC class II genes (DQB and DRB) in 32 S. chinensis from PRE collected between 2003 and 2011 were investigated. High ratios of non-synonymous to synonymous substitution rates, codon-based selection analysis, and trans-species polymorphism (TSP) support the hypothesis that balancing selection acted on S. chinensis MHC sequences. However, only 2 haplotypes were detected at either DQB or DRB loci. Moreover, the lack of deviation from the Hardy-Weinberg expectation at DRB locus combined with the relatively low heterozygosity at both DQB locus and microsatellite loci suggested that balancing selection might not be sufficient, which further suggested that genetic drift associated with historical bottlenecks was not mitigated by balancing selection in terms of the loss of MHC and neutral variation in S. chinensis. The combined results highlighted the importance of maintaining the genetic diversity of the endangered S. chinensis.

A review of molecular approaches for investigating patterns of coevolution in marine host-parasite relationships

Parasites and their relationships with hosts play a crucial role in the evolutionary pathways of every living organism. One method of investigating host-parasite systems is using a molecular approach. This is particularly important as analyses based solely on morphology or laboratory studies of parasites and their hosts do not take into account historical evolutionary interactions that can shape the distribution, abundance and population structure of parasites and their hosts. However, the predominant host-parasite coevolution literature has focused on terrestrial hosts and their parasites, and there still is a lack of studies in marine environments. Given that marine systems are generally more open than terrestrial ones, they provide fascinating opportunities for large-scale (as well as small-scale) geographic studies. Further, patterns and processes of genetic structuring and systematics are becoming more available across many different taxa (but especially fishes) in many marine systems, providing an excellent basis for examining whether parasites follow host population/species structure. In this chapter, we first highlight the factors and processes that challenge our ability to interpret evolutionary patterns of coevolution of hosts and their parasites in marine systems at different spatial, temporal and taxonomic scales. We then review the use of the most commonly utilized genetic markers in studying marine host-parasite systems. We give an overview and discuss which molecular methodologies resolve evolutionary relationships best and also discuss the applicability of new approaches, such as next-generation sequencing and studies utilizing functional markers to gain insights into more contemporary processes shaping host-parasite relationships.

Characterization of class I- and class II-like major histocompatibility complex loci in pedigrees of North Atlantic right whales

North Atlantic right whales have one of the lowest levels of genetic variation at minisatellite loci, microsatellite loci, and mitochondrial control region haplotypes among mammals. Here, adaptive variation at the peptide binding region of class I and class II DRB-like genes of the major histocompatibility complex was assessed. Amplification of a duplicated region in 222 individuals revealed at least 11 class II alleles. Six alleles were assigned to the locus Eugl-DRB1 and 5 alleles were assigned to the locus Eugl-DRB2 by assessing segregation patterns of alleles from 81 parent/offspring pedigrees. Pedigree analysis indicated that these alleles segregated into 12 distinct haplotypes. Genotyping a smaller subset of unrelated individuals (n = 5 and 10, respectively) using different primer sets revealed at least 2 class II pseudogenes (with ≥ 4 alleles) and at least 3 class I loci (with ≥ 6 alleles). Class II sequences were significantly different from neutrality at peptide binding sites suggesting loci may be under the influence of balancing selection. Trans-species sharing of alleles was apparent for class I and class II sequences. Characterization of class II loci represents the first step in determining the relationship between major histocompatibility complex variability and factors affecting health and reproduction in this species.

Comparison of 454 pyrosequencing methods for characterizing the major histocompatibility complex of nonmodel species and the advantages of ultra deep coverage

Characterization and population genetic analysis of multilocus genes, such as those found in the major histocompatibility complex (MHC) is challenging in nonmodel vertebrates. The traditional method of extensive cloning and Sanger sequencing is costly and time-intensive and indirect methods of assessment often underestimate total variation. Here, we explored the suitability of 454 pyrosequencing for characterizing multilocus genes for use in population genetic studies. We compared two sample tagging protocols and two bioinformatic procedures for 454 sequencing through characterization of a 185-bp fragment of MHC DRB exon 2 in wolverines (Gulo gulo) and further compared the results with those from cloning and Sanger sequencing. We found 10 putative DRB alleles in the 88 individuals screened with between two and four alleles per individual, suggesting amplification of a duplicated DRB gene. In addition to the putative alleles, all individuals possessed an easily identifiable pseudogene. In our system, sequence variants with a frequency below 6% in an individual sample were usually artefacts. However, we found that sample preparation and data processing procedures can greatly affect variant frequencies in addition to the complexity of the multilocus system. Therefore, we recommend determining a per-amplicon-variant frequency threshold for each unique system. The extremely deep coverage obtained in our study (approximately 5000×) coupled with the semi-quantitative nature of pyrosequencing enabled us to assign all putative alleles to the two DRB loci, which is generally not possible using traditional methods. Our method of obtaining locus-specific MHC genotypes will enhance population genetic analyses and studies on disease susceptibility in nonmodel wildlife species.

Evolution of the MHC-DQB exon 2 in marine and terrestrial mammals

On the basis of a general low polymorphism, several studies suggest that balancing selection in the class II major histocompatibility complex (MHC) is weaker in marine mammals as compared with terrestrial mammals. We investigated such differential selection among Cetacea, Artiodactyla, and Primates at exon 2 of MHC-DQB gene by contrasting indicators of molecular evolution such as occurrence of transpecific polymorphisms, patterns of phylogenetic branch lengths by codon position, rates of nonsynonymous and synonymous substitutions as well as accumulation of variable sites on the sampling of alleles. These indicators were compared between the DQB and the mitochondrial cytochrome b gene (cytb) as a reference of neutral expectations and differences between molecular clocks resulting from life history and historical demography. All indicators showed that the influence of balancing selection on the DQB is more variable and overall weaker for cetaceans. In our sampling, ziphiids, the sperm whale, monodontids and the finless porpoise formed a group with lower DQB polymorphism, while mysticetes exhibited a higher DQB variation similar to that of terrestrial mammals as well as higher occurrence of transpecific polymorphisms. Different dolphins appeared in the two groups. Larger variation of selection on the cetacean DQB could be related to greater stochasticity in their historical demography and thus, to a greater complexity of the general ecology and disease processes of these animals.

Considerable MHC diversity suggests that the functional extinction of baiji is not related to population genetic collapse

To further extend our understanding of the mechanism causing the current nearly extinct status of the baiji (Lipotes vexillifer), one of the most critically endangered species in the world, genetic diversity at the major histocompatibility complex (MHC) class II DRB locus was investigated in the baiji. Nine highly divergent DRB alleles were identified in 17 samples, with an average of 28.4 (13.2%) nucleotide difference and 16.7 (23.5%) amino acid difference between alleles. The unexpectedly high levels of DRB allelic diversity in the baiji may partly be attributable to its evolutionary adaptations to the freshwater environment which is regarded to have a higher parasite diversity compared to the marine environment. In addition, balancing selection was found to be the main mechanisms in generating sequence diversity at baiji DRB gene. Considerable sequence variation at the adaptive MHC genes despite of significant loss of neutral genetic variation in baiji genome might suggest that intense selection has overpowered random genetic drift as the main evolutionary forces, which further suggested that the critically endangered or nearly extinct status of the baiji is not an outcome of genetic collapse.

Complex interplay of evolutionary forces in the ladybird homeobox genes of Drosophila melanogaster

Tandemly arranged paralogous genes lbe and lbl are members of the Drosophila NK homeobox family. We analyzed population samples of Drosophila melanogaster from Africa, Europe, North and South America, and single strains of D. sechellia, D. simulans, and D. yakuba within two linked regions encompassing partial sequences of lbe and lbl. The evolution of lbe and lbl is highly constrained due to their important regulatory functions. Despite this, a variety of forces have shaped the patterns of variation in lb genes: recombination, intragenic gene conversion and natural selection strongly influence background variation created by linkage disequilibrium and dimorphic haplotype structure. The two genes exhibited similar levels of nucleotide diversity and positive selection was detected in the noncoding regions of both genes. However, synonymous variability was significantly higher for lbe: no nonsynonymous changes were observed in this gene. We argue that balancing selection impacts some synonymous sites of the lbe gene. Stability of mRNA secondary structure was significantly different between the lbe (but not lbl) haplotype groups and may represent a driving force of balancing selection in epistatically interacting synonymous sites. Balancing selection on synonymous sites may be the first, or one of a few such observations, in Drosophila. In contrast, recurrent positive selection on lbl at the protein level influenced evolution at three codon sites. Transcription factor binding-site profiles were different for lbe and lbl, suggesting that their developmental functions are not redundant. Combined with our previous results on nucleotide variation in esterase and other homeobox genes, these results suggest that interplay of balancing and directional selection may be a general feature of molecular evolution in Drosophila and other eukaryote genomes.

Sequence polymorphism and geographical variation at a positively selected MHC-DRB gene in the finless porpoise (Neophocaena phocaenoides): implication for recent differentiation of the Yangtze finless porpoise?

Sequence polymorphism at the MHC class II DRB locus was investigated in three finless porpoise (Neophocaena phocaenoides) populations in Chinese waters. Intragenic recombination and strong positive selection were the main forces in generating sequence diversity in the DRB gene. MHC sequence diversity changed significantly along the study period. Significant decrease in heterozygosity and lost alleles have been detected in the Yangtze River population and South China Sea population since 1990. Furthermore, there is a trend of increasing population differentiation over time. Especially, the genetic differentiation between the Yangtze River population and the Yellow Sea population was very low prior to 1990 (F (ST) = 0.036, P = 0.009), but became very significant after 1990 (F (ST) = 0.134, P < 0.001), suggesting a recent augmentation of genetic differentiation between both populations probably in a relatively short-term period. Porpoises from the Yangtze River displayed divergent frequencies of shared and private alleles from those displayed by two marine populations, which suggest that the former riverine population has been under a different selection regime (characteristic of a fresh water environment) than that of its marine counterparts.

Comparative genetic analysis of inflammatory bowel disease and type 1 diabetes implicates multiple loci with opposite effects

Inflammatory bowel disease, including Crohn's disease (CD) and ulcerative colitis (UC), and type 1 diabetes (T1D) are autoimmune diseases that may share common susceptibility pathways. We examined known susceptibility loci for these diseases in a cohort of 1689 CD cases, 777 UC cases, 989 T1D cases and 6197 shared control subjects of European ancestry, who were genotyped by the Illumina HumanHap550 SNP arrays. We identified multiple previously unreported or unconfirmed disease associations, including known CD loci (ICOSLG and TNFSF15) and T1D loci (TNFAIP3) that confer UC risk, known UC loci (HERC2 and IL26) that confer T1D risk and known UC loci (IL10 and CCNY) that confer CD risk. Additionally, we show that T1D risk alleles residing at the PTPN22, IL27, IL18RAP and IL10 loci protect against CD. Furthermore, the strongest risk alleles for T1D within the major histocompatibility complex (MHC) confer strong protection against CD and UC; however, given the multi-allelic nature of the MHC haplotypes, sequencing of the MHC locus will be required to interpret this observation. These results extend our current knowledge on genetic variants that predispose to autoimmunity, and suggest that many loci involved in autoimmunity may be under a balancing selection due to antagonistic pleiotropic effect. Our analysis implies that variants with opposite effects on different diseases may facilitate the maintenance of common susceptibility alleles in human populations, making autoimmune diseases especially amenable to genetic dissection by genome-wide association studies.