Spatio-temporal analysis of toxigenic genes expression in the growing Bufo gargarizans based on RNA sequencing data

BACKGROUND

Bufo gargarizans Cantor, a widely distributed amphibian species in Asia, produces and releases toxins through its retroauricular and granular glands. Although various tissues have been sequenced, the molecular mechanisms underlying the toxin production remain unclear. To elucidate these mechanisms, abdominal skin (non-toxic secretory glands) and retroauricular gland (toxic secreting glands) samples were collected at different time points (3, 6, 12, 24, and 36 months) for RNA sequencing (RNA-seq) and analysis.

RESULTS

In comparison to the S group during the same period, a total of 3053, 3026, 1516, 1028, and 2061 differentially expressed genes (DEGs) were identified across five developmental stages. Gene Ontology (GO) analysis revealed that DEGs were primarily enriched in biological processes including cellular processes, single-organism processes, metabolic processes, and biological regulation. In terms of cellular components, the DEGs were predominantly localized in the cell and cell parts, whereas molecular function indicated significant enrichment in binding and catalytic activity. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that the metabolism and synthesis of various substances, such as lipid metabolism, cofactor and vitamin metabolism, tryptophan metabolism, steroid biosynthesis, and primary bile acid biosynthesis, were accompanied by the development of toads. Additionally, using trend analysis, we discovered candidate genes that were upregulated in the retroauricular glands during development, and the abundance of these genes in the abdominal skin was extremely low. Finally, we identified 26 genes that are likely to be involved in toxin production and that are likely to be involved in toxin anabolism.

CONCLUSION

Overall, these results provide new insights into the genes involved in toxin production in B. gargarizans, which will improve our understanding of the molecular mechanisms underlying toxigenic gene expression.

An amphibian-derived cathelicidin accelerates cutaneous wound healing through its main regulatory effect on phagocytes

Cathelicidins are an important family of antimicrobial peptides (AMPs) involved in the innate immunity in vertebrates. The mammalian cathelicidins have been well characterized, but the relationship between structure and function in amphibian cathelicidins is still not well understood. In this study, a novel 29-residue cathelicidin antimicrobial peptide (BugaCATH) was identified from the skin of Bufo gargarizans. Unlike other AMPs, BugaCATH does not display any direct antimicrobial effects in vitro. However, it effectively promotes full-thickness wound repair in mice. Following injury, BugaCATH initiates and expedites the inflammatory stage by recruiting neutrophils and macrophages to the wound site. BugaCATH not only regulates neutrophil phagocytic activity but also stimulates the generation of cytokines (TNF-α, IL-6, and IL-1β) and chemokines (CXCL1, CXCL2, CCL2, and CCL3) in macrophages and in mice. Furthermore, it promotes macrophage M2 polarization that facilitates the conversion from a pro-inflammatory macrophage-dominated wound environment to an anti-inflammatory one during the mid to late stages, which is crucial for reducing inflammation and effective wound repair. The MAPK (ERK, JNK, and p38) and NF-κB-NLRP3 signaling pathways are involved in the activity. Moreover, BugaCATH directly enhances the migration of keratinocytes and vascular endothelial cells without affecting their proliferation. Notably, BugaCATH significantly improves the proliferation of keratinocytes and endothelial cells in the presence of macrophages. The current study revealed that in addition to proliferation of keratinocytes and endothelial cells, BugaCATH possesses the ability to modulate inflammatory processes during skin injury through its regulatory effect on phagocytes. The combination of these capabilities makes BugaCATH a potent candidate for skin wound therapy.

Unraveling the antimicrobial potential of Lactiplantibacillus plantarum strains TE0907 and TE1809 sourced from Bufo gargarizans: advancing the frontier of probiotic-based therapeutics

Introduction

In an era increasingly defined by the challenge of antibiotic resistance, this study offers groundbreaking insights into the antibacterial properties of two distinct Lactiplantibacillus plantarum strains, TE0907 and TE1809, hailing from the unique ecosystem of Bufo gargarizans. It uniquely focuses on elucidating the intricate components and mechanisms that empower these strains with their notable antibacterial capabilities.

Methods

The research employs a multi-omics approach, including agar diffusion tests to assess antibacterial efficacy and adhesion assays with HT-29 cells to understand the preliminary mechanisms. Additionally, gas chromatography-mass spectrometry (GC-MS) is employed to analyze the production of organic acids, notably acetic acid, and whole-genome sequencing is utilized to identify genes linked to the biosynthesis of antibiotics and bacteriocin-coding domains.

Results

The comparative analysis highlighted the exceptional antibacterial efficacy of strains TE0907 and TE1809, with mean inhibitory zones measured at 14.97 and 15.98 mm, respectively. A pivotal discovery was the significant synthesis of acetic acid in both strains, demonstrated by a robust correlation coefficient (cor ≥ 0.943), linking its abundance to their antimicrobial efficiency. Genomic exploration uncovered a diverse range of elements involved in the biosynthesis of antibiotics similar to tetracycline and vancomycin and potential regions encoding bacteriocins, including Enterolysin and Plantaricin.

Conclusion

This research illuminates the remarkable antibacterial efficacy and mechanisms intrinsic to L. plantarum strains TE0907 and TE1809, sourced from B. gargarizans. The findings underscore the strains' extensive biochemical and enzymatic armamentarium, offering valuable insights into their role in antagonizing enteric pathogens. These results lay down a comprehensive analytical foundation for the potential clinical deployment of these strains in safeguarding animal gut health, thereby enriching our understanding of the role of probiotic bacteria in the realm of antimicrobial interventions.

3D Morphological Scanning and Environmental Correlates of Bufo gargarizans in the Yellow River Basin

Morphology plays a crucial role in understanding the intricacies of biological forms. Traditional morphometric methods, focusing on one- or two-dimensional geometric levels, often fall short of accurately capturing the three-dimensional (3D) structure of organisms. The advent of 3D scanning techniques has revolutionized the study of organismal morphology, enabling comprehensive and accurate measurements. This study employs a 3D structured light scanning system to analyze the morphological variations in the Chinese toad (Bufo gargarizans Cantor, 1842) along the Yellow River Basin. The 3D digital model obtained from the scan was used to calculate various morphological parameters including body surface area, volume, fractal dimensions, and limb size. The research explores geographic variability patterns and identifies environmental drivers affecting the 3D phenotypic variation of B. gargarizans. Results reveal a bimodal pattern of variation in the toad population, with higher elevations exhibiting smaller body sizes, greater appendage proportions, and more complex body structures. Linear regression analyses highlight the influence of elevation and annual mean temperature on the morphological variation of B. gargarizans, with elevation playing a significant role. This study underscores the significance of 3D morphometric analysis in unraveling the intricacies of organismal morphology and understanding the adaptive strategies of species in diverse environments.

Elevational variation in metabolic rate, feeding capacity and their associations in the Asiatic toad Bufo gargarizans

Foraging behavior is known to place demands on the metabolic characteristics of anurans. Active foragers feeding on sedentary prey typically have high aerobic capacity and low anaerobic capacity, whereas sit-and-wait foragers feeding on active and mobile prey have the opposite pattern. Thus, the energetic demands of foraging may influence their metabolic adaptations to harsh environments, such as high elevations. Anurans that engage in active foraging have been found to increase maximum metabolic rate (MMR) and aerobic scope (AS, the difference between MMR and resting metabolic rate, RMR) at high elevations. However, data are lacking in amphibian ambush foragers. In this study, we examined the RMR, MMR, AS, and feeding capacity of a sit-and-wait forager ─the Asiatic toad (Bufo gargarizans), from two populations that are in close geographic proximity but differ by 1350 m in elevation. Our results show that there is no elevational variation in RMR and feeding capacity in either males or females. However, there are sex-specific variations in MMR and AS along an elevational gradient; females from high elevations have lower MMR and smaller net AS than their counterparts from low elevations while males maintain similar MMR and net AS across elevations. Furthermore, aerobic performances do not appear to be associated with feeding capacity at either the individual or population level. Our results support the hypothesis that sit-and-wait foragers may not increase their aerobic capacity as a strategy in hypoxic and low food availability environments and the role of sex in these adaptive adjustments should not be overlooked.

Sexual Dimorphism in the Limb Bones of Asiatic Toad (Bufo gargarizans) in Relation to Sexual Selection

Sexual dimorphism is often considered to be the result of differences in the intensity of sexual selection between sexes. From this point of view, the sexual dimorphism of the limb bones of the Bufo gargarizans in southwest China was studied. Results showed that the fore- and hindlimb skeletons of this species were sexually dimorphic in anatomy. The humerus, radioulna, and total lengths of the forelimb skeleton of males were substantially longer than those of females, but the hand length of males was smaller than that of females. Several other features of males, such as deltoid and medial crest areas and humerus and radioulnar weights, were also significantly larger than those of females. The femoris, tibiofibula, talus-calcaneus, and foot lengths; total hindlimb skeleton length; and femoral upper crest areas of males were significantly greater than those of females. However, no significant intersexual difference in femoris and tibiofibular weights was observed. These findings suggested that robust forelimb bones and long hindlimb bones could contribute to the mating success of males; if so, sexual selection promotes the evolution of sexual size and shape dimorphism in the limb bones of the B. gargarizans.

Identification of Protein Quality Markers in Toad Venom from Bufo gargarizans

Toad venom is a traditional Chinese medicine with high medicinal value. The existing quality evaluation standards of toad venom have obvious limitations because of the lack of research on proteins. Thus, it is necessary to screen suitable quality markers and establish appropriate quality evaluation methods for toad venom proteins to guarantee their safety and efficacy in clinical applications. SDS-PAGE, HPLC, and cytotoxicity assays were used to analyze differences in protein components of toad venom from different areas. Functional proteins were screened as potential quality markers by proteomic and bioinformatic analyses. The protein components and small molecular components of toad venom were not correlated in content. Additionally, the protein component had strong cytotoxicity. Proteomics analysis showed that 13 antimicrobial proteins, four anti-inflammatory and analgesic proteins, and 20 antitumor proteins were differentially expressed extracellular proteins. A candidate list of functional proteins was coded as potential quality markers. Moreover, Lysozyme C-1, which has antimicrobial activity, and Neuropeptide B (NPB), which has anti-inflammatory and analgesic activity, were identified as potential quality markers for toad venom proteins. Quality markers can be used as the basis of quality studies of toad venom proteins and help to construct and improve safe, scientific, and comprehensive quality evaluation methods.

Copper-mediated shifts in transcriptomic responses of intestines in Bufo gargarizans tadpoles to lead stress

The differential transcriptomic responses of intestines in Bufo gargarizans tadpoles to Pb alone or in the presence of Cu were evaluated. Tadpoles were exposed to 30 μg/L Pb individually and in combination with Cu at 16 or 64 μg/L from Gosner stage (Gs) 26 to Gs 38. After de novo assembly, 105,107 unigenes were generated. Compared to the control group, 7387, 6937, and 11139 differentially expressed genes (DEGs) were identified in the treatment of Pb + Cu0, Pb + Cu16, and Pb + Cu64, respectively. In addition, functional annotation and enrichment analysis of DEGs revealed substantial transcriptional reprogramming of diverse molecular and biological pathways were induced in all heavy metal treatments. The relative expression levels of genes associated with intestinal epithelial barrier and bile acids (BAs) metabolism, such as mucin2, claudin5, ZO-1, Asbt, and Ost-β, were validated by qPCR. This study demonstrated that Pb exposure induced transcriptional responses in tadpoles, and the responses could be modulated by Cu.

Comparison of the characteristics of gut microbiota response to lead in Bufo gargarizans tadpole at different developmental stages

In amphibians, lead (Pb) exposure could alter the composition and structure of gut microbiota, but changes involving microbiota of several successive phases following Pb exposure have been less studied. In the present study, we compared the effects of Pb exposure on morphological parameters and gut microbiota of Bufo gargarizans at Gosner stage (Gs) 33, Gs36, and Gs42. Our results showed that total length (TL), snout-vent length (SVL), and body wet weight (TW) of B. gargarizans at Gs33, as well as TL and SVL at Gs42, were significantly increased after Pb exposure. In addition, high-throughput sequencing analysis indicated that gut microbiota has distinct responses to Pb exposure at different developmental stages. The diversity of gut microbiota was significantly reduced under Pb exposure at Gs33, while it was significantly increased at Gs42. In terms of community composition, Spirochaetota, Armatimonadota, and Patescibacteria appeared in the control groups at Gs42, but not after Pb treatment. Furthermore, functional prediction indicated that the relative abundance of metabolism pathway was significantly decreased at Gs33 and Gs36, and significantly increased at Gs42. Our results fill an important knowledge gap and provide comparative information on the gut microbiota of tadpoles at different developmental stages following Pb exposure.

Effects of waterborne Pb/Cu mixture on Chinese toad, Bufo gargarizans tadpoles: morphological, histological, and intestinal microbiota assessment

Coexistence of heavy metals in aquatic environments exert complex effects on amphibians. Here, the adverse effects of Pb (0.14 μM) combined with Cu at concentrations of 0, 0.25, and 1.0 μM were investigated in Bufo gargarizans tadpoles. Tadpoles were chronically exposed from Gosner stage (Gs) 26 to Gs 38, and morphology of tadpoles as well as intestinal histology and bacterial community were assessed. Our results indicated that Pb+Cu1.0 exposure induced significant retardation of somatic mass, total length, intestine mass, and intestine length as well as intestinal histological alterations. Pb+Cu0.25 and Pb+Cu1.0 exposure were associated with the loss of gut bacterial diversity. Proteobacteria and Bacteroidetes were two dominant phyla in tadpoles independently of heavy metal exposure, but the abundance of Proteobacteria increased significantly in Pb+Cu1.0 group and Bacteroidetes decreased significantly in all treatment groups. Furthermore, functional prediction indicated that metabolic disorders were associated with Pb+Cu0.25 and Pb+Cu1.0 exposure. Overall, relative limited shifts in intestinal bacterial diversity, composition, and functionality caused by Pb+Cu0 exposure, while coexistence of Pb and Cu induced gut dysbiosis and might further cause disturbance of metabolic homeostasis. The findings of this study provide insights into the effects of Pb and Cu coexistence on the health of amphibians.

Intestinal changes associated with nitrite exposure in Bufo gargarizans larvae: Histological damage, immune response, and microbiota dysbiosis

Nitrite is a ubiquitous toxic compound in aquatic ecosystems and has negative effects on aquatic organisms. The intestine and the trillions of microbes that inhabit it, play an integral role in maintaining digestive and immune functions. However, the effects of nitrite on intestinal health and microflora have been poorly investigated. Therefore, the present study evaluated the response of intestinal histology, immunity, digestive enzyme activities and microbiota to nitrite exposure in Bufo gargarizans tadpoles. The results showed that nitrite caused damage to the intestine and impaired digestive performance. Significant changes in the transcriptional profiles of genes involved in oxidative stress (sod, gpx and hsp), inflammation, and immunity (socs3, il-27, il-1β and il-17d) were observed in the NO₂-N treatment groups. In addition, exposure to nitrite induced alterations of intestinal microbial diversity, structure and composition, suggesting that nitrite may lead to intestinal microbiota dysbiosis. It is noteworthy that probiotics (e.g., Bacteroidetes and Fusobacteria) were decreased after exposure to nitrite, whereas potentially opportunistic pathogens such as Proteobacteria and Enterobacteriaceae were elevated. Functional prediction and correlation analysis suggested that the above changes may interfere with metabolic function and trigger various diseases. Taken together, we concluded that nitrite exposure induced intestinal microbial dysbiosis, which may lead to immune dysfunction and metabolic disorder, and ultimately to histological damages in B. gargarizans. Further, this study will provide a scientific basis for further understanding the risk of nitrite pollution on the intestinal health of amphibians.

Effects of nitrite exposure on metamorphosis and skeletal development of Bufo gargarizans

Nitrite, as a part of nitrogen cycle, is one of the most common toxic compounds in aquatic ecosystems. Since skeletal development is an essential process during amphibian metamorphosis, exposure of larval amphibians to nitrite might disrupt skeletal development. To evaluate whether nitrite affects skeletal development of amphibian larvae, Bufo gargarizans larvae at Gs26 were exposed to 10, 100, 500 and 1000 μg/L nitrite-nitrogen (NO₂-N) in the present study. The metamorphosis rate, body weight, body length, forelimb length and hindlimb length of B. gargarizans exposed to NO₂-N were decreased. The microscopic structures of thyroid gland were altered under NO₂-N exposure at Gs42. The skeletal lengths of the humerus, femur and fibulare of tadpole at Gs42 were significantly reduced under 100, 500 and 1000 μg/L NO₂-N treatment groups, and the lengths of humerus, tibia-fibula and tibiale of tadpole at Gs46 were significantly reduced under 1000 μg/L NO₂-N treatment groups. In addition, the expression levels of thyroid hormone (TH) and endochondral ossification-related genes of tadpoles at Gs42 and Gs46 were tested by qRT-PCR. Overall, NO₂-N exposure could affect the expressions of these genes and then may influence the activity and function of thyroid gland, further disturbing the amphibian metamorphosis and skeletal development of amphibian larvae.

Effects of Separate and Combined Exposure of Cadmium and Lead on the Endochondral Ossification in Bufo gargarizans

Cadmium (Cd) and lead (Pb) are ubiquitous in aquatic environments and most studies have examined the potential effects of Cd or Pb alone on aquatic organisms. In the present study, chronic effects of Cd and Pb, alone and in combination, on Bufo gargarizans were investigated by exposing embryos to these contaminants throughout metamorphosis. Significant reductions in body mass and snout-to-vent length were observed in B. gargarizans at Gosner stage 42 (Gs 42) and Gs 46 exposed to a Cd/Pb mixture. Single and combined exposure with Cd and Pb induced histological alterations of the thyroid gland characterized by reduced colloid area and thickness of epithelial cells. There was a significant decrease in the maximum jump distance of froglets exposed to Cd alone and the Cd/Pb mixture, and the jumping capacity showed a positive correlation with hind limb length and tibia/fibula. Moreover, single metals and their mixture induced reduction of endochondral bone formation in B. gargarizans. Transcriptomic and real-time quantitative polymerase chain reaction results showed that genes involved in skeletal ossification (TRα, TRβ, Dio2, Dio3, MMP9, MMP13, Runx1, Runx2, and Runx3) were transcriptionally dysregulated by Cd and Pb exposure alone or in combination. Our results suggested that despite the low concentration tested, the Cd/Pb mixture induced more severe impacts on B. gargarizans. In addition, the Cd/Pb mixture might reduce chances of survival for B. gargarizans froglets by decreasing size at metamorphosis, impaired skeletal ossification, and reduction in jumping ability, which might result from dysregulation of genes involved in thyroid hormone action and endochondral ossification. The findings obtained could add a new dimension to understanding of the mechanisms underpinning skeletal ossification response to heavy metals in amphibians. Environ Toxicol Chem 2022;41:1228-1245. © 2022 SETAC.

Comparative transcriptomes of three different skin sites for the Asiatic toad (Bufo gargarizans)

Toads release toxic dry secretions from glands in their skin. Toxin possesses a wide range of biological effects, but little is known about its specific gene expression pattern and regulatory mechanisms. The Asiatic toad (Bufo gargarizans) is widely used to produce toxin. Here, we explored the gene expression of 30 tissue samples from three different skin sites (parotoid gland, dorsal skin, and abdomen skin) of B. gargarizans. After de novo assembly, 783,130 unigenes with an average length of 489 bp (N50 = 556 bp) were obtained. A total of 9,248 significant differentially expressed genes (DEGs) were detected. There were 8,819 DEGs between the parotoid gland and abdomen skin and 1,299 DEGs between the dorsal skin and abdomen skin, while only 1,283 DEGs were obtained between the parotoid gland and dorsal skin. Through enrichment analysis, it was found that the detected differential gene expressions corresponded to the different functions of different skin sites. Our key findings were the genetic expression of toxin secretion, the protection function of skin, and the related genes such as HSD3B, Cyp2c, and CAT, LGALS9. In conclusion, we provide useful transcript resources to study the gene expression and gene function of B. gargarizans and other amphibians. The detected DEGs between different sites of the skin provided better insights into the genetic mechanisms of toxin secretion and the protection function of skin for amphibians.

Effects of perchlorate and exogenous T4 exposures on development, metamorphosis and endochondral ossification in Bufo gargarizans larvae

Several endocrine-disrupting chemicals (EDCs) have been proven to interfere with the physiological function of thyroid hormone (TH), which affected growth and development. However, few studies have investigated the effects of EDCs on TH axis with consequence for skeletal development in amphibians. This study thus examined the potential role of perchlorate and T4 in growth, development and endochondral ossification during metamorphosis of Bufo gargarizans. Our studies showed that NaClO₄ treatment caused weight gain and delayed the developmental stage in B. gargarizans tadpoles, while T4 decreased body size and survival rate, accelerated metamorphic duration and increased the risk of early death. Histological sections suggested that NaClO₄ and T4 treatments caused damages to thyroid tissue, such as decreased thyroid gland size, follicle size, colloid area, the height of follicular epithelial cells and the number of follicles. In addition, the double skeletal staining and RT-qPCR showed that NaClO₄ and T4 treatments inhibited the endochondral ossification by regulating TH synthesis (TRs, Dios) and endochondral ossification-related genes (MMPs, Runxs, VEGFs and VEGFRs) expression levels, which might affect terrestrial locomotion and terrestrial life. Altogether, these thyroid injury and gene expression changes as caused by NaClO₄ and T4 may have an influence on development and endochondral ossification during the metamorphosis of amphibians.

Altered temperature affect body condition and endochondral ossification in Bufo gargarizans tadpoles

Bufo gargarizans is one kind of economic animals with higher medicinal value in China. In this study, B. gargarizans (Bufo gargarizans) tadpoles were reared at three different water temperature (15, 22 and 29 °C) from Gosner stages 28-46. We investigated the effects of temperature on growth, development, survival, metamorphic duration, size and skeletal ossification at Gosner stage 40, 42, and 46, as well as thyroid tissue reached metamorphic climax (Gs42). Besides, we examined the transcription levels of endochondral ossification-related genes in hind limb at metamorphic climax (Gs42). Our results showed that the growth and development of tadpoles conform to the temperature-size rule (TSR). While warm temperature resulted in the decrease in body size and hind limb length, and shorten larval period, cold temperature led to increase in body size and hind limb length but prolonged larval period. Histological examinations revealed that warm and cold temperatures caused damage to thyroid tissue. Also, warm and cold temperatures inhibited the degree of ossification with the double staining methodology. Additionally, the real-time PCR results suggested that warm and cold temperatures significantly up-regulated Runx2, VEGF and VEGFR mRNA levels, and down-regulated TRβ, MMP9, MMP13 and Runx3 mRNA levels. The up-regulation of Dio2 level and down-regulation of Dio3 level were observed in warm temperature. TRα mRNA level was significantly increased in warm temperature, but decreased in cold temperature. Collectively, these observations demonstrated that warm and cold temperatures affected endochondral ossification in B. gargarizans tadpoles, which might influence their capacity to terrestrial locomotion.

Metabolic cold adaptation in the Asiatic toad: intraspecific comparison along an altitudinal gradient

The metabolic cold adaptation (MCA) hypothesis predicts an increase in metabolic rate and thermal sensitivity of poikilotherms from cold environments as compared to those from warm environments, when measured under standardized conditions. This compensatory response is also expected to evolve in life history and behavioral traits if the reductions in these phenotypic traits at low temperature involves in a reduction in fitness. We investigated the extent to which the level of energy intake (measured as feeding rate), energy turnover (measured as standard metabolic rate, SMR) and the energy budget (energy allocation to growth and physical activity) are influenced by climatic conditions in three populations of the Asiatic toad (Bufo gargarizans) distributed across an altitudinal gradient of 1350 m in the Qionglai Mountains of Western China. We found a similar thermal reaction norm of SMR at both population and individual levels; therefore, the data did not support the MCA hypothesis. However, there was a co-gradient variation (CoGV) for mass change rate in which the high and medium altitudinal populations displayed slower mass change rates than their counterparts from low altitudes. Moreover, this CoGV pattern was accompanied by a low feeding rate and high physical activity for the high- and medium-altitude populations. Our results highlight that adjustments in energy intake and energy allocation to behaviors, but not energy allocation to metabolism of maintenance, could act as an energetic strategy to accommodate the varied growth efficiency in Asiatic toads along an altitudinal gradient.

Intestinal response characteristic and potential microbial dysbiosis in digestive tract of Bufo gargarizans after exposure to cadmium and lead, alone or combined

The gastrointestinal tract is the largest immune organ in the body and meanwhile, accommodates a large number of microorganisms. Heavy metals could disturb the intestinal homeostasis and change the gut microbial composition. However, the information regarding the links between dysbiosis of gut microbiota and imbalance of host intestinal homeostasis induced by the mixture of heavy metals is insufficient. The present study investigates the effects of Cd/Pb, both single and combination exposure, on the growth performance, intestinal histology, digestive enzymes activity, oxidative stress and immune parameters, and intestinal microbiota in Bufo gargarizans tadpoles. Our results revealed that co-exposure of Cd-Pb induced more severe impacts not only on the host, but the intestinal microbiota. On the one hand, co-exposure of Cd-Pb significantly induced growth retardation, intestinal histological injury, decreased activities of digestive enzymes. On the other hand, Cd and Pb exposure, especially in mixed form, changed the diversity and richness, structure of microbiota. Also, the intestinal microbial composition was altered by Cd/Pb exposure (alone and combination) both at the different levels. Proteobacteria, act as front-line responder, was significantly increased in tadpoles under the exposure of metals. Finally, the functional prediction revealed that the disorders of metabolism and immune responses of intestinal microbiota was increased in tadpoles exposed to Cd/Pb (especially the mixture of Cd and Pb). Our research complements the understanding of links between changes in host fitness loss and intestinal microbiota and will add a new dimension of knowledge to the ecological risks of mixed heavy metals in amphibian.

Physiological, developmental, and behavioral plasticity in response to thermal acclimation

Organisms with complex ecologies and life-cycle processes may shift physiologically (acclimation in tolerance), developmentally, and/or behaviorally (thermoregulation) in response to changes in climate. As such, climate change may trigger multiple, interacting phenotypic responses, which underscores the nuances of characterizing a species capacity to adapt and respond to climate change. In this study, we use a model frog species, Bufo gargarizans, to examine how three phenotypes, thermal tolerance limits (critical thermal minimum, CT_min and critical thermal maximum, CT_max), ontogeny, and behavioral preferences in temperature (T_pref) respond to different levels of thermal exposure (i.e., acclimation ranging from 10 °C to 30 °C). Acclimation temperature had little effect on T_pref of tadpoles, yet behaviorally they showed strong signs of thermal selection towards an optimum. Both CT_min and CT_max increased with acclimation temperature with an approximate 10% increase in tolerance limits per 1 °C increase in exposure. Development and body size both responded to acclimation temperature, both of which also influenced lower but not upper thermal limits. Our study highlights the idiosyncrasies of estimating climate vulnerability, where multiple phenotypes can respond to shifts in temperature-a complexity that is especially apparent in species with complex life-cycles.

Effects of temperature on growth, development and the leptin signaling pathway of Bufo gargarizans

Climate change is one of the most important causes of the decline in amphibians. Changes in temperature have an important effect on the growth and development and energy metabolism of amphibians. The aim of this study is to unravel the effects of temperature on the leptin signaling pathway of Bufo gargarizans and its molecular mechanisms. Our results showed that high temperature accelerated the development rate of tadpoles, but reduced body size and mass, while low temperature deferred the development of tadpoles, but increased size and mass. Both high temperature and low temperature exposure caused pathological damage of the liver in B. gargarizans. The results of RT-qPCR revealed that the high temperature treatment significantly upregulated the transcript levels of genes related to thyroid hormone (DIO2 (D2), Thyroid Hormone Receptor-α (TRα)) and the leptin signaling pathway (Leptin Receptor (LepR), Janus kinase 1 (JAK1), Janus kinase 2 (JAK2), Tyrosine kinase 2 (TYK2), Signal Transducer And Activator Of Transcription 3 (STAT3), Signal Transducer And Activator Of Transcription 3.1 (STAT3.1), and Signal Transducer And Activator Of Transcription 6 (STAT6)), while there was a decrease of mRNA expression of these genes (TRα, Thyroid Hormone Receptor-Beta (TRβ), LepR, JAK1, and TYK2) in the liver of tadpoles exposed to high temperature compared with the intermediate temperature treatment. Therefore, our results suggested that temperature extremes might interfere with the thyroid and leptin signaling pathways and affect the growth and development of B. gargarizans. Furthermore, tissue injury of the liver could occur due to exposure to temperature extremes. This work promotes public awareness of environmental protection and species conservation needs, also provides valuable experimental data and a theoretical basis for the protection of amphibians.

Effects of sodium perchlorate and exogenous L-thyroxine on growth, development and leptin signaling pathway of Bufo gargarizans tadpoles during metamorphosis

Sodium perchlorate (NaClO₄) and exogenous L-thyroxine (T4), two kinds of endocrine-disrupting chemicals (EDCs), mainly affect the circulating thyroid hormones, which regulate the initiation and rate of metamorphosis in amphibian. The aim of this study is to evaluate the potential role of EDCs in regulating the development of tadpoles and leptin signaling pathway of liver during the metamorphosis of Bufo gargarizans. There was completely opposite result of average development stage of tadpoles and morphological parameters between the NaClO₄ and T4 exposure groups. Histological analysis revealed that NaClO₄ and T4 exposure both caused liver injury, such as the decreased size of hepatocytes, atrophy of nucleus, increased melanomacrophage centres and disappearance of hepatocyte membranes. In addition, the results of RT-qPCR revealed that NaClO₄ treatment significantly inhibited the transcript levels of genes related to thyroid hormone (D2, TRα and TRβ) and leptin signaling pathway (LepR, JAK1, JAK2, and TYK2), while there was an increase of mRNA expression of these genes in the liver of tadpoles administrated with T4 compared with control. This work lays an important foundation for assessing the risk of EDCs in relation to amphibian development during metamorphosis.

Nitrate exposure induces intestinal microbiota dysbiosis and metabolism disorder in Bufo gargarizans tadpoles

Excess nitrate has been reported to be associated with many adverse effects in humans and experimental animals. However, there is a paucity of information of the effects of nitrate on intestinal microbial community. In this study, the effects of nitrate on development, intestinal microbial community, and metabolites of Bufo gargarizans tadpoles were investigated. B. gargarizans were exposed to control, 5, 20 and 100 mg/L nitrate-nitrogen (NO₃-N) from eggs to Gosner stage 38. Our data showed that the body size of tadpoles significantly decreased in the 20 and 100 mg/L NO₃-N treatment group when compared to control tadpoles. Exposure to 20 and 100 mg/L NO₃-N also caused indistinct cell boundaries and nuclear pyknosis of mucosal epithelial cells in intestine of tadpoles. In addition, exposure to NO₃-N significantly altered the intestinal microbiota diversity and structure. The facultative anaerobic Proteobacteria occupy the niche of the obligately anaerobic Bacteroidetes and Fusobacteria under the pressure of NO₃-N exposure. According to the results of functional prediction, NO₃-N exposure affected the fatty acid metabolism pathway and amino acid metabolism pathway. The whole-body fatty acid components were found to be changed after exposure to 100 mg/L NO₃-N. Therefore, we concluded that exposure to 20 and 100 mg/L NO₃-N could induce deficient nutrient absorption in intestine, resulting in malnutrition of B. gargarizans tadpoles. High levels of NO₃-N could also change the intestinal microbial communities, causing dysregulation of fatty acid metabolism and amino acid metabolism in B. gargarizans tadpoles.

Probing the effects of hexavalent chromium exposure on histology and fatty acid metabolism in liver of Bufo gargarizans tadpoles

Hexavalent chromium is one of the major detrimental heavy metal pollutants. B. gargarizans tadpoles were treated with different concentrations of Cr⁶⁺ (0, 13, 52, 104, 208, and 416 μg Cr⁶⁺ L^-1) from Gs 2 to Gs 42. The effect of Cr⁶⁺ on histopathological alterations and transcript levels of fatty acid metabolism-related genes as well as fatty acids composition and content in liver were examined. Histopathological changes were observed in liver at 52, 104, 208, and 416 μg Cr⁶⁺ L^-1. Moreover, RT-qPCR analyses showed the downregulated mRNA levels of the genes related to fatty acid synthesis (SCD, MECR, TECR and ELOVL1) and fatty acid β-oxidation (ACOT1, PPT1, HADH and ACAA2) at 416 μg Cr⁶⁺ L^-1. However, the mRNA expression of fatty acid β-oxidation-related genes (ECHS1, HADHA and ACAA1) were significantly upregulated at 13, 52, 104, 208 and 416 μg Cr⁶⁺ L^-1. In situ hybridization revealed BSEP was expressed in hepatocyte nucleus and plasma membrane, and HSD17B12 was abundantly expressed in the plasma membrane. The HSD17B12 mRNA levels were significantly upregulated in tadpoles exposed to all Cr⁶⁺ treatment groups, while the BSEP mRNA levels were downregulated at 104, 208 and 416 μg Cr⁶⁺ L^-1 groups compared to control. In addition, an increase in polyunsaturated fatty acids and a decrease in monounsaturated fatty acids were found in 52, 104 and 416 μg Cr⁶⁺ L^-1 groups. Overall, chronic exposure to Cr⁶⁺ may suppress fatty acid synthesis, disturb fatty acid β-oxidation, aggravate disorders of hepatic function and induce hepatic impairment in B. gargarizans tadpoles.

Changes in the community structure of the symbiotic microbes of wild amphibians from the eastern edge of the Tibetan Plateau

Environment has a potential effect on the animal symbiotic microbiome. Here, to study the potential relationship of the symbiotic microbiomes of wild amphibians with altitude, we collected the gut and skin samples from frogs (nine species) and the environmental samples (water and soil samples) from the Leshan Mountains (altitude: 360–410 m) and Gongga Mountains (altitude: 3340–3989 m) on the eastern edge of the Tibetan Plateau. Bufo gargarizans (Bg) samples were collected from both the Leshan and Gongga mountain regions (Bg was the only species sampled on both mountains). The DNA extracted from each sample was performed high‐throughput sequencing (MiSeq) of bacterial 16S rRNA gene amplicons. High relative abundance of Caulobacteraceae and Sphingomonadaceae was found in skin samples from both Bg and the other high‐altitude amphibians (nine species combined). High relative abundance of Coxiellaceae and Mycoplasmataceae was found in gut samples from both Bg and the other high‐altitude amphibians. Furthermore, the alpha and beta diversities of skin and gut samples from Bg and the other amphibian species (nine species combined) were similar. In terms of the symbiotic microbial community, the low‐altitude samples were less diverse and more similar to each other than the high‐altitude samples were. We speculated that extreme high‐altitude environments and host phylogeny may affect the amphibian microbiome. Despite the distinct microbial community differences between the skin and gut microbiomes, some functions were similar in the Bg and combined high‐altitude samples. The Bg and high‐altitude skin samples had higher oxidative stress tolerance and biofilm formation than the low‐altitude skin samples. However, the opposite results were observed for the Bg and high‐altitude gut samples. Further study is required to determine whether these characteristics favor high‐altitude amphibian adaptation to extreme environments.

Metabolites from Bufo gargarizans (Cantor, 1842): A review of traditional uses, pharmacological activity, toxicity and quality control

ETHNOPHARMACOLOGICAL RELEVANCE

Bufo gargarizans (Cantor, 1842) (BGC), a traditional medicinal animal distributed in many provinces of China, is well known for the pharmaceutical value of Chansu and Chanpi. As traditional Chinese medicines (TCMs), Chansu and Chanpi, with their broad-spectrum of therapeutic applications, have long been applied to detoxification, anti-inflammation, analgesia, etc. OVERARCHING OBJECTIVE: We critically analyzed the current evidence for the traditional uses, chemical profiles, pharmacological activity, toxicity and quality control of BGC (Bufonidae family) to provide a scientific basis for future in-depth studies and perspectives for the discovery of potential drug candidates.

METHODOLOGY

All of the available information on active constituents and TCMs derived from BGC was obtained using the keywords "Bufo gargarizans", "Chansu", "Chanpi", "Huachansu", or "Cinobufacini" through different electronic databases, including PubMed, Web of Science, Chinese National Knowledge Infrastructure (CNKI), the Wanfang Database, and Pharmacopoeia of China. In addition, Chinese medicine books from different times were used to elucidate the traditional uses of BGC. Electronic databases, including the "IUCN Red List of Threatened Species", "American Museum of Natural History" and "AmphibiaWeb Species Lists", were used to validate the scientific name of BGC.

RESULTS

To date, about 118 bufadienolide monomers and 11 indole alkaloids have been identified from BGC in total. The extracts and isolated compounds exhibit a wide range of in vitro and in vivo pharmacological effects. The literature search demonstrated that the ethnomedicinal uses of BGC, such as detoxification, anti-inflammation and the ability to reduce swelling and pain associated with infections, are correlated with its modern pharmacological activities, including antitumor, immunomodulation and attenuation of cancer-derived pain. Bufadienolides and indole alkaloids have been regarded as the main active substances in BGC, among which bufadienolides have significant antitumor activity. Furthermore, the cardiotoxicity of bufadienolides was discussed, and the main molecular mechanism involves in the inhibition of Na⁺/K⁺-ATPase. Besides, with the development of modern analytical techniques, the quality control methods of BGC-derived TCMs are being improved constantly.

CONCLUSIONS

An increasing number of reports suggest that BGC can be regarded as an excellent source for exploring the potential antitumor constituents. However, the future antitumor research of BGC needs to follow the standard pharmacology guidelines, so as to provide comprehensive pharmacological information and aid the reproducibility of the data. Besides, to ensure the efficacy and safety of BGC-derived TCMs, it is vital to construct a comprehensive quality evaluation model on the basis of clarifying pharmacodynamic-related and toxicity-related compositions.

Chronic effects of lead on metamorphosis, development of thyroid gland, and skeletal ossification in Bufo gargarizans

We examined the Pb²⁺ exposure on tadpoles of Bufo gargarizans from Gosner stage 26-42. Mortality, growth and development, time to metamorphosis, size, and skeletal ossification at metamorphic climax of Bufo gargarizans were examined. Also, histological characteristics of thyroid glands in tadpoles at Gosner stage 33, 38, and 42 as well as transcript levels of thyroid hormone-related genes in the hind-limb, tail, and liver of tadpoles at metamorphic climax were examined. Pb²⁺ exposure induced mortality in a concentration-dependent manner in Bufo gargarizans larvae. The significant increase in growth and development, percent metamorphosis, size at metamorphic climax, and skeletal ossification were observed at 50 μg Pb²⁺ L^-1; however, exposure to 1000 μg Pb²⁺ L^-1 resulted in the opposite effects in tadpoles. In addition, histological alterations of thyroid gland, such as follicular cell hyperplasia and colloid depletion could be found in 50-1000 μg Pb²⁺ L^-1 treatments. Furthermore, Pb²⁺ exposure at 1000 μg L^-1 resulted in significantly decreased transcript levels of Dio2, TRα and TRβ, and increased transcript levels of Dio3. In contrast, 50 μg Pb²⁺ L^-1 significantly upregulated the mRNA levels of Dio2, TRα, and TRβ, but it reduced the Dio3 expression. These results suggested that Pb²⁺ might disrupt TH homeostasis in tadpoles by histological alterations of thyroid gland and disturb the transcript levels of Dio2, Dio3, TRα, and TRβ, leading to altered growth and development, as well as percent metamorphosis and skeletal ossification. Further studies are needed to elucidate the underlying mechanisms of low-dose stimulation and high-dose inhibition effects.

Effects of fluoride on the histology, lipid metabolism, and bile acid secretion in liver of Bufo gargarizans larvae

In our study, Bufo gargarizans (B. gargarizans) larvae were exposed to control, 0.5, 5, 10 and 50 mg/L of NaF from Gs 26 to 42. At Gs 42, we evaluated the changes of liver histology and the mRNA levels of target genes in liver. In addition, we also examined the composition and content of fatty acids. Histological analysis revealed that fluoride caused liver injury, such as the increase of number of melanomacrophage centres, atrophy of nucleus, dilation of bile canaliculus, and decrease of quantity, degradation and deposition of lipid droplets. The results of RT-qPCR indicated that exposure to 5, 10 and 50 mg/L of NaF significantly decreased the transcript levels of genes related to fatty acid synthesis (FASN, FAE, MECR, KAR and TECR) in liver. Besides, mRNA expression of genes involved in fatty acid β-oxidation (ECHS1, HADHA, SCP2, CPT2, ACAA1 and ACAA2) and oxidative stress (SOD, GPx, MICU1 and HSP90) was significantly downregulated in 0.5, 5, 10 and 50 mg/L of NaF treatment groups. Also, in the relative expression of genes associated with synthesis and secretion of bile acid, BSEP significantly increased at 0.5, 5 and 50 mg/L of NaF while HSD3B7 significantly reduced in 0.5, 5, 10 and 50 mg/L of NaF. Finally, the fatty acid extraction and GC-MS analysis showed that the content of saturated fatty acids (SFAs) was decreased and the content of polyunsaturated fatty acids (PUFAs) was increased in all fluoride treatment groups. Taken together, the present results indicated that fluoride-induced the histological alterations of liver might be linked to the disorder of lipid metabolism, oxidative damage.

Effects of nitrate on development and thyroid hormone signaling pathway during Bufo gargarizans embryogenesis

Nitrate is known to disrupt the thyroid hormone, which is essential for the metamorphosis of amphibians. However, few studies are focused on the effects of nitrate on the maternal thyroid hormone in early amphibian embryos. We aimed to determine the impact of nitrate on maternal thyroid hormone signaling pathway in Bufo gargarizans embryos. B. gargarizans embryos were exposed to different concentrations of nitrate-nitrogen (NO3-N) for 7 days. High concentration of NO3-N (50, 100, and 200 mg/L) could induce embryonic malformation and influence the development of embryos. In addition, maternal T4 and components of the thyroid hormone (TH) signaling pathway were detected by ELISA and RNA-seq, respectively. The expression levels of mRNA related to thyroid hormone and oxidative stress were affected in the early developing embryos in all NO3-N treatment groups. However, the T4 levels and the spatial expression patterns of type II iodothyronine deiodinase (D2), type III iodothyronine deiodinase (D3), thyroid hormone receptor α (TRα), and thyroid hormone receptor β (TRβ) mRNA were not changed by nitrate. In conclusion, the results of our study highlight the crucial role of the maternal thyroid hormone signaling pathway in normal embryonic development, and the adverse effects of nitrate on the expression levels of mRNA related to thyroid hormone signaling pathway and oxidative stress in B. gargarizans embryos.

Immunohistochemical distribution of glucagon-like peptide 1 in the digestive system of different aquatic vertebrates

The distribution of glucagon-like peptide 1 (GLP-1)-positive cells in digestive tracts and pancreases of aquatic vertebrates was investigated by immunohistochemical staining method. The results suggested that GLP-1-positive cells were distributed in the columnar mucous epithelium and tubular glands of lamina propria in the digestive system. However, GLP-1-positive cells were also found in subepithelial lamina propria of the mucosae and muscularis in each segment of the digestive tract of Rana nigromaculata. The distribution densities of these cells reached peaks in the stomachs, and the middle or end segments of small intestines of Chinese softshell turtle, Bufo gargarizans, R. nigromaculata and catfish, and there was the third distribution density peak in the rectum of catfish. The total amount or overall density of GLP-1-positive cells varied a lot in the digestive tracts of different animal species. The distribution density was relatively low in the digestive tract of chub and reached the maximum in the digestive tracts of snakehead and catfish, but no GLP-1-positive cells were found in the digestive tract of bighead carp. GLP-1-positive cells were densely distributed in the pancreases of Chinese softshell turtle, B. gargarizans and R. nigromaculata. These cells spread over the superficial layers of islets or scattered in exocrine pancreas in the pancreas of B. gargarizans, spread in the endocrine cells or scattered in the pancreas of Chinese softshell turtle, scattered in the pancreas of R. nigromaculata and distributed in the superficial layers of islets in the pancreas of catfish.

Exposure to nitrate alters the histopathology and gene expression in the liver of Bufo gargarizans tadpoles

Nitrate is a ubiquitous component in aquatic environment and the concentrations of anthropogenic nitrate-nitrogen (NO₃N) can exceed 25 mg/L in surface waters and 100 mg/L in ground waters. The exceed nitrate has adverse effects on survival, development, and metamorphosis of amphibian. Liver is the hub of many biological processes, including lipid metabolism and bile salts secretion. However, there is little information about the effects of nitrate on the liver in amphibians during metamorphosis. In this study, B. gargarizans was exposed to different concentrations of nitrate from embryo to metamorphosis climax to investigate the effects of nitrate on the liver. The survival rate, metamorphosis percent, body mass, total length, and hind-limb length were measured. The histopathological changes and transcriptome responses in the liver of B. gargarizans to nitrate were examined. Results indicated exposure to 50 and 100 mg/L NO₃N delayed the metamorphosis and decreased the metamorphosis percent of B. gargarizans. The body size of B. gargarizans at 10 and 50 mg/L NO₃N groups were decreased while it was increased at 100 mg/L NO₃N group. In addition, exposure to 100 mg/L NO₃N caused severe histopathological changes, including cellular atrophy, increased intercellular areas, degraded lipid droplets, hepatic fibrosis, bile canaliculus contraction and degraded mitochondria in liver. The results of RNA-seq and qRT-PCR interpreted the molecular responses, which might be the factors to induce histopathological changes in the liver of B. gargarizans under the pressure of nitrate exposure.

Modulation of Gene Expression in Liver of Hibernating Asiatic Toads (Bufo gargarizans)

Hibernation is an effective energy conservation strategy that has been widely adopted by animals to cope with unpredictable environmental conditions. The liver, in particular, plays an important role in adaptive metabolic adjustment during hibernation. Mammalian studies have revealed that many genes involved in metabolism are differentially expressed during the hibernation period. However, the differentiation in global gene expression between active and torpid states in amphibians remains largely unknown. We analyzed gene expression in the liver of active and torpid Asiatic toads (Bufo gargarizans) using RNA-sequencing. In addition, we evaluated the differential expression of genes between females and males. A total of 1399 genes were identified as differentially expressed between active and torpid females. Of these, the expressions of 395 genes were significantly elevated in torpid females and involved genes responding to stresses, as well as contractile proteins. The expression of 1004 genes were significantly down-regulated in torpid females, most which were involved in metabolic depression and shifts in the energy utilization. Of the 715 differentially expressed genes between active and torpid males, 337 were up-regulated and 378 down-regulated. A total of 695 genes were differentially expressed between active females and males, of which 655 genes were significantly down-regulated in males. Similarly, 374 differentially expressed genes were identified between torpid females and males, with the expression of 252 genes (mostly contractile proteins) being significantly down-regulated in males. Our findings suggest that expression of many genes in the liver of B. gargarizans are down-regulated during hibernation. Furthermore, there are marked sex differences in the levels of gene expression, with females showing elevated levels of gene expression as compared to males, as well as more marked down-regulation of gene-expression in torpid males than females.

Transcriptome analyses reveal molecular mechanisms that regulate endochondral ossification in amphibian Bufo gargarizans during metamorphosis

BACKGROUND

A developmental transition from aquatic to terrestrial existence is one of the most important events in the evolution of terrestrial vertebrates. Amphibian metamorphosis is a classic model to study this transition. The development of the vertebrate skeleton can reflect its evolutionary history. Endochondral ossification serves a vital role in skeletal development. Thus, we sought to unravel molecular mechanisms that regulate endochondral ossification during Bufo gargarizans metamorphosis.

METHODS

The alizarin red-alcian blue double staining method was used to visualize the skeletal development of B. gargarizans during metamorphosis. RNA sequencing (RNA-seq) was used to explore the transcriptome of B. gargarizans in four key developmental stages during metamorphosis. Real-time quantitative PCR (RT-qPCR) was used to validate the expression patterns of endochondral ossification related genes.

RESULTS

Endochondral ossification increased gradually in skeletal system of B. gargarizans during metamorphosis. A total of 137,264 unigenes were assembled and 44,035 unigenes were annotated. 10,352 differentially expressed genes (DEGs) were further extracted among four key developmental stages. In addition, 28 endochondral ossification related genes were found by searching for DEG libraries in B. gargarizans. Of the 28 genes, 10 genes were validated using RT-qPCR.

CONCLUSIONS

The exquisite coordination of the 28 genes is essential for regulation of endochondral ossification during B. gargarizans metamorphosis.

GENERAL SIGNIFICANCE

The present study will not only provide an invaluable genomic resource and background for further research of endochondral ossification in amphibians but will also aid in enhancing our understanding of the evolution of terrestrial vertebrates.

Transcriptomics provides mechanistic indicators of fluoride toxicology on endochondral ossification in the hind limb of Bufo gargarizans

Endochondral ossification, the process by which most of the bone is formed, is regulated by many specific groups of molecules and extracellular matrix components. Hind limb of Bufo gargarizans is a model to study endochondral ossification during metamorphosis. Chinese toad (Bufo gargarizans) were exposed to different fluoride concentrations (0, 1, 5, 10 and 20 mg L^-1) from G3 to G42. The development of hind limb of B. gargarizans was observed using the double staining methodology. The transcriptome of hind limb of B. gargarizans was conducted using RNA-seq approach, and differentially expressed gene was also validated. In addition, the location of Sox9 and Ihh in the growth cartilage was determined using in situ hybridization. Our results showed that 5 mg L^-1 stimulated bone mineralization, while 10 and 20 mg L^-1 exposure could inhibit the tibio-fibula, tarsus and metacarpals ossification. Besides, 10 mg F/L treatment could down-regulate Ihh, Sox9, D2, D3, TRα, TRβ, Wnt10, FGF3 and BMP6 expression, while up-regulate ObRb and HHAT mRNA expression in the hind limb of B. gargarizans. Transcript level changes of Ihh, Sox9, D2, D3, TRα, TRβ, Wnt10, FGF3 and BMP6 were consistent with the results of RT-qPCR. In situ hybridization revealed that Ihh was expressed in prehypertrophic chondrocytes, while Sox9 was abundantly expressed in proliferous, prehypertrophic and hypertrophic chondrocytes. However, 10 mg F-/L did not cause any affect in the location of the Ihh and Sox9 mRNA. Therefore, high concentration of fluoride could affect the ossification-related genes mRNA expression and then inhibit the endochondral ossification. The present study thus will greatly contribute to our understanding of the effect of environmental contaminant on ossification in amphibian.

Changes in intestinal microbiota of Bufo gargarizans and its association with body weight during metamorphosis

The assembly of intestinal microbial communities can play major roles in animal development. We hypothesized that intestinal microbial communities could mirror the developmental programs of amphibian metamorphosis. Here, we surveyed the morphological parameters of the body and intestine of Bufo gargarizans at varying developmental stages and inventoried the intestinal microbial communities of B. gargarizans at four key developmental stages via 16S rDNA gene sequencing. Firstly, our survey showed that during metamorphosis, body weight and intestinal weight were reduced by 56.8 and 91.8%, respectively. Secondly, the gut bacterial diversity of B. gargarizans decreased with metamorphosis and the composition of the tadpoles' intestinal microbiota varied across metamorphosis. Compared to aquatic larvae, terrestrial juveniles showed major shifts in microbial composition, including reduction in Proteobacteria and Actinobacteria, increases in Bacteroidetes and Fusobacteria, and the appearance of Verrucomicrobia. Firmicutes in four developmental stages showed similar abundance at the phylum level, but in each stage was driven by distinct genera. Enterobacter, Aeromonas, Mucinivorans and Bacteroides also changed in abundance and were found to be significantly correlated with loss of body or intestinal tissue during metamorphosis. These results indicate a shift in intestinal microbial community composition throughout amphibian metamorphosis.

The morphological changes and molecular biomarker responses in the liver of fluoride-exposed Bufo gargarizans larvae

The goal of the current study was to evaluate the negative influences of fluoride on liver of Bufo gargarizans larvae. B. gargarizans larvae were treated with 42.4mgF^-/L for 0, 24, 48 and 72h at Gosner stage 37. The morphological changes and responses of molecular biomarkers involved in lipid metabolism, oxidative stress and apoptosis were examined in liver. Disappearance of cell boundaries, degeneration of hepatic parenchyma cells and significant increase in the number of melanomacrophage centres and the quantity of lipid droplets were found in the liver treated with 42.4mgF^-/L for 72h. In addition, in the relative expression of acetyl CoA carboxylase 1 (ACC-1), fatty acid elongase 1 (FAE-1), sterol carrier protein 2 (SCP-2), and carnitine palmitoyltransferase-1 (CPT-1), decrease was observed after 24, 48 and 72h of 42.4mgF^-/L exposure. Furthermore, the transcript levels of superoxide dismutase (SOD) and glutathione peroxidase (GPx) were downregulated in tadpoles exposed for 24, 48 and 72h to 42.4mgF^-/L, while the transcript level of heat shock protein 90 (HSP90) was upregulated at 42.4mgF^-/L for 72h. Also, mRNA expression of Bcl-2-associated transcription factor 1(BCLAF1) and thyroid hormone receptors (TRα and TRβ) was significantly upregulated in tadpoles treated with 42.4mgF^-/L for 72h. Therefore, our results suggested that the liver injury induced by fluoride might result from disruption of lipid metabolism, oxidative damage and apoptosis.

Effects of fluoride on morphology, growth, development, and thyroid hormone of Chinese toad (Bufo gargarizans) embryos

Excessive fluoride in natural water ecosystem has the potential to detrimentally affect amphibians, but little is known of such effects or underlying mechanisms in Bufo gargarizans embryos. In the present study, the effects of fluoride exposure on B. gargarizans embryos were investigated. First, fluoride teratogenic experiment showed that the 9 days EC₅₀ of fluoride on B. gargarizans embryos was 177.62 mg/L. Then, we studied the sublethal effects of fluoride on B. gargarizans embryos at control, 0.7, 4.1, 19.6, 41.9, and 62.7 mg/L fluoride concentration. Malformation, growth, and development of embryos were monitored, and type 2 and 3 iodothyronine deiodinase (Dio2 and Dio3), thyroid hormone receptors (TRα and TRβ) mRNA levels were measured. Our results showed the morphological malformations, such as tail curvature (lordosis), edema, cuticularized ciliated cells, and hyperplasia were occurred during fluoride exposure. Growth and development were all inhibited at 19.5, 41.9, and 62.7 mg/L fluoride-treated groups after 9 days' exposure. According to real-time PCR results, exposure to fluoride upregulated Dio3 and TRβ mRNA expression and downregulated Dio2 and TRα mRNA level. All above indicated that excessive fluoride could induce morphology malformations, inhibit embryonic growth and development, and disrupt the normal function of maternal thyroid hormone in B. gargarizans embryos. Environ. Mol. Mutagen. 59:123-133, 2018. © 2017 Wiley Periodicals, Inc.

Responses of growth, malformation, and thyroid hormone-dependent genes expression in Bufo gargarizans embryos following chronic exposure to Pb2

The aim of this study was to examine the adverse effects of lead (Pb) exposure on Bufo gargarizans embryos. The 96 h-LC₅₀ of Pb²⁺ for B. gargarizans embryos was determined to be 26.6 mg L^-1 after an acute test. In the chronic test, B. gargarizans embryos at Gosner stage 3 were exposed to 10~2000 μg Pb²⁺ L^-1 during embryogenesis. Total length, weight, developmental stage, and malformation were monitored. In addition, the transcript levels of type II and type III iodothyronine deiodinase (Dio2 and Dio3) and thyroid hormone receptors (TRα and TRβ) were determined to assess the thyroid-disrupting effects of Pb²⁺. Slightly increased growth and development of B. gargarizans embryos were observed at low concentrations of Pb²⁺ (10, 50, and 100 μg L^-1), while retarded growth and development were found at high concentrations of Pb²⁺ (1000 and 2000 μg L^-1). In addition, Pb²⁺ exposure induced morphological abnormalities, which were characterized by edema at tail, wavy fin, abdominal edema, stunted growth, hyperplasia, and axial flexures in B. gargarizans embryos. Furthermore, our results showed that exposure to 2000 μg Pb²⁺ L^-1 decreased the transcript levels of Dio2, TRα, and TRβ, but it increased Dio3 mRNA level. In contrast, exposure to 50 μg Pb²⁺ L^-1 increased TRα mRNA level and decreased Dio3 mRNA level. These results suggested that Pb²⁺ might have thyroid-disrupting effects, leading to the disruption of growth and development in B. gargarizans embryos.

Phylogeographic and population insights of the Asian common toad (Bufo gargarizans) in Korea and China: population isolation and expansions as response to the ice ages

The effects of ice ages on speciation have been well documented for many European and North American taxa. In contrast, very few studies have addressed the consequences of such environmental and topographical changes in North East Asian species. More precisely, the Korean Peninsula offers a unique model to assess patterns and processes of speciation as it hosts the northern- and eastern-most distribution limit of some widespread Asian taxa. Despite this, studies addressing phylogeographic patterns and population genetics in the peninsula and surrounding countries are few and studies for most families are lacking. Here we inferred the phylogenetic relationships of the common toad (Bufo gargarizans) from South Korea and their North East Asian counterpart populations, based on mitochondrial data. Korean B. gargarizans GenBank BLASTs matched few individuals from nearby China, but the presence of a Korean clade suggests isolation on the Korean Peninsula, previous to the last glacial maximum, linked to sea level resurgence. Molecular clock calibrations within this group were used to date the divergence between clades and their relationship to paleo-climatic events in the area. Lack of genetic structure among South Korean populations and strong homogeneity between the Korean and some Chinese localities suggest weak isolation and recent expansion. Geographical projection of continuous coalescent maximum-clade-credibility trees shows an original Chinese expansion towards the Korean Peninsula through the Yellow Sea circa two million years ago with colonisation events dating circa 800 thousand years ago (K. y. a.). Following this colonisation, the data point to outgoing Korean Peninsula dispersal events throughout different periods, towards the North through land, and West through land bridge formations over the Yellow Sea during sea level falls. In accordance, demographic analyses revealed a population expansion in the Koran Peninsula circa 300 K. y. a., likely attributed to glacial cycle fluctuations.

Toxic effects of NH4+-N on embryonic development of Bufo gargarizans and Rana chensinensis

Although nitrogen fertilizer is commonly used worldwide, little information is currently available about NH₄⁺-N toxicity on amphibians. This study determined the acute and chronic toxic effects of NH₄⁺-N on two native Chinese amphibian species (Bufo gargarizans and Rana chensinensis), and compared the negative sensitivity of different embryos to NH₄⁺-N. Static renewal aqueous exposures were performed using B. gargarizans and R. chensinensis embryos at Gosner stage 2 over 96 h. In terms of 96 h-LC₅₀, B. gargarizans and R. chensinensis embryos had significantly different responses to NH₄⁺-N, and the latter was more sensitive to NH₄⁺-N than the former. In the chronic toxicity test, exposure to 10 mg L^-1NH₄⁺-N or higher significantly decreased the hatching rate of embryos in both species. Significant increases in the abnormality rate of embryos at 50 mg L^-1NH₄⁺-N or higher were observed and morphological abnormalities were characterized by axial flexures, yolk sac edema, and hyperplasia in both species. Additionally, the total length of embryos decreased in a dose-dependent manner after exposure to NH₄⁺-N. The results indicate that NH₄⁺-N exposure can increase abnormality and inhibit the hatching and development of embryos in B. gargarizans and R. chensinensis.

Histopathological changes and lipid metabolism in the liver of Bufo gargarizans tadpoles exposed to Triclosan

In the current study, the adverse effects of TCS on liver health of B. gargarizans tadpoles were assessed. B. gargarizans larvae were exposed to TCS at 0, 10, 30, 60, and 150 μg L^-1 from Gosner stage 3 until metamorphic climax. The hepatosomatic index (HSI), hepatic histological and ultrastructural features, and transcript levels of genes associated with detoxification and oxidative stress as well as lipid metabolism in the livers were determined. Exposure to 150 μg L^-1 TCS resulted in increased HSI of tadpoles at metamorphic climax. Histological changes characterized by an increase in the number of melanomacrophage, nucleus pyknosis, and deposition of collagen fibers were observed in liver at 60 and 150 μg L^-1 TCS. Moreover, marked ultrastructural alterations including high electron dense in mitochondrial matrix and lipid accumulation were also observed. In addition, abundances of transcripts of Cu/Zn superoxide dismutase (SOD), phospholipid hydroperoxide glutathione peroxidase (PHGPx), and heat shock protein 90 (HSP₉₀) were decreased in larvae exposed to 60 and 150 μg L^-1 TCS, while transcript level of HSP₉₀ was increased at 30 μg L^-1 TCS. Also, abundances of transcripts of acetyl-CoA carboxylase (ACC), carnitine palmitoyltransferase 2 (CPT2), peroxisome proliferator-activated receptor alpha (PPARa), fatty acid elongase 1 (FAE), sterol carrier protein 2 (SCP) were significantly lesser in larvae exposed to 60 and 150 μg L^-1 TCS. Overall, TCS at high levels induced histopathological changes in the liver of B. gargarizans tadpoles. This might have been due to the alteration of oxidative stress-related genes and lipid metabolism-related genes expression levels.

Histological changes, lipid metabolism and oxidative stress in the liver of Bufo gargarizans exposed to cadmium concentrations

Chinese toad (Bufo gargarizans) were exposed to different concentrations of cadmium (5, 50, 100, 200 and 500 μg Cd L^-1) from Gosner stage 3-42. Metamorphosis rate, body weight, total length and body length were measured. Histological alterations in thyroid gland and liver were examined. Changes in hepatocyte were also examined using Transmission electron microscopic. In addition, the mRNA expression of several genes involved in lipid metabolism, oxidative stress and thyroid hormones signaling pathways were also measured. Our results showed that 200 and 500 μg Cd L^-1 decreased the metamorphosis rate and inhibited the body size of B. gargarizans larvae at G42. Moreover, histological examinations have clearly exhibited that cadmium caused liver damage. Ultrastructural examination revealed lipid accumulation and abnormal mitochondria. Exposure to 200 and 500 μg Cd L^-1 significantly up-regulated mRNA expression of D2, SOD, GPx, ACC and FAE, but down-regulated mRNA expression of TRα, TRβ, PPARα, ACOX, CPT and SCP. However, low Cd concentration (5, 50 and 100) exposure did not cause any effect in genes expression. Thus, we conclude that high Cd concentrations could affect the normal processes of lipid metabolism though increasing lipid synthesis and reducing the ability of fatty acid β-oxidation, and disturb thyroid hormone pathways in liver, and induced oxidative stress. In addition, lipid metabolism might be regulated by THs. To our knowledge, the present study is the first to report the influence of cadmium on hepatic lipid metabolism in B. gargarizans and will greatly provide new insights into cadmium hepatotoxicity in amphibian.

Immune responsiveness to phytohemagglutinin displays species but not sex differences in three anuran species

Phytohemagglutinin (PHA)-induced skin swelling response is widely used as a rough surrogate of integrative cell-mediated and innate immunity across multiple vertebrate taxa due to its simplification and feasibility. However, little is known whether there are sex and interspecific differences of immune responsiveness to PHA in ectotherms, especially for anurans. Therefore, we studied sex and species differences of PHA response in three anurans, Asiatic toads (Bufo gargarizans), Dark-spotted frogs (Pelophylax nigromaculatus) and Mongolian toads (Pseudepidalea raddei), captured in northern regions of Anhui Province (China). Footpad thickness was measured prior to (0 h) and after (6, 12, 24, 48 and 72 h) a PHA injection and normalized against saline injection in the opposite footpad. Body mass was recorded at the beginning (0 h) and end of each assay (72 h). Results showed effects of PHA assay, sex and taxa on body mass. Relative maximum swelling response (PHA_max) ranged from 18.58–29.75%, 9.77 to 20.56% and 21.97 to 31.78% and its occurrence over time was apparent 10.6–19.72 h , 7.74–14.01 h and 17.39–23.94 h postinjection for Asiatic toads, Dark-spotted frogs and Mongolian toads, respectively. Finally, the magnitude or timing of PHA_max in Dark-spotted frogs was significantly thinner and faster than in Mongolian toads, and Asiatic toads had an in-between value, not different from the other two species. The magnitude of PHA_max was significantly positively correlated with the timing of PHA_max considering individuals altogether, but not when analyzed within species. Our results indicate that male and female anuran species respond similarly to PHA antigen stimulation, but the magnitude and timing of PHA_max is species-specific. Briefly, we provide new evidence for the suitability of PHA assay in non-model anuran species with different body sizes, and exhort the need to further investigate the nature of PHA assay at the hematological and histological levels in order to extend its application in ecoimmunological studies of amphibians.

[ffects of Caulis Sinomenii decoction on electrophysiological character of sciatic nerve trunk of toad]

OBJECTIVE

To investigate the effects of Caulis Sinomenii (traditional Chinese medicine) on electrophysiological character of sciatic nerve trunk of toad.

METHODS

Forty toads were randomly divided into four groups:Ringer's solution (as a control group) and three experimental groups were treated with 0.20, 0.10, 0.05 g/ml Caulis Sinomenii decoction, each group was soaked for 15 and 30 min respectively. The conduction velocity, amplitude and threshold intensity of the sciatic nerves from each group were measured by RM6240C Physiological signals recording system.

RESULTS

Compared with the control group, the conduction velocity in high dose group was significantly decreased(P<0.01), the amplitude of the action potential in medium and high dose groups were lower(P<0.01); but the threshold intensity of the nerve trunk action potential in high dose group was increased significantly compared with that in control group (P<0.01).

CONCLUSIONS

The Caulis Sinomenii decoction can reduce the excitability sciatic nerve stem of toad and block the conduction of the action potential, which may play important roles in inhibiting of sciatica.

Chronic Effects of Fluoride Exposure on Growth, Metamorphosis, and Skeleton Development in Bufo gargarizans Larvae

Bufo gargarizans tadpoles were chronically exposed to waterborne fluoride at measured concentrations ranging from 0.4 to 61.2 mg F^-/L for 70 days from Gosner stage 26 to completion of metamorphosis. The chronic exposure caused a concentration-dependent mortality in all tested fluoride concentrations. Total length, snout-to-vent length (SVL), body mass, and developmental stage of tadpoles were significantly inhibited at 42.6 mg F^-/L. In addition, significant metamorphic delay and increase in size at completion of metamorphosis occurred after exposure to 19.8 mg F^-/L. Moreover, 19.8 mg F^-/L suppressed the bone mineralization of larvae at completion of metamorphosis. However, the bone mineralization could be enhanced by 4.1 mg F^-/L. In conclusion, our results suggested that the presence of high concentrations of fluoride could increase mortality risk, delay metamorphosis, and suppress skeletal ossification in B. gargarizans larvae.

Oxidative stress, endocrine disruption, and malformation of Bufo gargarizans embryo exposed to sub-lethal cadmium concentrations

Thyroid hormone (TH) is critical for vertebrate postembryonic development as well as embryonic development. Chinese toad (Bufo gargarizans) embryos were exposed to different concentrations of cadmium (5, 50, 100, 200 and 500μg Cd L-1) for 7days. Malformations were monitored daily, and growth and development of embryos were measured at day 4 and 7, and type 2 and 3 iodothyronine deiodinase (Dio2 and Dio3), thyroid hormone receptors (TRα and TRβ) mRNA levels were also measured to assess disruption of TH synthesis. In addition, superoxide dismutase (SOD), glutathione peroxidase (GPx) and heat shock proteins (HSPs) mRNA expression were examined to evaluate the ability of scavenging ROS. Our results demonstrated a bimodal inhibitory effect of Cd on the embryo growth and development of Bufo gargarizans. Reduced mean stage, total length and weight were observed at 5, 50, 200 and 500, but not at 100μg Cd L-1. Embryos malformation occurred in all cadmium treatments. Morphological abnormalities of embryos are characterized by axial flexures, abdominal edema, stunted growth and fin flexure. Real-time PCR results show that exposure to cadmium down-regulated TRα and Dio3 mRNA expression and up-regulated Dio2 mRNA level. SOD and GPx mRNA expression was significantly up-regulated after cadmium exposure. We concluded that cadmium could change mRNA expression of TRα, Dio2 and Dio3 leading the inhibition of growth and development of B. gargarizans embryo, which suggests that cadmium might have the endocrine-disrupting effect in embryos. Moreover, the reduced ability of scavenging ROS induced by cadmium might be responsible for the teratogenic effects of cadmium.

Chronic effects of triclosan on embryonic development of Chinese toad, Bufo gargarizans

Triclosan (TCS) is commonly used worldwide in a range of personal care and sanitizing products. The aim of this study was to evaluate potential effects of TCS exposure on embryonic development of Bufo gargarizans, an endemic frog species in China. Standard Gosner stage 3 B. gargarizans embryos were exposed to 10 ~ 150 μg/L TCS during embryogenesis. Survival, total length, weight, developmental stage, duration of different embryo stages, malformation, and type II and III deiodinase (D2 and D3) expression were measured. Inhibitory effects on embryo developmental stage, total length and weight were found at 30 ~ 150 μg/L TCS. Moreover, the duration of embryonic development was increased at gastrula, neural, circulation, and operculum development stage in TCS-treated groups. In addition, TCS exposure induced morphological malformations in B. gargarizans embryos, which are characterized by hyperplasia, abdominal edema, and axial flexures. Furthermore, our results showed that the expression of D2 in embryos was probably down-regulated at 60 and 150 μg/L TCS, but its spatial expression patterns was not affected by TCS. In summary, our study suggested that TCS exposure not only resulted in delayed growth and development but also caused teratogenic effects in B. gargarizans embryos, and the developmental effects of TCS at high concentrations may be associated with disruption of THs homeostasis. Although further studies are necessary, the present findings could provide a basis for understanding on harmful effects and the potential mechanisms of TCS in amphibian embryos.

[Effects of different stocking density and diet on development of Bufo gargarizans tadpoles]

The effects of four different stocking densities and five different diets on the growth of Bufo gargarizans tadpoles were studied to determine the optimum stocking density and diet. For stocking density experiment, the tadpoles were fed respectively at different density of 200, 500, 1 000 and 2 000 tadpoles per square meter. For diet experiment, the tadpoles were divided into five groups fed respectively with five different diets. The body weight, snout-vent length and tail length were measured every seven days, and mortality was recorded. The results showed that: the survival rates of tadpole before metamorphosis and after metamorphosis were from 68.7% to 96.3% and from 5.7% to 36.0%, respectively; the optimum stocking density is 1 000 tadpoles per square meter for the stocking density had no effect on the survival rate of tadpole before metamorphosis, and the tadpoles had the relative large body weight and survival numbers in 1000 tadpoles per square meter; the diet Ⅱ(37.9% crude protein and 5.7% crude fat), Ⅳ (25.1% crude protein and 4.0% crude fat), and Ⅴ (egg yolk) were all the optimum diets for the diet had no effect on the survival rate of tadpole before metamorphosis and the tadpoles fed with three kinds of diet above had relatively large body weight, and one of these three diets based on their availability and cost should be adopted during breeding period.

Different morphologic formation patterns of dark patches in the black-spotted frog (Pelophylax nigromaculata) and the Asiatic toad (Bufo gargarizans)

The black-spotted frog (Pelophylax nigromaculata) and Asiatic toad (Bufo gargarizans), two relatively distantly related species, live in different habitats with different adaptive dark patches. To explain the formation of dark patches, the distribution patterns of melanin granules were examined with light microscopy and transmission electron microscopy. Melanin granules were produced and gathered into the "cap" structures on top of the nuclei in most epidermal cells. The "cap" structures may play a role in forming the dorsal dark patches coupled with three-layer melanophores, which can give rise to three layers of interconnected melanin networks in the dorsal dermis in P. nigromaculata. Epidermal melanocytes are rare and do not have a definitive role in forming dorsal dark patches in either P. nigromaculata or B. gargarizans. In B. gargarizans, the dermal melanophores only give rise to a single-layered melanin network, which hardly results in dark patches in the dorsal skin. However, the dermal melanophores migrate twice and form into pseudostratified networks, leading to dark patch formation in the ventral skin in B. gargarizans. The melanin granules precisely coregulate dark patches in the dermis and/or epidermis in P. nigromaculata and B. gargarizans. The dark patch formation depends on melanin granules in the epidermis or/and dermis in P. nigromaculata and B. gargarizans.

Effects of copper on growth, metamorphosis and endocrine disruption of Bufo gargarizans larvae

Chinese toad (Bufo gargarizans) tadpoles were exposed to copper (1, 6.4, 32 and 64μgL(-1) copper) from the beginning of larval period through completion of metamorphosis. We examined the effects of chronic copper exposure on mortality, growth, time to metamorphosis, tail resorption time, body size at the metamorphic climax (Gs 42) and completion of metamorphosis (Gs 46) and thyroid gland histology. In addition, type 2 and 3 iodothyronine deiodinase (Dio2 and Dio3), thyroid hormone receptors (TRα and TRβ) mRNA levels were also measured to assess disruption of TH synthesis. Our result showed that 6.4-64μgL(-1) copper concentration increased the mortality and inhibited the growth of B. gargarizans tadpoles. In addition, significant reduction in size at Gs 42 and a time delay to Gs 42 were observed at 6.4-64μgL(-1) copper treatments. Moreover, histological examinations have clearly revealed that 64μgL(-1) copper caused follicular cell hyperplasia in thyroid gland. According to real-time PCR results, exposure to 32 and 64μgL(-1) copper significantly up-regulated mRNA expression of Dio3, but down-regulated mRNA expression of TRα and TRβ mRNA level. We concluded that copper delayed amphibian metamorphosis through changing mRNA expression of Dio3, TRα and TRβ, which suggests that copper might have the endocrine-disrupting effect.

Effects of nitrate on metamorphosis, thyroid and iodothyronine deiodinases expression in Bufo gargarizans larvae

Chinese toad (Bufo gargarizans) tadpoles were exposed to nitrate (10, 50 and 100mg/L NO3-N) from the beginning of the larval period through metamorphic climax. We examined the effects of chronic nitrate exposure on metamorphosis, mortality, body size and thyroid gland. In addition, thyroid hormone (TH) levels, type II iodothyronine deiodinase (Dio2) and type III iodothyronine deiodinase (Dio3) mRNA levels were also measured to assess disruption of TH synthesis. Results showed that significant metamorphic delay and mortality increased were caused in larvae exposed to 100mg/L NO3-N. The larvae exposed to 100mg/L NO3-N clearly exhibited a greater reduction in thyroxine (T4) and 3,5,3'-triiodothyronine (T3) levels. Moreover, treatment with NO3-N induced down-regulation of Dio2 mRNA levels and up-regulation of Dio3 mRNA levels, reflecting the disruption of thyroid endocrine. It seems that increased mass and body size may be correlated with prolonged metamorphosis. Interestingly, we observed an exception that exposure to 100mg/L NO3-N did not exhibit remarkable alterations of thyroid gland size. Compared with control groups, 100mg/L NO3-N caused partial colloid depletion in the thyroid gland follicles. These results suggest that nitrate can act as a chemical stressor inducing retardation in development and metamorphosis. Therefore, we concluded that the presence of high concentrations nitrate can influence the growth, decline the survival, impair TH synthesis and induce metamorphosis retardation of B. gargarizans larvae.