SMRT and Illumina RNA sequencing reveal the complexity of terpenoid biosynthesis in Zanthoxylum armatum

Comparative transcriptome analysis and HPLC reveal candidate genes associated with synthesis of bioactive constituents in Rheum palmatum complex

Content of bioactive constituents is one of the most important characteristics in Rheum palmatum complex. Increasing ingredient content through genetic breeding is an effective strategy to solve the contradiction between large market demand and resource depletion, but currently hampered by limited understanding of metabolite biosynthesis in rhubarb. In this study, deep transcriptome sequencing was performed to compare roots, stems, and leaves of two Rheum species (PL and ZK) that show different levels of anthraquinone contents. Approximately 0.52 billion clean reads were assembled into 58,782 unigenes, of which around 80% (46,550) were found to be functionally annotated in public databases. Expression patterns of differential unigenes between PL and ZK were thoroughly investigated in different tissues. This led to the identification of various differentially expressed genes (DEGs) involved in shikimate, MEP, MVA, and polyketide pathways, as well as those involved in catechin and gallic acid biosynthesis. Some structural enzyme genes were shown to be significantly up-regulated in roots of ZK with high anthraquinone content, implying potential central roles in anthraquinone synthesis. Taken together, our study provides insights for future functional studies to unravel the mechanisms underlying metabolite biosynthesis in rhubarb.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01492-z.

Genomic-wide identification and expression analysis of AP2/ERF transcription factors in Zanthoxylum armatum reveals the candidate genes for the biosynthesis of terpenoids

Terpenoids are the main active components in the Zanthoxylum armatum leaves, which have extensive medicinal value. The Z. armatum leaf is the main by-product in the Z. armatum industry. However, the transcription factors involved in the biosynthesis of terpenoids are rarely reported. This study was performed to identify and classify the APETALA2/ethylene-responsive factor (AP2/ERF) gene family of Z. armatum. The chromosome distribution, gene structure, conserved motifs, and cis-acting elements of the promoter of the species were also comprehensively analyzed. A total of 214 ZaAP2/ERFs were identified. From the obtained transcriptome and terpenoid content data, four candidate ZaAP2/ERFs involved in the biosynthesis of terpenoids were selected via correlation and weighted gene co-expression network analysis. A phylogenetic tree was constructed using 13 AP2/ERFs related to the biosynthesis of terpenoids in other plants. ZaERF063 and ZaERF166 showed close evolutionary relationships with the ERFs in other plant species and shared a high AP2-domain sequence similarity with the two closest AP2/ERF proteins, namelySmERF8 from Salvia miltiorrhiza and AaERF4 from Artemisia annua. Further investigation into the effects of methyl jasmonate (MeJA) treatment on the content of terpenoids in Z. armatum leaves revealed that MeJA significantly induced the upregulation of ZaERF166 and led to a significant increase in the terpenoids content in Z. armatum leaves, indicating that ZaERF166 might be involved in the accumulation of terpenoids of Z. armatum. Results will be beneficial for the functional characterization of AP2/ERFs in Z. armatum and establishment of the theoretical foundation to increase the production of terpenoids via the manipulation of the regulatory elements and strengthen the development and utilization of Z. armatum leaves.

The Classification, Molecular Structure and Biological Biosynthesis of Flavonoids, and Their Roles in Biotic and Abiotic Stresses

With the climate constantly changing, plants suffer more frequently from various abiotic and biotic stresses. However, they have evolved biosynthetic machinery to survive in stressful environmental conditions. Flavonoids are involved in a variety of biological activities in plants, which can protect plants from different biotic (plant-parasitic nematodes, fungi and bacteria) and abiotic stresses (salt stress, drought stress, UV, higher and lower temperatures). Flavonoids contain several subgroups, including anthocyanidins, flavonols, flavones, flavanols, flavanones, chalcones, dihydrochalcones and dihydroflavonols, which are widely distributed in various plants. As the pathway of flavonoid biosynthesis has been well studied, many researchers have applied transgenic technologies in order to explore the molecular mechanism of genes associated with flavonoid biosynthesis; as such, many transgenic plants have shown a higher stress tolerance through the regulation of flavonoid content. In the present review, the classification, molecular structure and biological biosynthesis of flavonoids were summarized, and the roles of flavonoids under various forms of biotic and abiotic stress in plants were also included. In addition, the effect of applying genes associated with flavonoid biosynthesis on the enhancement of plant tolerance under various biotic and abiotic stresses was also discussed.

Genome-wide identification of HD-ZIP gene family and screening of genes related to prickle development in Zanthoxylum armatum

Zanthoxylum armatum is an important cash crop for medicinal and food purposes in Asia. However, its stems and leaves are covered with a large number of prickles, which cause many problems in the production process. The homeodomain leucine zipper (HD-ZIP) gene family is a class of transcription factors unique to plants that play an important role in biological processes such as morphogenesis, signal transduction, and secondary metabolite synthesis. However, little is known about HD-ZIP gene information that may be involved in prickle development of Z. armatum. Here, we identified 76 ZaHDZ genes from the Z. armatum genome and classified them into four subfamilies (I-IV) based on phylogenetic analysis, a classification further supported by gene structure and conserved motif analysis. Seventy-six ZaHDZ genes were unevenly distributed on chromosomes. Evolutionary analysis revealed that the expansion of ZaHDZ genes mainly were due to whole-genome duplication (WGD) or segmental duplication, and they experienced strong purifying selection pressure in the process of evolution. A total of 47 cis-elements were identified in the promoter region of ZaHDZ genes. Quantitative real-time polymerase chain reaction analysis was performed on subfamily IV ZaHDZ gene expression levels in five tissues and under four hormone treatments. Finally, ZaHDZ16 was predicted to be the candidate gene most likely to be involved in prickle development of Z. armatum. These results contribute to a better understanding of the characteristics of HD-ZIP gene family and lay a foundation for further study on the function of genes related to prickle development of Z. armatum.

Full-length transcriptome sequencing reveals the molecular mechanism of monoterpene and sesquiterpene biosynthesis in Cinnamomum burmannii

Essential oil of Cinnamomum burmannii is rich in monoterpenes and sesquiterpenes and is widely used in cosmetics and medicines. Knowledge about the enzymes that catalyze the formation of monoterpenes and sesquiterpenes in C. burmannii is insufficient. Therefore, anatomy observation of C. burmannii at the four developmental stages (7 days, CBS1; 14 days, CBS2; 21 days, CBS3, and 28 days, CBS4) were conducted to elucidate the origins of essential oil production. Twelve full-length transcriptomes of C. burmannii leaves at the four stages were generated using Oxford Nanopore Technologies. GC-MS analysis revealed 15 monoterpene and sesquiterpenes dramatically increased from CBS1 to CBS4. A weighted correlation network analysis (WGCNA) in association and differentially expressed genes across four developmental stages were performed. A total of 44 differentially expressed genes (DEGs) were involved in terpenoid syntheses during leaf development. Among them, the DEGs of the mevalonate acid (MVA) pathway were predominantly expressed at CBS1, while those of the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway showed increased expression from CBS2 to CBS4. Besides, fourteen genes were associated with monoterpene synthesis and nine with sesquiterpene synthesis. Functions of these DEGs were further predicted with regard to gene expression profile and phylogenetic relationship with those characterized in previous studies. In addition, 922 long noncoding RNAs (lncRNAs) were detected, of which twelve were predicted to regulate monoterpene and sesquiterpene biosynthesis. The present study provided new insights the molecular mechanisms of monoterpenoid and sesquiterpenoid syntheses of C. burmannii.

Comparative Transcriptome Analyses of Different Rheum officinale Tissues Reveal Differentially Expressed Genes Associated with Anthraquinone, Catechin, and Gallic Acid Biosynthesis

Rheum officinale Baill. is an important traditional Chinese medicinal herb, its dried roots and rhizomes being widely utilized to cure diverse diseases. However, previous studies mainly focused on the active compounds and their pharmacological effects, and the molecular mechanism underlying the biosynthesis of these ingredients in R. officinale is still elusive. Here, we performed comparative transcriptome analyses to elucidate the differentially expressed genes (DEGs) in the root, stem, and leaf of R. officinale. A total of 236,031 unigenes with N50 of 769 bp was generated, 136,329 (57.76%) of which were annotated. A total of 5884 DEGs was identified after the comparative analyses of different tissues; 175 and 126 key enzyme genes with tissue-specific expression were found in the anthraquinone, catechin/gallic acid biosynthetic pathway, respectively, and some of these key enzyme genes were verified by qRT-PCR. The phylogeny of the PKS III family in Polygonaceae indicated that probably only PL_741 PKSIII1, PL_11549 PKSIII5, and PL_101745 PKSIII6 encoded PKSIII in the polyketide pathway. These results will shed light on the molecular basis of the tissue-specific accumulation and regulation of secondary metabolites in R. officinale, and lay a foundation for the future genetic diversity, molecular assisted breeding, and germplasm resource improvement of this essential medicinal plant.

Integrated transcriptome and plant growth substance profiles to identify the regulatory factors involved in floral sex differentiation in Zanthoxylum armatum DC

Zanthoxylum armatum is a prominent plant for food industries. Its male flowers often occur in gynogenesis plants; however, the potential mechanism remains poorly understood. Herein, a total of 26 floral sex differentiation stages were observed to select four vital phases to reveal key factors by using RNA-seq, phytohormones and carbohydrates investigation. The results showed that a selective abortion of stamen or pistil primordia could result in the floral sex differentiation in Z. armatum. Carbohydrates might collaborate with cytokinin to effect the male floral differentiation, whereas female floral differentiation was involved in SA, GA₁, and ABA biosynthesis and signal transduction pathways. Meanwhile, these endogenous regulators associated with reproductive growth might be integrated into ABCDE model to regulate the floral organ differentiation in Z. armatum. Furthermore, the 21 crucial candidates were identified in co-expression network, which would contribute to uncovering their roles in floral sex differentiation of Z. armatum in further studies. To the best of our knowledge, this study was the first comprehensive investigation to link floral sex differentiation with multi-level endogenous regulatory factors in Z. armatum. It also provided new insights to explore the regulatory mechanism of floral sex differentiation, which would be benefited to cultivate high-yield varieties in Z. armatum.

Alternative Splicing and Its Roles in Plant Metabolism

Plant metabolism, including primary metabolism such as tricarboxylic acid cycle, glycolysis, shikimate and amino acid pathways as well as specialized metabolism such as biosynthesis of phenolics, alkaloids and saponins, contributes to plant survival, growth, development and interactions with the environment. To this end, these metabolic processes are tightly and finely regulated transcriptionally, post-transcriptionally, translationally and post-translationally in response to different growth and developmental stages as well as the constantly changing environment. In this review, we summarize and describe the current knowledge of the regulation of plant metabolism by alternative splicing, a post-transcriptional regulatory mechanism that generates multiple protein isoforms from a single gene by using alternative splice sites during splicing. Numerous genes in plant metabolism have been shown to be alternatively spliced under different developmental stages and stress conditions. In particular, alternative splicing serves as a regulatory mechanism to fine-tune plant metabolism by altering biochemical activities, interaction and subcellular localization of proteins encoded by splice isoforms of various genes.