Integrated Analysis of mRNA and microRNA Elucidates the Regulation of Glycyrrhizic Acid Biosynthesis in Glycyrrhiza uralensis Fisch

Probing of plant transcriptomes reveals the hidden genetic diversity of the family Secoviridae

Secoviruses are single-stranded RNA viruses that infect plants. In the present study, we identified 61 putative novel secoviral genomes in various plant species by mining publicly available plant transcriptome data. These viral sequences represent the genomes of 13 monopartite and 48 bipartite secovirids. The genome sequences of 52 secovirids were coding-complete, and nine were partial. Except for small open reading frames (ORFs) determined in waikaviral genomes and RNA2 of torradoviruses, all of the recovered genomes/genome segments contained a large ORF encoding a polyprotein. Based on genome organization and phylogeny, all but three of the novel secoviruses were assigned to different genera. The genome organization of two identified waika-like viruses resembled that of the recently identified waika-like virus Triticum aestivum secovirus. Phylogenetic analysis revealed a pattern of host-virus co-evolution in a few waika- and waika-like viruses and increased phylogenetic diversity of nepoviruses. The study provides a basis for further investigation of the biological properties of these novel secoviruses.

MicroRNAs in Medicinal Plants

Medicinal plant microRNAs (miRNAs) are an endogenous class of small RNA central to the posttranscriptional regulation of gene expression. Biosynthetic research has shown that the mature miRNAs in medicinal plants can be produced from either the standard messenger RNA splicing mechanism or the pre-ribosomal RNA splicing process. The medicinal plant miRNA function is separated into two levels: (1) the cross-kingdom level, which is the regulation of disease-related genes in animal cells by oral intake, and (2) the intra-kingdom level, which is the participation of metabolism, development, and stress adaptation in homologous or heterologous plants. Increasing research continues to enrich the biosynthesis and function of medicinal plant miRNAs. In this review, peer-reviewed papers on medicinal plant miRNAs published on the Web of Science were discussed, covering a total of 78 species. The feasibility of the emerging role of medicinal plant miRNAs in regulating animal gene function was critically evaluated. Staged progress in intra-kingdom miRNA research has only been found in a few medicinal plants, which may be mainly inhibited by their long growth cycle, high demand for growth environment, immature genetic transformation, and difficult RNA extraction. The present review clarifies the research significance, opportunities, and challenges of medicinal plant miRNAs in drug development and agricultural production. The discussion of the latest results furthers the understanding of medicinal plant miRNAs and helps the rational design of the corresponding miRNA/target genes functional modules.

Comparative transcriptome analysis and identification of candidate genes involved in cucurbitacin IIa biosynthesis in Hemsleya macrosperma

Hemsleya macrosperma (H. macrosperma) is widely used in southwestern China as folk medicine with various bioactivities. Cucurbitacin IIa is the main active component in H. macrosperma and draws increased attention for its potential pharmacological activities. In order to reveal the mechanism of cucurbitacin IIa biosynthesis and regulation in H. macrosperma, transcriptome analysis was performed to compare differentially expressed genes in three tissues (root tuber, stem and leaf). A total of 47 946 unigenes were generated from these tissues and 55 unigenes were identified as candidate genes involved in triterpenoid backbone biosynthesis. Three homologous genes encoding squalene epoxidase (HmSE) were discovered and successfully expressed in a prokaryotic system. HmSE1 was found to be responsible for oxidization of squalene. In addition, several cytochrome P450s and transcription factors were predicted as candidates associated to cucurbitacin IIa biosynthesis. Notably, the expression profiles of those putative genes showed a positive correlation with elevated curcurbitacin IIa production in methyl jasmonate-elicited suspension cells of H. macrosperma., suggesting probable functions of the candidates on curcurbitacin IIa biosynthesis. These findings provide insights on cucurbitacin IIa biosynthesis and regulation in H. macrosperma.

GURFAP: A Platform for Gene Function Analysis in Glycyrrhiza Uralensis

Glycyrrhiza uralensis (Licorice), which belongs to Leguminosae, is famous for the function of pharmacologic action and natural sweetener with its dried roots and rhizomes. In recent years, the whole-genome sequence of G. uralensis has been completed, which will help to lay the foundation for the study of gene function. Here, we integrated the available genomic and transcriptomic data of G. uralensis and constructed the G. uralensis gene co-expression network. We then annotated gene functions of G. uralensis via aligning with public databases. Furthermore, gene families of G. uralensis were predicted by tools including iTAK (Plant Transcription factor and Protein kinase Identifier and Classifier), HMMER (hidden Markov models), InParanoid, and PfamScan. Finally, we constructed a platform for gene function analysis in G. uralensis (GURFAP, www.gzybioinfoormatics.cn/GURFAP). For analyzed and predicted gene function, we introduced various tools including BLAST (Basic local alignment search tool), GSEA (Gene set enrichment analysis), Motif, Heatmap, and JBrowse. Our analysis based on this platform indicated that the biosynthesis of glycyrrhizin might be regulated by MYB and bHLH. We also took CYP88D6, CYP72A154, and bAS gene in the synthesis pathway of glycyrrhizin as examples to demonstrate the reliability and availability of our platform. Our platform GURFAP will provide convenience for researchers to mine the gene function of G. uralensis and thus discover more key genes involved in the biosynthetic pathway of active ingredients.

Genome-Wide Analysis of Alternative Splicing and Non-Coding RNAs Reveal Complicated Transcriptional Regulation in Cannabis sativa L

It is of significance to mine the structural genes related to the biosynthetic pathway of fatty acid (FA) and cellulose as well as explore the regulatory mechanism of alternative splicing (AS), microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) in the biosynthesis of cannabinoids, FA and cellulose, which would enhance the knowledge of gene expression and regulation at post-transcriptional level in Cannabis sativa L. In this study, transcriptome, small RNA and degradome libraries of hemp 'Yunma No.1' were established, and comprehensive analysis was performed. As a result, a total of 154, 32 and 331 transcripts encoding key enzymes involved in the biosynthesis of cannabinoids, FA and cellulose were predicted, respectively, among which AS occurred in 368 transcripts. Moreover, 183 conserved miRNAs, 380 C. sativa-specific miRNAs and 7783 lncRNAs were predicted. Among them, 70 miRNAs and 17 lncRNAs potentially targeted 13 and 17 transcripts, respectively, encoding key enzymes or transporters involved in the biosynthesis of cannabinoids, cellulose or FA. Finally, the crosstalk between AS and miRNAs or lncRNAs involved in cannabinoids and cellulose was also predicted. In summary, all these results provided insights into the complicated network of gene expression and regulation in C. sativa.

Genome-Wide Characterization and Expression Analysis of bZIP Gene Family Under Abiotic Stress in Glycyrrhiza uralensis

bZIP gene family is one of the largest transcription factor families. It plays an important role in plant growth, metabolic, and environmental response. However, complete genome-wide investigation of bZIP gene family in Glycyrrhiza uralensis remains unexplained. In this study, 66 putative bZIP genes in the genome of G. uralensis were identified. And their evolutionary classification, physicochemical properties, conserved domain, functional differentiation, and the expression level under different stress conditions were further analyzed. All the members were clustered into 13 subfamilies (A–K, M, and S). A total of 10 conserved motifs were found in GubZIP proteins. Members from the same subfamily shared highly similar gene structures and conserved domains. Tandem duplication events acted as a major driving force for the evolution of bZIP gene family in G. uralensis. Cis-acting elements and protein–protein interaction networks showed that GubZIPs in one subfamily are involved in multiple functions, while some GubZIPs from different subfamilies may share the same functional category. The miRNA network targeting GubZIPs showed that the regulation at the transcriptional level may affect protein–protein interaction networks. We suspected that domain-mediated interactions may categorize a protein family into subfamilies in G. uralensis. Furthermore, the tissue-specific gene expression patterns of GubZIPs were analyzed using the public RNA-seq data. Moreover, gene expression level of 66 bZIP family members under abiotic stress treatments was quantified by using qRT-PCR. The results of this study may serve as potential candidates for functional characterization in the future.

Current advances in environmental stimuli regulating the glycyrrhizic acid biosynthesis pathway

Glycyrrhizic acid, the main active ingredient of licorice, has good antibacterial, anti-tumor, anti-viral, anti-inflammatory, and immunostimulatory activities. However, the content of glycyrrhizic acid fluctuates greatly in different licorice cultivars, and production depends on plant sources, which greatly limits its development and applications. Therefore, increasing glycyrrhizic acid content has become a research priority. In recent years, regulation of the glycyrrhizic acid biosynthesis pathway has been analyzed, the downstream synthesis pathway in licorice has been fully investigated, some key genes have been cloned, polymorphisms have been studied, and the content of glycyrrhizic acid was shown to be regulated by environmental stimuli. This work has provided a basis for studying the regulation mechanism of the glycyrrhizic acid synthesis pathway. This review summarizes and discusses relevant research to provide a current understanding of the glycyrrhizic acid synthesis pathway and its regulation in licorice.

Regulatory protein genes and microRNAs in response to selenium stimuli in Pueraria lobata (Willd.) Ohwi

Regulatory protein genes and microRNAs (miRNAs) play important roles in response to abiotic and biotic stress, and the biosynthesis of secondary metabolites in plants. However, their responses to selenium (Se) stimuli have not been comprehensively studied in Pueraria lobata (Willd.) Ohwi, a selenocompound-rich medicinal and edible plant. In this study, we identified a total of 436/556/1161/624 transcription factors, 134/157/308/172 transcriptional regulators, and 341/456/250/518 protein kinases, which were co-expressed with at least one selenocompound-related structural gene/sulfate transporter or phosphate transporter/reactive oxygen species (ROS) scavenging structural gene/isoflavone-related structural gene, respectively. Then, we identified a total of 87 expressed miRNAs by Se disposure, in which 11 miRNAs, including miR171f-3p, miR390b-3P, miR-N111b, miR-N118, miR-N30, miR-N38-3P, miR-N61a, miR-N61b, miR-N80-3p, miR-N84-3P, and miR-N90.2-3P, were significantly upregulated. We also identified a total of 1172 target genes for the 87 expressed miRNAs. Gene Ontology enrichment analysis of these target genes showed that regulation of transcription, DNA-templated, integral component of membrane, nucleus, ATP binding, and plasma membrane are the top five subclassifications. Finally, we revealed that 5 miRNAs targeted 10 regulatory protein genes, which are highly correlated with at least one selenocompound-related structural gene or transporter gene; 5 miRNAs targeted 10 regulatory protein genes, which are highly correlated with at least one ROS scavenging structural gene; and 5 miRNAs targeted 9 regulatory protein genes, which are potentially involved in the isoflavone biosynthesis. Overall, the study provides us the comprehensive insight into the roles of regulatory proteins and miRNAs in response to Se stimuli in P. lobata.