Transcriptome analysis reveals the genetic basis underlying the seasonal development of keratinized nuptial spines in Leptobrachium boringii

Integrated Morphological, Comparative Transcriptomic, and Metabolomic Analyses Reveal Mechanisms Underlying Seasonal Patterns of Variation in Spines of the Giant Spiny Frog (Quasipaa spinosa)

Quasipaa spinosa, commonly known as the spiny frog, is an economically valued amphibian in China prized for its tender meat and nutritional value. This species exhibits marked sexual dimorphism, most notably the prominent spiny structures on males that are pivotal for mating success and species identification. The spines of Q. spinosa exhibit strong seasonal variation, changing significantly with the reproductive cycle, which typically spans from April to October. Sexually mature males develop densely packed, irregularly arranged round papillae with black spines on their chests during the breeding season, which may then reduce or disappear afterward, while females have smooth chest skin. Despite their ecological importance, the developmental mechanisms and biological functions of these spines have been inadequately explored. This study integrates morphological, transcriptomic, and metabolomic analyses to elucidate the mechanisms underlying the seasonal variation in spine characteristics of Q. spinosa. Our results demonstrate that spine density inversely correlates with body size and that spine development is accompanied by significant changes in epidermal thickness and keratinization during the breeding season. Comparative transcriptomic analysis across different breeding stages revealed significant gene expression alterations in pathways related to extracellular matrix interactions, tyrosine metabolism, Wnt signaling, and melanogenesis. Metabolomic analysis further identified significant seasonal shifts in metabolites essential for energy metabolism and melanin synthesis, including notable increases in citric acid and β-alanine. These molecular changes are consistent with the observed morphological adaptations, suggesting a complex regulatory mechanism supporting spine development and functionality. This study provides novel insights into the molecular basis of spine morphogenesis and its seasonal dynamics in Q. spinosa, contributing valuable information for the species' conservation and aquaculture.

Investigation of the Alternations in the Muscle Quality of Swimming Crab (Ovalipes punctatus) during Cold-Chain Transportation Using Physicochemical and TMT-Based Quantitative Proteomic Analysis

Physicochemical properties and protein alterations in Ovalipes punctatus during cold-chain transportation were examined via sensory scores, water-holding capacity (WHC), glucose (GLU) content, catalase (CAT) activity, urea nitrogen (UN) content, and tandem mass tag (TMT)-based proteomic analysis. The results revealed that sensory characteristics and texture of crab muscle deteriorated during transportation. Proteomic analysis revealed 442 and 470 different expressed proteins (DEPs) in crabs after 18 h (FC) and 36 h (DC) of transportation compared with live crabs (LC). Proteins related to muscle structure and amino acid metabolism significantly changed, as evidenced by the decreased WHC and sensory scores of crab muscle. Glycolysis, calcium signaling, and peroxisome pathways were upregulated in the FC/LC comparison, aligning with the changes in GLU content and CAT activity, revealing the stress response of energy metabolism and immune response in crabs during 0-18 h of transportation. The downregulated tricarboxylic acid (TCA) cycle and carcinogenesis-reactive oxygen species pathways were correlated with the decreasing trend in CAT activity, suggesting a gradual retardation in both energy and antioxidant metabolism in crabs during 18-36 h of transportation. Furthermore, the regulated purine nucleoside metabolic and nucleoside diphosphate-related processes, with the increasing changes in UN content, revealed the accumulation of metabolites in crabs.

Comparative transcriptomic analysis reveals the gonadal development-related gene response to environmental temperature in Mauremys mutica

The Asian yellow pond turtle (Mauremys mutica) displays temperature-dependent sex determination (TSD), in which incubation temperature during embryonic development determines the sexual fate of the individual. However, the mechanism of the sex determination/differentiation of Mauremys mutica remains a mystery. Here, we first analyzed the temperature-specific gonadal transcriptomes of Mauremys mutica prior to gonad formation and gonads during the thermosensitive period. We uncovered a list of candidates that respond to temperature stimuli enriched in several categories, such as heat shock protein family members dnajb6a, dnaja4, hspa8 and hsp90aa1, temperature sensor genes mmp17 and mmp28, and putative novel temperature-responsive genes tmco6, gria3 and eif3f. Notably, striking differences were identified in the expression profiles of genes underlying sexual development, such as tex15, insr, igf1r, cirbp, esr1, dmrt2 and Serpinh1. Moreover, we analyzed the similarity and divergence of the timecourse of gene expression among Mauremys mutica and two other reported TSD turtles (Trachemys scripta and Chrysemys picta). The shared genes revealed the common gonad-specific regulatory mechanisms existing in these three TSD turtles that initiate their sexual development. Therefore, our findings could provide basic data to elucidate the mechanisms of sex determination/differentiation of M. mutica, even contributing to further understanding of these mechanisms in other TSD turtles.

Immunoreactivities of AR, ERα, ERβ and aromatase in the nuptial pad of Chinese brown frog (<em>Rana dybowskii</em>) during pre-hibernation and the breeding period

There is a prominent local raised pad called nuptial pad on the forelimb of Chinese brown frog (Rana dybowskii), which is hypothetically concluded as an enhancement of the grip and a spreader of pheromone during the amplexus. In this study, we investigated the immunolocalization and protein expression levels of androgen receptors (AR), estrogen receptor α (ERα), ERβ and aromatase in the nuptial pad of R. dybowskii during pre-hibernation and the breeding period. Histologically, the annual development of the nuptial pad in R. dybowskii is manifested as the larger area of specialized mucous gland and the longer length of papillary epidermal projection during the breeding period. AR, ERα, ERβ and aromatase are present in the stratum granulosum, stratum spinosum, stratum basale and the secretory portion of specialized mucous glands during both periods. Western blotting results confirmed that AR, ERα and ERβ protein levels are higher during pre-hibernation than those during the breeding season. These results suggest that nuptial pad is the direct target organ of androgen and estrogen. Androgen may participate in the regulation of annual development and glandular function of nuptial pad, and estrogen may play an endocrine, autocrine or paracrine role during pre-hibernation and the breeding period.

Hepatic transcriptome of the freeze-tolerant Cope's gray treefrog, Dryophytes chrysoscelis: responses to cold acclimation and freezing

Background

Cope’s gray treefrog, Dryophytes chrysoscelis, withstands the physiological challenges of corporeal freezing, partly by accumulating cryoprotective compounds of hepatic origin, including glycerol, urea, and glucose. We hypothesized that expression of genes related to cryoprotectant mobilization and stress tolerance would be differentially regulated in response to cold. Using high-throughput RNA sequencing (RNA-Seq), a hepatic transcriptome was generated for D. chrysoscelis, and gene expression was compared among frogs that were warm-acclimated, cold-acclimated, and frozen.

Results

A total of 159,556 transcripts were generated; 39% showed homology with known transcripts, and 34% of all transcripts were annotated. Gene-level analyses identified 34,936 genes, 85% of which were annotated. Cold acclimation induced differential expression both of genes and non-coding transcripts; freezing induced few additional changes. Transcript-level analysis followed by gene-level aggregation revealed 3582 differentially expressed genes, whereas analysis at the gene level revealed 1324 differentially regulated genes. Approximately 3.6% of differentially expressed sequences were non-coding and of no identifiable homology. Expression of several genes associated with cryoprotectant accumulation was altered during cold acclimation. Of note, glycerol kinase expression decreased with cold exposure, possibly promoting accumulation of glycerol, whereas glucose export was transcriptionally promoted by upregulation of glucose-6-phosphatase and downregulation of genes of various glycolytic enzymes. Several genes related to heat shock protein response, DNA repair, and the ubiquitin proteasome pathway were upregulated in cold and frozen frogs, whereas genes involved in responses to oxidative stress and anoxia, both potential sources of cellular damage during freezing, were downregulated or unchanged.

Conclusion

Our study is the first to report transcriptomic responses to low temperature exposure in a freeze-tolerant vertebrate. The hepatic transcriptome of Dryophytes chrysoscelis is responsive to cold and freezing. Transcriptomic regulation of genes related to particular pathways, such as glycerol biosynthesis, were not all regulated in parallel. The physiological demands associated with cold and freezing, as well as the transcriptomic responses observed in this study, are shared with several organisms that face similar ecophysiological challenges, suggesting common regulatory mechanisms. The role of transcriptional regulation relative to other cellular processes, and of non-coding transcripts as elements of those responses, deserve further study.

Genomic and transcriptomic insights into molecular basis of sexually dimorphic nuptial spines in Leptobrachium leishanense

Sexually dimorphic (SD) traits are important in sexual selection and species survival, yet the molecular basis remains elusive, especially in amphibians where SD traits have evolved repeatedly. We focus on the Leishan moustache toad (Leptobrachium leishanense), in which males develop nuptial spines on their maxillary skin. Here we report a 3.5 Gb genome assembly with a contig N50 of 1.93 Mb. We find a specific expansion of the intermediate filament gene family including numerous keratin genes. Within these genes, a cluster of duplicated hair keratin genes exhibits male-biased and maxillary skin-specific expression, suggesting a role in developing nuptial spines. We identify a module of coexpressed genes significantly associated with spine formation. In addition, we find several hormones likely to be involved in regulating spine development. This study not only presents a high-quality anuran genome but also provides a reference for studying skin-derived SD traits in amphibians.

The basis of sexual dimorphism in non-model species may be elusive, in part due to a lack of genomic and molecular resources. Here, Li et al. report a high-quality anuran genome and reveal candidate genes and pathways associated with shaping sexually dimorphic nuptial spines in a moustache toad.

Multi-Tissue Transcriptomes Yield Information on High-Altitude Adaptation and Sex-Determination in Scutiger cf. sikimmensis

The Himalayas are one of earth’s hotspots of biodiversity. Among its many cryptic and undiscovered organisms, including vertebrates, this complex high-mountain ecosystem is expected to harbour many species with adaptations to life in high altitudes. However, modern evolutionary genomic studies in Himalayan vertebrates are still at the beginning. Moreover, in organisms, like most amphibians with relatively high DNA content, whole genome sequencing remains bioinformatically challenging and no complete nuclear genomes are available for Himalayan amphibians. Here, we present the first well-annotated multi-tissue transcriptome of a Greater Himalayan species, the lazy toad Scutiger cf. sikimmensis (Anura: Megophryidae). Applying Illumina NextSeq 500 RNAseq to six tissues, we obtained 41.32 Gb of sequences, assembled to ~111,000 unigenes, translating into 54362 known genes as annotated in seven functional databases. We tested 19 genes, known to play roles in anuran and reptile adaptation to high elevations, and potentially detected diversifying selection for two (TGS1, SENP5) in Scutiger. Of a list of 37 genes, we also identify 27 candidate genes for sex determination or sexual development, all of which providing the first such data for this non-model megophryid species. These transcriptomes will serve as a valuable resource for further studies on amphibian evolution in the Greater Himalaya as a biodiversity hotspot.

Chromosome-level assembly of the mustache toad genome using third-generation DNA sequencing and Hi-C analysis

BACKGROUND

The mustache toad, Vibrissaphora ailaonica, is endemic to China and belongs to the Megophryidae family. Like other mustache toad species, V. ailaonica males temporarily develop keratinized nuptial spines on their upper jaw during each breeding season, which fall off at the end of the breeding season. This feature is likely result of the reversal of sexual dimorphism in body size, with males being larger than females. A high-quality reference genome for the mustache toad would be invaluable to investigate the genetic mechanism underlying these repeatedly developing keratinized spines.

FINDINGS

To construct the mustache toad genome, we generated 225 Gb of short reads and 277 Gb of long reads using Illumina and Pacific Biosciences (PacBio) sequencing technologies, respectively. Sequencing data were assembled into a 3.53-Gb genome assembly, with a contig N50 length of 821 kb. We also used high-throughput chromosome conformation capture (Hi-C) technology to identify contacts between contigs, then assembled contigs into scaffolds and assembled a genome with 13 chromosomes and a scaffold N50 length of 412.42 Mb. Based on the 26,227 protein-coding genes annotated in the genome, we analyzed phylogenetic relationships between the mustache toad and other chordate species. The mustache toad has a relatively higher evolutionary rate and separated from a common ancestor of the marine toad, bullfrog, and Tibetan frog 206.1 million years ago. Furthermore, we identified 201 expanded gene families in the mustache toad, which were mainly enriched in immune pathway, keratin filament, and metabolic processes.

CONCLUSIONS

Using Illumina, PacBio, and Hi-C technologies, we constructed the first high-quality chromosome-level mustache toad genome. This work not only offers a valuable reference genome for functional studies of mustache toad traits but also provides important chromosomal information for wider genome comparisons.