Abstract
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.
Introduction
Glycyrrhiza uralensis (Leguminosae family), an herb medicine, is distributed all over the world, especially in the northwest and northeast of China (). It is wildly used in traditional Chinese medicine with its dried roots and rhizomes for thousands of years (). Ancient records in Shennong Materia Medica demonstrated that G. uralensis was good for relieving cough and reducing sputum, regulating painful menstruation and toxin resolving (). For example, compound licorice tablets (CPLTs) contain various ingredients including G. uralensis and have antitussive and expectorant effects (). Modern pharmacological studies show that the active extract of G. uralensis has the function of anti-oxidant, anti-inflammatory, anti-tumor, anti-viral, liver protection, memory-enhancement, anti-aging, and so on (; ). Due to its antioxidant and anti-aging activities, G. uralensis is also widely used in cosmetics (; ). Some brands of cosmetics such as LANCOME, CHANEL, CLINIQUE, and Shiseido contain extracts of G. uralensis. In addition, G. uralensis is also wildly used in food processing as a natural non-caloric sweetener (). With the research development, many chemical components of G. uralensis are gradually discovered, such as flavonoids, triterpenoid saponins, coumarin, and polysaccharides (), which give it multiple pharmaceutical functions. Many drugs approved by China Food and Drug Administration (CFDA) contain the extract of G. uralensis, such as Licorice extract powder, Glycyrrhetinic acid and MgIG (). Glycyrrhizin is one of the important components of G. uralensis to exert pharmacological effects and is used for the treatment of chronic hepatitis in Asian countries (; ).
All the reasons above demonstrate that G. uralensis has good development opportunities and application prospects. However, increasing demand is easy to result in over-exploitation. Although the cultivated G. uralensis is wildly distributing in northwest China, the secondary metabolites had a lower level than that in wild samples (). Therefore, exploring the molecular mechanism of secondary metabolites biosynthesis will improve the accumulation of secondary metabolites in cultivated G. uralensis.
In recent years, with the wide application of high-throughput sequencing technology in the field of life science, several medicinal plant genomes have been decoded, such as Lonicera japonica (), Gastrodia elata (), Catharanthus roseus (), Eriobotrya japonica Lindl (), and Carthamus tinctorius (), which provides an important guarantee for protection and development medicinal plants genomes.To better share genome information, genome functional databases of many medical plants are constantly published (; ), which will be of great importance for gene functional research. She et al () collected genome sequence and 53 RNA-seq datasets of Catharanthus roseus from public databases and constructed the online croFGD database, then added miRNA-target pairs and several tools to predicted gene function. Yang et al () integrated genome, transcriptome and other relevant data of Gastrodia elata to annotate the gene function and identified gene families, and then constructed the online database GelFAP. At the same time, a variety of gene function analysis tools were introduced to GelFAP, which facilitated the research on gene function of G. elata. Xiao et al () established the functional genome database (LjaFGD) of Lonicera japonica. In this platform, the researchers collected 77 transcriptome data to construct the co-expression network of L. japonica, and then they found that MYB and WRKY transcription factor family co-expressed with key enzyme genes in the biosynthesis of chlorogenic acid and luteolin. This database will provide ideas for studying the biosynthesis mechanism of chlorogenic acid and luteolin in L. japonica. In addition to the research about medicine plant databases, other plant databases had also been reported, such as Coriander Genomics Database (CGDB) (), Kiwifruit Genome Database (KGD) (), Portal of Juglandaceae (PJU) (), Malvaceae plants (MaGenDB) ().
The whole-genome sequence of G. uralensis has been published (), which will provide more information for gene functional research. However, no comprehensive database of G. uralensis for gene functional mining and analysis is available for researchers. Therefore, combined with the genome and 35 transcriptome samples of G. uralensis from public platforms, we established an online gene functional platform, named GURFAP. In this platform, we constructed gene co-expression network and annotated gene function of G. uralensis by aligning with public platform. We also predicted the protein-protein interaction (PPI) network of G. uralensis based on the PPI network of A. thaliana. In addition, we introduced various analysis tools including basic local alignment search tool (BLAST), gene sets enrichment analysis (GSEA), motif analysis, and heatmap for researchers to explore gene functions. Based on GURFAP, we found that the biosynthesis of glycyrrhizin might be regulated by transcription factors such as MYB and bHLH, we also cited examples to demonstrate the reliability and availability of our platform. We hope this platform will provide convenience for the research about G. uralensis gene function.
Materials and Methods
Data Resource
Genome data was derived from G. uralensis genome database (http://ngs-data-archive.psc.riken.jp/Gur-genome/index.pl), including genome sequence, gene location file (gff3 file), transcript sequence, protein sequence. 35 transcriptome data samples of G. uralensis are from National Center for Biotechnology Information (NCBI, https://www.ncbi.nlm.nih.gov/) Sequence Read Archive (SRA, https://www.ncbi.nlm.nih.gov/sra/) and the detailed information is listed in Supplementary Table S1. Public protein sequences used for gene functional annotation are derived from non-redundant protein sequence database (NR, https://ftp.ncbi.nlm.nih.gov/blast/db/FASTA), cluster of orthologous groups of proteins (COG, https://www.ncbi.nlm.nih.gov/COG), the Arabidopsis information resource (TAIR, https://www.arabidopsis.org/) (), swissprot, and translated EMBL nucleotide sequence data library (TrEMBL) database (https://www.uniprot.org/downloads). Gene ontology (GO) functional annotation is based on AgriGOv2 (). Kyoto Encyclopedia of Genes and Genomes (KEGG) (https://www.kegg.jp/) annotation is based on GhostKOALA (https://www.kegg.jp/ghostkoala/) (). Pfam domain annotation information comes from the Pfam database (http://pfam.xfam.org/) ().
Gene Family Collection and Prediction
The Ethylene-responsive element binding factor-associated Amphiphilic Repression (EAR) motif-containing proteins were collected from PlantEAR (http://structuralbiology.cau.edu.cn/plantEAR) (). Transcription factors (TF), transcription regulators (TR), and protein kinase (PK) family members were predicted by iTAK software (http://itak.feilab.net/cgi-bin/itak/index.cgi) (). Ubiquitin proteases were predicted by HMMER software () based on the hidden Markov model files of conserved domains of ubiquitin proteases downloaded from ubiquitin and ubiquitin-like conjugation database (UUCD, http://uucd.biocuckoo.org/) (). For cytochrome P450 (CYP450), we firstly collected all the CYP450 protein from public databases (http://drnelson.utmem.edu/CytochromeP450.html) (), and then predicted the homology between CYP450 protein and G. uralensis protein by InParanoid (bootstrap>60%) (). Finally, we obtained the predicted CYP450 proteins of G. uralensis. To identify the carbohydrate-active enzymes (CAZy), we also collected the CAZy gene of A.thinana from the CAZy database (http://www.cazy.org/) (), and then matched the CAZy gene family to G. uralensis based on orthologous relationship.
Transcriptome Data Processing
RNA-seq transcriptomic data samples of G. uralensis were obtained from NCBI SRA database. We used Hisat2 software () to map clean reads to the reference genome and then calculated the fragments per kilobase of exon model per million reads mapped (FPKM) values by the Stringtie software (). In addition, We used the RseQC package () and WigToBigwig () software to generate wiggle files for further display and analysis.
Co-Expression network Construction
The gene expression matrix was constructed based on the FPKM expression values of each gene in each sample, and then the expression value of each gene was normalized by z-score method to calculate the Pearson Correlation Coefficient (PCC) between every two genes. Different PCC thresholds were set to construct co-expression networks. By evaluating the change of scale-free fit index and network density with the PCC value, the most appropriate PCC threshold was selected to construct the co-expression network, so that the constructed co-expression network was the most suitable for the scale-free network distribution with a relatively lower density.
Protein-Protein Interaction network Construction
PPIs of Arabidopsis we previously collected from the BioGRID (http://thebiogrid.org/) (), TAIR (), and BAR (http://bar.utoronto.ca/welcome.htm) () were used to predict the PPI network of G. uralensis. The orthologous proteins between Arabidopsis and G. uralensis were predicted by InPranoid software (bootstrap>0.6), and then the PPI network of Arabidopsis was mapped to G. uralensis to construct its protein interaction network.
Platform Construction and Visualization
We constructed the gene functional platform of G. uralensis based on classical LAMP (Linux, Apache, MySQL and PHP) architecture. Network was display by cytoscape.js (http://js.cytoscape.org/) () and heatmap was display by highcharts javascript (https://www.highcharts.com/).
Gene Set Enrichment Analysis Tool
Gene set enrichment analysis (GSEA) is an analysis method for gene functional annotation of gene sets. PlantGSEA provides gene set enrichment analysis tool for various plants. Here, we used the same method to construct GSEA of G. uralensis and the background gene sets including GO annotation, KEGG annotation, and gene family.
Motif Enrichment Analysis Tool
Cis-acting elements (Motif) enrichment analysis tool we used in the platform was as our previous research (; ). Motifs significantly enriched (p-value<0.05) were defined by calculating p-value and z-score as below formula:
‾X) Number of matches to a motif in the promoter (3-kb) of candidate genes; μ) Average number of matches to the same motif in the promoter of 1000 random lists of genes; σ) Standard deviation of the motif from 1000 random selection numbers; n) Gene number. Pnorm) Distributed function of the normal distribution.
Other Tools
In the platform, Online JBrowse () was used to display G. uralensis genes and transcriptome data. BLAST was used for sequence alignment. Heatmap analysis tool was used to display and compare the expression of genes in various transcriptome samples.
Results
Functional Annotation
We annotated 28454, 19286, 9940, 28481, 25903 genes respectively by comparison with protein bank of gene annotation database NR, Swissprot, COG, trEMBL, and TAIR () (Figure 1A). KEGG annotations of 3531 genes were predicted by KEGG database (Figure 1A). 21196 gene GO term and annotation of G. uralensis were obtained from AgriGOv2 (). Protein domains of 21960 genes were predicted by localized PfamScan software (Figure 1A).
FIGURE 1
Gene Family Prediction and Collection
In the research, 1823 transcription factors, 359 transcriptional regulators, and 1125 protein kinases were predicted by iTAK software. UUCD database () provided the HMM profile of conserved domains of ubiquitin protease family, we used it to predict 1091 ubiquitin proteases of G. uralensis. Orthologous relationships can be used to predict gene families in different species, so we collected CYP450 and CAZy sequences from the cytochrome P450 homepage () and CAZy database () respectively, and predicted 259 CYP450 and 1147 CAZy genes of G. uralensis via InParanoid softward (). In addition, 299 EAR motif-containing proteins were collected from PlantEAR (Figure 1B).
Analysis Tools
We firstly introduced BLAST analysis tool in GURFAP. BLAST analysis can help users to analyze homology between specific sequences and transcripts or proteins of G. uralensis. We also provided gene set enrichment analysis tools (GSEA) () and motif enrichment analysis tools. For GSEA analysis, we provided GO, KEGG, and gene family as background gene sets. For motif analysis, we took the motif we collected as the analysis background for users to conduct motif enrichment analysis. In addition, we also built a genome browser for G. uralensis based on JBrowse () to demonstrate gene structure and transcriptome expression level. Heatmap analysis was added to visually show the expression values of candidate gene lists in different samples.
Co-Expression Network
We first mapped reads from each RNA-seq sample to the genome of G. uralensis and screened the transcriptome with an overall mapping ratio more than 60% (Supplementary Table S1) and then calculated the expression level of each gene in all samples. The correlation between genes expression was calculated by the PCC algorithm to construct the gene co-expression network of G. uralensis. According to the PCC distribution result, we found that the correlation between most genes was not high, mainly concentrated in the middle part (Figure 1C), and with the increase or decrease of correlation, the number of gene pairs gradually decreased (Figure 1C). Biological networks are usually scale-free with low network density (; ; ). For the positive co-expression network, we evaluated and set the threshold of PCC in 0.6–0.95 to calculate the scale-free model fit index (R2) and network density respectively. We found that the scale-free fit index was the highest with PCC>0.85 and was the lowest with PCC < −0.65 (Figures 1D,E). Therefore, we chose PCC>0.85 and PCC < −0.65 to construct the positive and negative co-expression networks of G. uralensis (Figures 1D,E). Finally, we obtained the positive co-expression network including 28289 nodes and 1823961 edges, and the negative co-expression network including 15984 nodes and 215,334 edges (Figure 1F).
PPI network
We predicted the orthologous relationship of Arabidopsis and G. uralensis proteins using the Inparanoid software (bootstroop>0.6) (), and then map the PPI network of Arabidopsis to G. uralensis according to their orthologous relationship. Finally, we obtained PPI network of G. uralensis with 6562 nodes and 34652 edges (Figure 1F).
Expression View of Co-expression network
To make full use of the collected transcriptome data, we conducted a comparative analysis of transcriptome data and obtained differential expression genes (Supplementary Figure S1). SRP053019 is a set of data containing four RNA-seq samples, moderate drought stress, and control (repeat twice). Comparative analysis showed that 2803 genes were significantly up-regulated and 2595 genes were significantly down-regulated under moderate drought stress (Supplementary Figure S1A). We also conducted comparative analysis of SRP065514 (Supplementary Figure S1B), DRP000996 (Supplementary Figure S1C), SRP215420 (Supplementary Figure S1D) and SRP188776 (Supplementary Figure S1E). Combining with differential expression genes and network display, we constructed the expression view of co-expression network. Based on this expression view, users can obtain up-regulation and down-regulation genes in gene networks under specific conditions.
Platform Framework
Integrating relevant annotation information and analysis tools, we constructed a platform for gene functional analysis in G. uralensis. In the platform, seven sections are displayed with different functions (Figure 2). The home section mainly contains the introduction of G. uralensis and the platform. The network section contains the search of co-expression network and PPI network. The tools section contains search, GSEA, BLAST, Motif, Heatmap Analysis, and JBrowse. The pathway section contains genes of pathway annotated by KEGG. Gene family section contains information about transcription factors (TF), transcription regulators (TR), protein kinase (PK), EAR motif-containing proteins (EAR), cytochrome P450 (CYP450), Carbohydrate-Active enzymes (CAZy), and Ubiquitin proteases families. Download and Help section are individual part that provides users with download information and operation information.
FIGURE 2
Function Application
Co-Expression Network Analysis of Key Enzyme Genes in Glycyrrhizin Biosynthesis
Glycyrrhizin is one of the main pharmacological components in G. uralensis and its biosynthetic pathways have been relatively clear (). Based on the KEGG signaling pathway annotation, we constructed the glycyrrhizin synthesis pathway, including the key enzymes isopentenyl diphosphate isomerase (IDI), farnesyl pyrophosphate synthase (FPS), squalene synthase (SQS), squalene monooxygenase or epoxidase (SQE), bamyrin synthetase (bAS), cytochrome P450 monooxygenases (CYP450), and glucuronosyltransferase (UGAT) (Figure 3A, Supplementary Table S2). Here we obtained 1 IDI, 1 FPS, 5 SQE, 3 SQS, and 1 bAS based on KEGG annotation. In addition, Seki et al found that two CYP450 were involved in glycyrrhizic acid synthesis (; ). Xu et al found that one UGTA was involved in the biosynthesis of glycyrrhizin (). According to the co-expression network analysis of these key enzymes, we found that these key enzymes are co-expressed with many transcription factors, such as MYB, WKRY, bHLH, and ERF (Figure 3B). Therefore, these transcription factors may play an important role in the synthesis of glycyrrhizin in G. uralensis. In addition, we found co-expression relationships among these key enzymes, which were divided into four different subgroups (Figure 3C). The first subgroup is the co-expression relationship between 1 IDI and 1 SQE (green, #389339), the second subgroup represents the co-expression relationship between 2 SQE (orange, #f3951b), the third subgroup represents 1FPS, 2 SQS and 1bAS are co-expressed (purple, #674b9c), the fourth subgroup represent the co-expressed relationship in 3 CYP450 genes (blue, #204fa2) (Figure 3C). To explore the upstream regulatory factors of these key enzyme genes with co-expression relationships, we separately analyzed the four groups genes by motif enrichment analysis tool as described in the materials and methods section, and found that many transcription factors were enriched in their promoter regions. MYB, NAC, HSE, and AGL15 are enriched in the first groups (Figure 3D), MYB, NAC, and HSE are also enriched in second group (Figure 3E), bZIP, MYB, CCA1 are enriched in third group (Figure 3F) and MYB, HB etc., are existed in the last group (Figure 3G). This suggests that these transcription factors may play a regulatory role in the synthesis of glycyrrhizin.
FIGURE 3
Functional Analysis of Key Enzyme Genes in Glycyrrhizin Biosynthesis
CYP88D6, CYP72A154, and bAS are key enzymes involved in the biosynthesis of glycyrrhizin (; ; ; ). Here, we took these key enzymes as examples to introduce the application of the platform. Functional annotations for CYP88D6 exhibited on the gene detail page are including annotated information of public platform (Figure 4A), gene sequence information (Figure 4B), gene structure information (Figure 4C), gene family (Figure 4D), KEGG metabolic pathway information (Figure 4E), GO annotation information (Figure 4F), protein domain information (Figure 4G), gene expression profile information (Figure 4H). This gene was annotated as ent-kaurenoic acid hydroxylase or Cytochrome P450 family 88 protein. Ent-kaurenoic acid hydroxylase is a key enzyme in the synthesis pathway of terpenoids and glycyrrhizin is one of the triterpenoid saponins. Therefore, annotation analysis indicated that CYP88D6 may be involved in the synthesis process of glycyrrhizin. The structure information section provided the link of JBrowse, from which one can link to the JBrowse interface to check the gene structure and expression.
FIGURE 4
KEGG annotation information indicated that CYP88D6 might be involved in the triterpenoid synthesis. Related studies have shown that the accumulation of glycyrrhizin will be up-regulated under drought stress (), and the accumulation of glycyrrhizin is in the roots (). Next, we conducted further analysis of the CYP88D6 gene. The expression of this gene in different transcriptomes was analyzed and it was found that the expression of this gene in the lines with high glycyrrhizin yield was higher than that in the lines with low glycyrrhizin yield (Figure 5A). It was also found that the expression of this gene in roots was significantly higher than that in leaves (Figure 5B). Those results are consistent with the accumulation trend of glycyrrhizin. According to the co-expression network of CYP88D6, we found that this gene had a positive co-expression relationship with 46 genes and a negative co-expression relationship with 4 genes (Figure 4C). GSEA analysis of co-expressed genes shows that genesets related to Ent-kaurenoic acid hydroxylase, starch catabolic process significantly enriched (Supplementary Figure S2). Several genes in the CYP88D6 co-expression network were significantly up-regulated under drought stress, and also significantly up-regulated in the root and glycyrrhizin high producing lines (Figure 4C). The results are consistent with the accumulation trend of glycyrrhiza. Finally, we analyzed the heatmaps of CYP88D6 expression in different samples and found the same results with network analysis (Figure 4D). Therefore, our analysis results indicated that CYP88D6 might be involved in the biosynthesis of glycyrrhizin in roots.
FIGURE 5
We searched for CYP72A154 and bAS genes in the platform respectively and obtained annotated information including gene function, gene sequence, gene structure information, gene family, KEGG metabolic pathway information, GO annotation information, protein domain information, and gene expression profile information (Supplementary Figures S3, S4). JBrowse was used to analyze the expression profile and found that CYP72A154 and bAS genes were expressed higher in high-yield lines than that in low-yield lines (Supplementary Figures S5A, S6A), and they were also expressed higher in roots than that in leaves (Supplementary Figures S5B, S6B). These results were consistent with the accumulation trend of glycyrrhizin (). In addition, we used JBrowse to analyze the expression of CYP72A154 under drought stress and salt stress, and found no significant change in CYP72A154 expression compared with the control (Supplementary Figures S5C,D). We also analyzed the co-expression networks of CYP72A154 and bAS under different conditions, and found that the co-expression genes were significantly up-regulated in the high yield lines of glycyrrhizin, and significantly down-regulated in the leaves (Supplemenatry Figure S6C and Supplementary Figure S5E). However, the expression differences of CYP72A154 under drought and salt stress are not obvious (Supplementary Figure S5E).
Discussion
G. uralensis is a traditional Chinese medicinal herb, its flavonoids, triterpenoid saponins, coumarin, and polysaccharides are important functional components (). The release of the genomic data for G. uralensis provides us a considerable gene resource to study underlying biosynthesis mechanisms of bioactive components. However, there is no comprehensive platform for functional genetic analysis and mining in G. uralensis. To ease use by researchers, we developed a publicly available platform called GURFAP, and many tools have been introduced to analyze gene functions. We hope that researchers in the field can use the platform to achieve the new discovery of their study.
Our platform currently offers a searchable system for users, allowing a quick function search and browse of any genes. The obtained information includes structural annotations, functional annotations, gene expression profiles, and co-expression networks of genes, which can provide references for biologists to study the function of G. uralensis genes. All useful tools are publicly available, such as JBrowse (), heatmap analysis, and GSEA, and can also contribute to the functional analysis of genes. JBrowse can be used to display omics data. By integrating the JBrowse tool, users can view and analyze gene expression more intuitively. The platform provides the Heatmap tool to display genes in the gene co-expression network or user-defined gene sets. The GSEA tool can also analyze the gene set provided by the user. In addition to analyzing the genes that have a co-expression relationship with candidate genes, it can also be used to analyze user-defined gene sets, such as differentially expressed genes, which are helpful for researchers to carry out the mining of gene function. There is a certain gap between the information provided by the platform and other databases, and the platform data information and analysis tools still need to be continuously upgraded.
Gene co-expression networks of our platform are based on the correlation of gene expression. Co-expression networks in numerous plants has been developed to infer the gene co-functional relationships, such as studies concerning the co-expression network in Paeonia lactiflora have mining genes related to Oleate Desaturase and Photosynthesis (). Genes in a co-expression network have similar expression patterns, so we reasoned that they may be regulated by the same transcription factors. Our previous co-expression network and motif analysis discovered the biosynthesis pathways of chlorogenic acid and luteolin in L. japonica might be regulated by MYB, bHLH, WRKY (). Here, gene co-expression network and motif analysis show transcription factors might participate in the regulation of glycyrrhizin biosynthesis. Glycyrrhizin is a triterpene saponin, whereas MYB, bHLH, ERF, WKRY involved in the regulation of the biosynthesis of triterpenoids have already been reported (). This analysis method can provide a reference for users to study other secondary metabolite synthesis pathways.
On the other hand, three cases were used to present the usage of our platform. Comprehensive analysis of CYP88D6, CYP72A154, and bAS suggests that it might involve in regulation of glycyrrhizic acid biosysthesis, which had also been reported (). This indicates that our platform is of certain usability and reliability. Although we used CYP88D6, CYP72A154, and bAS to illustrate the usability and reliability of this platform, our data came from public platforms, and most of the information was obtained by computational methods, such as functional annotation information and GO annotation information, etc. Therefore, the platform information is only for reference. To explore the possible functions of genes, molecular biology experiments are still needed.
We believe that omics data of G. uralensis will keep increasing with the development of high-throughput sequencing technologies and the continuous reduction in their cost. Timely and effective collection and the process will assist researchers in their project. Our future work will prioritize effective integration of omics information, and update platform timely. Interesting users can obtain the URL at www.gzybioinformatics.cn/GURFAP.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author contributions
ZS and WX designed the study. JY and QX analyzed data, constructed platform, and wrote the manuscript. YL, LD, and HY participated in the construction of the platform.
Funding
This work was supported by the National Natural Science Foundation of China (NO.32160139) and Ph.D. Startup Foundation of Guizhou University of Traditional Chinese Medicine (2020)32 and (2019)141.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.823966/full#supplementary-material
References
1
BroidoA. D.ClausetA. (2019). Scale-Free Networks Are Rare. Nat. Commun.10 (1), 1017. 10.1038/s41467-019-08746-5
2
BuelsR.YaoE.DieshC. M.HayesR. D.Munoz-TorresM.HeltG.et al (2016). JBrowse: a Dynamic Web Platform for Genome Visualization and Analysis. Genome Biol.17, 66. 10.1186/s13059-016-0924-1
3
ChengM.ZhangJ.YangL.ShenS.LiP.YaoS.et al (2021). Recent Advances in Chemical Analysis of Licorice (Gan-Cao). Fitoterapia149, 104803. 10.1016/j.fitote.2020.104803
4
WangC.ChenL.XuC.ShiJ.ChenS.TanM.et al (2020). A Comprehensive Review for Phytochemical, Pharmacological, and Biosynthesis Studies onGlycyrrhizaspp. Am. J. Chin. Med.48 (1), 17–45. 10.1142/S0192415X20500020
5
WangD.FanW.GuoX.WuK.ZhouS.ChenZ.et al (2020). MaGenDB: a Functional Genomics Hub for Malvaceae Plants. Nucleic Acids Res.48 (D1), D1076–D1084. 10.1093/nar/gkz953
6
DuBoisG. E.PrakashI. (2012). Non-Caloric Sweeteners, Sweetness Modulators, and Sweetener Enhancers. Annu. Rev. Food Sci. Technol.3, 353–380. 10.1146/annurev-food-022811-101236
7
El-GebaliS.MistryJ.BatemanA.EddyS. R.LucianiA.PotterS. C.et al (2019). The Pfam Protein Families Database in 2019. Nucleic Acids Res.47 (D1), D427–D432. 10.1093/nar/gky995
8
FranzM.LopesC. T.HuckG.DongY.SumerO.BaderG. D. (2016). Cytoscape.js: a Graph Theory Library for Visualisation and Analysis. Bioinformatics32 (2), btv557–311. 10.1093/bioinformatics/btv557
9
GaoT.LiuZ.WangY.ChengH.YangQ.GuoA.et al (2013). UUCD: a Family-Based Database of Ubiquitin and Ubiquitin-like Conjugation. Nucleic Acids Res.41 (Database issue), D445–D451. 10.1093/nar/gks1103
10
GuoW.ChenJ.LiJ.HuangJ.WangZ.LimK.-J. (2020). Portal of Juglandaceae: A Comprehensive Platform for Juglandaceae Study. Hortic. Res.7, 35. 10.1038/s41438-020-0256-x
11
HayashiH.HirotaA.HiraokaN.IkeshiroY. (1999). Molecular Cloning and Characterization of Two cDNAs for Glycyrrhiza Glabra Squalene Synthase. Biol. Pharm. Bull.22 (9), 947–950. 10.1248/bpb.22.947
12
HouZ.SunG. (2021). Predictive Quality Control for Compound Liquorice Tablets by the Intelligent Mergence Fingerprint Method Combined with the Systematic Quantitative Fingerprint Method. Phytochem. Anal.32 (6), 1118–1130. 10.1002/pca.3053
13
HuanC.XuY.ZhangW.GuoT.PanH.GaoS. (2021). Research Progress on the Antiviral Activity of Glycyrrhizin and its Derivatives in Liquorice. Front. Pharmacol.12, 680674. 10.3389/fphar.2021.680674
14
KanehisaM.SatoY.MorishimaK. (2016). BlastKOALA and GhostKOALA: KEGG Tools for Functional Characterization of Genome and Metagenome Sequences. J. Mol. Biol.428 (4), 726–731. 10.1016/j.jmb.2015.11.006
15
KellnerF.KimJ.ClavijoB. J.HamiltonJ. P.ChildsK. L.VaillancourtB.et al (2015). Genome-guided Investigation of Plant Natural Product Biosynthesis. Plant J.82 (4), 680–692. 10.1111/tpj.12827
16
KhaninR.WitE. (2006). How Scale-free Are Biological Networks. J. Comput. Biol.13 (3), 810–818. 10.1089/cmb.2006.13.810
17
KimD.PaggiJ. M.ParkC.BennettC.SalzbergS. L. (2019). Graph-based Genome Alignment and Genotyping with HISAT2 and HISAT-Genotype. Nat. Biotechnol.37 (8), 907–915. 10.1038/s41587-019-0201-4
18
KwonY.-J.SonD.-H.ChungT.-H.LeeY.-J. (2020). A Review of the Pharmacological Efficacy and Safety of Licorice Root from Corroborative Clinical Trial Findings. J. Med. Food23 (1), 12–20. 10.1089/jmf.2019.4459
19
LameschP.BerardiniT. Z.LiD.SwarbreckD.WilksC.SasidharanR.et al (2012). The Arabidopsis Information Resource (TAIR): Improved Gene Annotation and New Tools. Nucleic Acids Res.40 (Database issue), D1202–D1210. 10.1093/nar/gkr1090
20
LeeC. M.BarberG. P.CasperJ.ClawsonH.DiekhansM.GonzalezJ. N.et al (2020). UCSC Genome Browser Enters 20th Year. Nucleic Acids Res.48 (D1), D756–D761. 10.1093/nar/gkz1012
21
LiW. D.HouJ. L.WangW. Q.TangX. M.LiuC. L.XingD. (2011). Effect of Water Deficit on Biomass Production and Accumulation of Secondary Metabolites in Roots of Glycyrrhiza Uralensis. Russ. J. Plant Physiol.58 (3), 538–542. 10.1134/S1021443711030101
22
LiY.ChenC.XieZ.XuJ.WuB.WangW. (2020). Integrated Analysis of mRNA and microRNA Elucidates the Regulation of Glycyrrhizic Acid Biosynthesis in Glycyrrhiza Uralensis Fisch. Ijms21 (9), 3101. 10.3390/ijms21093101
23
LombardV.Golaconda RamuluH.DrulaE.CoutinhoP. M.HenrissatB. (2014). The Carbohydrate-Active Enzymes Database (CAZy) in 2013. Nucl. Acids Res.42 (Database issue), D490–D495. 10.1093/nar/gkt1178
24
MochidaK.SakuraiT.SekiH.YoshidaT.TakahagiK.SawaiS.et al (2017). Draft Genome Assembly and Annotation of Glycyrrhiza Uralensis , a Medicinal Legume. Plant J.89 (2), 181–194. 10.1111/tpj.13385
25
NelsonD. R. (2009). The Cytochrome P450 Homepage. Hum. Genomics4 (1), 59–65. 10.1186/1479-7364-4-1-59
26
OughtredR.StarkC.BreitkreutzB.-J.RustJ.BoucherL.ChangC.et al (2019). The BioGRID Interaction Database: 2019 Update. Nucleic Acids Res.47 (D1), D529–D541. 10.1093/nar/gky1079
27
PerteaM.PerteaG. M.AntonescuC. M.ChangT. C.MendellJ. T.SalzbergS. L. (2015). StringTie Enables Improved Reconstruction of a Transcriptome From RNA-Seq Reads. Nat. Biotechnol.33, 290–295. 10.1038/nbt.3122
28
PotterS. C.LucianiA.EddyS. R.ParkY.LopezR.FinnR. D. (2018). HMMER Web Server: 2018 Update. Nucleic Acids Res.46 (W1), W200–W204. 10.1093/nar/gky448
29
PuX.LiZ.TianY.GaoR.HaoL.HuY.et al (2020). The Honeysuckle Genome Provides Insight into the Molecular Mechanism of Carotenoid Metabolism Underlying Dynamic Flower Coloration. New Phytol.227 (3), 930–943. 10.1111/nph.16552
30
ReiserL.SubramaniamS.LiD.HualaE. (2017). Using the Arabidopsis Information Resource (TAIR) to Find Information about Arabidopsis Genes. Curr. Protoc. Bioinformatics60, 1 11 11–11 11 45. 10.1002/cpbi.36
31
SekiH.OhyamaK.SawaiS.MizutaniM.OhnishiT.SudoH.et al (2008). Licorice β-amyrin 11-oxidase, a Cytochrome P450 with a Key Role in the Biosynthesis of the Triterpene Sweetener Glycyrrhizin. Proc. Natl. Acad. Sci. U.S.A.105 (37), 14204–14209. 10.1073/pnas.0803876105
32
SekiH.SawaiS.OhyamaK.MizutaniM.OhnishiT.SudoH.et al (2011). Triterpene Functional Genomics in Licorice for Identification of CYP72A154 Involved in the Biosynthesis of Glycyrrhizin. Plant Cell23 (11), 4112–4123. 10.1105/tpc.110.082685
33
SheJ.YanH.YangJ.XuW.SuZ. (2019). croFGD: Catharanthus Roseus Functional Genomics Database. Front. Genet.10, 238. 10.3389/fgene.2019.00238
34
ShengM.SheJ.XuW.HongY.SuZ.ZhangX. (2020). HpeNet: Co-expression Network Database for De Novo Transcriptome Assembly of Paeonia Lactiflora Pall. Front. Genet.11, 570138. 10.3389/fgene.2020.570138
35
ShibuyaM.KatsubeY.OtsukaM.ZhangH.TansakulP.XiangT.et al (2009). Identification of a Product Specific β-amyrin Synthase from Arabidopsis thaliana. Plant Physiol. Biochem.47 (1), 26–30. 10.1016/j.plaphy.2008.09.007
36
SongX.NieF.ChenW.MaX.GongK.YangQ.et al (2020). Coriander Genomics Database: a Genomic, Transcriptomic, and Metabolic Database for Coriander. Hortic. Res.7, 55. 10.1038/s41438-020-0261-0
37
SonnhammerE. L. L.ÖstlundG. (2015). InParanoid 8: Orthology Analysis between 273 Proteomes, Mostly Eukaryotic. Nucleic Acids Res.43 (Database issue), D234–D239. 10.1093/nar/gku1203
38
SuW.JingY.LinS.YueZ.YangX.XuJ.et al (2021). Polyploidy Underlies Co-option and Diversification of Biosynthetic Triterpene Pathways in the Apple Tribe. Proc. Natl. Acad. Sci. U.S.A.118 (20), e2101767118. 10.1073/pnas.2101767118
39
TangZ.-H.LiT.TongY.-G.ChenX.-J.ChenX.-P.WangY.-T.et al (2015). A Systematic Review of the Anticancer Properties of Compounds Isolated from Licorice (Gancao). Planta Med.81 (18), 1670–1687. 10.1055/s-0035-1558227
40
TianT.LiuY.YanH.YouQ.YiX.DuZ.et al (2017). agriGO v2.0: a GO Analysis Toolkit for the Agricultural Community, 2017 Update. Nucleic Acids Res.45 (W1), W122–W129. 10.1093/nar/gkx382
41
WaeseJ.ProvartN. J. (2017). The Bio-Analytic Resource for Plant Biology. Methods Mol. Biol.1533, 119–148. 10.1007/978-1-4939-6658-5_6
42
WangL.WangS.LiW. (2012). RSeQC: Quality Control of RNA-Seq Experiments. Bioinformatics28 (16), 2184–2185. 10.1093/bioinformatics/bts356
43
WangH.SongW.TaoW.ZhangJ.ZhangX.ZhaoJ.et al (2021). Identification Wild and Cultivated Licorice by Multidimensional Analysis. Food Chem.339, 128111. 10.1016/j.foodchem.2020.128111
44
WuZ.LiuH.ZhanW.YuZ.QinE.LiuS.et al (2021). The Chromosome‐scale Reference Genome of Safflower ( Carthamus tinctorius ) Provides Insights into Linoleic Acid and Flavonoid Biosynthesis. Plant Biotechnol. J.19, 1725–1742. 10.1111/pbi.13586
45
XiaoQ.LiZ.QuM.XuW.SuZ.YangJ. (2021). LjaFGD: Lonicera japonica Functional Genomics Database. J. Integr. Plant Biol.63 (8), 1422–1436. 10.1111/jipb.13112
46
XieW.HaoZ.ZhouX.JiangX.XuL.WuS.et al (2018). Arbuscular Mycorrhiza Facilitates the Accumulation of Glycyrrhizin and Liquiritin in Glycyrrhiza Uralensis under Drought Stress. Mycorrhiza28 (3), 285–300. 10.1007/s00572-018-0827-y
47
XuG.CaiW.GaoW.LiuC. (2016). A Novel Glucuronosyltransferase Has an Unprecedented Ability to Catalyse Continuous Two‐step Glucuronosylation of Glycyrrhetinic Acid to Yield Glycyrrhizin. New Phytol.212 (1), 123–135. 10.1111/nph.14039
48
YangR.WangL.-q.LiuY. (2014). Antitumor Activities of Widely-Used Chinese Herb-Licorice. Chin. Herbal Medicines6 (4), 274–281. 10.1016/s1674-6384(14)60042-3
49
YangR.YuanB.-C.MaY.-S.ZhouS.LiuY. (2017). The Anti-inflammatory Activity of Licorice, a Widely Used Chinese Herb. Pharm. Biol.55 (1), 5–18. 10.1080/13880209.2016.1225775
50
YangJ.LiuY.YanH.TianT.YouQ.ZhangL.et al (2018). PlantEAR: Functional Analysis Platform for Plant EAR Motif-Containing Proteins. Front. Genet.9, 590. 10.3389/fgene.2018.00590
51
YangJ.XiaoQ.XuJ.DaL.GuoL.HuangL.et al (2020). GelFAP: Gene Functional Analysis Platform for Gastrodia Elata. Front. Plant Sci.11, 563237. 10.3389/fpls.2020.563237
52
YaoL.LuJ.WangJ.GaoW.-Y. (2020). Advances in Biosynthesis of Triterpenoid Saponins in Medicinal Plants. Chin. J. Nat. Med.18 (6), 417–424. 10.1016/S1875-5364(20)30049-2
53
YiX.DuZ.SuZ. (2013). PlantGSEA: a Gene Set Enrichment Analysis Toolkit for Plant Community. Nucleic Acids Res.41 (Web Server issue), W98–W103. 10.1093/nar/gkt281
54
YiZ.ChenJ.SunH.RosliH. G.PomboM. A.ZhangP.et al (2016). iTAK: a Program for Genome-wide Prediction and Classification of Plant Transcription Factors, Transcriptional Regulators, and Protein Kinases. Mol. Plant9 (12), 1667–1670. 10.1016/j.molp.2016.09.014
55
YuP.LiQ.FengY.MaS.ChenY.LiG. (2021). Extraction and Analysis of Six Effective Components in Glycyrrhiza Uralensis Fisch by Deep Eutectic Solvents (DES) Combined with Quantitative Analysis of Multi-Components by Single Marker (QAMS) Method. Molecules26 (5), 1310. 10.3390/molecules26051310
56
YuanY.JinX.LiuJ.ZhaoX.ZhouJ.WangX.et al (2018). The Gastrodia Elata Genome Provides Insights into Plant Adaptation to Heterotrophy. Nat. Commun.9 (1), 1615. 10.1038/s41467-018-03423-5
57
YueJ.LiuJ.TangW.WuY. Q.TangX.LiW.et al (2020). Kiwifruit Genome Database (KGD): a Comprehensive Resource for Kiwifruit Genomics. Hortic. Res.7, 117. 10.1038/s41438-020-0338-9
Summary
Keywords
Glycyrrhiza uralensis, Platform, gene function analyses, co-expression network, glycyrrhizin
Citation
Yang J, Yan H, Liu Y, Da L, Xiao Q, Xu W and Su Z (2022) GURFAP: A Platform for Gene Function Analysis in Glycyrrhiza Uralensis. Front. Genet. 13:823966. doi: 10.3389/fgene.2022.823966
Received
28 November 2021
Accepted
08 March 2022
Published
12 April 2022
Volume
13 - 2022
Edited by
Jia Meng, Xi’an Jiaotong-Liverpool University, China
Updates
Copyright
© 2022 Yang, Yan, Liu, Da, Xiao, Xu and Su.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qiaoqiao Xiao, xqqiao2021@163.com; Wenying Xu, x_wenying@yahoo.com; Zhen Su, zhensu@cau.edu.cn
This article was submitted to Computational Genomics, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.