Search

The leptin receptor overlapping transcript (LepROT) has been suggested to play several roles in immunomodulatory mechanisms; however, the understanding of its role in Rana amurensis immunity is still very limited. Here, we performed hematoxylin-eosin staining, quantitative reverse-transcription polymerase chain reaction (qRT-PCR), immunofluorescence and western blotting to investigate the roles of LepROT in the immunomodulatory mechanism and the influence of its expression on the nuclear factor-kappa B (NF-κB) signaling pathway, such as the activation of IκB kinase and NF-кB, in amphibian resistance to infection with Aeromonas hydrophila ( Ah ). After Ah infection, the liver, lung, kidney, skin, muscle, and stomach of R. amurensis showed cell structure disturbance, bleeding, and texture abnormalities. In addition, the relative expression levels of LepROT, NF-кB, IKKα, and IKKβ were all upregulated after Ah infection; however, they showed time-dependent differential expression. The NF-кB signaling pathway exhibited robust expression levels, which might be explained by the positive feedback regulation function of LepROT. Overall, this study provides a basis for further assessment of the biological functions of LepROT and highlights its role in the regulation of immune mechanisms.

Green algae usually assigned to the genus Oophila are known to colonize egg capsules of amphibian egg masses across the Nearctic and Palearctic regions. We study the phylogenetic relationships of these algae using a phylotranscriptomic data set of 76 protein-coding single-copy nuclear genes. Our data set includes novel RNAseq data for six amphibian-associated and five free-living green algae, and draft genomes of two of the latter. Within the Oophila clade (nested within Moewusinia), we find samples from two European frogs (Rana dalmatina and R. temporaria) closely related to those of the North American frog R. aurora (Oophila subclade III). An isolate from the North American R. sylvatica (subclade IV) appears to be sister to the Japanese isolate from the salamander Hynobius nigrescens (subclade J1), and subclade I algae from Ambystoma maculatum are sister to all other lineages in the Oophila clade. Two free-living algae (Chlamydomonas nasuta and Cd. pseudogloeogama) are nested within the Oophila clade, and a strain of the type species of Chlorococcum (Cc. infusionum) is related to this assemblage. Our phylotranscriptomic tree suggests that recognition of different species within the Oophila clade ("clade B" of earlier studies) is warranted, and calls for a comprehensive taxonomic revision of Moewusinia., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Lead (Pb) is a ubiquitously detected heavy metal pollutant in aquatic ecosystems. Previous studies focused mainly on the response of gut microbiota to Pb stress, with less emphasis on gene expression in intestine, thereby limiting the information about impacts of Pb on intestinal homeostasis in amphibians. Here, microbial community and transcriptional response of intestines in Rana zhenhaiensis tadpoles to Pb exposure were evaluated. Our results showed that 10 μg/L Pb significantly decreased bacterial diversity compared to the controls by the Simpson index. Additionally, 1000 μg/L Pb exposure resulted in a significant reduction in the abundance of Fusobacteriota phylum and Cetobacterium genus but a significant expansion in Hafnia-Obesumbacterium genus. Moreover, transcriptome analysis revealed that about 90 % of the DEGs (8458 out of 9450 DEGs) were down-regulated in 1000 μg/L Pb group, mainly including genes annotated with biological functions in fatty acid degradation, and oxidative phosphorylation, while up-regulated DEGs involved in metabolism of xenobiotics by cytochrome P450. The expression of Gsto1, Gsta5, Gstt4, and Nadph showed strong correlation with the abundance of genera Serratia, Lactococcus, and Hafnia-Obesumbacterium. The findings of this study provide important insights into understanding the influence of Pb on intestinal homeostasis in amphibians., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Wood frogs are freeze-tolerant vertebrates that can endure weeks to months frozen during the winter without breathing and with as much as 65% of total body water frozen as extracellular ice. Underlying tolerances of anoxia and of cellular dehydration support whole body freezing. One pro-survival mechanism employed by these frogs is epigenetic modifications via DNA hypomethylation processes facilitating transcriptional repression or activation. These processes involve proteins such as DNA Methyltransferases (DNMTs), Methyl Binding Domain proteins (MBDs), Ten-Eleven Translocases (TETs), and Thymine Deglycosylase (TDG). The present study evaluates the responses of these proteins to dehydration and anoxia stresses in wood frog liver. DNMT relative protein expression was reduced in liver, but nuclear DNMT activity did not change significantly under anoxia stress. By contrast, liver DNMTs and nuclear DNMT activity were upregulated under dehydration stress. These stress-specific differences were speculated to arise from Post-Translational Modifications (PTMs). DNMT3A and DNMT3B showed increased relative protein expression during recovery from dehydration and anoxia. Further, MBD1 was elevated during both conditions suggesting transcriptional repression. TET proteins showed varying responses to anoxia likely due to the absence of oxygen, a main substrate required by TETs. Similarly, TDG, an enzyme that corrects DNA damage, was downregulated under anoxia potentially due to lower levels of reactive oxygen species that damage DNA, but levels returned to normal during reperfusion of oxygen. Our results indicate differential stress-specific responses that indicate the need for more research in the DNA hypomethylation mechanisms employed by the wood frog during stress., Competing Interests: Declaration of competing interest The authors declare no conflict of interest., (Copyright © 2023. Published by Elsevier Inc.)

The oviduct of the Chinese brown frog (Rana dybowskii) expands during pre-brumation rather than the breeding period, exhibiting a special physiological feature. Vitamin A is essential for the proper growth and development of many organisms, including the reproductive system such as ovary and oviduct. Vitamin A is metabolized into retinoic acid, which is crucial for oviduct formation. This study examined the relationship between oviducal expansion and vitamin A metabolism. We observed a significant increase in the weight and diameter of the oviduct in Rana dybowskii during pre-brumation. Vitamin A and its active metabolite, retinoic acid, notably increased during pre-brumation. The mRNA levels of retinol binding protein 4 (rbp4) and its receptor stra6 gene, involved in vitamin A transport, were elevated during pre-brumation compared to the breeding period. In the vitamin A metabolic pathway, the mRNA expression level of retinoic acid synthase aldh1a2 decreased significantly during pre-brumation, while the mRNA levels of retinoic acid α receptor (rarα) and the retinoic acid catabolic enzyme cyp26a1 increased significantly during pre-brumation, but not during the breeding period. Immunohistochemical results showed that Rbp4, Stra6, Aldh1a2, Rarα, and Cyp26a1 were expressed in ampulla region of the oviduct. Western blot results indicated that Aldh1a2 expression was lower, while Rbp4, Stra6, RARα, and Cyp26a1 were higher during pre-brumation compared to the breeding period. Transcriptome analyses further identified differential genes in the oviduct and found enrichment of differential genes in the vitamin A metabolism pathway, providing evidences for our study. These results suggest that the vitamin A metabolic pathway is more active during pre-brumation compared to the breeding period, and retinoic acid may regulate pre-brumation oviductal expansion through Rarα-mediated autocrine/paracrine modulation., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Gametogenesis produces gametes as a piece of genetic information transmitted to the offspring. While during sexual reproduction, progeny inherits a mix of genetic material from both parents, asexually reproducing organisms transfer a copy of maternal or paternal DNA to the progeny clonally. Parthenogenetic, gynogenetic and hybridogenetic animals have developed various mechanisms of gametogenesis, however, their inheritance is not fully understood. Here, we focused on the inheritance of asexual gametogenesis in hybrid Pelophylax esculentus (RL), emerging after crosses of P. lessonae (LL) and P. ridibundus (RR). To understand the mechanisms of gametogenesis in hybrids, we performed three-generation experiments of sexual P. ridibundus females and hybrids from all-male hybrid populations. Using fluorescent in situ hybridization, micronuclei analysis, flow cytometry and genotyping, we found that most adult hybrid males simultaneously produced two types of clonal sperm. Also, most male tadpole progeny in two successive backcrossed generations simultaneously eliminated L and R parental genomes, while some progeny produced only one type of sperm. We hypothesize that the reproductive variability of males producing two kinds of sperm is an adaptive mechanism to reproduce in mixed populations with P. ridibundus and may explain the extensive distribution of the all-male lineage across the European River Basin., (© 2024. The Author(s).)

Oviductus Ranae is the dried oviduct from Rana dybowskii, a forest frog species with medicinal, tonic, and cosmetic properties. Due to the high price and resource shortage, counterfeit varieties of Oviductus Ranae often appear in the market. However, traditional identification methods cannot accurately differentiate between Oviductus Ranae and its adulterants. In this study, a rapid molecular identification method has been established. The method involves extracting genomic DNA in just 30 s using filter paper purification, species-specific rapid polymerase chain reaction (PCR) amplification, and finally, fluorescence detection of the products. It can accurately identify Oviductus Ranae and its three common adulterants in about 30 min, making the process simple, fast, and highly specific., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Inc. All rights reserved.)

Epigenetic regulation, notably histone post-translational modification (PTM), has emerged as a major transcriptional control of gene expression during cellular stress adaptation. In the present study, we use an acid extraction method to isolate total histone protein and investigate dynamic changes in 23 well-characterized histone methylations/acetylations in the brains of wood frogs subject to 24-h freezing and subsequent 8-h thawed recovery conditions. Our results identify four histone PTMs (H2BK5ac, H3K14ac, H3K4me3, H3K9me2) and three histone proteins (H1.0, H2B, H4) that were significantly ( p < 0.05) responsive to freeze-thaw in freeze-tolerant R. sylvatica brains. Two other permissive modifications (H3R8me2a, H3K9ac) also trended downwards following freezing stress. Together, these data are strongly supportive of the proposed global transcriptional states of hypometabolic freeze tolerance and rebounded thawed recovery. Our findings shed light on the intricate interplay between epigenetic regulation, gene transcription and energy metabolism in wood frogs' adaptive response to freezing stress.

A common goal in evolutionary biology is to discern the mechanisms that produce the astounding diversity of morphologies seen across the tree of life. Aposematic species, those with a conspicuous phenotype coupled with some form of defence, are excellent models to understand the link between vivid colour pattern variations, the natural selection shaping it, and the underlying genetic mechanisms underpinning this variation. Mimicry systems in which species share a conspicuous phenotype can provide an even better model for understanding the mechanisms of colour production in aposematic species, especially if comimics have divergent evolutionary histories. Here we investigate the genetic mechanisms by which mimicry is produced in poison frogs. We assembled a 6.02-Gbp genome with a contig N50 of 310 Kbp, a scaffold N50 of 390 Kbp and 85% of expected tetrapod genes. We leveraged this genome to conduct gene expression analyses throughout development of four colour morphs of Ranitomeya imitator and two colour morphs from both R. fantastica and R. variabilis which R. imitator mimics. We identified a large number of pigmentation and patterning genes differentially expressed throughout development, many of them related to melanophores/melanin, iridophore development and guanine synthesis. We also identify the pteridine synthesis pathway (including genes such as qdpr and xdh) as a key driver of the variation in colour between morphs of these species, and identify several plausible candidates for colouration in vertebrates (e.g. cd36, ep-cadherin and perlwapin). Finally, we hypothesise that keratin genes (e.g. krt8) are important for producing different structural colours within these frogs., (© 2024 John Wiley & Sons Ltd.)

Global warming and environmental pollutants both pose a threat to the behavior and physiology of animals, but research on the combined effects of the two is limited. Atrazine, a widely used herbicide, has toxic effects on organisms. In this study, the effects of environmental concentrations of atrazine exposure (100 μg/L) for seven days on the movement, metabolism and gene expression related to motility of Pelophylax nigromaculatus larvae (GS8) were investigated under global warming. The results showed that compared to the optimal growth temperature (18 °C), atrazine treatment under global warming (21 °C) significantly increased the average speed (about 11.2 times) and maximum acceleration (about 1.98 times) of P. nigromaculatus larvae, altered the relative abundance of 539 metabolites, including Formyl-5-hydroxykynurenamine, 2,4-Dihydroxybenzophenone, and FAPy-adenine, and changed the nucleotide metabolism, pyrimidine metabolism, glycerophospholipid metabolism, and purine metabolism, as well as increased the gene expression of SPLA2 (about 6.46 times) and CHK (about 3.25 times). In summary, atrazine treatment under global warming caused metabolic disorders in amphibian larvae and increased the expression of some movement-related genes in the brain, resulting in abnormally active., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Inc. All rights reserved.)

Reintroduction efforts are increasingly used to mitigate biodiversity losses, but are frequently challenged by inadequate planning and uncertainty. High quality information about population status and threats can be used to prioritize reintroduction and restoration efforts and can transform ad hoc approaches into opportunities for improving conservation outcomes at a landscape scale. We conducted comprehensive environmental DNA (eDNA) and visual encounter surveys to determine the distribution of native and non-native aquatic species in two high-priority watersheds to address key uncertainties-such as the distribution of threats and the status of existing populations-inherent in restoration planning. We then used these occurrence data to develop a menu of potential conservation actions and a decision framework to benefit an endangered vertebrate (foothill yellow-legged frog, Rana boylii) in dynamic stream systems. Our framework combines the strengths of multiple methods, allowing managers and conservation scientists to incorporate conservation science and site-specific knowledge into the planning process to increase the likelihood of achieving conservation goals., (© 2024. The Author(s).)

Harmful cyanobacterial blooms are frequent and intense worldwide, creating hazards for aquatic biodiversity. The potential estrogen-like effect of Microcystin-LR (MC-LR) is a growing concern. In this study, we assessed the estrogenic potency of MC-LR in black-spotted frogs through combined field and laboratory approaches. In 13 bloom areas of Zhejiang province, China, the MC-LR concentrations in water ranged from 0.87 to 8.77 μg/L and were correlated with sex hormone profiles in frogs, suggesting possible estrogenic activity of MC-LR. Tadpoles exposed to 1 μg/L, an environmentally relevant concentration, displayed a female-biased sex ratio relative to controls. Transcriptomic results revealed that MC-LR induces numerous and complex effects on gene expression across multiple endocrine axes. In addition, exposure of male adults significantly increased the estradiol (E2)/testosterone (T) ratio by 3.5-fold relative to controls. Downregulation of genes related to male reproductive endocrine function was also identified. We also showed how MC-LR enhances the expression of specific estrogen receptor (ER) proteins, which induce estrogenic effects by activating the ER pathway and hypothalamic-pituitary-gonadal (HPG) axis. In aggregate, our results reveal multiple lines of evidence demonstrating that, for amphibians, MC-LR is an estrogenic endocrine disruptor at environmentally relevant concentrations. The data presented here support the need for a shift in the MC-LR risk assessment. While hepatoxicity has historically been the focus of MC-LR risk assessments, our data clearly demonstrate that estrogenicity is a major mode of toxicity at environmental levels and that estrogenic effects should be considered for risk assessments on MC-LR going forward.

This study compares the skin structures of Rana kukunoris with two different skin colors living in the same area of Haibei in the Northeastern Qinghai-Tibet Plateau. The skin thickness of the khaki R. kukunoris was significantly greater than that of the brown R. kukunoris (P < 0.01), and significantly more mucous and granular glands were present on the dorsal skin of the khaki frog (P < 0.05). Meanwhile, the melanocytes on the dorsal skin of the brown frog were significantly larger than those on the khaki one (P < 0.05). Morphological changes in the expansion and aggregation of melanocytes seemed to deepen the skin color of R. kukunoris. Moreover, transcriptome sequencing identified tyrosine metabolism, melanogenesis, and riboflavin metabolism as the main pathways involved in melanin formation and metabolism in brown R. kukunoris. TYR, MC1R was upregulated as the skin color of R. kukunoris was deepened and contributed to melanin production and metabolism. In contrast, the khaki frog had significantly more upregulated genes and metabolic pathways related to autoimmunity. The khaki frog appeared to defend against ultraviolet (UV) radiation-induced damage by secreting mucus and small molecular peptides, whereas the brown frog protected itself by distributing a large amount of melanin. Hence, the different skin colors of R. kukunoris might represent different adaptation strategies for survival in the intense UV radiation environment of the Qinghai-Tibet Plateau., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024. Published by Elsevier Inc.)

Background: Cold hardiness is fundamental for amphibians to survive during the extremely cold winter on the Qinghai-Tibet plateau. Exploring the gene regulation mechanism of freezing-tolerant Rana kukunoris could help us to understand how the frogs survive in winter., Results: Transcriptome of liver and muscle of R. kukunoris collected in hibernation and spring were assisted by single molecule real-time (SMRT) sequencing technology. A total of 10,062 unigenes of R. kukunoris were obtained, and 9,924 coding sequences (CDS) were successfully annotated. Our examination of the mRNA response to whole body freezing and recover in the frogs revealed key genes concerning underlying antifreeze proteins and cryoprotectants (glucose and urea). Functional pathway analyses revealed differential regulated pathways of ribosome, energy supply, and protein metabolism which displayed a freeze-induced response and damage recover. Genes related to energy supply in the muscle of winter frogs were up-regulated compared with the muscle of spring frogs. The liver of hibernating frogs maintained modest levels of protein synthesis in the winter. In contrast, the liver underwent intensive high levels of protein synthesis and lipid catabolism to produce substantial quantity of fresh proteins and energy in spring. Differences between hibernation and spring were smaller than that between tissues, yet the physiological traits of hibernation were nevertheless passed down to active state in spring., Conclusions: Based on our comparative transcriptomic analyses, we revealed the likely adaptive mechanisms of R. kukunoris. Ultimately, our study expands genetic resources for the freezing-tolerant frogs., (© 2024. The Author(s).)

Sex determination systems have greatly diversified between amphibians and reptiles, with such as the different sex chromosome compositions within a single species and transition between temperature-dependent sex determination (TSD) and genetic sex determination (GSD). In most sex chromosome studies on amphibians and reptiles, the whole-genome sequence of Xenopous tropicalis and chicken have been used as references to compare the chromosome homology of sex chromosomes among each of these taxonomic groups, respectively. In the present study, we reviewed existing reports on sex chromosomes, including karyotypes, in amphibians and reptiles. Furthermore, we compared the identified genetic linkages of sex chromosomes in amphibians and reptiles with the chicken genome as a reference, which is believed to resemble the ancestral tetrapod karyotype. Our findings revealed that sex chromosomes in amphibians are derived from genetic linkages homologous to various chicken chromosomes, even among several frogs within single families, such as Ranidae and Pipidae. In contrast, sex chromosomes in reptiles exhibit conserved genetic linkages with chicken chromosomes, not only across most species within a single family, but also within closely related families. The diversity of sex chromosomes in amphibians and reptiles may be attributed to the flexibility of their sex determination systems, including the ease of sex reversal in these animals., (© 2023 The Authors. Journal of Experimental Zoology Part A: Ecological and Integrative Physiology published by Wiley Periodicals LLC.)

Recently, populations of Chinese spiny frogs (Quasipaa spinosa), an important amphibian species in China, have decreased, mainly due to a disease caused by the gram-negative bacteria Proteus mirabilis. To elucidate the immune response of the frogs, this study aimed to identify novel candidate genes functionally associated with P. mirabilis infection-induced "rotting skin" disease. Chinese spiny frogs were infected with P. mirabilis, and the skin transcriptome was sequenced using the MGISEQ-2000 platform. A total of 233,965 unigenes were obtained by sequencing, of which 27.23 % were known genes. Screening of differentially expressed genes (DEGs) indicated 210 unigenes differentially expressed after P. mirabilis infection, of which 132 unigenes were up-regulated, and 78 unigenes were down-regulated. Using Kyoto Encyclopedia of Genes and Genomes enrichment analysis, DEGs were identified as enriched in signal pathways, such as oxidative phosphorylation, apoptosis, and the Janus kinase-signal transducer and activator of transcription pathway. Of the DEGs, there was a significant upregulation of the colony stimulating factor 2 receptor beta common subunit, interleukin 2 receptor subunit gamma, cathelicidin antimicrobial peptide, interleukin-17 receptor E, receptor-interacting serine/threonine-protein kinase 3, and pulmonary surfactant-associated protein D immune genes following P. mirabilis infection. Conversely, scavenger receptor cysteine-rich domain-containing group B protein, tumor protein p53 inducible nuclear protein 2, suppressor of cytokine signaling 2, and metalloreductase STEAP3 were significantly downregulated. In conclusion, the first skin transcriptome database of Chinese spiny frogs was established, and several immune genes were identified to elucidate the pathogenic mechanism of "skin rot" in Chinese spiny frogs and other cultured frogs., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Inc. All rights reserved.)

Sex chromosome turnover is the transition between sex chromosomes and autosomes. Although many cases have been reported in poikilothermic vertebrates, their evolutionary causes and genetic mechanisms remain unclear. In this study, we report multiple transitions between the Y chromosome and autosome in the Japanese Tago's brown frog complex. Using chromosome banding and molecular analyses (sex-linked and autosomal single nucleotide polymorphisms, SNPs, from the nuclear genome), we investigated the frogs of geographic populations ranging from northern to southern Japan of two species, Rana tagoi and Rana sakuraii (2n = 26). Particularly, the Chiba populations of East Japan and Akita populations of North Japan in R. tagoi have been, for the first time, investigated here. As a result, we identified three different sex chromosomes, namely chromosomes 3, 7, and 13, in the populations of the two species. Furthermore, we found that the transition between the Y chromosome (chromosome 7) and autosome was repeated through hybridization between two or three different populations belonging to the two species, followed by restricted chromosome introgression. These dynamic sex chromosome turnovers represent the first such findings in vertebrates and imply that speciation associated with inter- or intraspecific hybridization plays an important role in sex chromosome turnover in frogs.

Reproduction of water frog hybrids Pelophylax esculentus (Pelophylax ridibundus x Pelophylax lessonae) is associated with hemiclonal reproduction and backcrossing. The hemiclonal mode of reproduction occurs within P. esculentus allodiploids. In this case, the unrecombined genome of one parental species is transmitted to the offspring after premeiotic elimination of the chromosome set of the second parental species. Usually, the chromosome set of P. lessonae is eliminated, and the altered genome of P. ridibundus is passed on to the progeny. The hemiclonal inheritance within diploid Pelophylax esculentus hybrids may be accompanied by certain aberrations of premeiotic elimination. As a result, the formation of P. ridibundus specimens with introgressions of the P. lessonae genetic material, or the formation of recombinant hybrids occurs, depending on which of the parental species backcrossing takes place. The aim of our study is to describe the aberration of premeiotic elimination within the water frog P. esculentus complex detected by the nuclear gene Ldh-B inheritance, with an attempt to find out the causes of this phenomenon. It has been established that aberrations of premeiotic elimination are widespread, but only within populations of water frog from the river system of Upper Dnieper within Ukraine. The highest level of introgression takes place in the water frog populations within Kiev metropolis under conditions of expressed anthropogenization, while the maximum frequency of recombinants was detected within populations from the basin of Desna River, that has preserved native ecosystems. It was demonstrated that the frequency of premeiotic aberrations does not correlate with the intensity of interspecific water frog hybridization. Populations with introgressions are more common than populations with recombinants, however, within the latter, the frequency of recombination events is higher. The primary factor of gametogenesis aberrations, most likely, is the genetic characteristics of the local populations of parental species, since unambiguous explanations of this phenomenon based on the action of environmental stress (pollution of water systems) are not obvious., (© 2024. The Author(s), under exclusive licence to Springer Nature Switzerland AG.)

A peculiar physiological characteristic of the Chinese brown frog (Rana dybowskii) is that its oviduct dilates during pre-brumation rather than during the breeding season. This research aimed to examine the expression of genes connected with lipid synthesis and metabolism in the oviduct of R. dybowskii during both the breeding season and pre-brumation. We observed significant changes in the weight and size of the oviduct between the breeding season and pre-brumation. Furthermore, compared to the breeding season, pre-brumation exhibited significantly lower triglyceride content and a marked increase in free fatty acid content. Immunohistochemical results revealed the spatial distribution of triglyceride synthase (Dgat1), triglyceride hydrolase (Lpl and Hsl), fatty acid synthase (Fasn), and fatty acid oxidases (Cpt1a, Acadl, and Hadh) in oviductal glandular cells and epithelial cells during both the breeding season and pre-brumation. While the mRNA levels of triglycerides and free fatty acid synthesis genes (dgat1 and fasn) did not show a significant difference between the breeding season and pre-brumation, the mRNA levels of genes involved in triglycerides and free fatty acid metabolism (lpl, cpt1a, acadl, acox and hadh) were considerably higher during pre-brumation. Furthermore, the R. dybowskii oviduct's transcriptomic and metabolomic data confirmed differential expression of genes and metabolites enriched in lipid metabolism signaling pathways during both the breeding season and pre-brumation. Overall, these results suggest that alterations in lipid synthesis and metabolism during pre-brumation may potentially influence the expanding size of the oviduct, contributing to the successful overwintering of R. dybowskii.

The oviduct of the Chinese brown frog ( Rana dybowskii ) expands in prehibernation rather than in prespawning, which is one of the physiological phenomena that occur in the preparation for hibernation. Steroid hormones are known to regulate oviductal development. Cholesterol synthesis and steroidogenesis may play an important role in the expansion of the oviduct before hibernation. In this study, we investigated the expression patterns of the markers that are involved in the de novo steroid synthesis pathway in the oviduct of R. dybowskii during prespawning and prehibernation. According to histological analysis, the oviduct of R. dybowskii contains epithelial cells, glandular cells, and tubule lumens. During prehibernation, oviductal pipe diameter and weight were significantly larger than during prespawning. 3-Hydroxy-3-methylglutaryl CoA reductase (HMGCR), low-density lipoprotein receptor (LDLR), steroidogenic acute regulatory protein (StAR), cytochrome P450 cholesterol side-chain cleavage enzyme (P450scc), and steroidogenic factor 1 (SF-1) were detected in epithelial cells in prehibernation and glandular cells during prespawning. HMGCR, LDLR, StAR, and P450scc protein expression levels were higher in prehibernation than during prespawning, but the SF-1 protein expression level did not significantly differ. HMGCR, LDLR, StAR, P450scc ( CYP11A1 ), and SF-1 ( NR5A1 ) mRNA expression levels were significantly higher in prehibernation compared with prespawning. The transcriptome results showed that the steroid synthesis pathway was highly expressed during prehibernation. Existing results indicate that the oviduct is able to synthesize steroid hormones using cholesterol, and that steroid hormones may affect the oviductal functions of R. dybowskii.

Hybrid taxa from the genus Pelophylax can propagate themselves in a modified way of sexual reproduction called hybridogenesis ensuring the formation of clonal gametes containing the genome of only one parental (host) species. Pelophylax grafi from South-Western Europe is a hybrid composed of P. ridibundus and P. perezi genomes and it lives with a host species P. perezi (P-G system). Yet it is unknown, whether non-Mendelian inheritance is fully maintained in such populations. In this study, we characterize P. perezi and P. grafi somatic karyotypes by using comparative genomic hybridization, genomic in situ hybridization, fluorescent in situ hybridization, and actinomycin D-DAPI. Here, we show the homeology of P. perezi and P. grafi somatic karyotypes to other Pelophylax taxa with 2n = 26 and equal contribution of ridibundus and perezi chromosomes in P. grafi which supports F1 hybrid genome constitution as well as a hemiclonal genome inheritance. We show that ridibundus chromosomes have larger regions of interstitial (TTAGGG)n repeats flanking the nucleolus organizing region on chromosome no. 10 and a high quantity of AT pairs in the centromeric regions. In P. perezi, we found species-specific sequences in metaphase chromosomes and marker structures in lampbrush chromosomes. Pericentromeric RrS1 repeat sequence was present in perezi and ridibundus chromosomes, but the blocks were stronger in ridibundus. Various cytogenetic techniques applied to the P-G system provide genome discrimination between ridibundus and perezi chromosomal sets. They could be used in studies of germ-line cells to explain patterns of clonal gametogenesis in P. grafi and broaden the knowledge about reproductive strategies in hybrid animals., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

We describe a new species of ranid frog of the genus Nidirana from northeast India based on morphological, molecular and acoustic evidence. The new species is phenotypically distinct from its congeners by a combination of morphological characters: body robust with SVL 46.559.1 mm (n= 3) in adult males and SVL 60.666.0 mm (n= 2) in adult females; a pair of subgular vocal sacs and two patches of nuptial pad on the first finger in adult males; toe tips slightly dilated and oval; circum-marginal grooves present on all toes; dorsal skin with scattered small tubercles. A pale cream-coloured mid-dorsal line from the snout tip to the vent is present. Phylogenetically, the new species differs from its congeners by a genetic divergence of 3.48.0% and 7.712.4% in 16S and COI genes respectively. Furthermore, the new species can be differentiated from its congeners by its advertisement call, which consists of two different types of notes, call duration (0.580.92 s) and dominant frequency of the call (473.7 Hz). The discovery of a new species validates the presence of the genus Nidirana from India and emphasizes the importance of exploring specialized habitats such as marshlands, which are often overlooked.

Amphibian species have the largest genome size enriched with repetitive sequences and relatively similar karyotypes. Moreover, many amphibian species frequently hybridize causing nuclear and mitochondrial genome introgressions. In addition, hybridization in some amphibian species may lead to clonality and polyploidization. All such events were found in water frogs from the genus Pelophylax. Among the species within the genus Pelophylax, P. esculentus complex is the most widely distributed and well-studied. This complex includes two parental species, P. ridibundus and P. lessonae, and their hybrids, P. esculentus, reproducing hemiclonally. Parental species and their hybrids have similar but slightly polymorphic karyotypes, so their precise identification is still required. Here, we have developed a complete set of 13 chromosome painting probes for two parental species allowing the precise identification of all chromosomes. Applying chromosomal painting, we identified homologous chromosomes in both parental species and orthologous chromosomes in their diploid hemiclonal hybrids. Comparative painting did not reveal interchromosomal exchanges between the studied water frog species and their hybrids. Using cross-specific chromosome painting, we detected unequal distribution of the signals along chromosomes suggesting the presence of species-specific tandem repeats. Application of chromosomal paints to the karyotypes of hybrids revealed differences in the intensity of staining for P. ridibundus and P. lessonae chromosomes. Thus, both parental genomes have a divergence in unique sequences. Obtained chromosome probes may serve as a powerful tool to unravel chromosomal evolution in phylogenetically related species, identify individual chromosomes in different cell types, and investigate the elimination of chromosomes in hybrid water frogs., (© 2023. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Sex chromosomes are popularized as a special role in driving speciation. However, the empirical evidence from natural population processes has been limited to organisms with degenerated sex chromosomes, where hemizygosity is mainly considered to act as the driver of reproductive isolation. Here, we examined several hybrid zones of torrent frog Amolops mantzorum species complex, using an approach by mapping species-diagnostic loci onto the reference genome to compare sex-linked versus autosomal patterns of introgression. We find little support in sex-linked incompatibilities for large X-effects for these populations in hybrid zones with homomorphic sex chromosomes, due to the absence of the hemizygous effects. As expected, the large X-effects were not found in those with heteromorphic but newly evolved sex chromosomes, owing to the absence of strong genetic differences between X and Y chromosomes. The available data so far on amphibians suggest little role for sex-linked genes in speciation. The large X-effects in those with nascent sex chromosomes may not be as ubiquitous as presumed across the animal kingdom., (© 2023 John Wiley & Sons Ltd.)

Skin microbial communities are an essential part of host health and can play a role in mitigating disease. Host and environmental factors can shape and alter these microbial communities and, therefore, we need to understand to what extent these factors influence microbial communities and how this can impact disease dynamics. Microbial communities have been studied in amphibian systems due to skin microbial communities providing some resistance to the amphibian chytrid fungus, Batrachochytrium dendrobatidis . However, we are only starting to understand how host and environmental factors shape these communities for amphibians. In this study, we examined whether amphibian skin bacterial communities differ among host species, host infection status, host developmental stage, and host habitat. We collected skin swabs from tadpoles and adults of three Ranid frog species ( Lithobates spp.) at the Mianus River Gorge Preserve in Bedford, New York, USA, and used 16S rRNA gene amplicon sequencing to determine bacterial community composition. Our analysis suggests amphibian skin bacterial communities change across host developmental stages, as has been documented previously. Additionally, we found that skin bacterial communities differed among Ranid species, with skin communities on the host species captured in streams or bogs differing from the communities of the species captured on land. Thus, habitat use of different species may drive differences in host-associated microbial communities for closely-related host species., Competing Interests: Christopher Nagy is the Director of Research & Education at the Mianus River Gorge Preserve., (© 2023 Gajewski et al.)

Adiponectin receptor 1 (AdipoR1) and Adiponectin receptor 2 (AdipoR2) can bind to adiponectin (AdipoQ) secreted by adipose tissue to participate in various physiological functions of the body. In order to explore the role of AdipoR1 and AdipoR2 in amphibians infected by Aeromonas hydrophila ( Ah ), the genes adipor1 and adipor2 of Rana dybowskii were cloned by reverse transcription-polymerase chain reaction (RT-PCR) and analyzed by bioinformatics. The tissue expression difference of adipor1 and adipor2 was analyzed by real-time fluorescence quantitative polymerase chain reaction (qRT-PCR), and an inflammatory model of R . dybowskii infected by Ah was constructed. The histopathological changes were observed by hematoxylin-eosin staining (HE staining); the expression profiles of adipor1 and adipor2 after infection were dynamically detected by qRT-PCR and Western blotting. The results show that AdipoR1 and AdipoR2 are cell membrane proteins with seven transmembrane domains. Phylogenetic tree also shows that AdipoR1 and AdipoR2 cluster with the amphibians in the same branch. qRT-PCR and Western blotting results show that adipor1 and adipor2 were up-regulated at different levels of transcription and translation upon Ah infection, but the response time and level were different. It is speculated that AdipoR1 and AdipoR2 participate in the process of bacterial immune response, providing a basis for further exploring the biological functions of AdipoR1 and AdipoR2 in amphibians.

Although numerous investigations on the adverse impact of Cr and Pb have been performed, studies on intestinal homeostasis in amphibians are limited. Here, single and combined effects of Cr (104 μg/L) and Pb (50 μg/L) on morphological and histological features, bacterial community, digestive enzymes activities, as well as transcriptomic profile of intestines in Rana chensinensis tadpoles were assessed. Significant decrease in the relative intestine length (intestine length/snout-to-vent length, IL/SVL) was observed after exposure to Pb and Cr/Pb mixture. Intestinal histology and digestive enzymes activities were altered in metal treatment groups. In addition, treatment groups showed significantly increased bacterial richness and diversity. Tadpoles in treatment groups were observed to have differential gut bacterial composition from controls, especially for the abundance of phylum Proteobacteria, Firmicutes, Verrucomicrobia, Actinobacteria, and Fusobacteria as well as genus Citrobacter, Anaerotruncus, Akkermansia, and Alpinimonas. Moreover, transcriptomic analysis showed that the transcript expression profiles of GPx and SOD isoforms responded differently to Cr and/or Pb exposure. Besides, transcriptional activation of pro-apoptotic and glycolysis-related genes, such as Bax, Apaf 1, Caspase 3, PK, PGK, TPI, and GPI were detected in all treatment groups but downregulation of Bcl2 in Pb and Cr/Pb mixture groups. Collectively, these results suggested that Cr and Pb exposure at environmental relevant concentration, alone and in combination, could disrupt intestinal homeostasis of R. chensinensis tadpoles., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS) are widely detected in the environment and wild animals, thus posing a threat to wildlife and public health; however, knowledge about their immunotoxicity and the underlying mechanism remains limited. In the present study, male black-spotted frogs (Rana nigromaculata) were exposed to environmentally relevant concentrations (0, 1, and 10 μg/L) of PFOA or PFOS for 21 days; subsequently, biochemical analysis, molecular docking, and gene expression determination were conducted. The results indicated that exposure to 10 μg/L PFOA decreased the serum levels of immunoglobulin A. PFOS exposure significantly increased the hepatic levels of interleukin-1β, interleukin-6, tumor necrosis factor-α, interferon-γ, and nitric oxide; but PFOA significantly increased the levels of only tumor necrosis factor-α. Furthermore, PFOA and PFOS exposure significantly decreased the activity of inducible nitric oxide synthase and total nitric oxide synthase. IBRv2 analysis indicated that PFOA and PFOS had a similar effect on these immune indicators, but PFOS was more toxic than PFOA. Molecular docking revealed that PFOA and PFOS can bind to nuclear factor-κB (NF-κB) by forming stable hydrogen bonds. PFOA and PFOS exposure upregulated the gene expression of NF-κB and its downstream genes. Significant correlations between the expression of genes involved in the NF-κB pathway and immune-related indicators suggests that PFOA- and PFOS-induced immunotoxicity was associated with the activation of NF-κB. Our findings provide novel insights into the potential role of NF-κB in immunotoxicity induced by PFOA and PFOS in frogs., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Elsevier Ltd. All rights reserved.)

Genetic sex-determination features male (XX/XY) or female heterogamety (ZZ/ZW). To identify similarities and differences in the molecular evolution of sex-linked genes between these systems, we directly compared the sex chromosome systems existing in the frog Glandirana rugosa . The heteromorphic X/Y and Z/W sex chromosomes were derived from chromosomes 7 (2n = 26). RNA-Seq, de novo assembly, and BLASTP analyses identified 766 sex-linked genes. These genes were classified into three different clusters (XW/YZ, XY/ZW, and XZ/YW) based on sequence identities between the chromosomes, probably reflecting each step of the sex chromosome evolutionary history. The nucleotide substitution per site was significantly higher in the Y- and Z-genes than in the X- and W- genes, indicating male-driven mutation. The ratio of nonsynonymous to synonymous nucleotide substitution rates was higher in the X- and W-genes than in the Y- and Z-genes, with a female bias. Allelic expression in gonad, brain, and muscle was significantly higher in the Y- and W-genes than in the X- and Z-genes, favoring heterogametic sex. The same set of sex-linked genes showed parallel evolution across the two distinct systems. In contrast, the unique genomic region of the sex chromosomes demonstrated a difference between the two systems, with even and extremely high expression ratios of W/Z and Y/X, respectively.

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Poulsen, Jackman, Hansen and Helbing.)

Brevinins exhibit a wide range of structural features and strong biological activities. Brevinin-2, derived from several amphibians, has shown antimicrobial activities. However, little is known about the wound-healing activity of brevinin-2. In this study, brevinin-2 cDNA was identified from the skin transcriptome of the dark-spotted frog (Pelophylax nigromaculatus) and it comprises a signal peptide, a propeptide, and a mature peptide. Sequence alignment with brevinin-2 derived from other amphibians showed variability of the mature peptide, and the presence of a C-terminal cyclic heptapeptide domain (Cys-Lys-Xaa4-Cys) in the mature peptide. Dark-spotted frog brevinin-2 belonged to the brevinin-2 cluster and was closely related to brevinin-2HB1 from Pelophylax hubeiensis. Synthetic dark-spotted frog brevinin-2 mature peptide (brevinin-2PN) exhibited antibacterial activity against several pathogens by destroying cell membrane integrity and hydrolysis of genomic DNA. Brevinin-2PN exhibited significant wound-healing activity by accelerating the healing of human skin fibroblast cell scratches, influencing cell migration, and stimulating gene expression of growth factors., (Copyright © 2022 Elsevier Ltd. All rights reserved.)

The Heilongjiang brown frog (Rana amurensis) is widely used in traditional Chinese medicine. In particular, the oviduct and skin have been developed into various health products. However, limited numbers of complete genomes of amphibian species have been reported, excluding the Heilongjiang brown frog. Here, the transcriptomes of 45 samples from the liver, spleen, heart, ovaries, thigh muscles, skin, oviduct, stomach and intestine of five Heilongjiang brown frog were reassembled and analyzed. A total of 1,085,532 unigenes with an average length of 676.6 bp and N50 of 722 bp were obtained. Comparative transcriptomics of different tissues detected tissue-specific expression. There were 3248 differentially expressed genes (DEGs) in the ovary, and the number of unique DEGs between the ovary and spleen was the largest. The results of DEGs enrichment showed there were many pathways and items related to protein synthesis and metabolism in the oviduct. The DEGs of the skin were enriched with many bacterial defense items, indicating that there were a large number of antimicrobial peptides in the skin. Thus, these were suitable as biological sources for the development and extraction of antimicrobial peptides. Through the assembly of transcriptome sequencing data and functional annotation of the Heilongjiang brown frog genome, this study provides reference materials for further exploring and utilizing functional gene resources of frogs and lays a foundation for medical research and the development of new products., (© 2022. The Author(s).)

Pseudis tocantins is the only frog species of the hylid genus Pseudis that possesses highly heteromorphic sex chromosomes. Z and W chromosomes of Ps. tocantins differ in size, morphology, position of the nucleolar organizer region (NOR) and the amount and distribution of heterochromatin. A chromosomal inversion and heterochromatin amplification on the W chromosome were previously inferred to be involved in the evolution of this sex chromosome pair. Despite these findings, knowledge related to the molecular composition of the large heterochromatic band of this W chromosome is restricted to the PcP190 satellite DNA, and no data are available regarding the gene content of either the W or the Z chromosome of Ps. tocantins. Here, we sequenced microdissected Z and W chromosomes of this species to further resolve their molecular composition. Comparative genomic analysis suggests that Ps. tocantins sex chromosomes are likely homologous to chromosomes 4 and 10 of Xenopus tropicalis. Analyses of the repetitive DNA landscape in the Z and W assemblies allowed for the identification of several transposable elements and putative satellite DNA sequences. Finally, some transposable elements from the W assembly were found to be highly diverse and divergent from elements found elsewhere in the genome, suggesting a rapid amplification of these elements on the W chromosome., (© 2022 European Society for Evolutionary Biology.)

Freeze-tolerant species, such as wood frogs (Rana sylvatica), are susceptible to multiple co-occurring stresses that they must overcome to survive. Freezing is accompanied by mechanical stress and dehydration due to ice crystal formation in the extracellular space, ischemia/anoxia due to interruption in blood flood, and hyperglycemia due to cryoprotective measures. Wood frogs can survive dehydration, anoxia, and high glucose stress independently of freezing, thereby creating a multifactorial model for studying freeze-tolerance. Oxidative stress and high glucose levels favors the production of pro-oxidant molecules and advanced glycation end product (AGE) adducts that could cause substantial cellular damage. In this study, the involvement of the high mobility group box 1 (HMGB1)-AGE/RAGE (receptor for AGE) axis and the regulation of ETS1 and EGR1-mediated angiogenic responses were investigated in liver of wood frogs expose to freeze/thaw, anoxia/reoxygenation and dehydration/rehydration treatments. HMGB1 and not AGE-adducts are likely to induce the activation of ETS1 and EGR1 via the RAGE pathway. The increase in nuclear localization of both ETS1 and EGR1, but not DNA binding activity in response to stress hints to a potential spatial and temporal regulation in inducing angiogenic factors. Freeze/thaw and dehydration/rehydration treatments increase the levels of both pro- and anti-angiogenic factors, perhaps to prepare for the distribution of cryoprotectants or enable the repair of damaged capillaries and wounds when needed. Overall, wood frogs appear to anticipate the need for angiogenesis in response to freezing and dehydration but not anoxic treatments, probably due to mechanical stress associated with the two former conditions., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Delineating conservation units (CUs, e.g., evolutionarily significant units, ESUs, and management units, MUs) is critical to the recovery of declining species because CUs inform both listing status and management actions. Genomic data have strengths and limitations in informing CU delineation and related management questions in natural systems. We illustrate the value of using genomic data in combination with landscape, dispersal and occupancy data to inform CU delineation in Nevada populations of the Great Basin Distinct Population Segment of the Columbia spotted frog (Rana luteiventris). R. luteiventris occupies naturally fragmented aquatic habitats in this xeric region, but beaver removal, climate change and other factors have put many of these populations at high risk of extirpation without management intervention. We addressed three objectives: (i) assessing support for ESUs within Nevada; (ii) evaluating and revising, if warranted, the current delineation of MUs; and (iii) evaluating genetic diversity, effective population size, adaptive differentiation and functional connectivity to inform ongoing management actions. We found little support for ESUs within Nevada but did identify potential revisions to MUs based on unique landscape drivers of connectivity that distinguish these desert populations from those in the northern portion of the species range. Effective sizes were uniformly small, with low genetic diversity and weak signatures of adaptive differentiation. Our findings suggest that management actions, including translocations and genetic rescue, might be warranted. Our study illustrates how a carefully planned genetic study, designed to address priority management goals that include CU delineation, can provide multiple insights to inform conservation action., (© 2022 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

In sharp contrast to birds and mammals, in numerous cold-blooded vertebrates, sex chromosomes have been described as homomorphic. This sex chromosome homomorphy has been suggested to result from the high turnovers often observed across deeply diverged clades. However, little is known about the tempo and mode of sex chromosome evolution among the most closely related species. Here, we examined the evolution of sex chromosome among nine species of the torrent frog genus Amolops . We analyzed male and female GBS and RAD-seq from 182 individuals and performed PCR verification for 176 individuals. We identified signatures of sex chromosomes involving two pairs of chromosomes. We found that sex-chromosome homomorphy results from both turnover and X-Y recombination in the Amolops species, which simultaneously exhibits heterogeneous evolution on homologous and non-homologous sex chromosomes. A low turnover rate of non-homologous sex chromosomes exists in these torrent frogs. The ongoing X-Y recombination in homologous sex chromosomes will act as an indispensable force in preventing sex chromosomes from differentiating.

Glandirana rugosa is known to include several geographic groups differing in sex chromosomes, and has been proven to be paraphyletic in mitochondrial phylogeny with respect to G. susurra. By analyzing genetic and morphological variation in a large number of individuals of Glandirana, we studied their taxonomic relationships. A mitochondrial DNA phylogeny, with the G. tientaiensis as outgroup, revealed two major lineages containing respectively (1) the East group of G. rugosa, G. susurra, and the Central and Southeast-Kyushu groups of G. rugosa; and (2) G. emeljanovi, and the North and West groups of G. rugosa. In contrast, in a nuclear DNA phylogeny based on SNP data, lineages of (1) G. susurra and East group, and (2) the remaining groups of G. rugosa and G. emeljanovi, were split, indicating a distinct status of the East group among G. rugosa. In adult morphology, there were only minor differences between the East group and the remaining groups of G. rugosa, but in larvae, the East group had significantly more sparse skin glands than the others. The exact type locality of G. rugosa is most probably in western Japan, not including the range of the East group. From these results, we describe the East group as a new species, G. reliquia, distinct from the remaining groups of G. rugosa. The new species with sexually homomorphic chromosomes is thought to represent a basic stock of Japanese Glandirana, which existed far before G. rugosa originated.

Secondary sympatry amongst sister lineages is strongly associated with genetic and ecological divergence. This pattern suggests that for closely related species to coexist in secondary sympatry, they must accumulate differences in traits that mediate ecological and/or reproductive isolation. Here, we characterized inter- and intraspecific divergence in three giant tree frog species whose distributions stretch across West and Central Africa. Using genome-wide single-nucleotide polymorphism data, we demonstrated that species-level divergence coincides temporally and geographically with a period of large-scale forest fragmentation during the late Pliocene. Our environmental niche models further supported a dynamic history of climatic suitability and stability, and indicated that all three species occupy distinct environmental niches. We found modest morphological differentiation amongst the species with significant divergence in tympanum diameter and male advertisement call. In addition, we confirmed that two species occur in secondary sympatry in Central Africa but found no evidence of hybridization. These patterns support the hypothesis that cycles of genetic exchange and isolation across West and Central Africa have contributed to globally significant biodiversity. Furthermore, divergence in both ecology and reproductive traits appear to have played important roles in maintaining distinct lineages. At the intraspecific level, we found that climatic refugia, precipitation gradients, marine incursions, and potentially riverine barriers generated phylogeographic structure throughout the Pleistocene and into the Holocene. Further studies examining phenotypic divergence and secondary contact amongst these geographically structured populations may demonstrate how smaller scale and more recent biogeographic barriers contribute to regional diversification., (© 2021 John Wiley & Sons Ltd.)

Sex chromosomes in poikilothermal vertebrates are characterized by rapid and diverse evolution at the species or population level. Our previous study revealed that the Taiwanese frog Odorrana swinhoana (2n = 26) has a unique system of multiple sex chromosomes created by three sequential translocations among chromosomes 1, 3, and 7. To reveal the evolutionary history of sex chromosomes in the Odorrana species complex, we first identified the original, homomorphic sex chromosomes, prior to the occurrence of translocations, in the ancestral-type population of O. swinhoana. Then, we extended the investigation to a closely related Japanese species, Odorrana utsunomiyaorum, which is distributed on two small islands. We used a high-throughput nuclear genomic approach to analyze single-nucleotide polymorphisms and identify the sex-linked markers. Those isolated from the O. swinhoana ancestral-type population were found to be aligned to chromosome 1 and showed male heterogamety. In contrast, almost all the sex-linked markers isolated from O. utsunomiyaorum were heterozygous in females and homozygous in males and were aligned to chromosome 9. Morphologically, we confirmed chromosome 9 to be heteromorphic in females, showing a ZZ-ZW sex determination system, in which the W chromosomes were heterochromatinized in a stripe pattern along the chromosome axis. These results indicated that after divergence of the two species, the ancestral homomorphic sex chromosome 1 underwent highly rapid and diverse evolution, i.e., sequential translocations with two autosomes in O. swinhoana, and turnover to chromosome 9 in O. utsunomiyaorum, with a transition from XY to ZW heterogamety and change to heteromorphy., (© 2022 Japanese Society of Developmental Biologists.)

Background: Hybrid taxa exist in nature, but their fitness can vary greatly. Hybrids are usually thought to have lower viability and survival rate than parental species due to the occurrence of two different genomes and divergent evolution in each species. On the other hand, the hybrid vigour of the F1 generation may give hybrids an advantage in mixed populations where they have to live and compete with parental taxa. Post-zygotic selection with endogenous genetic mechanisms may be a significant evolutionary force in hybrid formation. Here we tested principles of post-zygotic reproductive dynamics in mixed populations of Pelophylax water frogs that would help us understand the origin and maintenance of such systems., Methods: Within experimental crosses, we combined various diploid Pelophylax genotypes resulting in 211 families. Statistical analysis of progeny was used to measure fertilization success, the rate of embryonic/tadpole mortality and the overall survival of the progeny till the time of metamorphosis. Using Generalized Estimating Equations models and variables defined by a mother/father included in mate pairs, we tested which factor best explains the successful embryonal development., Results: The development of Pelophylax offspring significantly varied in survival rate and morphological malformations. These post-zygotic reproductive dynamics were driven by parental combinations of species pairs. The best values in the proportion of developing eggs, embryos, tadpoles and overall survival showed progeny of homospecific P. lessonae crosses. Total survival rates were relatively similar between L-E and R-E population systems but much lower than homospecific crosses in parental taxa. However, once the early stages passed this period, tadpoles mostly of hybrid hemiclonal origin performed even better than pure P. ridibundus progeny. Hybrid × hybrid crosses showed the highest mortality values. Statistical testing revealed that high mortality affected paternal genetic input., Conclusions: Combined three water frog taxa and both sexes provided patterns of post-zygotic reproduction dynamics of early development in the widespread population systems in Central Europe. The results further showed high survival rates of hybrid F1s created de novo from parental species despite significant divergence between P. ridibundus and P. lessonae DNA. Potential conservation measures of sexual-asexual systems in natural populations are discussed., Competing Interests: The authors declare no conflict of interest., (© 2022 The Author(s). Published by IMR Press.)

The Black-spotted Pond Frog, Pelophylax nigromaculatus, is widely distributed across mainland China, Korean Peninsula, and Japan. The taxonomic relationships among P. n. nigromaculatus, Rana nigromaculata mongolia (sensu P. n. mongolicus), and P. tenggerensis have long been ambiguous. Here we examine the topotype specimens of P. tenggerensis and R. n. mongolia, and provide phylogenic analyses based on four mitochondrial DNA sequences. The combined evidences from morphology and molecular phylogeny have shown the distinct specific-level of P. n. mongolicus that distant from P. nigromaculatus, while indicating the homogeneity between P. n. mongolicus and P. tenggerensis. Thus, we suggest elevating P. n. mongolicus as a full species Pelophylax mongolicus comb. nov., and place P. tenggerensis to be a secondary synonym of P. mongolicus comb. nov.

To date, all captive breeding of the dusky gopher frog, Lithobates sevosus, a federally listed endangered species, has been accomplished using in vitro fertilization (IVF). Here, we describe multievent and highly fecund captive reproduction of dusky gopher frogs driven solely by natural environmental factors. Six pairs of L. sevosus were kept in a 3.7 × 4.4 m 2 outdoor enclosure designed to resemble their natural breeding habitat, which included a pool and three artificial burrows. Modifications to the enclosure that simulated temperatures and conditions within their natural range during winter were added in October and removed in late February. Following a warm, rainy period, five egg masses were laid between March 5 and 11, 2020. The number of oocytes per egg mass was 2300 ± 409 (range = 1341-3565), with the total across all five egg masses being 11,501. Of these oocytes, the hatching rate was 68.58 ± 10.05% (range = 37.53%-95.59%), with a total of 7887 successful hatchlings overall. Clutch sizes were similar to those in the wild and greater than those typically produced using IVF; thus, natural breeding can substantially increase the number of frogs available for reintroduction programs. Although assisted reproductive technologies such as IVF will continue to be useful for ensuring the success of strategic genetic pairings of captive L. sevosus, the new tool of nonassisted reproduction in specifically designed outdoor enclosures is an important advancement for the conservation and recovery of this endangered species., (© 2022 The Authors. Zoo Biology published by Wiley Periodicals LLC.)

Sex chromosomes constantly exist in a dynamic state of evolution: rapid turnover and change of heterogametic sex during homomorphic state, and often stepping out to a heteromorphic state followed by chromosomal decaying. However, the forces driving these different trajectories of sex chromosome evolution are still unclear. The Japanese frog Glandirana rugosa is one taxon well suited to the study on these driving forces. The species has two different heteromorphic sex chromosome systems, XX-XY and ZZ-ZW, which are separated in different geographic populations. Both XX-XY and ZZ-ZW sex chromosomes are represented by chromosome 7 (2n = 26). Phylogenetically, these two systems arose via hybridization between two ancestral lineages of West Japan and East Japan populations, of which sex chromosomes are homomorphic in both sexes and to date have not yet been identified. Identification of the sex chromosomes will give us important insight into the mechanisms of sex chromosome evolution in this species. Here, we used a high-throughput genomic approach to identify the homomorphic XX-XY sex chromosomes in both ancestral populations. Sex-linked DNA markers of West Japan were aligned to chromosome 1, whereas those of East Japan were aligned to chromosome 3. These results reveal that at least two turnovers across three different sex chromosomes 1, 3 and 7 occurred during evolution of this species. This finding raises the possibility that cohabitation of the two different sex chromosomes from ancestral lineages induced turnover to another new one in their hybrids, involving transition of heterogametic sex and evolution from homomorphy to heteromorphy., (© 2022 John Wiley & Sons Ltd.)

Frogs in the family Ranidae are diverse in Asia and are thought to have dispersed to the Sahul Shelf approximately 10 million years ago, where they radiated into more than a dozen species. Ranid species in the intervening oceanic islands of Wallacea, such as Hylarana florensis and H. elberti from the Lesser Sundas and H. moluccana from eastern Wallacea, are assumed to belong to the subgenus Papurana, yet this has not been confirmed with molecular data. We analyzed mitochondrial DNA of Hylarana species from five islands spanning the reported ranges of H. florensis and H. elberti and compared them to confirmed Papurana species and closely related subgenera within Hylarana. We find that the Lesser Sunda H. florensis and H. elberti form a clade that is sister to the rest of the Australo-Papuan Papurana assemblage. Species delimitation analyses and divergence time estimates suggest that populations of H. florensis on Lombok may be distinct from those on Flores at the species level. Likewise, populations of H. elberti on Sumba and Timor may be distinct from each other and from those on Wetar, tshe type locality of H. elberti. Samples from Babar Island thought to be members of H. elberti in fact belong to the wide-ranging H. daemeli, which occurs in northern Australia, across New Guinea, and on the neighboring island of Tanimbar. These results suggest that the Lesser Sundas may have served as a stepping-stone for colonization of the Sahul Shelf and that species diversity of Papurana frogs is underestimated in the Lesser Sundas.

The characteristics of acoustic signals co-evolve with preferences of the auditory sensory system. However, how the brain perceives call variations and whether it can reveal phylogenetic relationships among signalers remains poorly understood. Here, we recorded the neural signals from the Emei music frogs (Nidirana daunchina) in response to broadcasted calls of five different species of the same genus. We found that responses in terms of the different amplitudes of various event-related potential (ERP) components were correlated with diversification trends in acoustic signals, as well as phylogenetic relationships between N. daunchina and heterospecific callers. Specifically, P2 decreased gradually along the ordinal decline of similarities in acoustic characteristics of calls compared with those from conspecifics. Moreover, P3a amplitudes showed increasing trends in correspondence with callers' genetic distances from the subject species. These observations collectively support the view that neural activities in music frogs can reflect call variations and phylogenetic relationships within the genus Nidirana., (© 2022. The Author(s).)

The oviduct of Chinese brown frog (Rana dybowskii) displays seasonal morphological and functional changes, which expands specifically during pre-brumation. To uncover the molecular mechanism underlying this phenomenon, we firstly confirmed the increased weight and enlarged diameter of the oviduct in pre-brumation by morphological observation. Interestingly, the glycogen content in the oviduct increased significantly during pre-brumation, indicating Rana dybowskii stores energy in the oviduct before brumation. Transcriptome analysis further identified the differentially expressed genes in the synthesis and metabolism pathways of carbohydrates in the oviduct during pre-brumation. Based on that evidence, we focused on the mRNA and protein expression of glycogenic genes in the oviduct of Rana dybowskii. qPCR confirmed that the expression of glycolysis and glycogenesis-related genes were up-regulated while gluconeogenesis-related genes were down-regulated during pre-brumation. Western blot data showed that glucose transporter GLUT1 and glycogen synthesis-regulation proteins including GYS, and p-GSK-3β were highly expressed in the oviduct during pre-brumation. Moreover, immunohistochemical data showed that GLUT1, GYS, p-GYS, GSK-3β and p-GSK-3β were expressed regionally in the oviduct of Rana dybowskii. The data suggests that glycogen synthesis may be involved in the oviductal expansion of Rana dybowskii during the pre-brumation., Competing Interests: Declaration of competing interest The author(s) declare that they have no conflict of interest., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Leptin receptor overlapping transcript (LepROT) plays multiple roles in the regulation of immune systems. However, very little information is available about the anti-infectious mechanisms of amphibians LepROT. In this study, the cDNA sequence of the Rana dybowskii LepROT gene was determined by using RT-PCR and bioinformatics analysis. Then, the Aeromonas hydrophila ( Ah ) and lipopolysaccharides (LPS) infected models of R. dybowskii was constructed to obtain histopathological characteristics. Constitutive expression of LepROT mRNA and NF-κB signaling pathway were detected by real-time quantitative PCR. The full-length cDNA of LepROT gene was 396 bp and encoded 131 amino acids. Amino acid sequence analysis revealed LepROT shares 93.74% and 86.39% identity with homologues from other amphibians and mammals respectively, and the LepROT gene was quite conserved among different species. After infection, the relative expression levels of LepROT , NF- κ B , IKKα and IKKβ mRNA were all significantly upregulated ( P < 0.01), but showed a diverse temporal pattern of up-regulation in different tissues. Therefore, it was proposed that the LepROT gene of R . dybowskii might activate the NF-κB signaling pathway to exert anti-infectious effects, thus providing evidence for further extending the biological function of LepROT.

The Bowman-Birk protease inhibitor (BBI) family is a prototype group found mainly in plants, particularly grasses and legumes, which have been subjected to decades of study. Recently, the discovery of attenuated peptides containing the canonical Bowman-Birk protease inhibitory motif has been detected in the skin secretions of amphibians, mainly from Ranidae family members. The roles of these peptides in amphibian defense have been proposed to work cooperatively with antimicrobial peptides and reduce peptide degradation. A novel trypsin inhibitory peptide, named livisin, was found in the skin secretion of the green cascade frog, Odorrana livida . The cDNA encoding the precursor of livisin was cloned, and the predicted mature peptide was characterized. The mature peptide was found to act as a potent inhibitor against several serine proteases. A comparative activity study among the native peptide and its engineered analogs was performed, and the influence of the P 1 and P 2' positions, as well as the C-terminal amidation on the structure-activity relationship for livisin, was illustrated. The findings demonstrated that livisin might serve as a potential drug discovery/development tool.

Wildlife monitoring programs are instrumental for the assessment of species, habitat status, and for the management of factors affecting them. This is particularly important for species found in freshwater ecosystems, such as amphibians, as they have higher estimated extinction rates than terrestrial species. We developed and validated two species-specific environmental DNA (eDNA) protocols and applied them in the field to detect the Hazara Torrent Frog (Allopaa hazarensis) and Murree Hills Frog (Nanorana vicina). Additionally, we compared eDNA surveys with visual encounter surveys and estimated site occupancy. eDNA surveys resulted in higher occurrence probabilities for both A. hazarensis and N. vicina than for visual encounter surveys. Detection probability using eDNA was greater for both species, particularly for A. hazarensis. The top-ranked detection model for visual encounter surveys included effects of both year and temperature on both species, and the top-ranked occupancy model included effects of elevation and year. The top-ranked detection model for eDNA data was the null model, and the top-ranked occupancy model included effects of elevation, year, and wetland type. To our knowledge, this is the first time an eDNA survey has been used to monitor amphibian species in the Himalayan region., (© 2022. The Author(s).)

We used genotyping-by-sequencing (GBS) to identify sex-linked markers in 43 wild-collected spiny frog ( Quasipaa boulengeri ) adults from a single site. We identified a total of 1049 putatively sex-linked GBS-tags, 98% of which indicated an XX/XY system, and finally confirmed 574 XY-type sex-linked loci. The sex specificity of five markers was further validated by PCR amplification using a large number of additional individuals from 26 populations of this species. A total of 27 sex linkage markers matched with the Dmrt1 gene, showing a conserved role in sex determination and differentiation in different organisms from flies and nematodes to mammals. Chromosome 1, which harbors Dmrt1 , was considered as the most likely candidate sex chromosome in anurans. Five sex-linked SNP makers indicated sex reversals, which are sparsely present in wild amphibian populations, in three out of the one-hundred and thirty-three explored individuals. The variety of sex-linked markers identified could be used in population genetics analyses requiring information on individual sex or in investigations aimed at drawing inferences about sex determination and sex chromosome evolution.