Search

Comprising a major clade of Anura, toads produce and secrete numerous toxins from both the parotoid glands behind their eyes and their dorsal skin. These toxins, made of various proteins and compounds, possess pharmacological potential to be repurposed to benefit human health. However, the detailed genetic regulation of toad toxin production is still poorly understood. A recent publication uncovering the genome of the representative Asiatic toad (Bufo gargarizans) provides a good reference to resolve this issue. In the present study, we sequenced the transcriptomes of parotoid gland, dorsal skin and liver from the Asiatic toad. Combining our data with 35 previously published transcriptomes across eight different tissues from the same species but from different locations, we constructed a comprehensive gene co-expression network of the Asiatic toad with the assistance of the reference genome assembly. We identified 2,701 co-expressed genes in the toxin-producing tissues (including parotoid gland and dorsal skin). By comparative genomic analysis, we identified 599 expanded gene families with 2,720 genes. Through overlapping these co-expressed genes in the toad toxin-producing tissues, we observed that three cytochrome P450 (Cyp) family members (Cyp27a1, Cyp2c29, and Cyp2c39) were significantly enriched in pathways related to cholesterol metabolism. Cholesterol is a critical precursor to steroids, and the known main steroidal toxins of bufadienolides are considered as the major bioactive components in the parotoid glands of Asiatic toad. We found 3-hydroxy-methylglutaryl CoA reductase (hmgcr), encoding the major rate-limiting enzyme for cholesterol biosynthesis, appears with multiple copies in both Asiatic toad and common toad, possibly originating from a tandem duplication event. The five copies of hmgcr genes consistently displayed higher transcription levels in the parotoid gland when compared with the abdominal skin, suggesting it as a vital candidate gene in the involvement of toad toxin production. Taken together, our current study uncovers transcriptomic and gene-family dynamic evidence to reveal the vital role of both expanded gene copies and gene expression changes for production of toad toxins.

Midas cichlid (Amphilophus citrinellus), a popular aquarium fish, attracts extensive attention from worldwide biologists mainly due to its morphological polymorphism (dark versus gold). Continuous efforts have therefore been paid to address mechanisms of its coloration variants, while it is far away from the detailed illustration of a clear regulatory network. Some limits may come from the absence of a high-quality genome assembly and a relatively accurate gene set. In this study, we sequenced about 149 Gb of nucleotide sequences of Midas cichlid, generating a genome assembly with a total size of 933.5 Mb, which exhibits a good genome continuity with a contig N50 of 10.5 Mb. A total of 25,911 protein-coding genes were annotated and about 90% completeness was achieved, which helps to build a good gene pool for understanding expressional differences of color variation. With the assistance of the final gene set, we identified a total of 277 differential expressional genes (DEGs), of which 97 up- and 180 downregulated were determined in dark-vs-gold comparisons. Two protein-protein interaction (PPI) networks were constructed from these DEGs, and three key functional modules were classified. Hub genes within each module were evaluated, and we found that the third key module contains tyrp1b, oca2, pmela, tyr, and slc24a5, which were previously proven to be associated with melanin formation. Two downregulated DEGs (myl1 and pgam2) in the first key module may be involved in muscle movement and spermatogenesis, implying that certain side effects could result from the morphological polymorphism. The first key module, consisting of proteins encoded by upregulated DEGs that were associated with MAPK signaling, Toll-like receptor signaling, and gonadotropin-releasing hormone pathways, may contribute to a negative upstream regulation or downstream influence on melanin biosynthesis. Taken together, our new genome assembly and gene annotation of Midas cichlid provide a high-quality genetic resource for biological studies on this species, and the newly identified key networks and hub genes in dark-vs-gold comparisons enhance our understanding of the transcriptional regulatory mechanisms underlying coloration changes not only in Midas cichlid but also in other fishes from freshwater to marine ecosystems.

As an important economic and medicinal crop, Amomum tsao-ko is rich in volatile oils and widely used in food additives, essential oils, and traditional Chinese medicine. However, the lack of the genome remains a limiting factor for understanding its medicinal properties at the molecular level. Here, based on 288.72 Gb of PacBio long reads and 105.45 Gb of Illumina paired-end short reads, we assembled a draft genome for A. tsao-ko (2.70 Gb in size, contig N50 of 2.45 Mb). Approximately 90.07% of the predicted genes were annotated in public databases. Based on comparative genomic analysis, genes involved in secondary metabolite biosynthesis, flavonoid metabolism, and terpenoid biosynthesis showed significant expansion. Notably, the DXS, GGPPS, and CYP450 genes, which participate in rate-limiting steps for terpenoid backbone biosynthesis and modification, may form the genetic basis for essential oil formation in A. tsao-ko. The assembled A. tsao-ko draft genome provides a valuable genetic resource for understanding the unique features of this plant and for further evolutionary and agronomic studies of Zingiberaceae species.

Serotonin is important in vertebrates for its crucial roles in regulation of various physiological functions. Investigations on how the biosynthesizing enzymes mediate serotonin production and conversion during biological processes have been active in the past decades. However, a clear-cut picture of these enzymes in molecular evolution is very limited, particularly when the complexity is imaginable in fishes since teleosts had experienced additional whole genome duplication (WGD) event(s) than tetrapods. Since serotonin is the main intermediate product during melatonin biosynthesis from tryptophan, we therefore summarize an overview of recent discoveries about molecular evolution of the four melatonin biosynthesizing enzymes, especially the L-aromatic amino acid decarboxylase (AAAD) for serotonin production and aralkylamine N-acetyltransferase (AANAT) for serotonin conversion in vertebrates. Novel copies of these genes, possibly due to WGD, were discovered in fishes. Detailed sequence comparisons revealed various variant sites in these newly identified genes, suggesting functional changes from the conventional recognition of these enzymes. These interesting advances will benefit readers to obtain new insights into related genomic differences between mammals and fishes, with an emphasis on the potential specificity for AANAT in naturally cave-restricted and deep-sea fishes.

A Yunnan-Guizhou Plateau fish, the Kanglang white minnow (Anabarilius grahami), is a typical “3E” (Endangered, Endemic, and Economic) species in China. Its distribution is limited to Fuxian Lake, the nation’s second deepest lake, with a significant local economic value but a drastically declining wild population. This species has been evaluated as VU (Vulnerable) in the China Species Red List. As one of the “Four Famous Fish” in Yunnan province, the artificial breeding has been achieved since 2003. It has not only re-established its wild natural populations by reintroduction of the artificial breeding stocks, but also brought a wide and popular utilization of this species to the local fish farms. Anabarilius grahami has become one of the main native aquaculture species in Yunnan province, and the artificial production has been emerging in steady growth each year. To promote the conservation and sustainable utilization of this fish, we initiated its whole genome sequencing project using an Illumina Hiseq2500 platform. The assembled genome size of A. grahami is 1.006 Gb, accounting for 98.63% of the estimated genome size (1.020 Gb), with contig N50 and scaffold N50 values of 26.4 kb and 4.41 Mb, respectively. Approximately about 50.38% of the genome was repetitive. A total of 25,520 protein-coding genes were subsequently predicted. A phylogenetic tree based on 4,580 single-copy genes from A. grahami and 18 other cyprinids revealed three well-supported subclades within the Cyprinidae. This is the first inter-subfamily relationship of cyprinids at genome level, providing a simple yet useful framework for understanding the traditional but popular subfamily classification systems. Interestingly, a further population demography of A. grahami uncovered a historical relationship between this fish and Fuxian Lake, suggesting that range expansion or shrinkage of the habitat has had a remarkable impact on the population size of endemic plateau fishes. Additionally, a total of 33,836 simple sequence repeats (SSR) markers were identified, and 11 loci were evaluated for a preliminary genetic diversity analysis in this study, thus providing another useful genetic resource for studying this “3E” species.