Full-Length Transcriptome Analysis Reveals Candidate Genes Involved in Terpenoid Biosynthesis in Artemisia argyi

Integrated multi-omics analysis reveals genes involved in flavonoid biosynthesis and trichome development of Artemisia argyi

Artemisia argyi is an herbaceous plant of the genus Artemisia. Its young and mature leaves are used as food and medicine, respectively. Glandular trichomes (GTs) are distributed on the leaf surface in A. argyi and are generally considered the location of flavonoid biosynthesis and accumulation. However, the mechanism of flavonoid biosynthesis and accumulation in A. argyi remains unclear. In this study, the coregulatory genes involved in flavonoid biosynthesis and trichome development in this species were screened and evaluated, and the biosynthetic pathways for key flavonoids in A. argyi were uncovered. AaMYB1 and AaYABBY1 were screened using weighted gene co-expression network analysis, and both genes were then genetically transformed into Nicotiana tabacum L. cv. K326 (tobacco). Simultaneously, AaYABBY1 was also genetically transformed into Arabidopsis thaliana. The total flavonoid and rutin contents were increased in tobacco plants overexpressing AaMYB1 and AaYABBY1, and the expression levels of genes participating in the flavonoid synthesis pathway, such as PAL, FLS, and F3H, were significantly up-regulated in plants overexpressing these genes. These results indicated that AaMYB1 and AaYABBY1 promote flavonoid biosynthesis in tobacco. Furthermore, compared to that in the wild-type, the trichome density was significantly increased in tobacco and A. thaliana plants overexpressing AaYABBY1. These results confirm that AaYABBY1 might be involved in regulating trichome formation in A. argyi. This indicates the potential genes involved in and provides new insights into the development of trichome cellular factories based on the "development-metabolism" interaction network and the cultivation of high-quality A. argyi.

Full-Length Transcriptome Sequencing and RNA-Seq Analysis Offer Insights into Terpenoid Biosynthesis in Blumea balsamifera (L.) DC

Blumea balsamifera (L.) DC., an important economic and medicinal herb, has a long history of being used as a traditional Chinese medicine. Its leaves have always been used as a raw material for the extraction of essential oils, comprising large amounts of terpenoids, which have good therapeutic effects on many diseases, such as eczema, bacterial infection, and hypertension. However, the genetic basis of terpenoid biosynthesis in this plant is virtually unknown on account of the lack of genomic data. Here, a combination of next-generation sequencing (NGS) and full-length transcriptome sequencing was applied to identify genes involved in terpenoid biosynthesis at five developmental stages. Then, the main components of essential oils in B. balsamifera were identified using GC-MS. Overall, 16 monoterpenoids and 20 sesquiterpenoids were identified and 333,860 CCS reads were generated, yielding 65,045 non-redundant transcripts. Among these highly accurate transcripts, 59,958 (92.18%) transcripts were successfully annotated using NR, eggNOG, Swissprot, KEGG, KOG, COG, Pfam, and GO databases. Finally, a total of 56 differently expressed genes (DEGs) involved in terpenoid biosynthesis were identified, including 38 terpenoid backbone genes and 18 TPSs, which provide a significant amount of genetic information for B. balsamifera. These results build a basis for resource protection, molecular breeding, and the metabolic engineering of this plant.

Full-length transcriptome and RNA-Seq analyses reveal the resistance mechanism of sesame in response to Corynespora cassiicola

BACKGROUND

Corynespora leaf spot is a common leaf disease occurring in sesame, and the disease causes leaf yellowing and even shedding, which affects the growth quality of sesame. At present, the mechanism of sesame resistance to this disease is still unclear. Understanding the resistance mechanism of sesame to Corynespora leaf spot is highly important for the control of infection. In this study, the leaves of the sesame resistant variety (R) and the sesame susceptible variety (S) were collected at 0-48 hpi for transcriptome sequencing, and used a combined third-generation long-read and next-generation short-read technology approach to identify some key genes and main pathways related to resistance.

RESULTS

The gene expression levels of the two sesame varieties were significantly different at 0, 6, 12, 24, 36 and 48 hpi, indicating that the up-regulation of differentially expressed genes in the R might enhanced the resistance. Moreover, combined with the phenotypic observations of sesame leaves inoculated at different time points, we found that 12 hpi was the key time point leading to the resistance difference between the two sesame varieties at the molecular level. The WGCNA identified two modules significantly associated with disease resistance, and screened out 10 key genes that were highly expressed in R but low expressed in S, which belonged to transcription factors (WRKY, AP2/ERF-ERF, and NAC types) and protein kinases (RLK-Pelle_DLSV, RLK-Pelle_SD-2b, and RLK-Pelle_WAK types). These genes could be the key response factors in the response of sesame to infection by Corynespora cassiicola. GO and KEGG enrichment analysis showed that specific modules could be enriched, which manifested as enrichment in biologically important pathways, such as plant signalling hormone transduction, plant-pathogen interaction, carbon metabolism, phenylpropanoid biosynthesis, glutathione metabolism, MAPK and other stress-related pathways.

CONCLUSIONS

This study provides an important resource of genes contributing to disease resistance and will deepen our understanding of the regulation of disease resistance, paving the way for further molecular breeding of sesame.

Applications of some advanced sequencing, analytical, and computational approaches in medicinal plant research: a review

The potential active chemicals found in medicinal plants, which have long been employed as natural medicines, are abundant. Exploring the genes responsible for producing these compounds has given new insights into medicinal plant research. Previously, the authentication of medicinal plants was done via DNA marker sequencing. With the advancement of sequencing technology, several new techniques like next-generation sequencing, single molecule sequencing, and fourth-generation sequencing have emerged. These techniques enshrined the role of molecular approaches for medicinal plants because all the genes involved in the biosynthesis of medicinal compound(s) could be identified through RNA-seq analysis. In several research insights, transcriptome data have also been used for the identification of biosynthesis pathways. miRNAs in several medicinal plants and their role in the biosynthesis pathway as well as regulation of the disease-causing genes were also identified. In several research articles, an in silico study was also found to be effective in identifying the inhibitory effect of medicinal plant-based compounds against virus' gene(s). The use of advanced analytical methods like spectroscopy and chromatography in metabolite proofing of secondary metabolites has also been reported in several recent research findings. Furthermore, advancement in molecular and analytic methods will give new insight into studying the traditionally important medicinal plants that are still unexplored.

Transcriptome analysis of Artemisia argyi following methyl jasmonate (MeJA) treatment and the mining of genes related to the stress resistance pathway

Artemisia argyi Lev. et Vant. (A. argyi) is a perennial grass in the Artemisia family, the plant has a strong aroma. Methyl jasmonate (MeJA) is critical to plant growth and development, stress response, and secondary metabolic processes. The experimental material Artemisia argyi was utilized in this study to investigate the treatment of A. argyi with exogenous MeJA at concentrations of 100 and 200 μmol/L for durations of 9 and 24 h respectively. Transcriptome sequencing was conducted using the Illumina HiSeq platform to identify stress resistance-related candidate genes. Finally, a total of 102.43 Gb of data were obtained and 162,272 unigenes were identified. Differential analysis before and after MeJA treatment resulted in the screening of 20,776 differentially expressed genes. The GO classification revealed that the annotated unigenes were categorized into three distinct groups: cellular component, molecular function, and biological process. Notably, binding, metabolic process, and cellular process emerged as the most prevalent categories among them. The results of KEGG pathway statistical analysis revealed that plant hormone signal transduction, MAPK signaling pathway-plant, and plant-pathogen interaction were significant transduction pathways in A. argyi's response to exogenous MeJA-induced abiotic stress. With the alteration of exogenous MeJA concentration and duration, a significant upregulation was observed in the expression levels of calmodulin CaM4 (ID: EVM0136224) involved in MAPK signaling pathway-plant and auxin response factor ARF (ID: EVM0055178) associated with plant-pathogen interaction. The findings of this study establish a solid theoretical foundation for the future development of highly resistant varieties of A. argyi.

Novel insights on genes and pathways involved in Pinus elliottii response to resinosis

Pinus elliottii, an important coniferous timber species, has recently become one of the most popular sources of resin in China. Resinosis is a common disease that may negatively affect pine tree growth and production. In this study, we used single-molecule real-time sequencing and Illumina RNA sequencing to generate an accurate transcriptome for P. elliottii. The transcriptome included 90,026 transcripts, 5160 long non-coding RNAs and 7710 transcription factors. We then analyzed RNA-sequencing, small RNA-sequencing and degradome data to identify genes, miRNAs and key miRNA-target pairs involved in response to resinosis in P. elliottii. We identified 1305 genes and 1151 miRNAs exhibiting significant differential expression in response to resinosis. According to the degradome sequencing analysis, 318 differentially expressed transcripts were targets of 14 differentially expressed miRNAs. Our study has provided resources for further functional characterization of genes and miRNAs involved in resinosis in P. elliottii, which should aid the future disease-resistance breeding of this species.

Full-length transcriptome analysis of Zanthoxylum nitidum (Roxb.) DC

Zanthoxylum nitidum (Roxb.) DC. (Z. nitidum) is a type of Chinese Dao-di herb, also called Liangmianzhen, which is widely used to treat arthralgia, rheumatic arthralgia, and stomach pain. However, genomic resources for Z. nitidum are still scarce. This study provides transcriptomic resources for Z. nitidum by applying single-molecule real-time (SMRT) sequencing technology. In total, 456,109 circular consensus sequencing (CCS) reads were generated with a mean length of 2,216 bp from Z. nitidum roots, old stems, young branches, leaves, flowers, and fruits. Of these total reads, 353,932 were full-length nonchimeric (FLNC) reads with an average length of 1,996 bp. A total of 16,163 transcripts with a mean length of 1,171 bp were acquired. Of these transcripts, 14,231 (88%) were successfully annotated using public databases. Across all the 16,163 transcripts, we identified 6,255 long non-coding RNAs (lncRNAs) and 22,780 simple sequence repeats (SSRs). Furthermore, 3,482 transcription factors were identified. Among the SSR loci, 1-3 nucleotide repeats were dominant, occupying 99.36% of the total SSR loci, with mono-, di-, and tri-nucleotide repeats accounting for 61.80%, 19.89%, and 5.02% of the total SSR loci, respectively. A total of 36 out of 100 randomly selected primer pairs were verified to be positive, 20 of which showed polymorphism. These findings enrich the genetic resources available for facilitating future studies and research on relevant topics such as population genetics in Z. nitidum.

Full-Length Transcriptomic Sequencing and Temporal Transcriptome Expression Profiling Analyses Offer Insights into Terpenoid Biosynthesis in Artemisia argyi

Artemisiae argyi Folium is a traditional herbal medicine used for moxibustion heat therapy in China. The volatile oils in A.argyi leaves are closely related to its medicinal value. Records suggest that the levels of these terpenoids components within the leaves vary as a function of harvest time, with June being the optimal time for A. argyi harvesting, owing to the high levels of active ingredients during this month. However, the molecular mechanisms governing terpenoid biosynthesis and the time-dependent changes in this activity remain unclear. In this study, GC–MS analysis revealed that volatile oil levels varied across four different harvest months (April, May, June, and July) in A. argyi leaves, and the primarily terpenoids components (including both monoterpenes and sesquiterpenes) reached peak levels in early June. Through single-molecule real-time (SMRT) sequencing, corrected by Illumina RNA-sequencing (RNA-Seq), 44 full-length transcripts potentially involved in terpenoid biosynthesis were identified in this study. Differentially expressed genes (DEGs) exhibiting time-dependent expression patterns were divided into 12 coexpression clusters. Integrated chemical and transcriptomic analyses revealed distinct time-specific transcriptomic patterns associated with terpenoid biosynthesis. Subsequent hierarchical clustering and correlation analyses ultimately identified six transcripts that were closely linked to the production of these two types of terpenoid within A. argyi leaves, revealing that the structural diversity of terpenoid is related to the generation of the diverse terpene skeletons by prenyltransferase (TPS) family of enzymes. These findings can guide further studies of the molecular mechanisms underlying the quality of A. argyi leaves, aiding in the selection of optimal timing for harvests of A. argyi.

Sequence Characteristics and Phylogenetic Analysis of the Artemisia argyi Chloroplast Genome

Artemisia argyi Levl. et Van is an important Asteraceae species with a high medicinal value. There are abundant A. argyi germplasm resources in Asia, especially in China, but the evolutionary relationships of these varieties and the systematic localization of A. argyi in the family Asteraceae are still unclear. In this study, the chloroplast (cp) genomes of 72 A. argyi varieties were systematically analyzed. The 72 varieties originated from 47 regions in China at different longitudes, latitudes and altitudes, and included both wild and cultivated varieties. The A. argyi cp genome was found to be ∼151 kb in size and to contain 114 genes, including 82 protein-coding, 28 tRNA, and 4 rRNA genes. The number of short sequence repeats (SSRs) in A. argyi cp genomes ranged from 35 to 42, and most of them were mononucleotide A/T repeats. A total of 196 polymorphic sites were detected in the cp genomes of the 72 varieties. Phylogenetic analysis demonstrated that the genetic relationship between A. argyi varieties had a weak relationship with their geographical distribution. Furthermore, inverted repeat (IR) boundaries of 10 Artemisia species were found to be significantly different. A sequence divergence analysis of Asteraceae cp genomes showed that the variable regions were mostly located in single-copy (SC) regions and that the coding regions were more conserved than the non-coding regions. A phylogenetic tree was constructed using 43 protein-coding genes common to 67 Asteraceae species. The resulting tree was consistent with the traditional classification system; Artemisia species were clustered into one group, and A. argyi was shown to be closely related to Artemisia lactiflora and Artemisia montana. In summary, this study systematically analyzed the cp genome characteristics of A. argyi and compared cp genomes of Asteraceae species. The results provide valuable information for the definitive identification of A. argyi varieties and for the understanding of the evolutionary relationships between Asteraceae species.

Application of High-Throughput Sequencing on the Chinese Herbal Medicine for the Data-Mining of the Bioactive Compounds

The Chinese Herbal Medicine (CHM) has been used worldwide in clinic to treat the vast majority of human diseases, and the healing effect is remarkable. However, the functional components and the corresponding pharmacological mechanism of the herbs are unclear. As one of the main means, the high-throughput sequencing (HTS) technologies have been employed to discover and parse the active ingredients of CHM. Moreover, a tremendous amount of effort is made to uncover the pharmacodynamic genes associated with the synthesis of active substances. Here, based on the genome-assembly and the downstream bioinformatics analysis, we present a comprehensive summary of the application of HTS on CHM for the synthesis pathways of active ingredients from two aspects: active ingredient properties and disease classification, which are important for pharmacological, herb molecular breeding, and synthetic biology studies.