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

Skin microbiota during metamorphosis of Quasipaa spinosa: guidance for maintaining mucosal symbiotic microbial flora homeostasis in early life of frogs

The skin microbiota plays an essential role in helping the host adapt to different environments and maintain health. By examining the characteristics of amphibian skin flora alongside ontogenetic traits, we can gain insights into the adaptation mechanisms of amphibian skin flora to environmental changes during development. In this study, we analyzed the skin microbiota of Quasipaa spinosa during metamorphosis using Illumina sequencing. Venn diagrams and UpSet analysis revealed that the LTS (hindlimb tadpoles' skin, aquatic habitat) and FTS (forelimb tadpoles' skin, shift from aquatic to amphibious habitats) groups exhibited a higher number of unique amplicon sequence variants (ASVs), while the TS (tadpoles' skin, aquatic habitat) and LFS (land frogs' skin, amphibious habitats) groups displayed a lower abundance of ASVs. Diversity analysis indicated similarities in the microorganisms between the LTS and the FTS groups, with higher microbial diversity compared to the TS and the LFS groups. Additionally, microbial co-occurrence network analysis indicated a more stable microecology in the LTS group and FTS group. Proteobacteria, Firmicutes, and Bacteroidota were identified as the dominant phyla, although their relative abundances varied widely among groups. LEfSe (Linear discriminant analysis effect size) showed significant enrichment of beneficial bacteria at various developmental stages, including Bacteroides, Bacillus, and Lactobacillus. Functional prediction analysis shows significant differences in skin microorganism functions across various developmental stages, with a primary focus on metabolic functions. This study provides valuable insights into the compositional dynamics of skin microbiota in Q. spinosa at various developmental stages.

Assessment of the antimicrobial and immunomodulatory activity of QS-CATH, a promising therapeutic agent isolated from the Chinese spiny frogs (Quasipaa spinosa)

Cathelicidins are important antimicrobial peptides in various vertebrate species where they are crucial parts of the innate immune system. The current understanding of amphibian cathelicidins is limited, particularly with regard to their immunomodulatory effects. To address this knowledge gap, we produced the cDNA sequence of the cathelicidin gene from a skin transcriptome of the Chinese spiny frog Quasipaa spinosa. The amino acid sequence of the Quasipaa spinosa cathelicidin (QS-CATH) was predicted to consist of a signal peptide, a cathelin domain, and a mature peptide. Comparative analysis of the QS-CATH amino acid sequence with that of other amphibian cathelicidins revealed high variability in the functional mature peptide among amphibians, whereas the cathelin domain was conserved. The QS-CATH gene was expressed in several tissues, with the highest level of expression in the spleen. Upregulation of QS-CATH after Aeromonas hydrophila infection occurred in the kidney, gut, spleen, skin, and liver. Chemically synthesized QS-CATH exhibited pronounced antibacterial activity against Shigella flexneri, Staphylococcus warneri, Escherichia coli, Salmonella enterica, and Listeria monocytogenes. Furthermore, QS-CATH disrupted the cell membrane integrity of S. flexneri, as evidenced by a lactate dehydrogenase release assay, and it hydrolyzed the genomic DNA of S. flexneri. Additionally, QS-CATH elicited chemotaxis and modulated the expression of inflammatory cytokine genes in RAW264.7 mouse leukemic monocyte/macrophage cells. These findings confirm the antimicrobial effects of amphibian cathelicidin and its ability to influence immune cell function. This will expedite the potential utilization of amphibian antimicrobial peptides as therapeutic agents.

Antimicrobial and immunomodulatory activity of beta-defensin from the Chinese spiny frog (Quasipaa spinosa)

The β-defensins are important components of the vertebrate innate immune system. While mammalian β-defensins have wide-ranging antibacterial and immunomodulatory activities, those of amphibians remain largely uncharacterised. In this study, β-defensin cDNA was identified from the skin transcriptome of the Chinese spiny frog Quasipaa spinosa. This β-defensin (QS-BD) consists of a signal and a mature peptide. Sequence alignments with other amphibian β-defensins showed conservation of the functional mature peptide and that its closest relative is β-defensin from Zhangixalus puerensis. Synthetic QS-BD showed antibacterial activity against Vibrio vulnificus, Vibrio harveyi, Streptococcus iniae, and Aeromonas hydrophila. QS-BD showed bactericidal activity by destroying the cell membrane integrity, but did not hydrolyse genomic DNA. QS-BD treatment promoted respiratory bursts and upregulated the expression of interleukin-1β and tumour necrosis factor-α in the murine leukemic monocyte/macrophage cell line RAW264.7. This is the first demonstration of immunomodulatory activity by an amphibian β-defensin.

Functional differences of three CXCL10 homologues in the giant spiny frog Quasipaa spinosa

Chemokines are a superfamily of structurally related chemotactic cytokines exerting significant roles in acting as a bridge between the innate and adaptive immune responses. In this study, we identified three CXC motif chemokine 10 (CXCL10) homologues (QsCXCL10-1, QsCXCL10-2 and QsCXCL10-3) from giant spiny frog Quasipaa spinosa. All three deduced QsCXCL10 proteins contained four conserved cysteine residues as found in other known CXC chemokines. Phylogenetic analysis showed that QsCXCL10-1, 2, 3 and other CXCL10s in amphibian were grouped together to form a separate clade. These three QsCXCL10s were highly expressed in spleen and blood. Upon infection with Staphylococcus aureus or Aeromonas hydrophila, the expressions of QsCXCL10s were markedly increased in spleen and blood during biotic stresses. Meanwhile, the QsCXCL10s transcription in liver could also be up-regulated under abiotic stresses such as cold and heat stresses. The recombinant proteins of frog CXCL10 homologues were produced and purified in E. coli and possessed similar but differential bioactivities. Both rCXCL10-1 and rCXCL10-2 had strong effects on the up-regulation of pro-inflammatory cytokines (TNF-α, IL-1β and IL-8) in vivo, whereas rCXCL10-3 induced a weak expression of these cytokines. Moreover, the rCXCL10-1 and rCXCL10-2 could strongly promote splenocyte proliferation and induce lymphocytes migration, while rCXCL10-3 had limited effects on these biological processes. All three frog chemokines triggered their functional activities by engaging CXC motif chemokine receptor 3 (CXCR3). Taken together, these results revealed that the three QsCXCL10s had similar but differential functional activities in mediating immune responses and host defenses, which might contribute to a better understanding of the functional evolution of CXCL10 in vertebrates.

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.

Interspecies introgressive hybridization in spiny frogs Quasipaa (Family Dicroglossidae) revealed by analyses on multiple mitochondrial and nuclear genes

Introgression may lead to discordant patterns of variation among loci and traits. For example, previous phylogeographic studies on the genus Quasipaa detected signs of genetic introgression from genetically and morphologically divergent Quasipaa shini or Quasipaa spinosa. In this study, we used mitochondrial and nuclear DNA sequence data to verify the widespread introgressive hybridization in the closely related species of the genus Quasipaa, evaluate the level of genetic diversity, and reveal the formation mechanism of introgressive hybridization. In Longsheng, Guangxi Province, signs of asymmetrical nuclear introgression were detected between Quasipaa boulengeri and Q. shini. Unidirectional mitochondrial introgression was revealed from Q. spinosa to Q. shini. By contrast, bidirectional mitochondrial gene introgression was detected between Q. spinosa and Q. shini in Lushan, Jiangxi Province. Our study also detected ancient hybridizations between a female Q. spinosa and a male Q. jiulongensis in Zhejiang Province. Analyses on mitochondrial and nuclear genes verified three candidate cryptic species in Q. spinosa, and a cryptic species may also exist in Q. boulengeri. However, no evidence of introgressive hybridization was found between Q. spinosa and Q. boulengeri. Quasipaa exilispinosa from all the sampling localities appeared to be deeply divergent from other communities. Our results suggest widespread introgressive hybridization in closely related species of Quasipaa and provide a fundamental basis for illumination of the forming mechanism of introgressive hybridization, classification of species, and biodiversity assessment in Quasipaa.

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.

Transcriptome analysis of the endangered Chinese giant salamander (Andrias davidianus): Immune modulation in response to Aeromonas hydrophila infection

The endangered Chinese giant salamander (Andrias davidianus) is the largest extant amphibian species. Disease outbreaks represent one of the major factors threatening A. davidianus populations in the wild and the viability of artificial breeding programmes. Development of future immune therapies to eliminate infectious disease in A. davidianus is dependent on a thorough understanding of the immune mechanisms elicited by pathogen encounters. To this end we have undertaken, for the first time in amphibians, differential transcriptome analysis of the giant salamander response to Aeromonas hydrophila, one of the most devastating pathogens affecting amphibian populations. Out of 87,204 non-redundant consensus unigenes 19,216 were annotated, 6834 of which were upregulated and 906 down-regulated following bacterial infection. 2058 unigenes were involved with immune system processes, including 287 differentially expressed unigenes indicative of the impact of bacterial infection on several innate and adaptive immune pathways in the giant salamander. Other pathways not directly associated with immune-related activity were differentially expressed, including developmental, structural, molecular and growth processes. Overall, this work provides valuable insights into the underlying immune mechanisms elicited during bacterial infection in amphibians that may aid in the future development of disease control measures in protecting the Chinese giant salamander. With the unique position of amphibians in the transition of tetrapods from aquatic to terrestrial habitats, our study will also be invaluable towards the further understanding of the evolution of tetrapod immunity.