Analysis of evolutionary and functional features of the bullfrog SULT1 family.

We identified the bullfrog Rana catesbeiana sulfotransferase 1 (SULT1) family from the BLAST search tool of the public databases based on the SULT1 families of Nanorana parkeri, Xenopus laevis, and Xenopus tropicalis as queries, revealing the characteristics of the anuran SULT1 family. The results showed that the anuran SULT1 family comprises six subfamilies, four of which were related to the mammalian SULT1 subfamily. Additionally, the bullfrog has two SULT1Cc subfamily members that are consistent with the characteristics of the expanded Xenopus SULT1C subfamily. Several members of the bullfrog SULT1 family were suggested to play important roles in sulfation during metamorphosis. Among these, cDNAs encoding SULT1Cc1 and SULT1Y1 were cloned, and the sulfation activity was analyzed using recombinant proteins. The affinity for 2-naphthol and 3'-phosphoadenosine 5'-phosphosulfate (PAPS) and the enzymatic reaction rate were higher in SULT1Cc1 than in SULT1Y1. Both the enzymes showed inhibitory effect of many thyroid hormones (THs) analogs on the sulfation of 2-naphthol. The potency of sulfation activities of SULT1Cc1 and SULT1Y1 against T4 indicated their possible role in the intracellular T4 clearance during metamorphosis.

Seasonal acclimatization and thermal acclimation induce global histone epigenetic changes in liver of bullfrog (Lithobates catesbeianus) tadpole.

The American bullfrog (Lithobates catesbeianus) is a eurythermal amphibian that is naturally distributed from subarctic to subtropical areas. The tadpoles of this species overwinter, in water, in cold environments. Therefore, they may have adapted to a wide range of temperatures in an active state. To understand the adaptation mechanisms to cope with low or high temperatures, we investigated global epigenetic modifications, histone variants, transcript levels of related genes, and the cellular acetyl coenzyme A (acetyl-CoA) and free CoA (CoA-SH) levels, in the livers of tadpoles collected in summer and winter and of those acclimated to 4 °C and 21 °C. Among epigenetic marks tested, the levels of acetylated histones and the histone variant H2A.Z were influenced by different temperature conditions. Histone acetylation levels were higher in summer than in winter and increased within 3 days of warm acclimation, whereas histone H2A.Z levels were higher in winter than in summer and decreased within 2 weeks of warm acclimation. Transcript analysis revealed that decreased expression of histone H2A.Z in warm acclimation was regulated at the transcriptional level. Acetyl-CoA levels were not correlated with those of the acetylated histones, indicating that cellular acetyl-CoA levels may not directly influence the state of histone acetylation in the tadpole liver. Such epigenetic and metabolic changes in the tadpole liver may contribute to the maintenance of energy balance during seasonal acclimatization and thermal acclimation.

Analysis of global and gene-specific acetylation of histones in the liver of American bullfrog (Rana catesbeiana) tadpoles acclimated to low temperature.

Severe environmental stressors such as low temperatures can affect gene expression by changing epigenetic states. American bullfrog (Rana catesbeiana) can overwinter as tadpoles, which can be active even in winter. However, the molecular mechanisms of epigenetic controls by which the tadpoles acclimate to low temperature are still unclear. In this study, we aimed to clarify the molecular mechanisms of global and gene-specific epigenetic regulations of low-temperature acclimation. We found that the global acetylation was decreased in the liver of bullfrog tadpoles acclimated to low temperature. The amounts of transcripts for two histone acetyltransferases were higher in the liver of tadpoles acclimated to low temperature than in those acclimated to warm temperature, while we observed no significant differences in the amounts of transcripts for histone deacetylases. We also found that the amounts of transcripts and acetylated histones on the specific temperature-responsive genes scd and cyp7a1 whose transcripts were increased and decreased, respectively, in response to low temperature were positively correlated. Cellular acetyl-CoA levels were higher in the liver of tadpoles acclimated to low temperature than in those acclimated to warm temperature. These results contradicted the states of histone acetylation, suggesting that bullfrog tadpoles have different epigenetic mechanisms to modify the histones when compared with those of other organisms such as reptiles and mammals, even though the relationship between the transcript amount and the states of histone acetylation on temperature-responsive genes was similar to that of mammals.

Behavioral and molecular analyses of olfaction-mediated avoidance responses of Rana (Lithobates) catesbeiana tadpoles: Sensitivity to thyroid hormones, estrogen, and treated municipal wastewater effluent.

Olfaction is critical for survival, facilitating predator avoidance and food location. The nature of the olfactory system changes during amphibian metamorphosis as the aquatic herbivorous tadpole transitions to a terrestrial, carnivorous frog. Metamorphosis is principally dependent on the action of thyroid hormones (THs), l-thyroxine (T4) and 3,5,3'-triiodothyronine (T3), yet little is known about their influence on olfaction during this phase of postembryonic development. We exposed Taylor Kollros stage I-XIII Rana (Lithobates) catesbeiana tadpoles to physiological concentrations of T4, T3, or 17-beta-estradiol (E2) for 48h and evaluated a predator cue avoidance response. The avoidance response in T3-exposed tadpoles was abolished while T4- or E2-exposed tadpoles were unaffected compared to control tadpoles. qPCR analyses on classic TH-response gene transcripts (thra, thrb, and thibz) in the olfactory epithelium demonstrated that, while both THs produced molecular responses, T3 elicited greater responses than T4. Municipal wastewater feed stock was spiked with a defined pharmaceutical and personal care product (PPCP) cocktail and treated with an anaerobic membrane bioreactor (AnMBR). Despite substantially reduced PPCP levels, exposure to this effluent abolished avoidance behavior relative to AnMBR effluent whose feed stock was spiked with vehicle. Thibz transcript levels increased upon exposure to either effluent indicating TH mimic activity. The present work is the first to demonstrate differential TH responsiveness of the frog tadpole olfactory system with both behavioral and molecular alterations. A systems-based analysis is warranted to further elucidate the mechanism of action on the olfactory epithelium and identify further molecular bioindicators linked to behavioral response disruption.

Influence of Temperature on the Thyroidogenic Effects of Diuron and Its Metabolite 3,4-DCA in Tadpoles of the American Bullfrog (Lithobates catesbeianus).

Temperature is a key variable affecting the timing of amphibian metamorphosis from tadpoles to tetrapods, through the production and subsequent function of thyroid hormones (TH). Thyroid function can be impaired by environmental contaminants as well as temperature. Tadpoles can experience large temperature fluctuations in their habitats and many species are distributed in areas that may be impacted by agriculture. Diuron is a widely used herbicide detected in freshwater ecosystems and may impact endocrine function in aquatic organisms. We evaluated the influence of temperature (28 and 34 °C) on the action of diuron and its metabolite 3,4-dichloroaniline (3,4-DCA) on thyroid function and metamorphosis in tadpoles of Lithobates catesbeianus. Exposure to both compounds induced more pronounced changes in gene expression and plasma 3,3',5-triiodothyronine (T3) concentrations in tadpoles treated at higher temperature. T3 concentrations were increased in tadpoles exposed to 200 ng/L of diuron at 34 °C and an acceleration of metamorphosis was observed for the same group. Transcriptomic responses included alteration of thyroid hormone induced bZip protein (thibz), deiodinases (dio2, dio3), thyroid receptors (trα, trß) and Krüppel-like factor 9 (klf9), suggesting regulation by temperature on TH-gene expression. These results suggest that environmental temperature should be considered in risk assessments of environmental contaminants for amphibian species.

Effects of acute exposure to the non-steroidal anti-inflammatory drug ibuprofen on the developing North American Bullfrog (Rana catesbeiana) tadpole.

A variety of pharmaceutical chemicals can represent constituents of municipal effluent outflows that are dispersed into aquatic receiving environments worldwide. Increasingly, there is concern as to the potential of such bioactive substances to interact with wildlife species at sensitive life stages and affect their biology. Using a combination of DNA microarray, quantitative real-time polymerase chain reaction, and quantitative nuclease protection assays, we assessed the ability of sub-lethal and environmentally relevant concentrations of ibuprofen (IBF), a non-steroidal anti-inflammatory agent and prevalent environmental contaminant, to function as a disruptor of endocrine-mediated post-embryonic development of the frog. While the LC50 of IBF for pre-metamorphic Rana catesbeiana tadpoles is 41.5 mg/L (95% confidence interval: 32.3-53.5 mg/L), exposure to concentrations in the ppb range elicited molecular responses both in vivo and in organ culture. A nominal concentration of 15 µg/L IBF (actual = 13.7 µg/L) altered the abundance of 26 mRNA transcripts within the liver of exposed pre-metamorphic R. catesbeiana tadpoles within 6 d. IBF-treated animals demonstrated subsequent disruption of thyroid hormone-mediated reprogramming in the liver transcriptome affecting constituents of several metabolic, developmental, and signaling pathways. Cultured tadpole tail fin treated with IBF for 48 h also demonstrated altered mRNA levels at drug concentrations as low as 1.5 µg/L. These observations raise the possibility that IBF may alter the post-embryonic development of anuran species in freshwater environs, where IBF is a persistent or seasonal pollutant.

Utilization of temperature-mediated activation of thyroid hormone-induced molecular memory to evaluate early signaling events in the olfactory epithelium of Rana [Lithobates] catesbeiana tadpoles.

The amphibian olfactory system is highly distinct between aquatic tadpole and terrestrial frog life stages and therefore must remodel extensively during thyroid hormone (TH)-dependent metamorphosis. Developmentally appropriate functioning of the olfactory epithelium is critical for survival. Previous studies in other Rana [Lithobates] catesbeiana premetamorphic tadpole tissues showed that initiation of TH-induced metamorphosis can be uncoupled from execution of TH-dependent programs by holding tadpoles in the cold rather than at warmer permissive temperatures. TH-exposed tadpoles at the nonpermissive (5 °C) temperature do not undergo metamorphosis but retain a "molecular memory" of TH exposure that is activated upon shift to a permissive warm temperature. Herein, premetamorphic tadpoles were held at permissive (24 °C) or nonpermissive (5 °C) temperatures and injected with 10 pmoles/g body weight 3,5,3'-triiodothyronine (T3) or solvent control. Olfactory epithelium was collected at 48 h post-injection. RNA-sequencing (RNA-Seq) and reverse transcriptase quantitative real-time polymerase chain reaction (RT-qPCR) analyses generated differentially expressed transcript profiles of 4328 and 54 contigs for permissive and nonpermissive temperatures, respectively. Translation, rRNA, spliceosome, and proteolytic processes gene ontologies were enriched by T3 treatment at 24 °C while negative regulation of cell proliferation was enriched by T3 at 5 °C. Of note, as found in other tissues, TH-induced basic leucine zipper-containing protein-encoding transcript, thibz, was significantly induced by T3 at both temperatures, suggesting a role in the establishment of molecular memory in the olfactory epithelium. The current study provides critical insights by deconstructing early TH-induced induction of postembryonic processes that may be targets for disruption by environmental contaminants.

Horizontal plasmid transfer promotes antibiotic resistance in selected bacteria in Chinese frog farms.

The emergence and dissemination of antibiotic resistance genes (ARGs) in the ecosystem are global public health concerns. One Health emphasizes the interconnectivity between different habitats and seeks to optimize animal, human, and environmental health. However, information on the dissemination of antibiotic resistance genes (ARGs) within complex microbiomes in natural habitats is scarce. We investigated the prevalence of antibiotic resistant bacteria (ARB) and the spread of ARGs in intensive bullfrog (Rana catesbeiana) farms in the Shantou area of China. Antibiotic susceptibilities of 361 strains, combined with microbiome analyses, revealed Escherichia coli, Edwardsiella tarda, Citrobacter and Klebsiella sp. as prevalent multidrug resistant bacteria on these farms. Whole genome sequencing of 95 ARB identified 250 large plasmids that harbored a wide range of ARGs. Plasmid sequences and sediment metagenomes revealed an abundance of tetA, sul1, and aph(3″)-Ib ARGs. Notably, antibiotic resistance (against 15 antibiotics) highly correlated with plasmid-borne rather than chromosome-borne ARGs. Based on sequence similarities, most plasmids (62%) fell into 32 distinct groups, indicating a potential for horizontal plasmid transfer (HPT) within the frog farm microbiome. HPT was confirmed in inter- and intra-species conjugation experiments. Furthermore, identical mobile ARGs, flanked by mobile genetic elements (MGEs), were found in different locations on the same plasmid, or on different plasmids residing in the same or different hosts. Our results suggest a synergy between MGEs and HPT to facilitate ARGs dissemination in frog farms. Mining public databases retrieved similar plasmids from different bacterial species found in other environmental niches globally. Our findings underscore the importance of HPT in mediating the spread of ARGs in frog farms and other microbiomes of the ecosystem.

Transcriptomic profiling of Rana [Lithobates] catesbeiana back skin during natural and thyroid hormone-induced metamorphosis under different temperature regimes with particular emphasis on innate immune system components.

As amphibians undergo thyroid hormone (TH)-dependent metamorphosis from an aquatic tadpole to the terrestrial frog, their innate immune system must adapt to the new environment. Skin is a primary line of defense, yet this organ undergoes extensive remodelling during metamorphosis and how it responds to TH is poorly understood. Temperature modulation, which regulates metamorphic timing, is a unique way to uncover early TH-induced transcriptomic events. Metamorphosis of premetamorphic tadpoles is induced by exogenous TH administration at 24 °C but is paused at 5 °C. However, at 5 °C a "molecular memory" of TH exposure is retained that results in an accelerated metamorphosis upon shifting to 24 °C. We used RNA-sequencing to identify changes in Rana (Lithobates) catesbeiana back skin gene expression during natural and TH-induced metamorphosis. During natural metamorphosis, significant differential expression (DE) was observed in >6500 transcripts including classic TH-responsive transcripts (thrb and thibz), heat shock proteins, and innate immune system components: keratins, mucins, and antimicrobial peptides (AMPs). Premetamorphic tadpoles maintained at 5 °C showed 83 DE transcripts within 48 h after TH administration, including thibz which has previously been identified as a molecular memory component in other tissues. Over 3600 DE transcripts were detected in TH-treated tadpoles at 24 °C or when tadpoles held at 5 °C were shifted to 24 °C. Gene ontology (GO) terms related to transcription, RNA metabolic processes, and translation were enriched in both datasets and immune related GO terms were observed in the temperature-modulated experiment. Our findings have implications on survival as climate change affects amphibia worldwide.

Impact of wastewater treatment configuration and seasonal conditions on thyroid hormone disruption and stress effects in Rana catesbeiana tailfin.

Improved endocrine disrupting compound (EDC) removal is desirable in municipal wastewater treatment plants (MWWTPs) although increased removal does not always translate into reduced biological activity. Suitable methods for determining reduction in biological activity of effluents are needed. In order to determine which MWWTPs are the most effective at removing EDC activities, we operated three configurations of pilot sized biological reactors (conventional activated sludge, CAS; nitrifying activated sludge, NAS; and biological nutrient removal, BNR) receiving the same influent under simulated winter and summer conditions. As frogs are model organisms for the study of thyroid hormone (TH) action, we used the North American species Rana catesbeiana in a cultured tadpole tailfin (C-fin) assay to compare the effluents. TH-responsive (thyroid hormone receptors alpha (thra) and beta (thrb)) and stress-responsive (superoxide dismutase, catalase, and heat shock protein 30) mRNA transcript levels were examined. Effluents infrequently perturbed stress-responsive transcript abundance but thra/thrb levels were significantly altered. In winter conditions, CAS caused frequent TH perturbations while BNR caused none. In summer conditions, however, BNR caused substantial TH perturbations while CAS caused few. Our findings contrast other studies of seasonal variations of EDC removal and accentuate the importance of utilizing appropriate biological readouts for assessing EDC activities.

Molecular cloning and expression of ranalexin, a bioactive antimicrobial peptide from Rana catesbeiana in Escherichia coli and assessments of its biological activities.

The coding sequence, which corresponds to the mature antimicrobial peptide ranalexin from the frog Rana catesbeiana, was chemically synthesized with preferred codons for expression in Escherichia coli. It was cloned into the vector pET32c (+) to express a thioredoxin-ranalexin fusion protein which was produced in soluble form in E. coli BL21 (DE3) induced under optimized conditions. After two purification steps through affinity chromatography, about 1 mg of the recombinant ranalexin was obtained from 1 L of culture. Mass spectrometrical analysis of the purified recombinant ranalexin demonstrated its identity with ranalexin. The purified recombinant ranalexin is biologically active. It showed antibacterial activities similar to those of the native peptide against Staphylococcus aureus, Streptococcus pyogenes, E. coli, and multidrug-resistant strains of S. aureus with minimum inhibitory concentration values between 8 and 128 µg/ml. The recombinant ranalexin is also cytotoxic in HeLa and COS7 human cancer cells (IC50 = 13-15 µg/ml).

Use of Lithobates catesbeianus tadpoles in a multiple biomarker approach for the assessment of water quality of the Reconquista River (Argentina).

The water quality of the Reconquista River (Argentina) water was monitored between 2009 and 2010 by means of a multiparametric approach. This periurban river is affected by agricultural, urban, and industrial discharges. Water samples were collected at a dam located in the headwaters and at 18 km downstream (M). Physicochemical profile and two water-quality indices (WQIs) were determined. Laboratory bioassays were performed by exposing Lithobates catesbeianus tadpoles to environmental samples for 96 h and determining the following parameters: (1) brain: acetylcholinesterase (AChE) activity; (2) gill: catalase and glutathione-S-transferase (GST) activities and glutathione (GSH) content; (3) liver: CAT and GST activities, superoxide dismutase, lipid peroxidation, and GSH content; (4) condition factor and hepatosomatic index; and (5) genotoxicity [micronucleus (MN) test in erythrocytes]. Physicochemical profile and WQIs corresponded with extensive pollution in M. Important temporal and spatial variability in biomarkers of tadpoles exposed to samples was found. Multivariate analyses showed that AChE in brain, MN frequency, liver and gill GST activities, and GSH content were key biomarkers.

Models and data of AMPlify: a deep learning tool for antimicrobial peptide prediction.

OBJECTIVES: Antibiotic resistance is a rising global threat to human health and is prompting researchers to seek effective alternatives to conventional antibiotics, which include antimicrobial peptides (AMPs). Recently, we have reported AMPlify, an attentive deep learning model for predicting AMPs in databases of peptide sequences. In our tests, AMPlify outperformed the state-of-the-art. We have illustrated its use on data describing the American bullfrog (Rana [Lithobates] catesbeiana) genome. Here we present the model files and training/test data sets we used in that study. The original model (the balanced model) was trained on a balanced set of AMP and non-AMP sequences curated from public databases. In this data note, we additionally provide a model trained on an imbalanced set, in which non-AMP sequences far outnumber AMP sequences. We note that the balanced and imbalanced models would serve different use cases, and both would serve the research community, facilitating the discovery and development of novel AMPs. DATA DESCRIPTION: This data note provides two sets of models, as well as two AMP and four non-AMP sequence sets for training and testing the balanced and imbalanced models. Each model set includes five single sub-models that form an ensemble model. The first model set corresponds to the original model trained on a balanced training set that has been described in the original AMPlify manuscript, while the second model set was trained on an imbalanced training set.

Natural occurrences and characterization of Elizabethkingia miricola infection in cultured bullfrogs (Rana catesbeiana).

Introduction: The bacterium Elizabethkingia miricola is a multispecies pathogen associated with meningitis-like disease that has been isolated from several amphibian species, including the bullfrog, but this is the first isolation in Guangxi. In the present study, the dominant bacteria were isolated from the brains of five bullfrogs with meningitis-like disease on a South China farm in Guangxi. Methods: The NFEM01 isolate was identified by Gram staining; morphological observations; 16S rRNA, rpoB, and mutT-based phylogenetic tree analysis; and physiochemical characterization and was subjected to drug sensitivity and artificial infection testing. Results and discussion: As a result of identification, the NFEM01 strain was found to be E. miricola. An artificial infection experiment revealed that NFEM01 infected bullfrogs and could cause symptoms of typical meningitis-like disease. As a result of the bacterial drug sensitivity test, NFEM01 is highly sensitive to mequindox, rifampicin, enrofloxacin, nitrofural, and oxytetracycline and there was strong resistance to gentamicin, florfenicol, neomycin, penicillin, amoxicillin, doxycycline, and sulfamonomethoxine. This study provides a reference to further study the pathogenesis mechanism of E. miricola-induced bullfrog meningitislike disease and its prevention and treatment.

RNA polymerase II phosphorylation at serine 2 and histone H3 tri-methylation at lysine 36 are key steps for thyroid hormone receptor ß gene activation by thyroid hormone in Rana catesbeiana tadpole liver.

Thyroid hormone (TH) is essential for amphibian metamorphosis, during which the expression of many genes is controlled directly or indirectly through TH receptors (TRs). Thyroid hormone binding to TRs induces coregulator switching on regulatory regions of TH-inducible genes: corepressors complexed with unliganded TRs are replaced by coactivators complexed with liganded TR resulting in transcriptionally active states. The coregulator switching is linked to histone acetylation. In our study, we have investigated the acetylation and methylation states of histones H3 and H4 using chromatin immunoprecipitation (ChIP) assays on the 5' coding region of the TRß gene, a primary TH-response gene, in the liver from Rana catesbeiana tadpoles either treated with or not treated with 3,3',5-triiodothyronine (T3). 3,3',5-Triiodothyronine treatment for 3 days increased the amount of TRß transcript by 19-fold. This increase was associated with increases in the acetylation of histone H4 and lysine 9 in histone H3 (H3-K9), and tri-methylation of lysine 36 in histone H3 (H3-K36). In addition, the amounts of RNA polymerase II (PolII) and serine 2 phosphorylation in PolII (PolII-S2) increased. These results suggest that T3 treatment enhances the elongation activity of PolII on the TRß gene in the liver by increasing histone H3-K36 tri-methylation through PolII-S2 phosphorylation.

Gene expression profile in the liver of Rana catesbeiana tadpoles exposed to low temperature in the presence of thyroid hormone.

Amphibian metamorphosis, which is controlled by thyroid hormone (TH), is highly temperature-sensitive. Using real-time PCR, we investigated the gene expression profile in the liver of Rana catesbeiana tadpoles kept at 28 and 4 °C and treated with 5 nM 3,3',5-triiodothyronine (T3). Out of the 48 genes tested, 12 were up-regulated at 4 °C in T3-treated or untreated tadpoles. These included genes involved in energy metabolism, transcription, and translation. Four TH-response genes, including TH receptor ß (TRß) gene, showed no response to T3 at 4 °C. Deiodinase III was the only gene down-regulated at 4 °C. Chromatin immunoprecipitation assays showed that CCAAT/enhancer-binding protein 2 gene activation by cold exposure was associated with an increase in the acetylation of histone H3 at lysine 9, whereas TRß gene activation by T3 at 28 °C was associated with an increase in the methylation of histone H3 at lysine 36 with no T3-dependent changes in methylation states on cold exposure. Our results suggest that the transfer of TH signal to chromatin modifications on a primary early TH response gene was specifically blocked by exposure to cold.

Genetic structure of American bullfrog populations in Brazil.

Non-native species are a major problem affecting numerous biomes around the globe. Information on their population genetics is crucial for understanding their invasion history and dynamics. We evaluated the population structure of the non-native American bullfrog, Aquarana catesbeiana, in Brazil on the basis of 324 samples collected from feral and captive groups at 38 sites in seven of the nine states where feral populations occur. We genotyped all samples using previously developed, highly polymorphic microsatellite loci and performed a discriminant analysis of principal components together with Jost's D index to quantify pairwise differentiation between populations. We then amplified 1,047 base pairs of the mitochondrial cytochrome b (cytb) gene from the most divergent samples from each genetic population and calculated their pairwise differences. Both the microsatellite and cytb data indicated that bullfrogs comprise two populations. Population grouping 1 is widespread and possesses two cytb haplotypes. Population grouping 2 is restricted to only one state and possesses only one of the haplotypes from Population grouping 1. We show that there were two imports of bullfrogs to Brazil and that there is low genetic exchange between population groupings. Also, we find that there is no genetic divergence among feral and captive populations suggesting continuous releases. The limited genetic variability present in the country is associated to the small number of introductions and founders. Feral bullfrogs are highly associated to leaks from farms, and control measures should focus on preventing escapes using other resources than genetics, as feral and captive populations do not differ.

Invasive Bullfrogs Maintain MHC Polymorphism Including Alleles Associated with Chytrid Fungal Infection.

Maintenance of genetic diversity at adaptive loci may facilitate invasions by non-native species by allowing populations to adapt to novel environments, despite the loss of diversity at neutral loci that typically occurs during founder events. To evaluate this prediction, we compared genetic diversity at major histocompatibility complex (MHC) and cytochrome b (cytb) loci from 20 populations of the American bullfrog (Rana catesbeiana) across theinvasive and native ranges in North America and quantified the presence of the pathogen Batrachochytrium dendrobatidis (Bd). Compared to native populations, invasive populations had significantly higher Bd prevalence and intensity, significantly higher pairwise MHC and cytb FST, and significantly lower cytb diversity, but maintained similar levels of MHC diversity. The two most common MHC alleles (LiCA_B and Rapi_33) were associated with a significant decreased risk of Bd infection, and we detected positive selection acting on four peptide binding residues. Phylogenetic analysis suggested invasive populations likely arose from a single founding population in the American Midwest with a possible subsequent invasion in the northwest. Overall, our study suggests that the maintenance of diversity at adaptive loci may contribute to invasion success and highlights the importance of quantifying diversity at functional loci to assess the evolutionary potential of invasive populations.

Relationship between serum thyroid hormones and their associated metabolites, and gene expression bioindicators in the back skin of Rana [Lithobates] catesbeiana tadpoles and frogs during metamorphosis.

Anuran metamorphosis is characterized by profound morphological changes including remodeling of tissues and organs. This transition is initiated by thyroid hormones (THs). However, the current knowledge of changing levels of THs during metamorphosis relies on pooled samples using methods known for high variability with sparse reporting of measured variation. Moreover, establishing a clear linkage between key gene expression bioindicators and TH levels throughout the metamorphic process is needed. Using state-of-the-art ultra-high performance liquid chromatography isotope-dilution tandem mass spectrometry, we targeted 12 THs and metabolites in the serum of Rana [Lithobates] catesbeiana (n=5-10) across seven distinct postembryonic stages beginning with premetamorphic tadpoles (Gosner stage 31-33) and continuing through metamorphosis to a juvenile frog (Gosner stage 46). TH levels were related to TH-relevant gene transcripts (thra, thrb, and thibz) in back skin of the same individual animals. Significant increases from basal levels were observed for thyroxine (T4) and 3,3',5-triiodothyronine (T3) at Gosner stage 41, reaching maximal levels at Gosner stage 44 (28 ± 10 and 2.3 ± 0.5 ng/mL, respectively), and decreasing to basal levels in juvenile frogs. In contrast, 3,5-diiodothyronine (T2) increased significantly at Gosner stage 40 and was maintained elevated until stage 44. While thra transcript levels remained constant and then decreased at the end of metamorphic climax, thrb and thibz were induced to maximal levels at Gosner stage 41, followed by a decrease to basal levels in the froglet. This exemplifies the exquisite timing of events during metamorphosis as classic early response genes are transcribed in anticipation of peak TH concentrations. The distinct T2 concentration profile suggests a biological role of this biomolecule in anuran postembryonic development and an additional aspect that may be a target of anthropogenic chemicals that can disrupt anuran metamorphosis and TH signalling. Hence, as a second aim of the study, we set out to find additional bioindicators of metamorphosis, which can aid future investigations of developmental disruption. Using a sensitive nanoLC-Orbitrap system an untargeted analysis workflow was applied. Among 6,062 endogenous metabolites, 421 showed metamorphosis-dependent concentration dynamics. These potential bioindicators included several carnitines, prostaglandins and some steroid hormones.

Microcystin-leucine arginine exposure induced intestinal lipid accumulation and MC-LR efflux disorder in Lithobates catesbeianus tadpoles.

Few studies exist on the toxic effects of chronic exposure to microcystins (MCs) on amphibian intestines, and the toxicity mechanisms are unclear. Here, we evaluated the impact of subchronic exposure (30 days) to environmentally realistic microcystin-leucine arginine (MC-LR) concentrations (0 µg/L, 0.5 µg/L and 2 µg/L) on tadpole (Lithobates catesbeianus) intestines by analyzing the histopathological and subcellular microstructural damage, the antioxidative and oxidative enzyme activities, and the transcriptome levels. Histopathological results showed severe damage accompanied by inflammation to the intestinal tissues as the MC-LR exposure concentration increased from 0.5 µg/L to 2 µg/L. RNA-sequencing analysis identified 634 and 1,147 differentially expressed genes (DEGs) after exposure to 0.5 µg/L and 2 µg/L MC-LR, respectively, compared with those of the control group (0 µg/L). Biosynthesis of unsaturated fatty acids and the peroxisome proliferator-activated receptor (PPAR) signaling pathway were upregulated in the intestinal tissues of the exposed groups, with many lipid droplets being observed on transmission electron microscopy, implying that MC-LR may induce lipid accumulation in frog intestines. Moreover, 2 µg/L of MC-LR exposure inhibited the xenobiotic and toxicant biodegradation related to detoxification, implying that the tadpoles' intestinal detoxification ability was weakened after exposure to 2 µg/L MC-LR, which may aggravate intestinal toxicity. Lipid accumulation and toxin efflux disorder may be caused by MC-LR-induced endoplasmic reticular stress. This study presents new evidence that MC-LR harms amphibians by impairing intestinal lipid metabolism and toxin efflux, providing a theoretical basis for evaluating the health risks of MC-LR to amphibians.