Next Article in Journal
Expanded Ligands Based upon Iron(II) Coordination Compounds of Asymmetrical Bis(terpyridine) Domains
Previous Article in Journal
Fractionation Coupled to Molecular Networking: Towards Identification of Anthelmintic Molecules in Terminalia leiocarpa (DC.) Baill
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Essay

Flavonoid Metabolism in Tetrastigma hemsleyanum Diels et Gilg Based on Metabolome Analysis and Transcriptome Sequencing

1
Zhejiang Provincial Key Laboratory of Resources Protection and Innovation of Traditional Chinese Medicine, Hangzhou 311300, China
2
College of Food and Health, Department of Traditional Chinese medicine, Zhejiang Agriculture & Forestry University, Hangzhou 311300, China
3
State Key Laboratory of Subtropical Silviculture, Zhejiang Agriculture and Forestry University, Hangzhou 311300, China
4
Botanical Garden, Zhejiang Agricultural and Forestry University, Hangzhou 311300, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Molecules 2023, 28(1), 83; https://doi.org/10.3390/molecules28010083
Submission received: 27 October 2022 / Revised: 13 December 2022 / Accepted: 14 December 2022 / Published: 22 December 2022

Abstract

:
Tetrastigma hemsleyanum Diels et Gilg, known as a “plant antibiotic”, possesses several attractive properties including anti-inflammatory, anti-tumor, and antioxidant effects, with its efficacy being attributed to flavonoids. However, the flavonoid biosynthesis of T. hemsleyanum has rarely been studied. In this study, we investigated the flavonoid metabolism of T. hemsleyanum through metabolome analysis and transcriptome sequencing. The metabolomic results showed differences in the flavonoids of the leaves and root tubers of T. hemsleyanum. A total of 22 flavonoids was detected, and the concentrations of most flavonoids in the leaves were higher than those in the root tubers. Transcriptome analysis revealed that differentially expressed genes (DEGs) in the leaves and root tubers were enriched in photosynthesis-antenna proteins. Pearson correlation analysis indicated that the expression levels of chalcone isomerase (CHI) and UDP-glycose flavonoid glycosyltransferase (UFGT) were highly correlated with the concentrations of most flavonoids. Further, this study found that the photosynthesis-antenna proteins essentially contributed to the difference in the flavonoids in T. hemsleyanum. The gene expressions and concentrations of the total flavonoids of leaves and root tubers in Hangzhou, Jinhua, Lishui, and Taizhou in Zhejiang Province, China, showed that CHI (CL6715.Contig1_All, Unigene19431_All, CL921.Contig4_All) and UFGT (CL11556.Contig3_All, CL11775.Contig1_All) were the potential key genes of accumulation of most flavonoids in T. hemsleyanum.

1. Introduction

Tetrastigma hemsleyanum Diels et Gilg is a rare Chinese herbal medicine that is also known as a “plant antibiotic”. It is a folk medicine mainly produced in Zhejiang Province, and is widely distributed in Hangzhou, Jinhua, Lishui, Taizhou, and Ningbo, among other regions. Traditionally, it can be used to cure fever, detoxify, and promote blood circulation [1]. T. hemsleyanum has gradually become widely accepted due to its anti-tumor, anti-inflammatory, and antioxidant effects. The Hua Shi Xuan Fei mixture, a prescription that makes a difference in the clinical treatment of COVID-19 (approval number of Zhejiang medicine, Z20200026000), is mainly composed of T. hemsleyanum [2]. Considering its use in medicine, T. hemsleyanum has been excessively developed for the multi-application of its root tubers over the last 3–5 years, which has led to a sharp decline in its natural resources [3,4]. Further, its chaotic quality management has resulted in the medicine being characterized by inconsistent curative effects. Consequently, we have focused our efforts toward flavonoid biosynthesis research in order to improve the quality of T. hemsleyanum.
Analysis of T. hemsleyanum’s metabolic components shows numerous chemical components in the whole plant, which mainly include flavonoids, triterpenoids, organic acids, and fatty acids [5]. Flavonoids, a type of defense substance, are produced when T. hemsleyanum is subjected to environmental stresses. Flavonoids are generally considered to be the main active constituent of T. hemsleyanum. However, reports on T. hemsleyanum in recent years have mainly focused on variety optimization, cultivation technology management [6,7], active ingredient extraction, and related pharmacological experiments. Information on the flavonoid composition and biosynthesis genes of T. hemsleyanum is limited. Phenylalanine ammonia lyase (PAL), a key bridge connecting primary metabolism, is the rate-limiting enzyme in the first step of phenylpropanoid metabolism [8]. The biosynthesis of flavonoids is branched from the general phenylpropanoid pathway via the catalysis of chalcone synthase (CHS) and chalcone isomerase (CHI) [9]. CHS is the first rate-limiting enzyme in the synthetic pathway and belongs to the plant polyketide synthase superfamily [10]. It is mainly responsible for the binding of one molecule of 4-coumaroyl-CoA to three molecules of malonyl-CoA. CHI, which is the second rate-limiting enzyme involved in the biosynthetic pathway of flavonoids, catalyzes the intramolecular cyclization reaction and converts bicyclic chalcone into tricyclic (2S)-flavanone [11]. In addition to PAL, CHS, and CHI, many enzymes are involved in the flavonoid synthesis pathway, such as flavonol synthase (FLS), flavone synthase (FNS), and UDP-glycose flavonoid glycosyltransferase (UFGT). FLS is a key enzyme in the synthesis of flavonols, and it synthesizes target products using dihydroflavonols as substrates in Arabidopsis thaliana (L.) Heynh [12]. FNS is a key enzyme in the synthesis of flavones; it introduces a double bond between the C2 and C3 positions of flavanones to convert flavanones into flavones [13]. The UFGT enzyme is the final gene of the anthocyanin synthesis pathway, and it transfers the glucosyl moiety from UDP-glucose to the 3-hydroxyl group of anthocyanidins [14].
The co-application of transcriptome and metabolome analyses has become an effective method to investigate biosynthetic pathways; it can identify novel functional genes and key genes and classify the regulatory mechanisms of pathways [15]. In this study, we utilized the abovementioned method to explore the mechanism of flavonoid accumulation in the leaves and root tubers of T. hemsleyanum, as well as to explore the mechanism of T. hemsleyanum’s flavonoid metabolic pathway and its core genes. We also clarified the synthesis mechanism of the flavonoids in the leaves and root tubers of T. hemsleyanum and discovered the key genes involved in the biosynthesis of these flavonoids.

2. Results

2.1. Flavonoid Components Differ in Leaves and Root Tubers of T. hemsleyanum

A total of 22 flavonoid monomers was identified in T. hemsleyanum, and their accumulation was obviously different in its leaves and root tubers. The principal component analysis (PCA) score plot showed that the two organs were markedly separated, and the metabolome sequencing analysis had good repeatability. Indeed, PC1 and PC2 (82.2% and 8.1%, respectively) were clearly separated among the root tubers and leaves (Figure 1A).
The cluster heatmap of the 22 metabolites (Figure 1B) showed that the accumulation levels of various components including kaempferol, eriodictyol, astragaline, apigenin, and isorhamnetin were clustered together. In the metabolome analysis, luteolin, luteoloside, hesperidin, rutin, lonicerin, baicalin, and apigenin-7-O-glycoside were found to only be contained in the leaves, and kaempferol was only found to be contained in the root tubers. Concentrations of vitexin, orientin, lonicerin, rutin, isoquercitrin, luteolin, luteoloside, hesperidin, and naringin in the leaves were higher than those in the root tubers. For instance, the concentration of vitexin, which reached 418.828 μg·g−1, was higher in the leaves than in the root tubers. The concentration of orientin in the leaves was 3283.919 times higher than that in the root tubers. By contrast, the astragaline content in the root tubers was 27.636 times higher than that in the leaves (Figure 1C,D and Table S1).

2.2. Transcriptome Sequencing of Leaves and Root Tubers of T. hemsleyanum

A total of 47,471,677 bp non-repetitive sequences was detected in the leaves, and 53,338,272 bp were found in the root tubers. The average lengths of the transcripts were 895.03 and 875.07 bp. The longest non-repetitive gene was 15,654 bp, and the shortest was 301 bp. Among them, the length of the N50 of the leaves was 1275 bp, and the length of the N50 of the root tubers was 1236 bp. These results indicated a high-quality assembly (Table 1).
To understand their relevant biological significance, the assembled unigenes were searched against seven functional databases (NR, GO, KOG, KEGG, and SwissProt) and annotated. As shown in Figure 2, 61,820 functional genes were annotated after searching the transcriptome sequences in KOG, NR, SwissProt, KEGG, and GO. Of the 61,820 annotated genes, 28,757 were found in leaves and 33,063 were observed in the root tubers. A total of 28,757 functional genes in the leaves were jointly annotated in five functional databases (Figure 2A). In the root tubers, 33,063 functional genes were annotated in all databases (Figure 2B).
The assembled unigenes were statistically annotated for three biological functions in GO using Blast2GO. The study showed 22 groups in biological processes, 25,003 in cellular components, and 303 in molecular function (Figure 3A). The unigenes that matched with the KEGG database were divided into seven branches, namely, cellular processes, environmental information processing, genetic information processing, human diseases, metabolism, organic systems, and drug development. In addition, 1751 genes were related to secondary metabolism, and 2517 genes were related to environmental adaptation (Figure 3B).

2.3. Differentially Expressed Genes of Leaves and Root Tubers of T. hemsleyanum

Through our further study of the various differentially expressed genes (DEGs) in the leaves and root tubers of T. hemsleyanum, a result similar to that of the differential metabolite analysis emerged. Compared with the leaves of T. hemsleyanum, 67,345 DEGs were identified in the root tubers of T. hemsleyanum, including 6876 up-regulated genes, 14,277 down-regulated genes, and 46,192 non-DEGs (Figure 4A). When they were divided into the seven KEGG pathways used for the assembled unigenes, the majority of the DEGs were obviously enriched in the metabolism-related pathways including cellular processes, environmental information processing, genetic information processing, human diseases, metabolism, organismal systems, and drug development (Figure 4B). The 20 most significant pathways in the leaves versus the roots are listed in Figure 4C. Additionally, the DEGs in the two different medicinal portions were enriched in various pathways, including the photosynthesis, photosynthesis-antenna proteins, flavone and flavonol biosynthesis, phenylpropanoid biosynthesis, anthocyanin biosynthesis, isoflavonoid biosynthesis, biosynthesis of the secondary metabolites, plant mitogen-activated protein kinase (MAPK) signaling, flavonoid biosynthesis, and plant hormone signal transduction pathways. Among the 20 most significant pathways, the biosynthesis of the secondary metabolite pathway was enriched with the highest number of DEGs, reaching 2250 DEGs. Additionally, DEG enrichment was most significant in the plant hormone signal transduction pathway. By contrast, DEG enrichment was the most insignificant in alpha-linolenic acid metabolism. The greatest enrichment of DEGs was found in the photosynthesis-antenna proteins pathway (Figure 4C). A total of 25 DEGs belonging to the chlorophyll a/b-binding protein family was enriched in the photosynthesis-antenna proteins pathway (Figure 4D).

2.4. Analysis Correlation between Flavonoid Biosynthetic Gene Expression and Flavonoid Concentration

In this part, the detected flavonoid monomers were divided into six classes of flavonoids according to their basic skeletons (Figure 5). Flavones, flavone glycosides, flavonols, flavonol glycosides, flavanones, and flavanone glycosides constituted the major fractions of the flavonoids in T. hemsleyanum. As illustrated in Figure 5, the flavones included vitexin, orientin, luteolin, apigenin, diosmetin, and chrysoeriol. The flavone glycosides were composed of lonicerin, luteoloside, baicalin, and apigenin-7-O-glycoside. The flavonols were composed of myricetin, kaempferol, and isorhamnetin. The flavonol glycosides consisted of isoquercitrin, rutin, astragalin, and isohamnetin-3-O-glucoside. The flavonones included naringenin and eriodictyol. Finally, the flavanone glycosides consisted of naringin and hesperidin.
A rough flavonoid metabolism pathway of T. hemsleyanum is listed in Figure 6; it contained 11 key genes related to flavonoid synthesis and 22 flavonoid metabolites. The heatmaps of the expression of synthetases revealed that the expression levels of CHI, UFGT, and FNS in the leaves were higher than those in the roots. For instance, the log10 (FPKM) value of CHI1 (CL6715.Contig1_All) in the leaves was higher than that in the root tubers by 0.213 units, and the log10 (FPKM) value of UFGT1 (CL11556.Contig3_All) in leaves was 1.746 times higher than that in root tubers. The expression levels of C4H, 4CL, F3′H, and F3′5′H in the leaves and root tubers were similar. The expression levels of PAL, CHS, F3H, and FNS in the roots were higher than those in the leaves. For instance, the log10 (FPKM) value of CHS (CL1413.Contig5_All) in the root tubers was 1.554 times higher than that in the leaves.
Next, Pearson correlation analysis of the relationship between flavonoid monomer concentrations and the typical FPKM values of the key genes involved in flavonoid biosynthesis were carried out. Interestingly, 15 flavonoids, namely vitexin, orientin, lonicerin, rutin, isoquercitrin, luteolin, luteoloside, hesperidin, naringin, eriodictyol, myricetin, baicalin, chrysoeriol, diosmetin, and quercitrin, were highly correlated with the gene expression of UFGT (Figure 7). Moreover, 12 flavonoids, including apigenin-7-glycoside, quercitrin, eriodictyol, luteolin, isoquercitrin, hesperidin, chrysoeriol, and naringin, were highly correlated with the gene expression of CHI. Three flavonoids, namely, astragaline, naringenin, and kaempferol, were highly correlated with the gene expression levels of CHI, PAL, and C4H. Orientin, isorhamnetin-3-O-glucoside, myricetin, baicalin, and diosmetin were moderately correlated with CHI. Astragaline, naringenin, and kaempferol were moderately correlated with F3H, F3′H, and FLS. Isorhamnetin-3-O-glucoside, apigenin, and apigenin-7-O-glycoside were moderately correlated with UFGT. A low correlation between the concentrations of astragaline, naringenin, and kaempferol and the gene expression of FNS was also found.

2.5. Validation of the DEGs

Next, we used quantitative real-time PCR (qRT-PCR) to determine the expression of 12 DEGs in the root tubers and leaves of T. hemsleyanum (Figure 8A). Consistent with the transcriptome analysis results, the flavonoid-synthesis-related key genes consisted of three CHIs (CL6715.Contig1_All, Unigene19431_All, and CL921.Contig4_All) and two UFGTs (CL11556.Contig3_All and CL11775.Contig1_All) whose gene expression in the leaves were much higher than those in the root tubers. The gene expression of four CHSs (CL1413.Contig5_All, CL1413.Contig12_All, CL1413.Contig4_All, and CL1413.Contig11_All), one FNS (CL21915.Contig3_All) and two PALs (CL2379.Contig9_All and CL2379.Contig4_All) in the root tubers was much higher than that in the leaves, a finding which concurred with the RNA-seq analysis results. In Figure 8B, the scatter plot of the correlation between RNA-seq and relative gene expression indicated a positive correlation (R2 = 0.8171). This result showed that the RNA-seq data were reliable.

2.6. Key Gene Expression and Total Flavonoid Content Trends for T. hemsleyanum from Different Locations

T. hemsleyaum specimens from Hangzhou, Jinhua, Lishui, and Taizhou were collected to detect their concentration of total flavonoids and expression levels of key genes including CHI 1 (CL6715.Contig1_All), CHI 2 (Unigene19431_All), CHI 3 (CL921.Contig4_All), UFGT 1 (CL11556.Contig3_All), and UFGT 2 (CL11775.Contig1_All). The flavonoid synthesis key gene expression results found that these key genes’ relative expressions in the leaves and root tubers showed similar trends to those obtained from the RNA-seq data. Furthermore, the maximum difference was observed in the expression of CHI 2 in the leaves and root tubers from Hangzhou, Jinhua, and Taizhou (Figure 9A). As shown in Figure 9B, across the four locations, the concentrations of the total flavonoids in the leaves were significantly higher than those in the root tubers. Further, a comparison of the total flavonoid content of the specimens from the four locations showed that the leaves and root tubers of Lishui T. hemsleyanum were of a better quality than the others, followed by Taizhou.

3. Discussion

A unified standard for the medicinal part of T. hemsleyanum has been lacking for several years in China. Flora of China records that the whole plant can be used for medicine, whereas Zhejiang Province, as the main production area of T. hemsleyanum, issued Standards of Processing Chinese Crude Drugs, which classifies the root tubers as drug-parts [22,23]. In the present study, metabolomic analysis of the whole plant indicated that the species and concentration of flavonoids in T. hemsleyanum’s leaves were higher than those in the root tubers. Interestingly, a great deal of anti-inflammatory and anti-tumor flavonoid metabolites were found in the leaves. Thus, we suggest that T. hemsleyanum’s leaves should also be used for medicine. Transcriptome analysis revealed that the DEGs in the photosynthesis-antenna proteins played a role in flavonoid metabolism in the leaves and root tubers. To date, few studies have thoroughly analyzed the flavonoid biosynthesis of T. hemsleyanum. In the present study, we explored flavonoid metabolism, finding key genes called CHI and UFGT via metabolomic analysis and transcriptome sequencing.

3.1. Key Role of DEGs in Photosynthesis-Antenna Proteins in Flavonoid Metabolism in Leaves and Root Tubers

In this study, metabolomic analysis showed that the concentrations of most flavonoids in the leaves were higher than those in the root tubers. Compared with the flavonoids in the root tubers, flavones and flavone glycosides were mostly produced in the leaves. Zhang et al. also found identical occurrences [24]. Further, many anti-inflammatory and anti-tumor ingredients were detected in the leaves, such as luteolin, luteoloside, hesperidin, baicalin, chrysoeriol, diosmetin, apigenin, vitexin, orientin, lonicerin, myricetin, eriodictyol, quercitrin, and isoquercitrin. Analysis of metabolites also indicated that the major flavonoids in the leaves of T. hemsleyanum were vitexin and orientin.
KEGG enrichment analysis indicated that the flavonoids in the leaves and root tubers were enriched in the photosynthesis, photosynthesis-antenna proteins, plant MAPK signaling, plant hormone signal transduction, phenylpropanoid biosynthesis, biosynthesis of the secondary metabolite, flavone and flavonol biosynthesis, and flavonoid biosynthesis pathways. The DEGs in these pathways may jointly play a role in regulating key enzymes and may ultimately lead to differences in the concentrations of flavonoid metabolites between the leaves and root tubers. Photosynthesis and photosynthesis-antenna proteins were found to be the key pathways involved in flavonoid synthesis. About 95% of plant biomass accumulation was attributed to photosynthesis, and the regulation of photosynthesis-antenna proteins was found to regulate the ability of the leaves to capture solar energy and transfer energy to reaction centers [25]. Located on the thylakoid membrane of chloroplasts, photosynthesis-antenna proteins were found to have a high expression in the leaves and low expression in the root tubers [26]. Different expression of DEGs enriched in the photosynthesis-antenna protein pathway between the root tubers and leaves resulted in varying levels of light energy being captured, which affected the synthesis of the C skeleton [27]. The C skeleton plays a key role in flavonoid synthesis. Wang et al. [28] found that the down-regulation of photosynthesis-antenna protein-related genes in citrus infected with Candidatus Liberibacter asiaticus resulted in a significant decrease in flavonoids and phenylalanine derivatives. Additionally, the leaves and root tubers of T. hemsleyanum are exposed to ultraviolet (UV) radiation with distinct intensities, also affecting photosynthesis-antenna proteins’ roles in energy capture and flavonoid biosynthesis. Hong et al. indicated that chlorophyll a/b-binding proteins are initially downregulated after shading, which is consistent with the expression pattern of genes including CHI and 3GT in the anthocyanin synthesis pathway [29]. Our preliminary study found that appropriate doses of UV-B and UV-C radiation (30 min–3 h) induced active stress and contributed to the accumulation of flavonoids [30]. Furthermore, we previously conducted different shading experiments on T. hemsleyanum and found that light intensity affected its flavonoid accumulation [15]. In the two different medicinal portions studied in this work, the DEGs in the photosynthesis-antenna proteins pathway affected the expression of DEGs in the secondary metabolite biosynthesis, flavonoid biosynthesis, and flavone and flavonol biosynthesis pathways. This phenomenon led to the discrepancy in flavonoids observed in the two studied medicinal portions.
According to the literature, MAPKs have important functions in abiotic stresses including salt stress, drought stress, and UV light [31,32,33]. MAPKs are also involved in the response to plant hormones [34]. Plant hormones regulate flavonoid metabolism by directly regulating the expression of key genes related to flavonoid biosynthesis [35,36]. However, whether MAPKs and phytohormones are involved in flavonoid biosynthesis in root tubers and leaves of T. hemsleyanum requires further study.

3.2. CHI and UFGT Are Critical Genes for Flavonoid Biosynthesis in T. hemsleyanum

Flavonoid metabolites may play a critical role in plant growth, as well as plant resistance to biological and abiotic stresses, such as extreme temperature, drought, UV light, nutrient deficiencies, pathogens, and herbivores; however, their biosynthetic pathways and core genes in T. hemsleyanum are yet to be uncovered [37].
Following Peng’s research [19], we plotted the framework of the flavonoid metabolic pathway of T. hemsleyanum and classified its flavonoid metabolites. We concretized the metabolic pathway reference to what is known of flavonoid metabolic pathways in plants such as citrus, celery, and buckwheat. PAL, C4H, and 4CL are jointly involved in the synthesis of P-coumaroyl-coA. CHS contributes to the production of chalcone. Chalcone can transform into naringenin by catalyzing the cyclization reaction with CHI. Naringenin is glucosylated and rhamnosylated by 7-glucosyltransferase and 1-2-rhamnosyltransferase to naringin [16], and hesperidin is also derived from naringenin. Naringenin’s B-3′ position is hydroxylated by F3′H to accumulate eriodictyol. F3′5′H can use eriodictyol to generate flavanones. FNS catalyzes the conversion of flavanones to various flavones including vitexin, orientin, apigenin, and chrysoeriol. The production of lonicerin, baicalin, and luteoloside undergoes glycosylation. Apigenin is methylated and glycosylated to diosmetin and apigenin-7-O-glycosides by OMT and UFGT, respectively [21]. By catalyzing with F3′H or F3′5′H, apigenin transforms into luteolin [21]. In another synthesis branch, F3H catalyzes flavanones at the C-3 position to accumulate dihydroflavonol. The conversion of dihydroflavonol to anthocyanidins with a stable structure is catalyzed by DFR, LDOX, and UFGT. FLS can also use dihydroflavonol as a substrate to produce various flavonols including myricetin, kaempferol, and quercitrin. Quercitrin is methylated and glycosylated to isoquercitrin and isorhamnetin, respectively [17]. By catalyzing with 1,6-Rhamnosyltransferase, isorhamnetin transforms into rutin [18]. Subsequently, flavonol is glycosylated to flavonol glycosides such as isorhamnetin-3-O-glucoside and astragalin.
In our study, some key genes related to the accumulation of flavonoids were found using Pearson correlation analysis. The total flavonoid and highest flavonoid concentration trends were similar to those of the gene expressions of CHI (CL6715.Contig1_All, Unigene19431_All, and CL921.Contig4_All) and UFGT (CL11556.Contig3_All and CL11775.Contig1_All) in T. hemsleyanum. Regarding the key gene expressions and total flavonoid content of T. hemsleyanum from different locations, the results indicated that CHI and UFGT were the candidate genes for the flavonoid biosynthesis pathway. Nam Il Park et al. also demonstrated that transgenic hairy roots with high SbCHI expression mediated by Agrobacterium rhizogenes can synthesize many flavones [38]. Although UFGT is a downstream enzyme that is generally considered a key gene related to the synthesis of anthocyanin, some studies have indicated that UFGT also plays an important role in the synthesis of other flavonoids [39,40]. The catalytic rate of the enzyme and the consumption rate of the substrate may lead to this phenomenon.

4. Materials and Methods

4.1. Plant Materials

Fresh and healthy three-year-old T. hemsleyanum plotted plants from Hangzhou, Jinhua, Lishui, and Taizhou were collected in summer from Zhejiang Agriculture and Forestry University (Hangzhou Zhejiang, 30°15′32″ N, 119°43′23.14″ E) in summer. Sample collection was repeated three times. The collected root tubers and leaves were immediately frozen in liquid nitrogen and then stored at −80 °C.

4.2. Determination of Flavonoid Metabolomes

Approximately 0.1 g of root tubers or leaf powder was placed accurately in a 2.0 mL EP tube, and to this was added 1.5 mL of 80% methanol (v/v). Vitexin, orientin, afzelin, rutin, isoquercitrin, luteolin, luteoloside, isoorientin, naringin, isorhamnetin-3-O-glucoside, isorhamnetin, astragaline, eriodictyol, taxifolin, isovitexin, calycosin-7-O-glucoside, hispidulin, naringenin, apigenin, quercitrin, and apigenin-7-O-glycoside were used as reference substances, and their fabrication processes were consistent with the abovementioned process. The samples were sonicated for 30 min, and the extractions were filtered through a 0.22 μm PTFE membrane to obtain the filtrate. The obtained filtrate was sent to Kyushu Kexin (Beijing, China) Testing Technology Co., Ltd. for metabolome determination. The analytic conditions were as follows:
Chromatographic conditions: The testing instrument was an API 6500 quadrupole-linear ion trap mass spectrometer (Framingham, MA, USA). An ACQUITYUPLCBEH C18 column (2.1 mm × 100 mm, 1.8 μm) was used with 0.1% formic acid-water (A)-0.1% formic acid-acetonitrile (B) as moving phase, gradient elution (0 min, 5% B; 1 min, 25% B; 3.5 min, 40% B; 4.5 min, 60% B), and flow rate of 0.6 mL/min. The column temperature was 40 °C, and the injection volume was 1 μL.
Ion source: Turbo V; ionization mode: ESI-; gas curtain gas (CUR) volume flow rate: 30 L/min; spray voltage (IS): −4500 V; nebulizing gas (GS1) volume flow rate: 55 L/min; heating assist gas (GS2) volume flow rate: 55 L·min−1.
Acquisition method: multiple reaction monitoring mode (MRM), ionization temperature (TEMP): 550 °C.
Investigation of linear relationship: The peak area was measured according to the above conditions, and the regression equation was obtained using the injection concentration X (μg/mL) as the horizontal coordinate and the chromatographic peak area Y as the vertical coordinate as shown in Table 2.

4.3. RNA Extraction and Transcriptome Sequencing

Total RNA was extracted using the plant total RNA extraction kit (Vazyme). We rapidly ground 50 mg of fresh leaves and root tubers in liquid nitrogen and then added buffer PSL for lysis. The RNA of T. hemsleyanum was extracted according to the manufacturer’s instruction. RNA purity, concentration, and integrity were tested with BERTHOLD Colibri LB 915 Microvolume (Berthold, Germany) and Bio-Rad PowerPac (Hercules, CA, USA). In addition, 7μL of RNA was used for reverse transcription into cDNA reference in accordance with the TaKaRa reverse transcription kit. The extracted cDNA samples of T. hemsleyanum of Lishui were sent to Sangon Biotech (Xuhui, SHH, China) Co., Ltd., for transcriptome sequencing. Libaries were constructed using an Agilent 2100 Bioanalyzer and ABI StepOnePlus Real-Time PCR System for quality control. The qualified libraries were sequenced on the Illumina HiSeq platform. After the raw data of low quality, adapter contamination, and high unknown N concentration were removed, clean reads were obtained.

4.4. De Novo Assembly and Sequence Clustering

Trinity v2.5.1 was used for de novo assembly, and Tgicl was used for clustering the assembled transcripts to obtain unigenes. The clean reads were segmented into a large number of individual de Bruijn plots from which full-length transcript shear isoforms were extracted. The transcripts were clustered and redundantly processed with Tgicl. For multiple samples, Tgicl was used again for each unigene of each sample to obtain the final unigene for subsequent analysis.

4.5. FPKM Calculation and DEGs Screening

The FPKM value calculation formula is as follows:
FPKM = 1,000,000 × C N × L × 1000
C is the number of fragments compared with the transcript, N is the total number of reads of all transcripts, and L is the number of transcript bases.
Differential expression of genes was detected using Audic’s [41] sequencing-based method. In this research, false discovery rate (FDR) ≤0.001 and fold change ≥2 were used as criteria for the determination of DEGs.

4.6. Functional Annotation and Classification of Genes

NR, GO, KOG, KEGG, and SwissProt belong to functional databases. In this section, the functional annotation and classification of unigenes were performed by Blast v2.2.23.
Differentially expressed genes were annotated and classified into the KEGG pathways. According to the annotation and classification results, DEGs were subjected to biological pathway classification. In the meantime, phyper in the R package was used for enrichment analysis. FDR ≤0.01 was used as a criterion for the rectification of p-values.

4.7. Quantitative Real-Time PCR

The primer sequences are listed in Table 3. These primers were designed using an online site (https://www.genscript.com/tools/real-time-pcr-taqman-primer-design-tool?act=register_success, accessed on 7 June 2022) according to the CDS sequences and synthesized by SUNYA. MDH was used as a reference gene [42]. A quantitative real-time PCR experiment using specialized reagents (TaKaRa) followed the CFX96 Real-Time PCR Detection System instructions, and three technical replicates were prepared for each sample. The 25 μL system contained 2 μL of cDNA, 12.5 μL of TB Green Premix Ex Taq Ⅱ(Tli RNaseH Plus), 8.5 μL of sterile water, and 2 μL of forward and reverse primers. The two-step PCR reaction was incubated at 95 °C for 30 s for 40 cycles and then incubated at 95 °C for 5 s and 60 °C for 30 s. Gene expression was calculated by 2−ΔΔCt, and special amplification of primers was judged with melting curve analysis.

4.8. Total Flavonoid Determination

After being freeze-dried, the root tubers and leaves of T. hemsleyanum were ground into a powder that could be passed through an 80 mesh sieve net. The total flavonoid content was determined by UV spectrophotometry according to a method that improved on that of Han et al. [6].
Linear relationship investigation: With rutin as the standard, the test scheme was consistent with the above procedure. The absorbance was measured at 501 nm (the maximum absorption wavelength obtained by scanning the field). With absorbance A as the ordinate and reference mass concentration C as the abscissa, the linear regression equation was A = 12.045 C + 0.0019, r2 = 0.99997.

4.9. Statistical Analysis

Metabolome and q-PCR results were calculated using Excel 2020 (Redmond, WA, USA), and data were expressed as mean ± standard error (SE). One-way analysis of variance (ANOVA) and t-test with IBM SPSS statistical software 24.0 (Armonk, NY, USA) were used for statistical significance analysis. As for RNA-seq results, the R package was used for analysis and drawing. PCR analysis, heatmaps, bar graphs, scatter plots, and Pearson correlation analysis were performed using Origin 2021 (Northampton, MA, USA).

5. Conclusions

Metabolite analysis and transcriptome sequencing were used to explore the mechanism of flavonoid accumulation in the leaves and roots and the flavonoid metabolic pathways of T. hemsleyanum, as well as to explore the core genes of flavonoid biosynthesis. Given that the active ingredient levels of flavonoids in the leaves were higher than those in the root tubers, we suggest that the leaves should be given more attention in future research on T. hemsleyanum. We found regulatory proteins called photosynthesis-antenna proteins that play a key role in the flavonoid biosynthesis pathway of T. hemsleyanum, which researchers have previously ignored. This study also illustrated that CHI (CL6715.Contig1_All, Unigene19431_All, and CL921.Contig4_All) and UFGT (CL11556.Contig3_All and CL11775.Contig1_All) were critical genes in most of the flavonoid metabolism pathways of T. hemsleyanum. This study provides a theoretical basis for the development and utilization of the leaves and the improvement of the flavonoids of T. hemsleyanum.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28010083/s1, Table S1: Concentration of flavonoid Active ingredients in different parts of T. hemsleyanum.

Author Contributions

Conceptualization, Y.B.; Methodology, Z.L., X.H. and S.L.; Software, L.J.; Validation, Y.B., L.J. and F.G.; Formal analysis, L.J.; Investigation, Y.B.; Resources, Y.B.; Data curation, L.J.; Writing—original draft preparation, L.J.; Writing—review & editing, Y.B., F.G. and S.L.; Visualization, L.J.; Supervision, Y.B.; Project administration, Y.B.; Funding acquisition, Y.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Key Research and Development Program of Zhejiang Province (grant number 2021C02043), Zhejiang Provincial Department of Education project (grant number Y202147328), Forestry special fund of Zhejiang Provincial Forestry Bureau (grant number Lqly2020-03), Forestry Science and Technology Extension Foundation of Central Government Financial Project (grant number 2019TS08), and Zhejiang Agriculture and Forestry University School Development Fund (grant number 2022LFR111).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All sequencing data are available for download at the NCBI Sequence Read Archive (SRA) under the accession number PRJNA826914.

Acknowledgments

J.Y. (YANG Jian, China Academy of Chinese Medical Sciences) and X.F. (FAN Xiao-hui, Institute of Genetics and Developmental Biology) helped determine the concentrations of flavonoid compounds.

Conflicts of Interest

The authors declare that they have no competing interest.

Abbreviations

T. hemsleyanumTetrastigma hemsleyanum Diels et Gilg
PALPhenylalanine ammonia lyase
C4HCinnamate 4-hydroxylase
4CL4-coumarate-coenzyme A
CHIChalcone isomerase
CHSChalcone synthase
F3HFlavanone-3-hydroxylase
F3′HFlavanone-3′-hydroxylase
F3′5′HFlavonoid-3′,5′-hydroxylase
FLSFlavonol synthase
FNSflavone synthase
DFRDihydroflavonol 4-reductase
LDOXLeucoanthocyanidin dioxygenase
3-GlctFlavonoid-3-O-glucosyltransferase
7-GlctFlavonoid-3-O-glucosyltransferase
1,2-Rhat1,2-rhamnosyltransferase
1,6-Rhat1,6-rhamnosyltransferase
UGTsUDP-Glycosyltransferases
OMTO-methyltransferase
UFGTUDP-glycose flavonoid glycosyltransferase
MAPKMitogen-activated protein kinase.
KOGClusters of orthologous groups for eukaryotic complete genomes
NRNon-Redundant Protein Sequence Database
GOGene Ontology
KEGGKyoto Encyclopedia of Genes and Genomes

References

  1. Gong, W.; Liu, T.; Zhou, Z.; Wu, D.; Shu, X.; Xiong, H.; Gong, W.; Liu, T.; Zhou, Z.; Wu, D.; et al. Physicochemical characterizations of starches isolated from Tetrastigma hemsleyanum Diels et Gilg. Int. J. Biol. Macromol. 2021, 183, 1540–1547. [Google Scholar] [CrossRef] [PubMed]
  2. Ji, T.; Ji, W.W.; Wang, J.; Chen, H.J.; Peng, X.; Cheng, K.J.; Qiu, D.; Yang, W.J. A comprehensive review on traditional uses, chemical compositions, pharmacology properties and toxicology of Tetrastigma hemsleyanum. J. Ethnopharmacol. 2021, 264, 1–19. [Google Scholar] [CrossRef] [PubMed]
  3. Song, W.; Liu, Y.X.; Fanady, B.; Han, Y.F.; Xie, L.; Chen, Z.Y.; Yu, K.B.; Peng, X.; Zhang, X.L.; Ge, Z.Y. Ultra-flexible light-permeable organic solar cells for the herbal photosynthetic growth. Nano Energy 2021, 86, 1–11. [Google Scholar] [CrossRef]
  4. Xiang, Q.; Hu, S.; Ligaba-Osena, A.; Yang, J.; Tong, F.; Guo, W. Seasonal variation in transcriptomic profiling of Tetrastigma hemsleyanum fully developed tuberous roots enriches candidate genes in essential metabolic pathways and phytohormone signaling. Front. Plant Sci. 2021, 12, 1–20. [Google Scholar] [CrossRef] [PubMed]
  5. Jiang, Y.L.; Xu, Z.J.; Cao, Y.F.; Wang, F.; Chu, C.; Zhang, C.; Tao, Y.; Wang, P. HPLC fingerprinting-based multivariate analysis of chemical components in Tetrastigma Hemsleyanum Diels et Gilg: Correlation to their antioxidant and neuraminidase inhibition activities. J. Pharm. Biomed. Anal. 2021, 205, 1–11. [Google Scholar] [CrossRef]
  6. Han, M.Q.; Xu, L.Y.; Bai, Y.; Hu, X.T.; Wu, X.Q.; Zheng, B.S. Effect of light quality on photosynthetic characteristics and total flavonoid content in Tetrastigma hemsleyanum. J. Plant Physiol. 2019, 55, 883–890. [Google Scholar]
  7. Li, Z.; Zhou, C.H.; Du, X.B.; Liu, S.Z.; Bai, Y.; Jiang, L.T.; Wang, J.Y.; Pan, J.Y.; Lou, Y.Y.; Zheng, B.S. Current situation and countermeasures on research and development of Tetrastigma hemsleyanum. J. Zhejiang Sci. Tech. 2020, 40, 87–92. [Google Scholar] [CrossRef]
  8. Williams, J.S.; Thomas, M.; Clarke, D.J. The gene stlA encodes a phenylalanine ammonia-lyase that is involved in the production of a stilbene antibiotic in Photorhabdus luminescens TT01. Microbiology 2005, 151, 2543–2550. [Google Scholar] [CrossRef] [Green Version]
  9. Cheng, A.X.; Zhang, X.; Han, X.J.; Zhang, Y.Y.; Gao, S.; Liu, C.J.; Lou, H.X. Identification of chalcone isomerase in the basal land plants reveals an ancient evolution of enzymatic cyclization activity for synthesis of flavonoids. New Phytol. 2018, 217, 909–924. [Google Scholar] [CrossRef] [Green Version]
  10. Liu, W.; Feng, Y.; Yu, S.; Fan, Z.; Li, X.; Li, J.; Yin, H. The Flavonoid Biosynthesis Network in Plants. Int. J. Mol. Sci. 2021, 22, 12824. [Google Scholar] [CrossRef]
  11. Sun, W.; Shen, H.; Xu, H.; Tang, X.; Tang, M.; Ju, Z.; Yi, Y. Chalcone Isomerase a Key Enzyme for Anthocyanin Biosynthesis in Ophiorrhiza japonica. Front. Plant Sci. 2019, 10, 1–12. [Google Scholar] [CrossRef] [PubMed]
  12. Owens, D.K.; Alerding, A.B.; Crosby, K.C.; Bandara, A.B.; Westwood, J.H.; Winkel, B.S. Functional analysis of a predicted flavonol synthase gene family in Arabidopsis. Plant Physiol. 2008, 147, 1046–1061. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Deng, Y.; Li, C.; Li, H.; Lu, S. Identification and characterization of flavonoid biosynthetic enzyme genes in Salvia miltiorrhiza (Lamiaceae). Molecules 2018, 23, 1467. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Zhao, Z.C.; Hu, G.B.; Hu, F.C.; Wang, H.C.; Yang, Z.Y.; Lai, B. The UDP glucose: Flavonoid-3-O-glucosyltransferase (UFGT) gene regulates anthocyanin biosynthesis in litchi (Litchi chinesis Sonn.) during fruit coloration. Mol. Biol. Rep. 2012, 39, 6409–6415. [Google Scholar] [CrossRef] [PubMed]
  15. He, Y.; Chen, H.; Zhao, J.; Yang, Y.; Yang, B.; Feng, L.; Zhang, Y.; Wei, P.; Hou, D.; Zhao, J.; et al. Transcriptome and metabolome analysis to reveal major genes of saikosaponin biosynthesis in Bupleurum chinense. BMC Genom. 2021, 22, 1–14. [Google Scholar] [CrossRef]
  16. Bar-Peled, M.; Fluhr, R.; Gressel, J. Juvenile-specific localization and accumulation of a Rhamnosyltransferase and its bitter flavonoid in foliage, flowers, and young Citrus fruits. Plant Physiol. 1993, 103, 1377–1384. [Google Scholar] [CrossRef] [Green Version]
  17. Gou, J.Y.; Felippes, F.F.; Liu, C.J.; Weigel, D.; Wang, J.W. Negative regulation of anthocyanin biosynthesis in Arabidopsis by a miR156-targeted SPL transcription factor. Plant Cell 2011, 23, 1512–1522. [Google Scholar] [CrossRef] [Green Version]
  18. Koja, E.; Ohata, S.; Maruyama, Y.; Suzuki, H.; Shimosaka, M.; Taguchi, G. Identification and characterization of a rhamnosyltransferase involved in rutin biosynthesis in Fagopyrum esculentum (common buckwheat). Biosci. Biotechnol. Biochem. 2018, 82, 1790–1802. [Google Scholar] [CrossRef] [Green Version]
  19. Peng, X.; Wu, H.; Chen, H.; Zhang, Y.; Qiu, D.; Zhang, Z. Transcriptome profiling reveals candidate flavonol-related genes of Tetrastigma hemsleyanum under cold stress. BMC Genom. 2019, 20, 1–15. [Google Scholar] [CrossRef] [Green Version]
  20. Guo, Q.S. Medicinal Plant Cultivation, 1st ed.; Higher Education Press: Beijing, China, 2009; p. 47. [Google Scholar]
  21. Benavente-Garcia, O.; Castillo, J.; del Rio Conesa, J.A. Changes in neodiosmin levels during the development of Citrus aurantium leaves and fruits. postulation of a neodiosmin biosynthetic pathway. J. Agric. Food Chem. 1993, 41, 1916–1919. [Google Scholar] [CrossRef]
  22. Flora of China Editorial Committee. Flora of China; Science Press: Beijing, China, 1998; p. 120. [Google Scholar]
  23. Zhejiang Food and Medical Products Administration. Zhejiang Provincial Standards of Processing Chinese Crud Drugs, 7th ed.; China Medical Science Press: Beijing, China, 2015; pp. 5–6. [Google Scholar]
  24. Zhang, J.J.; Zhang, Z.F.; Lv, Y.G.; Lu, M.D.; Fang, L.L. Comparison of antioxidant activities of aerial part and underground part of T. hemsleyanum from Zhejiang Province. J. Zhejiang Agric. Sci. 2021, 62, 709–711. [Google Scholar]
  25. Su, Q.; Sun, Z.; Liu, Y.; Lei, J.; Zhu, W.; Nanyan, L. Physiological and comparative transcriptome analysis of the response and adaptation mechanism of the photosynthetic function of mulberry (Morus alba L.) leaves to flooding stress. Plant Signal. 2022, 17, e2094619. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, D.; Gao, Y.; Sun, S.; Li, L.; Wang, K. Expression profiles and characteristics of apple lncRNAs in roots, phloem, leaves, flowers, and fruit. Int. J. Mol. Sci 2022, 23, 5931. [Google Scholar] [CrossRef] [PubMed]
  27. Chen, J.; Zeng, H.; Zhang, X. Integrative transcriptomic and metabolomic analysis of D-leaf of seven pineapple varieties differing in N-P-K% contents. BMC Plant Biol. 2021, 21, 1–19. [Google Scholar] [CrossRef]
  28. Wang, F.; Wu, Y.; Wu, W.; Huang, Y.; Zhu, C.; Zhang, R.; Chen, J.; Zeng, J. Integrative analysis of metabolome and transcriptome profiles provides insight into the fruit pericarp pigmentation disorder caused by ’Candidatus Liberibacter asiaticus’ infection. BMC Plant Biol. 2021, 21, 1–16. [Google Scholar] [CrossRef]
  29. Hong, Y.; Li, M.; Dai, S. iTRAQ-Based protein profiling provides insights into the mechanism of light-induced anthocyanin biosynthesis in Chrysanthemum (Chrysanthemum × morifolium). Genes 2019, 10, 1024. [Google Scholar] [CrossRef] [Green Version]
  30. Bai, Y.; Gu, Y.W.; Liu, Z.S.; Jiang, L.T.; Han, M.Q.; Geng, D.J. Flavonoids metabolism and physiological response to ultraviolet treatments in Tetrastigma hemsleyanum Diels et Gilg. Front. Plant Sci. 2022, 13, 926197. [Google Scholar] [CrossRef]
  31. Ma, H.; Chen, J.; Zhang, Z.; Ma, L.; Yang, Z.; Zhang, Q.; Li, X.; Xiao, J.; Wang, S. MAPK kinase 10.2 promotes disease resistance and drought tolerance by activating different MAPKs in rice. Plant J. 2017, 92, 557–570. [Google Scholar] [CrossRef] [Green Version]
  32. Ryu, T.K.; Roh, E.; Shin, H.S.; Kim, J.E. Inhibitory effect of lotusine on solar UV-induced matrix metalloproteinase-1 expression. Plants 2022, 11, 773. [Google Scholar] [CrossRef]
  33. Zhou, H.; Ren, S.; Han, Y.; Zhang, Q.; Qin, L.; Xing, Y. Identification and analysis of mitogen-activated protein kinase (MAPK) cascades in Fragaria vesca. Int J Mol Sci 2017, 18, 1766. [Google Scholar] [CrossRef] [Green Version]
  34. Yao, Y.; Zhao, H.; Sun, L.; Wu, W.; Li, C.; Wu, Q. Genome-wide identification of MAPK gene family members in Fagopyrum tataricum and their expression during development and stress responses. BMC Genom. 2022, 23, 1–18. [Google Scholar] [CrossRef] [PubMed]
  35. Li, X.; Zhang, L.P.; Zhang, L.; Yan, P.; Ahammed, G.J.; Han, W.Y. Methyl salicylate enhances flavonoid biosynthesis in tea leaves by stimulating the Phenylpropanoid Pathway. Molecules 2019, 24, 362. [Google Scholar] [CrossRef] [PubMed]
  36. Liang, D.; Shen, Y.; Ni, Z.; Wang, Q.; Lei, Z.; Xu, N.; Deng, Q.; Lin, L.; Wang, J.; Lv, X.; et al. Exogenous melatonin application delays senescence of kiwifruit leaves by regulating the antioxidant capacity and biosynthesis of flavonoids. Front. Plant Sci. 2018, 9, 1–14. [Google Scholar] [CrossRef] [PubMed]
  37. Isah, T. Stress and defense responses in plant secondary metabolites production. Biol. Res. 2019, 52, 1–25. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Park, N.I.; Xu, H.; Li, X.; Kim, S.J.; Park, S.U. Enhancement of flavone levels through overexpression of chalcone isomerase in hairy root cultures of Scutellaria baicalensis. Funct. Integr. Genom. 2011, 11, 491–496. [Google Scholar] [CrossRef] [PubMed]
  39. Chen, Z.; Liu, G.; Tang, N.; Li, Z. Transcriptome analysis reveals molecular signatures of luteoloside accumulation in senescing leaves of Lonicera macranthoides. Int. J. Mol. Sci. 2018, 19, 1012. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Yu, X.; Zhu, Y.; Fan, J.; Wang, D.; Gong, X.; Ouyang, Z. Accumulation of flavonoid glycosides and UFGT gene expression in Mulberry Leaves (Morus alba L.) before and after frost. Chem Biodivers 2017, 14, e1600496. [Google Scholar] [CrossRef]
  41. Audic, S.; Claverie, J.M. The significance of digital gene expression profiles. Genome Res. 1997, 7, 986–995. [Google Scholar] [CrossRef]
  42. Liu, Y.; Pan, J.; Ni, S.; Xing, B.; Cheng, K.; Peng, X. Transcriptome and Metabonomics Combined Analysis Revealed the Defense Mechanism Involved in Hydrogen-Rich Water-Regulated Cold Stress Response of Tetrastigma hemsleyanum. Front. Plant Sci. 2022, 13, 1–17. [Google Scholar] [CrossRef]
Figure 1. Metabolomic sequencing analysis of flavonoids. (A) PCA score plot of leaves and root tubers of T. hemsleyanum; (B) Cluster contour heatmap of 22 metabolites; (C) Concentration of flavonoid active ingredients in different parts of T. hemsleyanum; (D) Flavonoids only existing in leaves. (Concentration of flavonoids in leaves or root tubers more than 0.1 μg·g−1; different symbol means, * represents p ≤ 0.05, ** represents p ≤ 0.01.)
Figure 1. Metabolomic sequencing analysis of flavonoids. (A) PCA score plot of leaves and root tubers of T. hemsleyanum; (B) Cluster contour heatmap of 22 metabolites; (C) Concentration of flavonoid active ingredients in different parts of T. hemsleyanum; (D) Flavonoids only existing in leaves. (Concentration of flavonoids in leaves or root tubers more than 0.1 μg·g−1; different symbol means, * represents p ≤ 0.05, ** represents p ≤ 0.01.)
Molecules 28 00083 g001
Figure 2. Venn diagrams of unigene annotations. (A) Venn diagram of KOG, NR, SwissProt, KEGG, and GO gene database annotation from leaves of T. hemsleyanum; (B) Venn diagram of NT, NR, SwissProt, KEGG, and GO gene database annotation from root tubers.
Figure 2. Venn diagrams of unigene annotations. (A) Venn diagram of KOG, NR, SwissProt, KEGG, and GO gene database annotation from leaves of T. hemsleyanum; (B) Venn diagram of NT, NR, SwissProt, KEGG, and GO gene database annotation from root tubers.
Molecules 28 00083 g002
Figure 3. Classification of the assembled unigenes. (A) GO classification of unigenes; (B) KEGG classification of unigenes.
Figure 3. Classification of the assembled unigenes. (A) GO classification of unigenes; (B) KEGG classification of unigenes.
Molecules 28 00083 g003
Figure 4. Results of DEGs in different parts of T. hemsleyanum. (A) Volcano plots for DEGs from leaves vs. roots; (B) KEGG classification of DEGs from leaves vs. roots; (C) Scatter plot of differentially expressed T. hemsleyanum genes in the 20 most enriched KEGG pathways (the y-axis shows KEGG pathways, and the x-axis shows rich factor; the size of a point represents the number of DEGs, and the color of a point represents its q-value); (D) Network of differentially expressed T. hemsleyanum genes enriched in the photosynthesis-antenna proteins.
Figure 4. Results of DEGs in different parts of T. hemsleyanum. (A) Volcano plots for DEGs from leaves vs. roots; (B) KEGG classification of DEGs from leaves vs. roots; (C) Scatter plot of differentially expressed T. hemsleyanum genes in the 20 most enriched KEGG pathways (the y-axis shows KEGG pathways, and the x-axis shows rich factor; the size of a point represents the number of DEGs, and the color of a point represents its q-value); (D) Network of differentially expressed T. hemsleyanum genes enriched in the photosynthesis-antenna proteins.
Molecules 28 00083 g004
Figure 5. Classification map of flavonoids.
Figure 5. Classification map of flavonoids.
Molecules 28 00083 g005
Figure 6. Biosynthetic pathway of flavonoids [16,17,18,19,20,21]. This pathway contains the major flavonoid metabolites and the expression of synthetases. The ingredients in the dotted box are the flavonoids found in the roots and leaves; the enzyme near the red dashed arrow was not involved in our experiment. PAL, phenylalanine ammonia lyase; C4H, Cinnamate 4-hydroxylase; 4CL, 4-coumarate-coenzyme A; CHI, chalcone isomerase; CHS, chalcone synthase; F3H, Flavanone-3-hydroxylase; F3′H, Flavanone-3′-hydroxylase; F3′5′H, Flavonoid-3′,5′-hydroxylase; FLS, Flavonol synthase; FNS, Flavone synthase; DFR, Dihydroflavonol 4-reductase; LDOX, Leucoanthocyanidin dioxygenase; 3-Glct, Flavonoid-3-O-glucosyltransferase; 7-Glct, Flavonoid-3-O-glucosyltransferase; 1,2-Rhat, 1,2-rhamnosyltransferase; 1,6-Rhat, 1,6-rhamnosyltransferase; UGTs, UDP-Glycosyltransferases; OMT, O-methyltransferase; UFGT, UDP-glycose flavonoid glycosyltransferase. The heatmaps represent log10 (FPKM) values of these enzymes in leaves and root tubers.
Figure 6. Biosynthetic pathway of flavonoids [16,17,18,19,20,21]. This pathway contains the major flavonoid metabolites and the expression of synthetases. The ingredients in the dotted box are the flavonoids found in the roots and leaves; the enzyme near the red dashed arrow was not involved in our experiment. PAL, phenylalanine ammonia lyase; C4H, Cinnamate 4-hydroxylase; 4CL, 4-coumarate-coenzyme A; CHI, chalcone isomerase; CHS, chalcone synthase; F3H, Flavanone-3-hydroxylase; F3′H, Flavanone-3′-hydroxylase; F3′5′H, Flavonoid-3′,5′-hydroxylase; FLS, Flavonol synthase; FNS, Flavone synthase; DFR, Dihydroflavonol 4-reductase; LDOX, Leucoanthocyanidin dioxygenase; 3-Glct, Flavonoid-3-O-glucosyltransferase; 7-Glct, Flavonoid-3-O-glucosyltransferase; 1,2-Rhat, 1,2-rhamnosyltransferase; 1,6-Rhat, 1,6-rhamnosyltransferase; UGTs, UDP-Glycosyltransferases; OMT, O-methyltransferase; UFGT, UDP-glycose flavonoid glycosyltransferase. The heatmaps represent log10 (FPKM) values of these enzymes in leaves and root tubers.
Molecules 28 00083 g006
Figure 7. Correlations between flavonoid monomers and gene expression of related enzymes.
Figure 7. Correlations between flavonoid monomers and gene expression of related enzymes.
Molecules 28 00083 g007
Figure 8. Validation of the differentially expressed genes. (A) Relative expressions of flavonoid-synthesis-related genes; (B) Correlation between RNA-seq and relative gene expression (the x-axis represents the log2 fold change value of the RNA-seq results, and the y-axis represents the log2 fold change value of the relative expression of related genes; all data were collected from three biological replicates and three technical replicates).
Figure 8. Validation of the differentially expressed genes. (A) Relative expressions of flavonoid-synthesis-related genes; (B) Correlation between RNA-seq and relative gene expression (the x-axis represents the log2 fold change value of the RNA-seq results, and the y-axis represents the log2 fold change value of the relative expression of related genes; all data were collected from three biological replicates and three technical replicates).
Molecules 28 00083 g008
Figure 9. Expressions of flavonoid synthesis key genes and total flavonoid content of Zhejiang T. hemsleyanum from different locations. (A) Expression of flavonoid synthesis key genes of Zhejiang T. hemsleyanum from different locations (CHI 1, CL6715.Contig1_All; CHI 2, Unigene19431_All; CHI 3, CL921.Contig4_All; UFGT 1, CL11556.Contig3_All; UFGT 2, CL11775.Contig1_All.). (B) The concentration of total flavonoids of Zhejiang T. hemsleyanum from different locations (different symbols mean a significant difference, ** represents p ≤ 0.01, *** represents p ≤ 0.001).
Figure 9. Expressions of flavonoid synthesis key genes and total flavonoid content of Zhejiang T. hemsleyanum from different locations. (A) Expression of flavonoid synthesis key genes of Zhejiang T. hemsleyanum from different locations (CHI 1, CL6715.Contig1_All; CHI 2, Unigene19431_All; CHI 3, CL921.Contig4_All; UFGT 1, CL11556.Contig3_All; UFGT 2, CL11775.Contig1_All.). (B) The concentration of total flavonoids of Zhejiang T. hemsleyanum from different locations (different symbols mean a significant difference, ** represents p ≤ 0.01, *** represents p ≤ 0.001).
Molecules 28 00083 g009
Table 1. Length distribution of assembled unigenes.
Table 1. Length distribution of assembled unigenes.
TermAll > 300 bp≥500 bp≥1000 bpN50Total LengthMax LengthMin LengthAverage Length
Y-Unigene53,03929,97314,690127547,471,67715,654301895.03
G-Unigene60,95333,60716,023123653,338,27215,654301875.07
Table 2. Linear relationship investigation results.
Table 2. Linear relationship investigation results.
CompoundRegression EquationrLinear Range
VitexinY = 8,542,235.9 X + 18,454,664.60.998720.000–2000.000
OrientinY = 9,071,337.1 X + 33,093.60.99961.000–100.000
AfzelinY = 9,658,807.6 X + 92,259.20.99930.001–0.050
RutinY = 4,950,860.8 X + 231,551.20.99940.010–0.500
IsoquercitrinY = 3,356,305.1 X + 593,847.80.99890.020–1.000
LuteolinY = 9,109,434.3 X + 35,359.90.99980.020–1.000
LuteolosideY = 3,168,802.6 X + 166,068.30.99940.100–5.000
IsoorientinY = 4,845,698.3 X + 4,222,807.80.99871.000–100.000
NaringinY = 3,972,039.5 X + 64,183.60.99760.010–0.500
Isorhamnetin-3-O-GlucosideY = 8,988,517.2 X + 5,992,231.10.99820.020–1.000
IsorhamnetinY = 1,476,424.2 X + 2,111,521.60.99860.001–0.050
AstragalineY = 2,210,508.5 X + 2,507,321.40.99860.010–0.500
EriodictyolY = 691,932.1 X + 84659.50.99910.010–0.500
TaxifolinY = 3,575,802.9 X + 284,658.70.99920.001–0.050
IsovitexinY = 8,542,235.9 X + 18,454,664.60.998720.000–2000.000
Calycosin-7-O-glucosideY = 22,042.7 X + 2463.80.99860.020–1.000
HispidulinY = 207,676.0 X + 105,248.90.99930.020–1.000
NaringeninY = 12,452,147.9 X +74,631.90.99920.001–0.050
ApigeninY = 2,885,631.5 X + 238,287.70.99890.001–0.050
QuercitrinY = 12,206,204.2 X + 1,425,428.80.99920.001–0.050
Apigenin-7-O-GlycosideY = 1,806,515.0 X + 16,658.00.99960.010–0.500
Table 3. Primer sequences of flavonoid-related genes in different parts of T. hemsleyanum.
Table 3. Primer sequences of flavonoid-related genes in different parts of T. hemsleyanum.
GeneGene IDTm (°C)Primer F (5′-3′)Primer R (5′-3′)
MDH\60TGTTGCTACGACTGATGTCCTGAGACTTGTAGATGGAA
PALCL2379.Contig9_All56.02ACCAACCATGTCCAAAGTGCCGCCACAAGGTATGTGGAAG
CL2379.Contig4_All56.06TGGTGACCTCACCTTCTCACGAAGCCCATCCATTCCGAAC
CHSCL1413.Contig5_All55.9GGAACTGTCCTTCGAACTGCGTCCGCGGAATGTAACAACA
CL1413.Contig12_All55.99TAGCACGTTGAGCGTTTCTGAGAGTTGGTGGCATAAGGCT
CL1413.Contig4_All55.71CCAACTTGTCTCAGCAGCTCAATCAAAGTGGGCACAGTGG
CL1413.Contig11_All55.89TGTACCATCAAGGGTGCCATGGGCTTGGCCAACTAAAGAG
CHICL6715.Contig1_All56.01GTGCAGGGTGTGAAGTTTGTTTTGGCTTCCTTCCAACAGC
Unigene19431_All56.47ACGCCATGGATAGAGAGCAACCCATGGTTGAGGATTCGGA
CL921.Contig4_All55.96ACTGAACACCATCCACGACTGACCAACGAATGCCTCGAAA
FNSUnigene21915_All55.99AAAGACGAACTCTCCACCGTGATGTTCGCAACTCCGTTCA
UFGTCL11556.Contig3_All56.05TTGACTTGCCTGAGTGTCCTAACTCGGATGCTGAGTTGGA
CL11775.Contig1_All55.90CAACGGCGGAATGAGCTAAATCTTGTGGTCCCTTCGTCAA
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bai, Y.; Jiang, L.; Li, Z.; Liu, S.; Hu, X.; Gao, F. Flavonoid Metabolism in Tetrastigma hemsleyanum Diels et Gilg Based on Metabolome Analysis and Transcriptome Sequencing. Molecules 2023, 28, 83. https://doi.org/10.3390/molecules28010083

AMA Style

Bai Y, Jiang L, Li Z, Liu S, Hu X, Gao F. Flavonoid Metabolism in Tetrastigma hemsleyanum Diels et Gilg Based on Metabolome Analysis and Transcriptome Sequencing. Molecules. 2023; 28(1):83. https://doi.org/10.3390/molecules28010083

Chicago/Turabian Style

Bai, Yan, Lingtai Jiang, Zhe Li, Shouzan Liu, Xiaotian Hu, and Fei Gao. 2023. "Flavonoid Metabolism in Tetrastigma hemsleyanum Diels et Gilg Based on Metabolome Analysis and Transcriptome Sequencing" Molecules 28, no. 1: 83. https://doi.org/10.3390/molecules28010083

APA Style

Bai, Y., Jiang, L., Li, Z., Liu, S., Hu, X., & Gao, F. (2023). Flavonoid Metabolism in Tetrastigma hemsleyanum Diels et Gilg Based on Metabolome Analysis and Transcriptome Sequencing. Molecules, 28(1), 83. https://doi.org/10.3390/molecules28010083

Article Metrics

Back to TopTop