Full-length transcriptome profiling of Gentiana straminea Maxim. provides new insights into iridoid biosynthesis pathway

TL;DR

This study uses transcriptome sequencing to explore the genetic basis of iridoid biosynthesis in the medicinal plant Gentiana straminea.

Contribution

The first full-length transcriptome analysis of G. straminea, identifying key genes and enzymes involved in iridoid biosynthesis.

Findings

A total of 117 isoforms encoding 19 key enzymes related to the iridoid synthesis pathway were identified.

Tissue-specific expression patterns of key genes like AACT, IDI, ISPH, and GCPE were highest in leaves.

Abstract

Gentiana straminea Maxim. is a traditional Chinese medicinal plant celebrated for its diverse array of bioactive compounds, particularly iridoid glycosides. These compounds are recognized as the main components that exert therapeutic effects against rheumatism, osteoarthritis, hepatitis, gastritis, and cholecystitis. Consequently, research on G. straminea has attracted significant interest, yet the genetic factors underlying the production and diversification of its secondary metabolites remain poorly understood, especially the pathway of iridoid biosynthesis. In the present study, a full-length transcriptome analysis based on PacBio Sequel II platform and Illumina sequencing was performed to identify genes that were differentially expressed in five G. straminea tissues, and proteins catalyzing iridoid biosynthesis were characterized. After sequence clustering and redundancy removal, a total of 32,776 isoforms were identified in PacBio sequencing, with an average length of 2,589.14 bp, an N50 value of 2,767 bp, and a guanine-cytosine (GC) content of 41.43%. Results of Illumina sequencing unveiled that a total of 31,330 genes were found in common in all the five tissues. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that the DEGs were mainly enriched in terms related to biosynthesis of secondary metabolites, metabolic pathways, mitogen-activated protein kinase (MAPK) signaling pathway, etc. A total of 117 isoforms encoding 19 key enzymes related to the iridoid synthesis pathway were identified, including one geranyl diphosphate synthase (GPPS) and two geranylgeranyl diphosphate synthases (GGPPS). A phylogenetic analysis further classified plant G(G)PPSs into three distinct branches. Profiling tissue-specific expression of key genes involved in iridoid synthesis revealed that the Quantitative reverse transcription polymerase chain reaction (RT-qPCR) results demonstrated the consistent trend with the FPKM (Fragments Per Kilobase of transcript per Million mapped reads) values in the root, stem, leaf, flower, ovary, non-embryonic callus (NEC), and embryonic callus (EC). Among them, AACT, IDI, ISPH, and GCPE had the highest expression levels in leaves, whereas DXS and GPPS had the highest expression levels in stems. This work provides the first transcriptomic analysis of G. straminea, serving as a foundational resource for elucidating mechanisms of bioactive compound biosynthesis, facilitating molecular studies, and supporting genomic exploration of this medicinal species.

Figures