Genomic differentiation with isolation by distance along a latitudinal gradient in the spotted-leg treefrog Polypedates megacephalus.

The patterns of isolation by distance (IBD) entailing increased genetic differentiation among populations have aroused extensive concerns for evolutionary biologists. Although the IBD may act on spatial processes contributing to the genetic differentiation among populations in anuran species, the factors shaping the IBD of frogs among populations in natural systems are largely unknown. Here, we studied the genetic differentiation among six populations with 24 individuals of the spotted-leg treefrog along a latitudinal gradient (1860.31 km) based on 1020 single nucleotide polymorphisms from restriction site-associated DNA sequencing. The results showed that the genetic diversity differed significantly among populations and that the insular populations had higher genetic diversity than the mainland populations. Furthermore, we also found a significant genetic differentiation among populations (FST = 0.277) and no sign of inbreeding (FUNI = -0.145). The IBD was detected for all populations, and a higher degree of the IBD was indicated when controlling for the effects of the isolation between Hainan and mainland populations caused by the Qiongzhou Strait. Our findings suggest that the form of the Qiongzhou Strait plays a key role in shaping the genetic diversity and population differentiation in treefrogs.

Two-Fold ND5 Genes, Three-Fold Control Regions, lncRNA, and the "Missing" ATP8 Found in the Mitogenomes of Polypedates megacephalus (Rhacophridae: Polypedates).

In prior research on the mitochondrial genome (mitogenome) of Polypedates megacephalus, the one copy of ND5 gene was translocated to the control region (CR) and the ATP8 gene was not found. Gene loss is uncommon among vertebrates. However, in this study, we resequenced the mitogenomes of P. megacephalus from different regions using a "primer bridging" approach with Sanger sequencing technologies, which revealed the "missing" ATP8 gene in P. megacephalus as well as three other previously published Polypedates. The mitogenome of this species was found to contain two copies of the ND5 genes and three copies of the control regions. Furthermore, multiple tandem repeats were identified in the control regions. Notably, we observed that there was no correlation between genetic divergence and geographic distance. However, using the mitogenome, gene expression analysis was performed via RT-qPCR of liver samples and it was thus determined that COIII, ND2, ND4, and ND6 were reduced to 0.64 ± 0.24, 0.55 ± 0.34, 0.44 ± 0.21 and 0.65 ± 0.17, respectively, under low-temperature stress (8 °C) as compared with controls (p < 0.05). Remarkably, the transcript of long non-coding RNA (lncRNA) between positions 8029 and 8612 decreased significantly with exposure to low-temperature stress (8 °C). Antisense ND6 gene expression showed a downward trend, but this was not significant. These results reveal that modulations of protein-coding mitochondrial genes and lncRNAs of P. megacephalus play a crucial role in the molecular response to cold stress.

Comprehensive assessment of the ecological risk of exposure to triphenyl phosphate in a bioindicator tadpole.

The toxicity of triphenyl phosphate (TPhP) to aquatic organisms in surface waters has been demonstrated; However, an understanding of toxicity profiles of TPhP in amphibians is limited. Therefore, the adverse effects and threshold concentrations of TPhP on metamorphosis, growth, locomotion, and hepatic antioxidants of Gosner stage 25 Polypedates megacephalus tadpoles under long-term (35 d) exposure to six TPhP concentrations until complete metamorphosis were assessed. Additionally, the overall effect of using integrated multiple biomarkers were determined to demonstrate the potential ecological risks of waterborne TPhP at environmentally relevant concentrations in amphibian tadpoles. With increasing TPhP concentrations, physical parameters (snout-vent length, body mass, condition factor, and hepatic somatic index), jumping distance, hepatic catalase, and superoxide dismutase activities decreased, whereas metamorphosis time and malondialdehyde content increased. The threshold concentration of TPhP that affected the tadpole biomarker, except for metamorphosis rate and jumping distance, was 50-400 µg/L. Furthermore, the standardized scores of the examined integrated biomarkers in the six TPhP concentrations were visualized using radar plots and calculated as the integrated biomarker responses (IBRs). The varying TPhP concentrations had different scores in the radar plots, and the threshold for affecting the IBR value was 10 µg/L, which was close to the TPhP concentration in surface waters. Additionally, IBR values were strongly positively correlated with the TPhP concentrations. These findings indicate that environmentally relevant exposure to waterborne TPhP can pose an ecological risk to amphibian tadpoles. This study can serve as a reference and assist in the formulation of relevant policies and strategies to control TPhP pollution in water bodies.

The revised complete mitogenome sequence of the tree frog Polypedatesmegacephalus (Anura, Rhacophoridae) by next-generation sequencing and phylogenetic analysis.

The mitochondrial genome (mitogenome) sequence of the tree frog Polypedates megacephalus (16,473 bp) was previously reported as having the unusual characteristic of lacking the ND5 gene. In this study, a new mitogenome of P. megacephalus (19,952 bp) was resequenced using the next-generation sequencing (NGS) and standard Sanger sequencing technologies. It was discovered that the ND5 gene was not lost but translocated to the control region (CR) from its canonical location between the ND4 and ND6 genes. In addition, a duplicated control region was found in the new mitogenome of this species. Conservative region identification of the ND5 gene and phylogenetic analysis confirmed that the ND5 gene was located between two control regions. The phylogenetic relationship among 20 related species of anura revealed a rearrangement of the ND5 gene during the evolutionary process. These results also highlighted the advantages of next-generation sequencing. It will not only decrease the time and cost of sequencing, but also will eliminate the errors in published mitogenome databases.

Inferring Microbial Interactions in the Gut of the Hong Kong Whipping Frog (Polypedates megacephalus) and a Validation Using Probiotics.

The concerted activity of intestinal microbes is crucial to the health and development of their host organisms. Investigation of microbial interactions in the gut should deepen our understanding of how these micro-ecosystems function. Due to advances in Next Generation Sequencing (NGS) technologies, various bioinformatic strategies have been proposed to investigate these microbial interactions. However, due to the complexity of the intestinal microbial community and difficulties in monitoring their interactions, at present there is a gap between the theory and biological application. In order to construct and validate microbial relationships, we first induce a community shift from simple to complex by manipulating artificial hibernation (AH) in the treefrog Polypedates megacephalus. To monitor community growth and microbial interactions, we further performed a time-course screen using a 16S rRNA amplicon approach and a Lotka-Volterra model. Lotka-Volterra models, also known as predator-prey equations, predict the dynamics of microbial communities and how communities are structured and sustained. An interaction network of gut microbiota at the genus level in the treefrog was constructed using Metagenomic Microbial Interaction Simulator (MetaMIS) package. The interaction network obtained had 1,568 commensal, 1,737 amensal, 3,777 mutual, and 3,232 competitive relationships, e.g., Lactococcus garvieae has a commensal relationship with Corynebacterium variabile. To validate the interacting relationships, the gut microbe composition was analyzed after probiotic trials using single strain (L. garvieae, C. variabile, and Bacillus coagulans, respectively) and a combination of L. garvieae, C. variabile, and B. coagulans, because of the cooperative relationship among their respective genera identified in the interaction network. After a 2 week trial, we found via 16S rRNA amplicon analysis that the combination of cooperative microbes yielded significantly higher probiotic concentrations than single strains, and the immune response (interleukin-10 expression) also significantly changed in a manner consistent with improved probiotic effects. By taking advantage of microbial community shift from simple to complex, we thus constructed a reliable microbial interaction network, and validated it using probiotic strains as a test system.