Patterns of reticulate evolution for the classical class I and II HLA loci

Evolution of major histocompatibility complex class I genes in the sable Martes zibellina (Carnivora, Mustelidae)

The molecules encoded by major histocompatibility complex (MHC) genes play an essential role in the adaptive immune response among vertebrates. We investigated the molecular evolution of MHC class I genes in the sable Martes zibellina. We isolated 26 MHC class I sequences, including 12 putatively functional sequences and 14 pseudogene sequences, from 24 individuals from two geographic areas of northeast China. The number of putatively functional sequences found in a single individual ranged from one to five, which might be at least 1-3 loci. We found that both balancing selection and recombination contribute to evolution of MHC class I genes in M. zibellina. In addition, we identified a candidate nonclassical MHC class I lineage in Carnivora, which may have preceded the divergence (about 52-57 Mya) of Caniformia and Feliformia. This may contribute to further understanding of the origin and evolution of nonclassical MHC class I genes. Our study provides important immune information of MHC for M. zibellina, as well as other carnivores.

The MHC Class Ia Genes in Chenfu's Treefrog (Zhangixalus chenfui) Evolved via Gene Duplication, Recombination, and Selection

Simple Summary

Amphibians, the first terrestrial vertebrates, provide materials for adaptive evolutionary studies, such as the evolution of the major histocompatibility complex (MHC). To date, various MHC evolutionary mechanisms have been identified in frogs, but more research is needed to determine the evolutionary mechanisms of the frog MHC. The main purpose of this study was to evaluate polymorphisms in the MHC class Ia genes of the Chenfu’s Treefrog. The MHC class Ia genes of the Chenfu’s Treefrog have high polymorphism. The mechanisms responsible for the formation of the polymorphisms include gene duplication, recombination, and selection.

Abstract

The molecular mechanisms underlying the evolution of adaptive immunity-related proteins can be deduced by a thorough examination of the major histocompatibility complex (MHC). Currently, in vertebrates, there is a relatively large amount of research on MHCs in mammals and birds. However, research related to amphibian MHC genes and knowledge about the evolutionary patterns is limited. This study aimed to isolate the MHC class I genes from Chenfu’s Treefrog (Zhangixalus chenfui) and reveal the underlying evolutionary processes. A total of 23 alleles spanning the coding region of MHC class Ia genes were identified in 13 individual samples. Multiple approaches were used to test and identify recombination from the 23 alleles. Amphibian MHC class Ia alleles, from NCBI, were used to construct the phylogenetic relationships in MEGA. Additionally, the partition strategy was adopted to construct phylogenetic relationships using MrBayes and MEGA. The sites of positive selection were identified by FEL, PAML, and MEME. In Chenfu’s Treefrog, we found that: (1) recombination usually takes place between whole exons of MHC class Ia genes; (2) there are at least 3 loci for MHC class Ia, and (3) the diversity of genes in MHC class Ia can be attributed to recombination, gene duplication, and positive selection. We characterized the evolutionary mechanisms underlying MHC class Ia genes in Chenfu’s Treefrog, and in so doing, broadened the knowledge of amphibian MHC systems.

Discovery of gorilla MHC-C expressing C1 ligand for KIR

In comparison to humans and chimpanzees, gorillas show low diversity at MHC class I genes (Gogo), as reflected by an overall reduced level of allelic variation as well as the absence of a functionally important sequence motif that interacts with killer cell immunoglobulin-like receptors (KIR). Here, we use recently generated large-scale genomic sequence data for a reassessment of allelic diversity at Gogo-C, the gorilla orthologue of HLA-C. Through the combination of long-range amplifications and long-read sequencing technology, we obtained, among the 35 gorillas reanalyzed, three novel full-length genomic sequences including a coding region sequence that has not been previously described. The newly identified Gogo-C*03:01 allele has a divergent recombinant structure that sets it apart from other Gogo-C alleles. Domain-by-domain phylogenetic analysis shows that Gogo-C*03:01 has segments in common with Gogo-B*07, the additional B-like gene that is present on some gorilla MHC haplotypes. Identified in ~ 50% of the gorillas analyzed, the Gogo-C*03:01 allele exclusively encodes the C1 epitope among Gogo-C allotypes, indicating its important function in controlling natural killer cell (NK cell) responses via KIR. We further explored the hypothesis whether gorillas experienced a selective sweep which may have resulted in a general reduction of the gorilla MHC class I repertoire. Our results provide little support for a selective sweep but rather suggest that the overall low Gogo class I diversity can be best explained by drastic demographic changes gorillas experienced in the ancient and recent past.

Characterization of the peripheral blood transcriptome and adaptive evolution of the MHC I and TLR gene families in the wolf (Canis lupus)

Background

The wolf (Canis lupus) is one of the most widely distributed terrestrial mammals, because it is well adapted to various ecological niches and their corresponding pathogen environments. Immunological competence is a crucial factor involved in adapting to a changing environment and fighting pathogen infection in animals. In this study, the peripheral blood transcriptome of wolves was generated via RNA-seq to advance understanding of the wolf immunome, with a special focus on the major histocompatibility complex class I (MHC I) and toll-like receptor (TLR) gene families, which are involved in pathogen recognition and defense.

Results

The blood transcriptomic libraries of eight wolves originating from Tibet and Inner Mongolia were sequenced, and approximately 383 million reads were generated. Using a genome-guided assembly strategy, we obtained 123,851 unigenes, with a mean length of 845 bp and an N50 length of 1121 bp. On the basis of BLAST searches against the NCBI non-redundant protein database (Nr), a total of 36,192 (29.22%) unigenes were annotated. For functional classification, 24,663 unigenes were assigned to 13,016 Gene Ontology (GO) terms belonging to 51 sub-categories of the three main GO categories. Additionally, 7682 unigenes were classified into 6 Kyoto Encyclopedia of Genes and Genomes (KEGG) categories, in which the most represented functional sub-categories were signal transduction and the immune system, and 16,238 unigenes were functionally classified into 25 Eukaryotic Orthologous Groups (KOG) categories. We observed an overall higher ω (d_N/d_S) value at antigen-binding sites (ABSs) than at non-ABS regions as well as clear evidence of intergenic/intragenic recombination events at wolf MHC I loci. Additionally, our analysis revealed that carnivorous TLRs were dominated by purifying selection, with mean ω values at each TLR locus ranging from 0.173 to 0.527. However, we also found significant instances of positive selection that acted on several codons in pathogen recognition domains and were linked to species-specific differences in pathogen recognition.

Conclusions

This study represents the first attempt to characterize the blood transcriptome of the wolf and to highlight the value of investigating the immune system. Balancing selection and recombination have contributed to the historical evolution of wolf MHC I genes. Moreover, TLRs in carnivores have undergone adaptive evolution against the background of purifying selection, and a high level of adaptive evolution was detected in the wolf TLR system.

Electronic supplementary material

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

Gorilla MHC class I gene and sequence variation in a comparative context

Comparisons of MHC gene content and diversity among closely related species can provide insights into the evolutionary mechanisms shaping immune system variation. After chimpanzees and bonobos, gorillas are humans’ closest living relatives; but in contrast, relatively little is known about the structure and variation of gorilla MHC class I genes (Gogo). Here, we combined long-range amplifications and long-read sequencing technology to analyze full-length MHC class I genes in 35 gorillas. We obtained 50 full-length genomic sequences corresponding to 15 Gogo-A alleles, 4 Gogo-Oko alleles, 21 Gogo-B alleles, and 10 Gogo-C alleles including 19 novel coding region sequences. We identified two previously undetected MHC class I genes related to Gogo-A and Gogo-B, respectively, thereby illustrating the potential of this approach for efficient and highly accurate MHC genotyping. Consistent with their phylogenetic position within the hominid family, individual gorilla MHC haplotypes share characteristics with humans and chimpanzees as well as orangutans suggesting a complex history of the MHC class I genes in humans and the great apes. However, the overall MHC class I diversity appears to be low further supporting the hypothesis that gorillas might have experienced a reduction of their MHC repertoire.

Evolution by selection, recombination, and gene duplication in MHC class I genes of two Rhacophoridae species

BACKGROUND

Comparison of major histocompatibility complex (MHC) genes across vertebrate species can reveal molecular mechanisms underlying the evolution of adaptive immunity-related proteins. As the first terrestrial tetrapods, amphibians deserve special attention because of their exposure to probably increased spectrum of microorganisms compared with ancestral aquatic fishes. Knowledge regarding the evolutionary patterns and mechanisms associated with amphibian MHC genes remains limited. The goal of the present study was to isolate MHC class I genes from two Rhacophoridae species (Rhacophorus omeimontis and Polypedates megacephalus) and examine their evolution.

RESULTS

We identified 27 MHC class I alleles spanning the region from exon 2 to 4 in 38 tree frogs. The available evidence suggests that these 27 sequences all belong to classical MHC class I (MHC Ia) genes. Although several anuran species only display one MHC class Ia locus, at least two or three loci were observed in P. megacephalus and R. omeimontis, indicating that the number of MHC class Ia loci varies among anuran species. Recombination events, which mainly involve the entire exons, played an important role in shaping the genetic diversity of the 27 MHC class Ia alleles. In addition, signals of positive selection were found in Rhacophoridae MHC class Ia genes. Amino acid sites strongly suggested by program to be under positive selection basically accorded with the putative antigen binding sites deduced from crystal structure of human HLA. Phylogenetic relationships among MHC class I alleles revealed the presence of trans-species polymorphisms.

CONCLUSIONS

In the two Rhacophoridae species (1) there are two or three MHC class Ia loci; (2) recombination mainly occurs between the entire exons of MHC class Ia genes; (3) balancing selection, gene duplication and recombination all contribute to the diversity of MHC class Ia genes. These findings broaden our knowledge on the evolution of amphibian MHC systems.

Evolution of MHC class I genes in the European badger (Meles meles)

The major histocompatibility complex (MHC) plays a central role in the adaptive immune system and provides a good model with which to understand the evolutionary processes underlying functional genes. Trans-species polymorphism and orthology are both commonly found in MHC genes; however, mammalian MHC class I genes tend to cluster by species. Concerted evolution has the potential to homogenize different loci, whereas birth-and-death evolution can lead to the loss of orthologs; both processes result in monophyletic groups within species. Studies investigating the evolution of MHC class I genes have been biased toward a few particular taxa and model species. We present the first study of MHC class I genes in a species from the superfamily Musteloidea. The European badger (Meles meles) exhibits moderate variation in MHC class I sequences when compared to other carnivores. We identified seven putatively functional sequences and nine pseudogenes from genomic (gDNA) and complementary (cDNA) DNA, signifying at least two functional class I loci. We found evidence for separate evolutionary histories of the α1 and α2/α3 domains. In the α1 domain, several sequences from different species were more closely related to each other than to sequences from the same species, resembling orthology or trans-species polymorphism. Balancing selection and probable recombination maintain genetic diversity in the α1 domain, evidenced by the detection of positive selection and a recombination event. By comparison, two recombination breakpoints indicate that the α2/α3 domains have most likely undergone concerted evolution, where recombination has homogenized the α2/α3 domains between genes, leading to species-specific clusters of sequences. Our findings highlight the importance of analyzing MHC domains separately.

HLA DNA sequence variation among human populations: molecular signatures of demographic and selective events

Molecular differences between HLA alleles vary up to 57 nucleotides within the peptide binding coding region of human Major Histocompatibility Complex (MHC) genes, but it is still unclear whether this variation results from a stochastic process or from selective constraints related to functional differences among HLA molecules. Although HLA alleles are generally treated as equidistant molecular units in population genetic studies, DNA sequence diversity among populations is also crucial to interpret the observed HLA polymorphism. In this study, we used a large dataset of 2,062 DNA sequences defined for the different HLA alleles to analyze nucleotide diversity of seven HLA genes in 23,500 individuals of about 200 populations spread worldwide. We first analyzed the HLA molecular structure and diversity of these populations in relation to geographic variation and we further investigated possible departures from selective neutrality through Tajima's tests and mismatch distributions. All results were compared to those obtained by classical approaches applied to HLA allele frequencies.Our study shows that the global patterns of HLA nucleotide diversity among populations are significantly correlated to geography, although in some specific cases the molecular information reveals unexpected genetic relationships. At all loci except HLA-DPB1, populations have accumulated a high proportion of very divergent alleles, suggesting an advantage of heterozygotes expressing molecularly distant HLA molecules (asymmetric overdominant selection model). However, both different intensities of selection and unequal levels of gene conversion may explain the heterogeneous mismatch distributions observed among the loci. Also, distinctive patterns of sequence divergence observed at the HLA-DPB1 locus suggest current neutrality but old selective pressures on this gene. We conclude that HLA DNA sequences advantageously complement HLA allele frequencies as a source of data used to explore the genetic history of human populations, and that their analysis allows a more thorough investigation of human MHC molecular evolution.

Evidence for widespread reticulate evolution within human duplicons

Approximately 5% of the human genome consists of segmental duplications that can cause genomic mutations and may play a role in gene innovation. Reticulate evolutionary processes, such as unequal crossing-over and gene conversion, are known to occur within specific duplicon families, but the broader contribution of these processes to the evolution of human duplications remains poorly characterized. Here, we use phylogenetic profiling to analyze multiple alignments of 24 human duplicon families that span >8 Mb of DNA. Our results indicate that none of them are evolving independently, with all alignments showing sharp discontinuities in phylogenetic signal consistent with reticulation. To analyze these results in more detail, we have developed a quartet method that estimates the relative contribution of nucleotide substitution and reticulate processes to sequence evolution. Our data indicate that most of the duplications show a highly significant excess of sites consistent with reticulate evolution, compared with the number expected by nucleotide substitution alone, with 15 of 30 alignments showing a >20-fold excess over that expected. Using permutation tests, we also show that at least 5% of the total sequence shares 100% sequence identity because of reticulation, a figure that includes 74 independent tracts of perfect identity >2 kb in length. Furthermore, analysis of a subset of alignments indicates that the density of reticulation events is as high as 1 every 4 kb. These results indicate that phylogenetic relationships within recently duplicated human DNA can be rapidly disrupted by reticulate evolution. This finding has important implications for efforts to finish the human genome sequence, complicates comparative sequence analysis of duplicon families, and could profoundly influence the tempo of gene-family evolution.

Molecular characterization of major histocompatibility complex class II alleles in wild tiger salamanders (Ambystoma tigrinum)

Major histocompatibility complex (MHC) class II genes are usually among the most polymorphic in vertebrate genomes because of their critical role (antigen presentation) in immune response. Prior to this study, the MHC was poorly characterized in tiger salamanders (Ambystoma tigrinum), but the congeneric axolotl (Ambystoma mexicanum) is thought to have an unusual MHC. Most notably, axolotl class II genes lack allelic variation and possess a splice variant without a full peptide binding region (PBR). The axolotl is considered immunodeficient, but it is unclear how or to what extent MHC genetics and immunodeficiency are interrelated. To study the evolution of MHC genes in urodele amphibians, we describe for the first time an expressed polymorphic class II gene in wild tiger salamanders. We sequenced the PBR of a class II gene from wild A. tigrinum (n=33) and identified nine distinct alleles. Observed heterozygosity was 73%, and there were a total of 46 polymorphic sites, most of which correspond to amino acid positions that bind peptides. Patterns of nucleotide substitutions exhibit the signature of diversifying selection, but no recombination was detected. Not surprisingly, trans-species evolution of tiger salamander and axolotl class II alleles was apparent. We have no direct data on the immunodeficiency of tiger salamanders, but the levels of polymorphism in our study population should suffice to bind a variety of foreign peptides (unlike axolotls). Our tiger salamander data suggest that the monomorphism and immunodeficiencies associated with axolotl class II genes is a relic of their unique historical demography, not their phylogenetic legacy.

Evolution by recombination and transspecies polymorphism in the MHC class I gene of Xenopus laevis

The patterns of major histocompatibility complex (MHC) evolution involve duplications, deletions, and independent divergence of loci during episodes punctuated by natural selection. Major differences in MHC evolution among taxa have previously been attributed to variation in linkage patterns of class I and class II MHC genes. Here we characterize patterns of evolution in the MHC class Ia gene of Xenopus laevis in terms of polymorphism, recombination, and extent of transspecies polymorphism. We also compare these patterns to see if a correlation exists with linkage or separation of the MHC class I and class II regions as seen in amphibians and teleost fishes. In X. laevis, we find high levels of polymorphism. Also, genetic exchange is relatively frequent and occurs in intron II, reshuffling allelic forms of exons 2 and 3. Evolutionary relationships among class I alleles show an intermingling of alleles from divergent Xenopus species rather than a species-specific clustering. Results indicate that the patterns of evolution are similar to those found in salmonid fishes and are different from the mode of evolution seen in primates. Similar patterns of class Ia evolution in salmonid fishes and X. laevis suggest that nonlinkage of class I and class II regions alone is insufficient to explain some patterns of MHC evolution in salmonids.

Geographic patterns of functional categories of HLA-DRB1 alleles: a new approach to analyse associations between HLA-DRB1 and disease

Because specific amino acids found within the peptide-binding cleft of human leukocyte antigen (HLA) molecules have been implicated in HLA/disease associations, an approach which consists in grouping the alleles according to their functional properties at the protein level may enable us to better understand HLA associations than the conventional allelic classification. In this study, we applied this methodology to investigate the associations between HLA-DRB1 and rheumatoid arthritis. The alleles were first classified into seven functional categories [restrictive supertype patterns (RSPs)], among which three were known to be significantly associated with susceptibility (one category) or resistance (two categories) to rheumatoid arthritis. The frequencies of these categories were then estimated in 104 population samples previously tested for HLA-DRB1, and their variability was analysed spatially on a worldwide scale by applying an original methodology for detecting discontinuities in geographically patterned data. RSP frequencies were also compared to known values of rheumatoid arthritis prevalence in some populations. The results indicated that the three RSP frequency distributions were geographically structured, and that these patterns could generally be explained by the history of human migrations. However, the peculiar pattern observed for RSP 'A' (conferring susceptibility to rheumatoid arthritis) indicated a possible association with some latitude-dependent disease. Furthermore, the very high correlation coefficient found between RSP 'A' frequencies and rheumatoid arthritis prevalence confirmed the significant disease association of this functional category. In contrast, the putative protective effect of the other RSPs ('De' and 'Q') was not detectable at the worldwide level, but may be significant in specific geographic areas. This study shows that population genetic diversity analyses based on a functional grouping of HLA alleles provide an efficient way to explore the mutual influence of HLA genetic variation and disease.

Identification of a novel DPB1 allele, DPB1*9301, by sequence-based typing in a Lahu ethnic minority of China

A novel DPB1 allele (DPB1*9301) has been identified in an individual from the Lahu ethnic minority of China. Its sequence was confirmed by sequencing the PCR products and clones. This allele differs from DPB1*0601, *1701 and *3701 by two nucleotides in each, and results in amino acids substitution.

Strong positive selection and habitat-specific amino acid substitution patterns in MHC from an estuarine fish under intense pollution stress

Population-level studies using the major histocompatibility complex (Mhc) have linked specific alleles with specific diseases, but data requirements are high and the power to detect disease association is low. A novel use of Mhc population surveys involves mapping allelic substitutions onto the inferred structural molecular model to show functional differentiation related to local selective pressures. In the estuarine fish Fundulus heteroclitus, populations experiencing strong differences in antigenic challenges show significant differences in amino acid substitution patterns that are reflected as variation in the structural location of changes between populations. Fish from a population genetically adapted to severe chemical pollution also show novel patterns of DNA substitution at a highly variable Mhc class II B locus including strong signals of positive selection at inferred antigen-binding sites and population-specific signatures of amino acid substitution. Heavily parasitized fish from an extreme PCB-contaminated (U.S. Environmental Protection Agency Superfund) site show enhanced population-specific substitutions in the a-helix portion of the inferred antigen-binding region. In contrast, fish from an unpolluted site show a significantly different pattern focused on the first strand of the B-pleated sheet. Whether Mhc population profile differences represent the direct effects of chemical toxicants or indirect parasite-mediated selection, the result is a composite habitat-specific signature of strong selection and evolution affecting the genetic repertoire of the major histocompatibility complex.

Modes of salmonid MHC class I and II evolution differ from the primate paradigm

Rainbow trout (Oncorhynchus mykiss) and brown trout (Salmo trutta) represent two salmonid genera separated for 15--20 million years. cDNA sequences were determined for the classical MHC class I heavy chain gene UBA and the MHC class II beta-chain gene DAB from 15 rainbow and 10 brown trout. Both genes are highly polymorphic in both species and diploid in expression. The MHC class I alleles comprise several highly divergent lineages that are represented in both species and predate genera separation. The class II alleles are less divergent, highly species specific, and probably arose after genera separation. The striking difference in salmonid MHC class I and class II evolution contrasts with the situation in primates, where lineages of class II alleles have been sustained over longer periods of time relative to class I lineages. The difference may arise because salmonid MHC class I and II genes are not linked, whereas in mammals they are closely linked. A prevalent mechanism for evolving new MHC class I alleles in salmonids is recombination in intron II that shuffles alpha 1 and alpha 2 domains into different combinations.

Evolution of HLA class II molecules: Allelic and amino acid site variability across populations

Analysis of the highly polymorphic beta1 domains of the HLA class II molecules encoded by the DRB1, DQB1, and DPB1 loci reveals contrasting levels of diversity at the allele and amino acid site levels. Statistics of allele frequency distributions, based on Watterson's homozygosity statistic F, reveal distinct evolutionary patterns for these loci in ethnically diverse samples (26 populations for DQB1 and DRB1 and 14 for DPB1). When examined over all populations, the DQB1 locus allelic variation exhibits striking balanced polymorphism (P < 10(-4)), DRB1 shows some evidence of balancing selection (P < 0.06), and while there is overall very little evidence for selection of DPB1 allele frequencies, there is a trend in the direction of balancing selection (P < 0.08). In contrast, at the amino acid level all three loci show strong evidence of balancing selection at some sites. Averaged over polymorphic amino acid sites, DQB1 and DPB1 show similar deviation from neutrality expectations, and both exhibit more balanced polymorphic amino acid sites than DRB1. Across ethnic groups, polymorphisms at many codons show evidence for balancing selection, yet data consistent with directional selection were observed at other codons. Both antigen-binding pocket- and non-pocket-forming amino acid sites show overall deviation from neutrality for all three loci. Only in the case of DRB1 was there a significant difference between pocket- and non-pocket-forming amino acid sites. Our findings indicate that balancing selection at the MHC occurs at the level of polymorphic amino acid residues, and that in many cases this selection is consistent across populations.