Multiple major histocompatibility complex class I genes in Asian anurans: Ontogeny and phylogeny

The Amphibian Major Histocompatibility Complex-A Review and Future Outlook

The major histocompatibility complex (MHC) is a cluster of functionally related genes encoding proteins which, among other functions, mediate immune system activation. While the MHC of many vertebrates has been extensively studied, less is known about the amphibian MHC. This represents an important knowledge gap because amphibians mark the evolutionary transition from an aquatic to a terrestrial lifestyle and often maintain a biphasic lifestyle. Hence, they tend to be exposed to both aquatic and terrestrial pathogen communities, providing opportunities to gain fundamental insights into how the immune system responds to different environmental challenges. Moreover, amphibians are globally threatened by invasive pathogens and the MHC may play a role in combating population decline. In this review, we summarize the current state of knowledge regarding the amphibian MHC and identify the major differences with other vertebrates. We also review how the number of MHC gene copies varies across amphibian groups and how MHC-based variation relates to amphibian ontogeny, behaviour, disease, and phylogeography. We conclude by identifying knowledge gaps and proposing priorities for future research.

Conserved Evolution of MHC Supertypes among Japanese Frogs Suggests Selection for Bd Resistance

The chytrid fungus Batrachochytrium dendrobatidis (Bd) is a major threat to amphibians, yet there are no reports of major disease impacts in East Asian frogs. Genetic variation of the major histocompatibility complex (MHC) has been associated with resistance to Bd in frogs from East Asia and worldwide. Using transcriptomic data collated from 11 Japanese frog species (one individual per species), we isolated MHC class I and IIb sequences and validated using molecular cloning. We then compared MHC from Japanese frogs and other species worldwide, with varying Bd susceptibility. Supertyping analysis, which groups MHC alleles based on physicochemical properties of peptide binding sites, identified that all examined East Asian frogs contained at least one MHC-IIb allele belonging to supertype ST-1. This indicates that, despite the large divergence times between some Japanese frogs (up to 145 million years), particular functional properties in the peptide binding sites of MHC-II are conserved among East Asian frogs. Furthermore, preliminary analysis using NetMHCIIpan-4.0, which predicts potential Bd-peptide binding ability, suggests that MHC-IIb ST-1 and ST-2 have higher overall peptide binding ability than other supertypes, irrespective of whether the peptides are derived from Bd, other fungi, or bacteria. Our findings suggest that MHC-IIb among East Asian frogs may have co-evolved under the same selective pressure. Given that Bd originated in this region, it may be a major driver of MHC evolution in East Asian frogs.

Expression Changes of MHC and Other Immune Genes in Frog Skin during Ontogeny

Simple Summary

Tadpoles undergo many changes in physiology and immunology until metamorphosis into adult frogs. Major histocompatibility complex (MHC) molecules are an important part of vertebrate adaptive immunity, and our study measured the expression of two MHC genes (MHC class I and II) in skin during six tadpole stages of the Montane Brown frog (Rana ornativentris). First, using a qPCR method, we found that both MHC class I and II expression significantly increased between stage 24/25 (‘early’) and stage 28 (‘mid’) tadpole skin. Then, we conducted next-generation sequencing for ‘early’, ‘mid’ and ‘late’ stage tadpole skin mRNA of both R. ornativentris and a model species, Xenopus tropicalis, and confirmed that MHC expression increased from the ‘mid’ stage. We also performed further analyses of transcriptome data and found that several immune-related gene ontology terms were upregulated from the ‘mid’ tadpole stage. Our findings probably support that both MHC class I and II have a functional role during tadpole development.

Abstract

Anuran amphibians undergo major physiological and immunological changes following metamorphosis. Genes of the major histocompatibility complex (MHC) code for receptors important for vertebrate adaptive immunity. We used qPCR to measure skin MHC expression in six different ontological stages of Rana ornativentris (n = 10 per stage); normalized MHC class I and II expression at the mRNA level was significantly higher in stage 28 (mid-larval) compared to stages 24/25 (early-larval) tadpoles. Subsequent transcriptomic analyses of three tadpole (early-, mid-, and late-larval) stages of R. ornativentris and model species Xenopus tropicalis focused on mRNA expression of immune-related genes in the skin. Normalized expression of most MHC class I and II transcripts in both species were significantly higher in mid- and late-larval stages compared to early-larval stage. In addition, gene ontology (GO) analyses of differentially expressed transcripts revealed several immune-related GO terms that were significantly upregulated from the mid-larval stage. Our study provides evidence that both MHC class I and II is expressed during development in both R. ornativentris and X. tropicalis.

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.

Origin and evolution of the major histocompatibility complex class I region in eutherian mammals

Major histocompatibility complex (MHC) genes in vertebrates are vital in defending against pathogenic infections. To gain new insights into the evolution of MHC Class I (MHCI) genes and test competing hypotheses on the origin of the MHCI region in eutherian mammals, we studied available genome assemblies of nine species in Afrotheria, Xenarthra, and Laurasiatheria, and successfully characterized the MHCI region in six species. The following numbers of putatively functional genes were detected: in the elephant, four, one, and eight in the extended class I region, and κ and β duplication blocks, respectively; in the tenrec, one in the κ duplication block; and in the four bat species, one or two in the β duplication block. Our results indicate that MHCI genes in the κ and β duplication blocks may have originated in the common ancestor of eutherian mammals. In the elephant, tenrec, and all four bats, some MHCI genes occurred outside the MHCI region, suggesting that eutherians may have a more complex MHCI genomic organization than previously thought. Bat-specific three- or five-amino-acid insertions were detected in the MHCI α1 domain in all four bats studied, suggesting that pathogen defense in bats relies on MHCIs having a wider peptide-binding groove, as previously assayed by a bat MHCI gene with a three-amino-acid insertion showing a larger peptide repertoire than in other mammals. Our study adds to knowledge on the diversity of eutherian MHCI genes, which may have been shaped in a taxon-specific manner.

Genetic variation and selection of MHC class I loci differ in two congeneric frogs

Major histocompatibility complex (MHC) genes encode proteins in the acquired immune response pathway that often show distinctive selection-driven patterns in wild vertebrate populations. We examined genetic variation and signatures of selection in the MHC class I alpha 1 (A1)- and alpha 2 (A2)-domain encoding exons of two frog congeners [Agalychnis callidryas (n = 20) and A. lemur (n = 20)] from a single locality in Panama. We also investigated how historical demographic processes may have impacted MHC genetic diversity by analyzing a neutral mitochondrial marker. We found that both MHC domains were highly variable in both species, with both species likely expressing three loci. Our analyses revealed different signatures of selection between the two species, most notably that the A. callidryas A2 domain had experienced positive selection while the A2 domain of A. lemur had not. Diversifying selection acted on the same number of A1 and A2 allelic lineages, but on a higher percentage of A1 sites compared to A2 sites. Neutrality tests of mitochondrial haplotypes predominately indicated that the two species were at genetic equilibrium when the samples were collected. In addition, two historical tests of demography indicated both species have had relatively stable population sizes over the past 100,000 years; thus large population size changes are unlikely to have greatly influenced MHC diversity in either species during this time period. In conclusion, our results suggest that the impact of selection on MHC diversity varied between these two closely related species, likely due to a combination of distinct ecological conditions and past pathogenic pressures.

Protective Immunity Induced by DNA Vaccination against Ranavirus Infection in Chinese Giant Salamander Andrias davidianus

Andrias davidianus ranavirus (ADRV) is an emerging viral pathogen that causes severe systemic hemorrhagic disease in Chinese giant salamanders. There is an urgent need for developing an effective vaccine against this fatal disease. In this study, DNA vaccines containing the ADRV 2L gene (pcDNA-2L) and the 58L gene (pcDNA-58L) were respectively constructed, and their immune protective effects were evaluated in Chinese giant salamanders. In vitro and in vivo expression of the vaccine plasmids were confirmed in transfected cells and muscle tissues of vaccinated Chinese giant salamanders by using immunoblot analysis or RT-PCR. Following ADRV challenge, the Chinese giant salamanders vaccinated with pcDNA-2L showed a relative percent survival (RPS) of 66.7%, which was significant higher than that in Chinese giant salamanders immunized with pcDNA-58L (RPS of 3.3%). Moreover, the specific antibody against ADRV was detected in Chinese giant salamanders vaccinated with pcDNA-2L at 14 and 21 days post-vaccination by indirect enzyme-linked immunosorbent assay (ELISA). Transcriptional analysis revealed that the expression levels of immune-related genes including type I interferon (IFN), myxovirus resistance (Mx), major histocompatibility complex class IA (MHCIA), and immunoglobulin M (IgM) were strongly up-regulated after vaccination with pcDNA-2L. Furthermore, vaccination with pcDNA-2L significantly suppressed the virus replication, which was seen by a low viral load in the spleen of Chinese giant salamander survivals after ADRV challenge. These results indicated that pcDNA-2L could induce a significant innate immune response and an adaptive immune response involving both humoral and cell-mediated immunity that conferred effective protection against ADRV infection, and might be a potential vaccine candidate for controlling ADRV disease in Chinese giant salamanders.

Transcriptome analyses of immune tissues from three Japanese frogs (genus Rana ) reveals their utility in characterizing major histocompatibility complex class II

Background

In Japan and East Asia, endemic frogs appear to be tolerant or not susceptible to chytridiomycosis, a deadly amphibian disease caused by the chytrid fungus Batrachochytridium dendrobatidis (Bd). Japanese frogs may have evolved mechanisms of immune resistance to pathogens such as Bd. This study characterizes immune genes expressed in various tissues of healthy Japanese Rana frogs.

Results

We generated transcriptome data sets of skin, spleen and blood from three adult Japanese Ranidae frogs (Japanese brown frog Rana japonica, the montane brown frog Rana ornativentris, and Tago’s brown frog Rana tagoi tagoi) as well as whole body of R. japonica and R. ornativentris tadpoles. From this, we identified tissue- and stage-specific differentially expressed genes; in particular, the spleen was most enriched for immune-related genes. A specific immune gene, major histocompatibility complex class IIB (MHC-IIB), was further characterized due to its role in pathogen recognition. We identified a total of 33 MHC-IIB variants from the three focal species (n = 7 individuals each), which displayed evolutionary signatures related to increased MHC variation, including balancing selection. Our supertyping analyses of MHC-IIB variants from Japanese frogs and previously studied frog species identified potential physiochemical properties of MHC-II that may be important for recognizing and binding chytrid-related antigens.

Conclusions

This is one of the first studies to generate transcriptomic resources for Japanese frogs, and contributes to further understanding the immunogenetic factors associated with resistance to infectious diseases in amphibians such as chytridiomycosis. Notably, MHC-IIB supertyping analyses identified unique functional properties of specific MHC-IIB alleles that may partially contribute to Bd resistance, and such properties provide a springboard for future experimental validation.

Electronic supplementary material

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