Comparative transcriptomic analysis delineates adaptation strategies of Rana kukunoris toward cold stress on the Qinghai-Tibet Plateau.

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.

Multiple Transitions between Y Chromosome and Autosome in Tago's Brown Frog Species Complex.

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.

Nuclear gene introgressions in hybrid populations of water frog Pelophylax esculentus complex: geographical analysis of the phenomenon and its interpretation.

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.

Transcriptomic analysis of skin immunity genes in the Chinese spiny frog (Quasipaa spinosa) after Proteus mirabilis infection.

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.

Unleashing diversity through flexibility: The evolutionary journey of sex chromosomes in amphibians and reptiles.

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.

Investigation of seasonal changes in lipid synthesis and metabolism-related genes in the oviduct of Chinese brown frog (Rana dybowskii).

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.

Cytogenetics of the Hybridogenetic Frog Pelophylax grafi and Its Parental Species Pelophylax perezi.

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.

Karyotypes of water frogs from the Pelophylax esculentus complex: results of cross-species chromosomal painting.

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.

Expression of cholesterol synthesis and steroidogenic markers in females of the Chinese brown frog (Rana dybowskii) during prespawning and prehibernation.

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.

Large X-effects are absent in torrent frogs with nascent sex chromosomes.

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.

Skin bacterial community differences among three species of co-occurring Ranid frogs.

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.

[Cloning of adipor1 and adipor2 genes in Rana dybowskii and its expression pattern upon infection].

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.

Intestinal response of Rana chensinensis larvae exposed to Cr and Pb, alone and in combination.

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.

Parallel Evolution of Sex-Linked Genes across XX/XY and ZZ/ZW Sex Chromosome Systems in the Frog Glandirana rugosa.

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.

Perfluorooctanoic acid and perfluorooctanesulfonic acid induce immunotoxicity through the NF-κB pathway in black-spotted frog (Rana nigromaculata).

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.

Cloning of adipor1 and adipor2 genes in Rana dybowskii and its expression pattern upon infection / 生物工程学报

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.

Musicians in the marsh: a new species of music frog (Anura: Ranidae: Nidirana) from Arunachal Pradesh, India.

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.

Comparative transcriptomes of nine tissues for the Heilongjiang brown frog (Rana amurensis).

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.

Heterogeneous Evolution of Sex Chromosomes in the Torrent Frog Genus Amolops.

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.

Taxonomic status of Rana nigromaculata mongolia and the validity of Pelophylax tenggerensis (Anura, Ranidae).

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.