An endemic frog harbors multiple expression loci with different patterns of variation in the MHC class II B gene

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.

The succession of gut microbiota in the concave-eared torrent frog (Odorrana tormota) throughout developmental history

The gut microbiota of amphibians plays a crucial role in maintaining health and adapting to various developmental stages. The composition of gut microbial community is influenced by the phylogeny, habitat, diet, and developmental stage of the host. The present study analyzed the microbiota in the intestine of O. tormota at 11 developmental stages (from the tadpole at Gosner stage 24 to the 3-year-old adult) using high-throughput 16S rRNA sequencing. Alpha diversity index analysis of the microbiota revealed that the index decreased from tadpole at Gosner stage 24 to adult frog stage, remained stable during the adult frog stages, but increased significantly at the early metamorphosis and hibernation preparation stages. The gut microbiota structure is similar in adult frogs but differs significantly in other developmental stages. Furthermore, the dominant phyla of gut microbiota in tadpoles were Proteobacteria, Firmicutes, Actinobacteria, and Bacteroidetes, whereas those in adult frogs were Proteobacteria, Firmicutes, Bacteroidetes, and Verrucomicrobia. Host and environmental factors jointly affected the gut microbial diversity and community composition of O. tormota, but developmental stage, feeding habit, and habitat type had a more significant influence. The microbial community in the gut varies with the developmental stage of the host and constantly adapts to the survival requirements of the host. These findings advance our understanding of the evolutionary mechanism of amphibian gut microbiota in maintaining health homeostasis and adaptation.

Dynamic transcriptome and histomorphology analysis of developmental traits of hindlimb thigh muscle from Odorrana tormota and its adaptability to different life history stages

Background

Systematic studies on the development and adaptation of hindlimb muscles in anura amphibians are rare. Here, we integrated analysis of transcriptome and histomorphological data for the hindlimb thigh muscle of Odorrana tormota (concave-eared torrent frog) at different developmental stages, to uncover the developmental traits of hindlimb thigh muscle from O. tormota and its adaptability to different life history stages.

Results

The development of hindlimb thigh muscle from O. tormota has the following characteristics. Before metamorphosis, myogenous cells proliferate and differentiate into myotubes, and form 11 muscle groups at G41; Primary myofibers and secondary myofibers appeared during metamorphosis; 11 muscle groups differentiated continuously to form myofibers, accompanied by myofibers hypertrophy after metamorphosis; During the growth process of O. tormota from G42 to G46, there were differences between the sexes in the muscle groups that differentiate into muscle fibers, indicating that there was sexual dimorphism in the hindlimb thigh muscles of O. tormota at the metamorphosis stages. Some genes and pathways related to growth, development, and movement ability of O. tormota at different developmental stages were obtained. In addition, some pathways associated with adaptation to metamorphosis and hibernation also were enriched. Furthermore, integrated analysis of the number of myofibers and transcriptome data suggested that myofibers of specific muscle groups in the hindlimbs may be degraded through lysosome and ubiquitin pathways to transform into energy metabolism and other energy-related substances to meet the physiological needs of hibernation.

Conclusions

These results provide further understanding the hindlimb thigh muscle development pattern of frogs and their adaption to life history stages.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07677-0.

Limited polymorphism of the functional MHC class II B gene in the black-spotted frog (Pelophylax nigromaculatus) identified by locus-specific genotyping

Amphibians can be more vulnerable to environmental changes and diseases than other species because of their complex life cycle and physiological requirements. Therefore, understanding the adaptation of amphibians to environmental changes is crucial for their conservation. Major histocompatibility complex (MHC) presents an excellent tool for the investigation of adaptive variations and the assessment of adaptive potential because it can be under strong diversifying selection. Here, we isolated the MHC class II B (MHCIIB) gene from cDNA sequences of the black-spotted frog (Pelophylax nigromaculatus), a widespread amphibian species in China, and designed locus-specific primers to characterize adaptive variability of this amphibian. Ten alleles were identified from 67 individual frogs of three populations and no more than two alleles were present in each individual animal. Furthermore, none of the sequences had indels or/and stop codons, which is in good agreement with locus-specific amplification of a functional gene. However, we found low polymorphism at both nucleotide and amino acid levels, even in the antigen-binding region. Purifying selection acting at this locus was supported by the findings that the d_N/d_S ratio across all alleles was below 1 and that negatively selected sites were detected by different methods. Allele frequency distributions were significantly different among geographic populations, indicating that physiographic factors may have strong effect on the genetic structure of the black-spotted frog. This study revealed limited polymorphism of three adjacent black-spotted frog populations at the functional MHCIIB locus, which may be attributed to region-specific differences. The locus-specific genotyping technique developed in this study would provide a foundation for future studies on adaptive divergence among different frog populations.

Evolutionary dynamics of an expressed MHC class IIβ locus in the Ranidae (Anura) uncovered by genome walking and high-throughput amplicon sequencing

The Major Histocompatibility Complex (MHC) is a genomic region encoding immune loci that are important and frequently used markers in studies of adaptive genetic variation and disease resistance. Given the primary role of infectious diseases in contributing to global amphibian declines, we characterized the hypervariable exon 2 and flanking introns of the MHC Class IIβ chain for 17 species of frogs in the Ranidae, a speciose and cosmopolitan family facing widespread pathogen infections and declines. We find high levels of genetic variation concentrated in the Peptide Binding Region (PBR) of the exon. Ten codons are under positive selection, nine of which are located in the mammal-defined PBR. We hypothesize that the tenth codon (residue 21) is an amphibian-specific PBR site that may be important in disease resistance. Trans-species and trans-generic polymorphisms are evident from exon-based genealogies, and co-phylogenetic analyses between intron, exon and mitochondrial based reconstructions reveal incongruent topologies, likely due to different locus histories. We developed two sets of barcoded adapters that reliably amplify a single and likely functional locus in all screened species using both 454 and Illumina based sequencing methods. These primers provide a resource for multiplexing and directly sequencing hundreds of samples in a single sequencing run, avoiding the labour and chimeric sequences associated with cloning, and enabling MHC population genetic analyses. Although the primers are currently limited to the 17 species we tested, these sequences and protocols provide a useful genetic resource and can serve as a starting point for future disease, adaptation and conservation studies across a range of anuran taxa.

Genetic polymorphism of major histocompatibility complex class IIB alleles and pathogen resistance in the giant spiny frog Quasipaa spinosa

Major histocompatibility complex (MHC) genes are candidates for determining disease susceptibility due to their pivotal role in both innate and adaptive immune responses. Accordingly, the association between the genetic variation of MHC genes and the pathogen resistance has been investigated in numerous vertebrates. To date, however, little is reported in amphibians. In this study, we investigate the genetic variation at the MHC class IIB gene in the giant spiny frog Quasipaa spinosa, which has high commercial value in China. The full length of MHC class IIB cDNA was cloned from Q. spinosa by homology cloning and rapid amplification of cDNA end-polymerase chain reaction (RACE-PCR). Two MHC class IIB loci were identified in Q. spinosa. We also developed PCR primers for a portion of the second exon of the MHC class IIB gene. A total of 26 MHC class IIB alleles were identified. The dN rate was significantly higher than the dS rate in the putative peptide-binding region, thereby proving the positive selection hypothesis. In addition, individuals intraperitoneally injected with Aeromonas hydrophila were used to study the association between MHC class IIB alleles and pathogen resistance/susceptibility, to explore the specific alleles in balancing selection. Eighty frogs were used after exposure to A. hydrophila infection. Nine alleles were used to study the association between the alleles and disease resistance. Two alleles, namely, Pasa-DAB(∗)1301 and Pasa-DAB(∗)0901, were significantly associated with resistance against A. hydrophila. This study provides valuable information on the structure of the MHC class IIB gene and confirms the association between MHC class IIB gene alleles and disease resistance to bacterial infection in Q. spinosa. Moreover, pathogen resistance-related MHC markers can be used for the selective breeding of the giant spiny frog.