Next Article in Journal
Effect of n-3 (Omega-3) Polyunsaturated Fatty Acid Supplementation on Metabolic and Inflammatory Biomarkers and Body Weight in Patients with Type 2 Diabetes Mellitus: A Systematic Review and Meta-Analysis of RCTs
Previous Article in Journal
Cord Blood Metabolome and BMI Trajectory from Birth to Adolescence: A Prospective Birth Cohort Study on Early Life Biomarkers of Persistent Obesity
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Metabolites
Volume 11
Issue 11
10.3390/metabo11110741
metabolites-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of Flavonoids and Transcripts Involved in Their Biosynthesis in Different Organs of Cissus rotundifolia Lam

by
Duncan Kiragu Gichuki
1,2,3,4,
Qingyun Li
1,2,4,
Yujun Hou
1,2,4,
Yuanshuang Liu
1,2,4,
Mengxue Ma
1,2,
Huimin Zhou
1,2,4,
Chen Xu
1,2,
Zhenfei Zhu
1,2,4,
Lina Wang
1,2,
Fredrick Mutie Musila
5,
Qingfeng Wang
1,2,3 and
Haiping Xin
1,2,3,*
1
Core Botanical Gardens/Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan 430074, China
2
Key Laboratory of Plant Germplasm Enhancement and Specialty Agriculture, Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan 430074, China
3
Sino-Africa Joint Research Center, Chinese Academy of Sciences, Wuhan 430074, China
4
University of Chinese Academy of Sciences, Beijing 100049, China
5
School of Biological and Life Sciences, Technical University of Kenya, Nairobi 52428-00200, Kenya
*
Author to whom correspondence should be addressed.
Metabolites 2021, 11(11), 741; https://0-doi-org.brum.beds.ac.uk/10.3390/metabo11110741
Submission received: 7 October 2021 / Revised: 25 October 2021 / Accepted: 26 October 2021 / Published: 28 October 2021
(This article belongs to the Section Plant Metabolism)
Browse Figures
Versions Notes

Abstract

:
Cissus rotundifolia Lam. is used as a medicinal herb and vegetable. Flavonoids are the major components for the therapeutic effects. However, flavonoids constituents and expression profiles of related genes in C. rotundifolia organs are unknown. Colorimetric assay showed the highest flavonoid concentration in roots compared to the stem and leaf. Widely target-based metabolome analysis allowed tentative identification of 199 compounds in three organs. Flavonols and flavones were the dominant flavonoids subclasses. Among the metabolites, 171 were common in the three organs. Unique accumulation profile was observed in the root while the stem and leaf exhibited relatively similar patterns. In the root, six unique compounds (jaceosidin, licoagrochalcone D, 8-prenylkaempferol, hesperetin 7-O-(6″malonyl) glucoside, aureusidin, apigenin-4′-O-rhamnoside) that are used for medicinal purposes were detected. In total, 18,427 expressed genes were identified from transcriptome of the three organs covering about 60% of annotated genes in C. rotundifolia genome. Fourteen gene families, including 52 members involved in the main pathway of flavonoids biosynthesis, were identified. Their expression could be found in at least one organ. Most of the genes were highly expressed in roots compared to other organs, coinciding with the metabolites profile. The findings provide fundamental data for exploration of metabolites biosynthesis in C. rotundifolia and diversification of parts used for medicinal purposes.

1. Introduction

Cissus rotundifolia Lam., a species of Cissus genus within the grape family, native to Africa, is widely used as a vegetable and medicinal herb [1,2,3,4]. As the leaves are rich in proteins, fatty acids, crude fibers, and minerals [5], C. rotundifolia is regarded as a probable source of healthy food. Due to its anti-diabetic and anti-parasitic properties, preparations derived from, leaves, stems, or a whole plant of Cissus rotundifolia Lam. are used as conventional remedies for diabetes, obesity, malaria, allergies, and bacterial infections [6,7,8].
To further understand its medicinal active ingredients and toxicity, the chemical constituents of C. rotundifolia, especially phenolic components, were extracted and studied [9,10]. Said et al. [10] found that methanolic extract of C. rotundifolia had significant central and peripheral analgesic effects, and their inhibitory effect on paw edema was better than indomethacin, while the total flavonoids in the extract accounted for 42.5% of the total phenolic content. Flavonoids possess a wide range of pharmaceutical properties, such as anti-tumor, antioxidant, anti-inflammatory, and anti-viral properties, as well as inhibitory properties against blood clots [11,12,13,14,15,16]. In medicinal plants, flavonoids subclasses, such as flavones, have been isolated from variable tissues and their medicinal traits have been reported. In Dracaena cambodiana and Dracaena cochinchinensis, which are used in Chinese traditional medicine, flavonoids and related metabolites have been isolated, exhibiting antibacterial and growth inhibitory properties [17,18,19]. Additionally, a flavanone isolated from Bauhinia variegata Linn. was reported to be effective against human cancer cell lines [20]. Recently, Alqahtani [21] pioneered and characterized three infrequent C-glycosyl flavones in C. rotundifolia and determined that 1-O-(4-coumaroyl)-β-D-glucopyranose was chiefly responsible for the glucose uptake stimulation. Therefore, it can be speculated that C. rotundifolia flavonoids are among the metabolites leading to their use in the pharmaceutical field. However, to our knowledge, there is limited information available about flavonoids and their biosynthesis process in C. rotundifolia.
As a consequence of the wide range of importance in biological systems and the medical field, the flavonoids biosynthesis pathway in plants has been widely explored [22]. Although there are several important features regarding modifications or decorative reactions of flavonoids still unrevealed, the main trunk biosynthesis pathway is, by and large, conserved across plant species [23]. Flavonoids are synthesized from phenylanine in the phenylpropanoid pathway [24,25,26]. This pathway and flavonoids diversification are regulated by different transcription factors, such as MYBs, bHLH (basic helix-loop-helix), WD40 proteins, and WRKYs [27,28], through regulation of expression for genes involved in this metabolic pathway [29], including phenylalanine ammonia-lyase (PAL), cinnamate-4-hydroxylase (C4H), 4-coumarate: CoA ligase (4CL), chalcone synthase (CHS), chalcone isomerase (CHI), dihydroflavonol 4-reductase (DFR), flavonoid 3′-hydroxylase (F3′H), isoflavone synthase (IFS), flavonoid 3′,5′-hydroxylase (F3′5′H), flavonol synthase (FLS), anthocyanidin synthase (ANS), anthocyanidin reductase (ANR), and UDP-glucose: flavonoid 3-O-glucosyltransferase (UFGT). Among these genes, the FLS gene family (FLS1 and FLS2), encoding for key branching enzymes, were characterized in O. caudatum [30]. Their functionality was reported in activating the conversion of dihydroflavonols to flavonols, as well as in the hydroxylation of flavanones to dihydroflavonols.
The development and integration of contemporary -omic technologies, including proteomics, transcriptomics, and metabolomics have enhanced understanding of metabolites biosynthesis mechanism at the molecular level [31,32]. Metabolomics represents the physiological events at the cellular level through the exploration of cellular metabolites and has been applied in the detection of low molecular weight metabolites, such as flavonoids in model plants, crops, and fruits [33,34,35,36,37]. However, associating the metabolome to the genome is challenging, even in model plants with plentiful genomic resources [38].
Flavonoids are a diverse group of plant secondary metabolites and have been widely characterized. However, molecular characterization of flavonoids in the Cissus genus is limited despite the wide range of medicinal applications of its members. For instance, C. quandrangularis is widely used in the treatment of bone fractures and body weight management, C. hypoglauca for sore throats, C. assamica to neutralize snake venoms, C. rubiginosa for anti-diarrhea, and C. rotundifolia for blood sugar management [4]. In the current work, through the integration of metabolomics and transcriptomic analysis, elucidation of flavonoids components, associated variations in accumulation and expression of corresponding genes, were explored in three organs (root, leaf, and stem) of C. rotundifolia. The expression patterns for flavonoids-related genes were also examined in the organs. This study aims to reveal metabolic variations across organs of C. rotundifolia providing a valuable foundation for further exploration of the species and other members in the genus in modern pharmaceuticals.

2. Results

2.1. Total Flavonoids Content Estimation

To determine the accumulation of flavonoids across C. rotundifolia organs, the concentration of total flavonoids in leaf, stem, and root was measured by colorimetric methods. The results (Figure 1) showed that root had the highest concentration of total flavonoids, up to 88.11 mg (RE)/g (DW), followed by the stem (24.82 mg (RE)/g (DW) and leaf (15.39 mg (RE)/g (DW).

2.2. Flavonoids Profiling in C. rotundifolia

A total of 199 compounds were tentatively identified in the 3 organs, including 50 flavonols, 42 flavones, 32 flavone c-glycosides, 17 anthocyanidins, 17 flavanones, 11 flavanols, 8 flavanonols, 6 flavanone c-glycosides, 5 isoflavones, 5 chalcones, 1 aurone, and 5 phenolic acids (Figure 2A). Organ-specific identification of the metabolites was also carried out identifying metabolites that were shared across organs as well as those unique to specific organs. In total, 177 metabolites were shared across the three organs with the leaf having no unique metabolite (Figure 2B). In stem, the two unique metabolites were tricetin, and homoeriodictyol. On the other hand, six root-specific compounds were detected, including jaceosidin, licoagrochalcone D, 8-prenylkaempferol, aureusidin, apigenin-4′-o-rhamnoside, and hesperetin-7-o-(6″-malonyl)glucoside. A comparison of metabolites abundance across the organs revealed a unique pattern, in which metabolites with higher abundance in root had variable abundance in stem and leaf tissues. On the other hand, stem and leaf seemingly shared similar patterns in metabolites, abundance (Figure 3A). Details of the detected metabolites including the ion abundance and retention time, are highlighted (Table 1 and Supplementary file S1).
Using Pearson’s correlation coefficient, the repeatability among the intragroup samples was evaluated (Supplementary Figure S1). The rate of contribution of the first two primary principal components was evaluated in PCA. As expected, the three tissues were separated into distinct clusters. The results suggest that standard reproducibility for the tissues and methods used was acceptable allowing further qualitative and quantitative analyses.

2.3. Identification of Differentially Accumulated Metabolites (DAMs)

The metabolites with fold change ≥2 or fold change ≤0.5 between different organs were selected. In instances where biological duplication was detected in the sample grouping, metabolites with variable importance in the projection (VIP) ≥1 were selected. VIP value indicates the influence intensity of the difference between the corresponding metabolites in the classification of samples in each group in the model. Generally speaking, the metabolites with VIP ≥ −1 were significantly different. The comparison was carried out across the three C. rotundifolia organs.
In total 72, 83, and 86 differentially accumulated metabolites were identified for leaf–vs–stem (L–vs–S), root–vs–leaf (R–vs–L), and root–vs–stem (R–vs–S) respectively (Supplementary Figure S2). In L–vs–S DAMs, flavones (18) were the majority followed by flavone c–glycosides (16), and flavonols (11) subclasses. In R-vs-L, flavone c–glycosides (22) followed by 21 flavonols were detected and there were 15 flavones among others. Flavonols (24) and flavone c–glycoside (20) were the subclasses with large representation in root–vs–stem DAMs (Figure 3B).

2.4. Transcript Sequencing and Mapping

Using the Illumina sequencing platform, mRNA libraries were generated for leaf, stem, and root tissues of C. rotundifolia. The quality of the reads was checked, eliminating low-quality reads and the adapter sequences. A summary of the sequencing statistics has been highlighted in supplementary Table S1. The Q-20 values averaged about 96%, while the GC content ranged between 44% and 47%. Moreover, read mapping ranged between 84.9–86.6% for leaf, 70–70.5% for root, and 86.5–86.9% for stem (Supplementary Table S2). A total of 18,427 expressed genes were identified from the transcriptome data of the three tissues, which covered about 60% of the total annotated genes in the C. rotundifolia genome (from our research group; the Cissus rotundifolia genome project was deposited at the National Genomics Data Center (https://ngdc.cncb.ac.cn/; to be released on 1 January 2022, (Accessed on: 25 October 2021)) under the BioProject number PRJCA005006).

2.5. Functional Annotation of Identified Genes

Using the eggNOG platform, functional annotation of identified genes was carried out. A total of 17,282 (93.8%) C. rotundifolia transcriptome-expressed genes were annotated to COG, KEGG, and GO functional categories. In total, 48.9% (9010 genes) of all genes were characterized into GO terms in the three main ontologies, with 17,109 GO functional terms (Supplementary Figure S3a). The main classes in the major ontologies are listed in supplementary file S2.
Based on the clusters of orthologous groups of proteins database (COG), a total of 17,282 C. rotundifolia genes were annotated into 25 COG functional groups. Among the 25 COG classes, most of the genes were classified under function unknown (S) (Supplementary Figure S3b). As was suggested by Galperin et al. [39], the category for uncharacterized proteins (S) is an important indicator of the progress in the integration of experimental characterization and digital profiling of protein families. To identify the active biochemical pathways in C. rotundifolia and improve the understanding of biological functions and gene interactions, KEGG analysis was carried out. A total of 7452 genes were assigned to 233 KEGG pathways, including metabolism (A091000), which was the dominant pathway (Supplementary file S3).

2.6. Candidate Genes Involved in Flavonoids Biosynthesis

An analysis of C. rotundifolia transcriptome revealed multiple transcripts that have been identified to encode enzymes involved in flavonoids metabolism. Most of these genes have higher expression in root compared to leaf and stem (Figure 4B), which is consistent with total flavonoids content. A brief schematic chart was developed (Figure 4A) using the KEGG database and modifying the flavonoids biosynthesis pathway previously described [40]. Generally, the initial process of flavonoids metabolism involves the conversion of phenylalanine through coumaroyl-CoA to chalcones/naringenin by the activation of several enzymes that include phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H), 4-coumarate CoA ligase (4CL), and chalcone synthase (CHS) through the phenylpropanoid pathway. Through the action of shikimate o-hydroxycinnamoyltransferase (HCT), among other enzymes, coumaroyl-CoA can also be converted to eriodictyol. Another enzyme, chalcone isomerase (CHI), catalyzes the cyclization of chalcone naringenin to naringenins or flavanones, while flavanone 3-hydroxylase (F3H), flavonoid 3′-hydroxylase (F3′H), and flavonoid 3′5′-hydroxylase (F3′5′H) are involved in the hydroxylation of flavanones to various flavonol classes. Additionally, flavonol synthase (FLS) also plays a key role in the conversion of flavanones to the respective flavonols. The number of key genes encoding for the identified enzymes has been highlighted in Table 2 and compared with those identified in D. cambodiana. Similar to D. cambodiana [40], no transcripts encoding for flavone synthase (FNS) were detected from our RNA sequencing analysis.
In addition to the phenylalanine branch of flavonoids biosynthesis, an anthocyanin branch through the action of dihydroflavonol 4-reductase (DFR) catalyzes the conversion of both dihydroquercetins and dihydromyricetin to either leucocyanidins or leucodelphinidins. In our study, enzymes that are involved in this branch of biosynthesis were identified and include dihydroflavanol 4-reductase (DFR), anthocyanidin reductase (ANR), and leucoanthocyanidin dioxygenase (ANS/LDOX).

2.7. Validation of RNA-Seq Data Using qPCR

The results from RNA-seq were validated using qPCR to determine the correlation in expression from the two techniques. In this study, 12 genes involved in flavonoid biosynthesis were selected for validation (Supplementary Table S3). The primers used were designed using NCBI Primer-BLAST (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/tools/primer-blast/ (Accessed on: 25 October 2021)). The primer details have been highlighted in supplementary Table S4. The qPCR analysis was performed using a 7500 Fast Real-Time PCR system (Applied Biosystems, Waltham, MA, USA) in a total of 10 µL reaction volume. The Ct value was determined using the instrument’s software, and the relative quantification of gene expression was monitored after normalization using actin (CRGY0219113) as the internal standard. A comparison of the genes expression patterns obtained from RNA-seq and qPCR was carried out, and the reliability of the RNA-seq was confirmed by the consistency in the expression trends detected by the two data sets (Supplementary Figure S4).

2.8. Candidate Transcription Factors Related to Flavonoids Biosynthesis

Several enzymatic and regulatory proteins related to flavonoids biosynthesis were identified and characterized. They included MYB proteins that were involved in the earlier steps in the pathway, regulating flavonol biosynthesis. In the late stages of the pathway leading to the production of anthocyanins and proanthocyanins, a complex of MYB, bHLH, and WD40 proteins (MYB-bHLH-WD40) activated the related genes. Our analysis identified 262 MYB, 70 bHLH, and 169 WD40 protein transcripts that were expressed in the three organs. Their expression profiles were examined, and distinct patterns were observed. Higher expression profiles were observed in roots when compared with both stem and root tissues for the three protein families (Supplementary Figure S5). The expression profiles for the transcription factors exhibited a pattern similar to that of flavonoids and the related genes in the three C. rotundifolia organs.

2.9. Differential Gene Expression between Tissues of C. rotundifolia

To understand the differences in expression of genes across three organs from C. rotundifolia, we used the FPKM (fragments per kb per million fragments) method to digitally profile the expression of genes between leaf, stem, and root [41]. The profiling was carried out for root–vs–stem (R–vs–S), leaf–vs–stem (L–vs–S), and root–vs–leaf (R–vs–L). Differentially expressed genes (DEGs) were reported as those with more than at least a two-fold- change between organs and with a p-value ≤ 0.05. Most of the DEGs were identified between root and leaf while a comparison between leaf and stem indicated the least number of DEGs (Figure 5).

2.10. GO Enrichment and KEGG Pathway Analysis for Differentially Expressed Genes

The identified DEGs were annotated into the three main GO ontologies. The GO annotation details for the DEGs are highlighted in supplementary Table S5. KEGG analysis for DEGs identified between organs was carried out. Many of the DEGs were associated with metabolic pathways, including secondary metabolites (B09110), which are an important part of medicinal plants (Supplementary Figure S6).
Functional characterization of DEGs into GO terms and KEGG pathways improved our inference about gene expression patterns across the tissues. Basic functions were the dominant GO terms identified across the board when gene expression levels were compared between organs. Organ-specific KEGG enrichment was observed. DEGs in leaf were annotated into photosynthesis-related pathways, such as photosynthesis and phenylanine metabolism, and carotenoid biosynthesis, among others. DEGs in root were highly enriched in secondary metabolite pathways, including flavonoids (00941), isoflavonoids (00943), and phenylpropanoid (00940), suggesting more flavonoids concentration in the root.

3. Discussion

Although C. rotundifolia is widely used as traditional medicine, and its extract has also been shown to have anti-diabetic and antioxidant properties, there are few studies on the specific active ingredients of its extract. The extract of C. rotundifolia leaves using different solvents, including methanol, acetone, and ethanol, have shown high phenolic content and antioxidant activity [9]. By analyzing the correlation of antioxidant activities and total phenols, Al-Mamary [42] postulated that the antioxidant activities of C. rotundifolia were not only related to total phenol content but also related to structures of phenolic compounds and primarily related to their hydroxylation and methylation patterns. Using methanolic extracts from aerial non-flowering tissues of C. rotundifolia, Said et al. [43] identified 27 compounds that were dominated by 16 phenolics. Among the identified phenolics, flavanols were the majority. The total flavonoids content from their study using above-ground tissues was 1.35 mg (QE)/g (DW). In our study, the total flavonoids contents were 24.82 mg (RE)/g (DW) and 15.93 mg (RE)/g (DW) for stem and leaf, respectively. The roots had the highest concentration of total flavonoids, suggesting higher potential compared to the other organs. Our results for leaf and stem were in the same range as those reported in Cissus quadrangularis and Cissus adnata aerial parts [44,45,46]. On the other hand, seven metabolites were identified by Alqahtani et al. [21], including cissoic acid, which belongs to a rare class of secondary metabolites, and cissuxinoside, which they characterized as a new sucrose diester. Further, their study identified and characterized three uncommon C-glycosyl flavones. In this study, 194 flavonoids were tentatively identified in C. rotundifolia, which was more than four times the number of metabolites previously identified in leaves [9]. We explored flavonoids in C. rotundifolia tissues, including roots that recorded the highest flavonoids concentration and found 8 tissue-specific metabolites (Figure 2B). In general, above-ground tissues are used for medicinal preparation in C. rotundifolia, reporting analgesic, anti-ulcerative, and anti-inflammatory properties [10]. However, unique metabolites identified from roots in our study indicated potential use, for example jaceosidin [47] for antidiabetic properties, 8-prenylkaempferol for osteogenesis properties [48], and aureusidin for anti-inflammatory properties, among others. Roots and other underground plant parts have been examined in other medicinal plants and compounds with therapeutic activities reported [49]. Therefore, we postulate that root tissues of C. rotundifolia can also be sources for medicinal preparations.
Flavonols, flavones, and flavone C-glycosides were the subclasses with most of the compounds detected in all the organs in this study (Figure 2B). Isoflavones have been associated with pharmacological effects in human medicine, including in weakening menopausal hot flashes [50]. Generally, isoflavones are mainly available in legumes conferring protective roles to plants and nodules. However, from our analysis, five isoflavones with a wide range of pharmaceutical properties, ranging from antioxidants anti-diabetic to anti-mutagenic, were detected, including, calycosin [51], formononetin 7-O-glucoside (Ononin) [52], and pratensein [53].
Cissus rotundifolia is a valuable medicinal plant. From our analysis, the species is enriched with flavonoids. However, limited information has hindered its exploration at the molecular level. In our study, we provide a transcriptome assembly for C. rotundifolia and examine the genes related to flavonoids biosynthesis. Digital expression patterns for candidate genes involved in flavonoids biosynthesis were studied using FPKM (Figure 4B). Flavonoids biosynthesis-related genes showed unique expression patterns in leaves, stems, and roots.
Anthocyanidins form an important component of metabolites, playing a key role in plants response to abiotic stress, as well as in human health, and anthocyanidin reductase (ANR) is a vital enzyme in their biosynthesis [54]. High expression levels for genes encoding for enzymes that participate in the flavonoids biosynthesis, such as DFR and ANS, have been associated with a higher accumulation of anthocyanins [55]. Our analysis detected two genes encoding for ANS and one for DFR. The three genes were highly expressed in root compared to other tissues, suggesting their activity in the root, which had the highest flavonoids concentration among the organs. In Dendrobium officinale, two CHSs were expressed in all of the tissues, but the expressions were especially high in leaves, while five DFRs were highly expressed in stems, one in leaves and one in roots, and one of the two LARs was specifically expressed in stems, and the other one was expressed in leaves [56]. In grapevine, the mRNA of CHS3 accumulated primarily within the berry skin of red cultivars throughout coloration, whereas those of CHS1 and CHS2 accumulated within the leaves and berry skin of both the white and red cultivars [57].
Through a reaction associated with both CHS and CHI, the genes encoding these enzymes were reported to be critical in flavonoids biosynthesis and correlated with flavonoids accumulation [58,59]. A single CHS encoding gene was identified in our study with the highest expression levels in roots. Moreover, seven CHI encoding genes were identified, and their expression levels were generally high in the root tissues. The variation in transcription profile for genes encoding these enzymes across tissues of C. rotundifolia could therefore imply their important role in flavonoids metabolism. Additionally, FLS is a key enzyme in the conversion of dihydroflavonol to flavonol, thereby associated with the accumulation of flavonols and their composition [60]. Expression levels for FLS in our study were variable, with three genes exhibiting higher levels in the leaf while root had two genes had higher expression levels for FLS encoding genes. The variable expression levels of FLSs across the three tissues could be suggested to have contributed to the higher abundance or the diverse flavonol compounds (Figure 2A). Expression levels of other genes involved in flavonoids biosynthesis have also been linked with flavonoids accumulation levels [61]. From our analysis, no gene encoding FNS was detected. This is possibly due to a lack of similarity sequences or due to small-sized transcript fragments, which could not be detected. Generally, from our study, flavonoids biosynthesis-related genes were highly expressed in root and lowly expressed in leaf. The combined expression of genes encoding for these enzymes may partly be associated with the variable concentration of flavonoids and related metabolites in the tissues of C. rotundifolia.
Transcription factors modulate the expression of genes related to biosynthesis, as well as accumulation of secondary metabolites, and are therefore critical in the molecular examination of metabolites accumulation and synthesis [62]. Flavonoids biosynthesis has been reported to be regulated by a complex of transcription factors comprising MYB, bHLH, and WD40 families of transcription factors [63,64]. Examining flavonoids biosynthesis and accumulation across tissues of Anoectochilus roxburghii, MYB transcription factors encoding genes were identified and correlated to the observed differential expression and accumulation of flavonoids-related genes [65]. In E. konishii, MYB unigenes were identified, and their expression patterns were related to the observed accumulation of metabolites in different tissues. Specifically, consistent expression patterns were reported for a MYB and FLS gene, suggesting their role in the observed tissue-specific accumulation of rutin [66]. In our analysis, we identified transcription factor families, including bHLH, MYB, and WD40. In the root tissues, the three transcription factors were generally highly expressed when compared with either leaf or stem tissues (Supplementary Figure S5). It can therefore be suggested that these transcription factors also could have contributed to the observed flavonoids accumulation patterns across tissues of C. rotundifolia.

4. Materials and Methods

4.1. Plant Materials

Mature stem cuttings of C. rotundifolia were collected from Kenya around the Cherangani hills forest reserve. The duplicate voucher specimens SAJIT Z0041 were deposited in the Herbarium of Wuhan Botanical Garden, CAS (HIB) and in the Herbarium of the National museums of Kenya (NMK). The stem cuttings were propagated in pots in the greenhouse for 60 days, with adequate watering and nutrients under natural lighting conditions. Mature tissues were collected, rapidly cleaned, and immediately frozen in liquid nitrogen, after which they were stored at −80 °C until use for metabolome analysis and RNA extraction. For total flavonoids content determination, after cleaning, the samples were sliced into smaller sections and dried in the oven. The dry material was milled and sifted through a 40 mesh filter and stored in tubes for use in total flavonoids content determination. All materials for the study of C. rotundifolia were obtained from a single individual main rootstock. Three replicates were collected from each tissue for RNA extraction and metabolomics.

4.2. Metabolome

Sample preparation for metabolomics, extraction, identification as well as quantification of the compounds was carried out following conventional procedures developed by Wuhan Metware Biotechnology Co., Ltd. (https://www.metware.cn/ (Accessed on: 25 October 2021)). The repeatability and reliability of the extraction and detection methods were evaluated by analyzing the overlapping of total ion current (TIC) by using quality control samples. The quality control sample was prepared by combining all sample extracts into a combined sample and was injected after every 10 experimental samples.

4.2.1. Sample Preparation and Extraction

The frozen samples were freeze-dried before grinding using a freeze-dryer (SCIENTZ-100F/A; Ningbo Scientz Biotechnology Co., Ltd. Ningbo, China). Freeze-dried samples were crushed using a mixer mill (MM 400, Retsch) with a zirconia bead for 1.5 min at 30 Hz. In total, 100 mg powder was weighed and extracted overnight at 4 °C with 1.0 mL 70% aqueous methanol. Following centrifugation at 10,000× g for 10 min, the extracts were absorbed, (CNWBOND Carbon-GCB SPE Cartridge, 250 mg, 3 mL; ANPEL, Shanghai, China, https://www.anpel.com.cn/ (Accessed on: 25 October 2021)) eluted, and filtrated (SCAA-104, 0.22μm pore size; ANPEL, Shanghai, China, https://www.anpel.com.cn/ (Accessed on: 25 October 2021)) before LC-MS analysis.

4.2.2. HPLC Conditions

The sample extracts were analyzed using a LC-ESI-MS/MS system (HPLC, Shim-pack UFLC SHIMADZU CBM30A system, https://www.shimadzu.com.cn/ MS, Applied Biosystems 4500 Q TRAP, www.appliedbiosystems.com.cn/ (Accessed on: 25 October 2021)). The analytical conditions were as follows: HPLC: column, waters ACQUITY UPLC HSS T3 C18 (1.8 µm, 2.1 mm × 100 mm); solvent system, solvent A (water, 0.04% acetic acid) solvent B (acetonitrile, 0.04% acetic acid); gradient program, 100:0 V(A)/V(B) at 0 min, 5:95 V(A)/V(B) at 11.0min, 5:95 V(A)/V(B) at 12.0 min, 95:5 V(A)/V(B) at 12.1 min, 95:5 V(A)/V(B) at 15.0 min; flow rate, 0.40 mL/min; temperature, 40 °C; injection volume: 5 μL. The effluent was alternatively connected to an ESI-triple quadrupole-linear ion trap (Q TRAP)-MS.

4.2.3. ESI-Q TRAP-MS/MS

Linear ion trap (LIT) and triple quadrupole (QQQ) scans were acquired on a triple quadrupole-linear ion trap mass spectrometer (Q TRAP; API 4500 Q TRAP LC/MS/MS System), equipped with an electrospray ionization (ESI) Turbo Ion-Spray interface, operating in a positive and negative ion mode and controlled by Analyst 1.6.3 software (AB Sciex). The ESI source operation parameters were as follows: an ion source, turbo spray; source temperature 550 °C; ion spray voltage (IS) 5500 V; ion source gas I (GSI), gas II (GSII), curtain gas (CUR) was set at 55, 60, and 25.0 psi, respectively; the collision gas (CAD) was high. Instrument tuning and mass calibration were performed with 10 and 100 μmol/L polypropylene glycol solutions in QQQ and LIT modes, respectively. QQQ scans were acquired as multiple reaction monitor (MRM) experiments with collision gas (nitrogen) set to 5 psi. Declustering potential (DP) and collision energy (CE) for individual MRM transitions were done with further DP and CE optimization. A specific set of MRM transitions were monitored for each period according to the metabolites eluted within this period.

4.2.4. Qualitative and Quantitative Analysis of Metabolites

Qualitative analysis of primary and secondary MS data was carried out by comparison of the accurate precursor ions (Q1), product ions (Q3) values, the retention time (RT), and the fragmentation patterns with those obtained by injecting standards using the same conditions if the standards were available (Sigma-Aldrich, St. Louis, MO, USA, http://www.sigmaaldrich.com/united-states.html (Accessed on: 25 October 2021)) or conducted using a self-compiled database MWDB (MetWare biological science and Technology Co., Ltd. Wuhan, China) and publicly available metabolite databases if the standards were unavailable. Repeated signals of K+, Na+, NH4+, and other large molecular weight substances were eliminated during identification. The quantitative analysis of metabolites was based on the multiple reaction monitoring (MRM) mode. In the quadrupole (Q Trap), the precursor ions (parent ions) of the target compound were first selected. To eliminate the interference by non-target substances, the precursor ions were ionized by the collision chamber forming other fragment ions. Fragment ions were screened through the triple quadrupole, to select the specific fragment ion while eliminating the interference of the non-target ions. The characteristic ions of each metabolite were screened through the QQQ mass spectrometer to obtain the signal strengths. Integration and correction of chromatographic peaks for similar metabolites in different samples were performed using MultiQuant version 3.0.2 (AB SCIEX, Concord, ON, Canada). The corresponding relative metabolite contents were represented as chromatographic peak area integrals. The analysis was carried out in triplicates for each set sample.

4.3. Determination of Total Flavonoids Content (TFC)

4.3.1. TFC Extraction

Flavonoids extraction was carried out using the methods described by AL-Bukhaiti et al. [9] for phenolic extraction with minor modifications. The method is based on the spectrometric determination of the complex formed when flavonoids react with aluminium chloride for quantification and have widely been accepted. In summary, 1 g ground material was mixed with 50 mL methanol (90%), soaked for 3 h, and subjected to ultrasonic-assisted extraction. Following ultrasonic treatment, centrifugation was carried out at 11,000× g for 10 min and the supernatant was obtained. The extraction process was repeated three times, and the collected solution was combined. Filtration of the supernatant was carried out using Whatmann filter paper, and the resulting extract was stored for use in TFC determination.

4.3.2. Determination of TFC

The total flavonoids content (TFC) was determined using the colorimetric method as described by Zuo et al. [67], with some modifications. Briefly, 80 µL of twice diluted crude extract was mixed with NaNO2 (80 µL 5% W/V) solution and then shaken for 6 min. AlCl3 (80 µL 10% W/V) was added and allowed to stand for 6 min. Then, NaOH (400 µL 4% W/V) solution was added and allowed to react for 15 min. Afterward, the absorbance of the reaction mixture was read at 510 nm with UV/VIS spectrophotometer with methanol used as the blank. TFC of each sample was determined from a rutin standard curve, and the results were expressed in mg of rutin in 1 g dry material.

4.4. Transcriptomics

4.4.1. Total RNA Extraction, RNA Library Construction, and Sequencing

Using the three tissues from C. rotundifolia, RNA was separately extracted using a general plant total RNA extraction kit (BioTeke Corporation., Ltd. cat. NO RP3301, Wuxi, China) following the methods by Chomczynski and Sacchi [68]. The extracted RNA samples were quantified using a NanoDropTM OneC spectrophotometer (Thermo Fisher Scientific Inc. Waltham, USA.) and the quality was confirmed by agarose gel electrophoresis. The RNA integrity number (RIN) for the samples was 7.8 for stem and leaf and 7.6 for root tissue. To generate libraries for Illumina sequencing, Oligo(dT)-attached magnetic beads were used to purify mRNA. Purified mRNA was fragmented into small pieces with fragment buffer at an appropriate temperature. First-strand cDNA was generated using random hexamer-primed reverse transcription, followed by second-strand cDNA synthesis. Afterward, A-Tailing Mix and RNA Index adapters were added by incubating to end repair. The cDNA fragments obtained from the previous step were amplified by PCR, and products purified by Ampure XP Beads, then dissolved in EB solution. The product was validated on the Agilent Technologies 2100 bioanalyzer for quality control. The double-stranded PCR products from the previous step were heat-denatured and circularized by the splint oligo sequence to get the final library. The single-strand circular DNA (ssCir DNA) was formatted as the final library. The final library was amplified with phi29 to make a DNA nanoball (DNB) which had more than 300 copies of one molecular. DNBs were loaded into the patterned nanoarray and paired-end 150 base reads were generated on the MGISEQ-2000 platform (BGI-Shenzhen, China). Raw data was submitted to the NCBI sequence read archive (SRA) database (https://dataview.ncbi.nlm.nih.gov/object/PRJNA728209?reviewer=oeqs1aq76g906vji4doppnf225 (Accessed on: 25 October 2021)).

4.4.2. RNA-Seq Data Analysis and Functional Annotation

Raw sequencing data were first filtered by Trimmomatic (version 0.38), low-quality reads were discarded, and the reads contaminated with adaptor sequences were trimmed. Clean reads were then mapped to the reference genome of C. rotundifolia (deposited at the Genomic Data Center, BioProject PRJCA005006) using TopHat software (version 2.1.1) with default parameters. Reads mapped to the exon regions of each gene were counted by Cufflinks (version 2.2.1) and then fragments per kilobase per million fragments (FPKM) were calculated.
Genes differentially expressed between groups were identified using the Cuffdiff (version 2.2.1). A p-value cutoff of 0.05 and a fold-change cutoff of 2 were used to judge the statistical significance of gene expression differences. Gene ontology (GO) (https://wego.genomics.cn/ (Accessed on: 25 October 2021)), Clusters of Orthologous Groups of proteins (COG), and Kyoto encyclopedia of genes and genomes (KEGG) annotation for expressed genes were implemented by eggNOG software. The visualization of GO annotation was by WEGO 2.0 (https://wego.genomics.cn/ (Accessed on: 25 October 2021)) while COG functional groups were visualized using GraphPad Prism version 8.0.0 for Windows (GraphPad Software, San Diego, CA, USA).

4.4.3. RNA-Seq Data Validation Using qPCR

Total RNA was extracted using a general plant total RNA extraction kit (BioTeke Corporation Co., Ltd. Cat. # RP3301, Wuxi, China). DNA was removed from the RNA sample using the RQ1 RNase-Free DNase kit (Promega Cat. # M6101, Beijing, China) following the manufacturer’s instructions. First-strand cDNA synthesis was carried out using HiScript® III 1st Strand cDNA Synthesis Kit (Vazyme Biotech Co., Ltd. Cat # R312-01, Nanjing, China) following the manufacturer’s instructions in a 20 µL total reaction volume. The qPCR analysis was performed using a 7500 Fast Real-Time PCR system (Applied Biosystems, MA, USA) in a total of 10 µL reaction volume using ChamQ universal SYBR qPCR Master Mix (Vazyme Biotech Co., Ltd., Vazyme code: Q711-02, Nanjing, China). The amplification conditions were 50 °C for 20 s, 95 °C for 10 min, 40 cycles of 95 °C for 15 s, and 60 °C for 1 min. Three biological and two technical replicates per sample were carried out, and Actin was used as the internal standard. The Ct value was determined using the instrument’s software, and the relative quantification of gene expression was monitored after normalization using Actin. The relative transcription levels were calculated using the ΔΔCt method [69], and leaf was considered as the control tissue for normalized relative expression.

5. Conclusions

In this study through metabolomic approaches, 194 flavonoids were tentatively identified in C. rotundifolia. Additionally, 18,427 expressed genes from leaf, stem, and root transcriptome were identified and mapped onto the C. rotundifolia genome. Regulatory mechanisms involved in flavonoids accumulation were also explored across the tissues through comparative analysis of metabolite accumulation and the expression profile for flavonoids-related biosynthesis genes. In addition to the enzymes involved in the central flavonoids biosynthesis pathway, which we suggest may have influenced flavonoids accumulation in different tissues of C. rotundifolia, we hypothesize that transcription factors may also have contributed to the variable flavonoid concentration across tissues. To our knowledge, no prior report has been made on differences in flavonoids accumulation among different tissues of C. rotundifolia that are used variably for medicinal value. Our study provides valuable information about flavonoids metabolites and contributes to molecular research in Cissus rotundifolia, as well as other members in this genus and facilitates exploration of their medicinal uses.

Supplementary Materials

S1: The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/metabo11110741/s1. Table S1: Flavonoids metabolites in the tissues of C. rotundifolia. The table represents the compound identification details including retention time, molecular ions and fragmentation. Differentially accumulated metabolites are also included, Table S2: GO details for major sub-ontologies for expressed genes in C. rotundifolia transcriptome, Table S3: KEGG details for expressed genes for C. rotundifolia, Table S4: qPCR primers. The primers were designed using NCBI Primer-BLAST. PAL, phenylalanine ammonialyase; C4H, Cinnamate 4-hydroxylase; CHS, chalcone synthase; CHI, chalcone isomerase; F3H, flavanone 3-hydroxylase; 4CL, 4-Coumaric acid: CoA ligase; FLS, flavonol synthase. Table S5: Statistics for GO annotation of DEGs identified in C. rotundifolia, Figure S1: Overall qualitative and quantitative analysis of metabolomic data. PCA analysis for the 3 C. rotundifolia tissues. The x-axis represents the first primary principal component while the y-axis represents the second primary principal component. Distinct patterns were observed within different tissues, Figure S2: Differentially accumulated metabolites (DAMs) among leaf, stem, and root for C. rotundifolia. (a) Volcano plot representing DAMs in L-vs-S; (b) DAMs in R-vs-L; (c) DAMs in R-vs-S. The spots represent the DAMs; red for up-accumulated, green for down-accumulated while black for those not significantly changed, Figure S3: Functional annotation of C. rotundifolia transcriptome. (a) Gene Ontology (GO) Classification. The three main categories were identified (cellular components, molecular function, and biological process). The left y-axis represents the gene percentage while the right y-axis indicates the number of genes in the categories. (b) COG terms. The genes were classified into 25 functional categories. The letters represent respective functional categories, Figure S4: RNA-seq validation by qPCR. The histograms indicate the qPCR results for 12 selected genes involved in flavonoid biosynthesis in 3 organs of C. rotundifolia. The error bars represent the mean SD of three biological replicates, Figure S5: Transcription factors expression profile. The expression major transcription factor families involved in flavonoid biosynthesis were analysed. (a) bHLH transcription factors, (b) MYB transcription factors, (c) WD40 transcription factors, Figure S6: Classification of C. rotundifolia identified DEGs to KEGG pathways. (a) Leaf-vs-Stem; (b); Root-vs-Leaf (c) Root-vs-Stem.

Author Contributions

D.K.G., H.X., and Q.W conceived and designed the experiments. D.K.G., Q.L., Y.L., and Y.H. carried out experiments and analyzed the data, H.Z., and M.M. collected samples. C.X., and Z.Z. prepared the figures and tables. L.W., F.M.M., and Q.W. contributed to the discussion of the project at different stages. D.K.G., H.X., and Q.L. wrote the initial manuscript. All authors revised and contributed to the final version of the text. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the International Cooperation of National Natural Science Foundation of China (NSFC No. 31961143026).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Transcriptome raw data has been deposited in the SRA database (to be released; https://dataview.ncbi.nlm.nih.gov/object/PRJNA728209?reviewer=oeqs1aq76g906vji4doppnf225 (Accessed on: 25 October 2021)). The genome data used in this study were deposited on National Genomics Data Center (NGDC, https://ngdc.cncb.ac.cn/ (Accessed on: 25 October 2021)) with the project numbers PRJCA005006. All data can be obtained from the corresponding author upon reasonable request.

Acknowledgments

We thank Gunlin Chen, Wuhan Botanical Garden, Chinese Academy of Sciences (CAS) for his support during our project. We are grateful to the Kenya Forest Service (KFS) for issuing fieldwork permits (permit number: RESEA/1/KFS/VOL. III (48)) to conduct the field investigations in the Cherangani forest reserve.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Alzoreky, N.S.; Nakahara, K. Antibacterial activity of extracts from some edible plants commonly consumed in Asia. Int. J. Food Microbiol. 2003, 80, 223–230. [Google Scholar] [CrossRef]
  2. Onyechi, U.A.; Judd, P.A.; Ellis, P.R. African plant foods rich in non-starch polysaccharides reduce postprandial blood glucose and insulin concentrations in healthy human subjects. Br. J. Nutr. 1998, 80, 419–428. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Alshehri, S.A. Antidiabetic Activity of Cissus Rotundifolia Plant Growing in Saudi Arabia; South Dakota State University: Brookings, SD, USA, 2020. [Google Scholar]
  4. Fernandes, G. Medicinal properties of plants from the genus Cissus: A Review. J. Med. Plants Res. 2012, 30, 6. [Google Scholar] [CrossRef]
  5. Korish, M. Nutritional evaluation of wild plant Cissus rotundifolia. Ital. J. Food Sci. 2016, 28, 43–49. [Google Scholar]
  6. Onyechi, U.A.; Ibeanu, V.N. Effects of diets containing Cissus rotundifolia flour on lipid profile of rats and postprandial glucose levels of normoglycemic human adults. Afr. J. Biotechnol. 2016, 15, 557–564. [Google Scholar]
  7. Al-Mehdar, A.A.; Al-Battah, A.M. Evaluation of hypoglycemic activity of Boswellia carterii and Cissus rotundifolia in streptozotocin/nicotinamide-induced diabetic rats. Yemeni J. Med. Sci. 2016, 10, 30–38. [Google Scholar] [CrossRef]
  8. Ali, A.A.N.; Al-Rahwi, K.; Lindequist, U. Some medicinal plants used in Yemeni Herbal Medicine to treat malaria. Afr. J. Tradit. Complement Altern. Med. 2004, 1, 72–76. [Google Scholar] [CrossRef]
  9. AL-Bukhaiti, W.Q.; Noman, A.; Mahdi, A.A.; Ali, A.H.; Abed, S.M.; Rashed, M.M.A.; Wang, H. Profiling of phenolic compounds and antioxidant activities of Cissus rotundifolia (Forssk.) as influenced by ultrasonic-assisted extraction conditions. J. Food Meas. Charact. 2019, 13, 634–647. [Google Scholar] [CrossRef]
  10. Said, A.A.; Aboutabl, E.A.; El Awdan, S.A.; Raslan, M.A. Proximate analysis, phytochemical screening, and bioactivities evaluation of Cissus rotundifolia (Forssk.) Vahl. (Fam. Vitaceae) and Sansevieria cylindrica Bojer ex Hook. (Fam. Dracaenaceae) growing in Egypt. Egypt Pharm. J. 2015, 14, 180. [Google Scholar] [CrossRef]
  11. Stavric, B. Antimutagens and anticarcinogens in foods. Food Chem. Toxicol. 1994, 32, 79–90. [Google Scholar] [CrossRef]
  12. Barreca, D.; Bellocco, E.; Caristi, C.; Leuzzi, U.; Gattuso, G. Elucidation of the flavonoid and furocoumarin composition and radical-scavenging activity of green and ripe chinotto (Citrus myrtifolia Raf.) fruit tissues, leaves and seeds. Food Chem. 2011, 129, 1504–1512. [Google Scholar] [CrossRef]
  13. Cushnie, T.P.T.; Lamb, A.J. Antimicrobial activity of flavonoids. Int. J. Antimicrob. Agents 2005, 26, 343–356. [Google Scholar] [CrossRef]
  14. García-Salas, P.; Gómez-Caravaca, A.M.; Arráez-Román, D.; Segura-Carretero, A.; Guerra-Hernández, E.; García-Villanova, B.; Fernández-Gutiérrez, A. Influence of technological processes on phenolic compounds, organic acids, furanic derivatives, and antioxidant activity of whole-lemon powder. Food Chem. 2013, 141, 869–878. [Google Scholar] [CrossRef]
  15. Heim, K.E.; Tagliaferro, A.R.; Bobilya, D.J. Flavonoid antioxidants: Chemistry, metabolism and structure-activity relationships. J. Nutr. Biochem. 2002, 13, 572–584. [Google Scholar] [CrossRef]
  16. Middleton, E.; Kandaswami, C. Effects of flavonoids on immune and inflammatory cell functions. Biochem. Pharmacol. 1992, 43, 1167–1179. [Google Scholar] [CrossRef]
  17. Chen, H.Q.; Zuo, W.J.; Wang, H.; Shen, H.Y.; Luo, Y.; Dai, H.F.; Mei, W.L. Two new antimicrobial flavanes from dragon’s blood of Dracaena cambodiana. J. Asian Nat. Prod. Res. 2012, 14, 436–440. [Google Scholar] [CrossRef]
  18. Luo, Y.; Wang, H.; Zhao, Y.X.; Zeng, Y.-B.; Shen, H.Y.; Dai, H.F.; Mei, W.L. Cytotoxic and antibacterial flavonoids from dragon’s blood of Dracaena cambodiana. Planta Med. 2011, 77, 2053–2056. [Google Scholar] [CrossRef] [Green Version]
  19. Zheng, Q.A.; Xu, M.; Yang, C.R.; Wang, D.; Li, H.Z.; Zhu, H.T.; Zhang, Y.J. Flavonoid oligomers from Chinese dragon’s blood, the red resins of Dracaena cochinchinensis. Nat. Prod. Bioprospect. 2012, 2, 111–116. [Google Scholar] [CrossRef] [Green Version]
  20. Rajkapoor, B.; Murugesh, N.; Rama Krishna, D. Cytotoxic activity of a flavanone from the stem of Bauhinia variegata Linn. Nat. Prod. Res. 2009, 23, 1384–1389. [Google Scholar] [CrossRef]
  21. Alqahtani, J.; Formisano, C.; Chianese, G.; Luciano, P.; Stornaiuolo, M.; Perveen, S.; Taglialatela-Scafati, O. Glycosylated phenols and an unprecedented diacid from the Saudi plant Cissus rotundifolia. J. Nat. Prod. 2020, 83, 3298–3304. [Google Scholar] [CrossRef]
  22. Nabavi, S.M.; Šamec, D.; Tomczyk, M.; Milella, L.; Russo, D.; Habtemariam, S.; Suntar, I.; Daglia, M.; Xu, S.; Rastrelli, L.; et al. Flavonoid biosynthetic pathways in plants: Versatile targets for metabolic engineering. Biotechnol. Adv. 2020, 38, 107316. [Google Scholar] [CrossRef]
  23. Tohge, T.; de Souza, L.P.; Fernie, A.R. Current understanding of the pathways of flavonoid biosynthesis in model and crop plants. J. Exp. Bot. 2017, 68, 4013–4028. [Google Scholar] [CrossRef]
  24. Davies, K.M.; Schwinn, K.E. Flavonoids. Chemistry, Biochemistry and Applications; Andersen, O.M., Markham, K.R., Eds.; Taylor & Francis Group: Abingdon, UK, 2006; pp. 149–178. [Google Scholar]
  25. Tripoli, E.; La Guardia, M.; Giammanco, S.; Majo Di, D.; Giammanco, M. Citrus flavonoids: Molecular structure, biological activity and nutritional properties: A review. Food Chem. 2007, 104, 466–479. [Google Scholar] [CrossRef]
  26. Wang, J.L.; Li, X.; Jiang, D.F.; Ma, P.; Yang, C.R. Chemical constituents of dragons blood resin from Dracaena cochinchinensis in Yunnan and their antifungal activity. Plant Divers 1995, 17, 1–3. [Google Scholar]
  27. Lloyd, A.; Brockman, A.; Aguirre, L.; Campbell, A.; Bean, A.; Cantero, A.; Gonzalez, A. Advances in the MYB-bHLH-WD Repeat (MBW) pigment regulatory model: Addition of a WRKY Factor and co-option of an anthocyanin MYB for betalain regulation. Plant Cell Physiol. 2017, 58, 1431–1441. [Google Scholar] [CrossRef] [Green Version]
  28. Amato, A.; Cavallini, E.; Walker, A.R.; Pezzotti, M.; Bliek, M.; Quattrocchio, F.; Koes, R.; Ruperti, B.; Bertini, E.; Zenoni, S.; et al. The MYB5-driven MBW complex recruits a WRKY factor to enhance the expression of targets involved in vacuolar hyper-acidification and trafficking in grapevine. Plant J. 2019, 99, 1220–1241. [Google Scholar] [CrossRef]
  29. Winkel, B. The Biosynthesis of Flavonoids. In The Science of Flavonoids; Springer: New York, NY, USA, 2006; pp. 71–95. [Google Scholar] [CrossRef]
  30. Sun, Y.J.; He, J.M.; Kong, J.Q. Characterization of two flavonol synthases with iron-independent flavanone 3-hydroxylase activity from Ornithogalum caudatum Jacq. BMC Plant Biol. 2019, 19, 195. [Google Scholar] [CrossRef]
  31. Matsuda, F.; Hirai, M.Y.; Sasaki, E.; Akiyama, K.; Yonekura-Sakakibara, K.; Provart, N.J.; Sakurai, T.; Shimada, Y.; Saito, K. AtMetExpress development: A phytochemical atlas of Arabidopsis development. Plant Physiol. 2010, 152, 566–578. [Google Scholar] [CrossRef] [Green Version]
  32. Albinsky, D.; Kusano, M.; Higuchi, M.; Hayashi, N.; Kobayashi, M.; Fukushima, A.; Masaki, M.; Takanari, I.; Keiko, M.; Yoko, H.; et al. Metabolomic screening applied to rice FOX Arabidopsis lines leads to the identification of a gene-changing nitrogen metabolism. Mol. Plant. 2010, 3, 125–142. [Google Scholar] [CrossRef]
  33. Abad-García, B.; Garmón-Lobato, S.; Berrueta, L.A.; Gallo, B.; Vicente, F. Online characterization of 58 phenolic compounds in Citrus fruit juices from Spanish cultivars by high-performance liquid chromatography with photodiode-array detection coupled to electrospray ionization triple quadrupole mass spectrometry. Talanta 2012, 99, 213–224. [Google Scholar] [CrossRef]
  34. Chen, W.; Gong, L.; Guo, Z.; Wang, W.; Zhang, H.; Liu, X.; Yu, S.; Xiong, L.; Luo, J. A novel integrated method for large-scale detection, identification, and quantification of widely targeted metabolites: Application in the study of rice metabolomics. Mol. Plant. 2013, 6, 1769–1780. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Dong, X.; Chen, W.; Wang, W.; Zhang, H.; Liu, X.; Luo, J. Comprehensive profiling and natural variation of flavonoids in rice. J. Integr. Plant Biol. 2014, 56, 876–886. [Google Scholar] [CrossRef] [PubMed]
  36. Gong, L.; Chen, W.; Gao, Y.; Liu, X.; Zhang, H.; Xu, C.; Yu, S.; Zhang, Q.; Luo, J. Genetic analysis of the metabolome exemplified using a rice population. Proc. Natl. Acad. Sci. USA 2013, 110, 20320–20325. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Liu, Q.; Wang, X.; Tzin, V.; Romeis, J.; Peng, Y.; Li, Y. Combined transcriptome and metabolome analyses to understand the dynamic responses of rice plants to attack by the rice stem borer Chilo suppressalis (Lepidoptera: Crambidae). BMC Plant Biol. 2016, 16, 259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Saito, K.; Yonekura-Sakakibara, K.; Nakabayashi, R.; Higashi, Y.; Yamazaki, M.; Tohge, T.; Fernie, A.R. The flavonoid biosynthetic pathway in Arabidopsis: Structural and genetic diversity. Plant Physiol. Biochem. 2013, 72, 21–34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Galperin, M.Y.; Kristensen, D.M.; Makarova, K.S.; Wolf, Y.I.; Koonin, E.V. Microbial genome analysis: The COG approach. Brief Bioinform. 2019, 20, 1063–1070. [Google Scholar] [CrossRef]
  40. Zhu, J.H.; Cao, T.J.; Dai, H.F.; Li, H.L.; Guo, D.; Mei, W.L.; Peng, S.Q. De Novo transcriptome characterization of Dracaena cambodiana and analysis of genes involved in flavonoid accumulation during formation of dragon’s blood. Sci. Rep. 2016, 6, 38315. [Google Scholar] [CrossRef] [Green Version]
  41. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Fan, L.; Zeng, Q.; Chen, Z.; Freidman, N.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [Green Version]
  42. Al-Mamary, M.A.J. Antioxidant activity of commonly consumed vegetables in Yemen. Malays J. Nutr. 2002, 8, 179–189. [Google Scholar]
  43. Said, A.; Aboutabl, E.A.; Melek, F.R.; Jaleel, G.A.R.B.; Raslan, M. Phytoconstituents profiling of Cissus rotundifolia (Forssk.) Vahl. by HPLC-MS/MS, and evaluation of its free radical scavenging activity (DPPH) and cytotoxicity. Trends Phytochem. Res. 2018, 2, 65–74. [Google Scholar]
  44. Kaewnarin, K.; Niamsup, H.; Shank, L.; Rakariyatham, N. Antioxidant and antiglycation activities of some edible and medicinal plants. Chiang Mai J. Sci. 2014, 41, 105–116. [Google Scholar]
  45. Vijayalakshmi, A.; Kumar, P.R.; Sakthi Priyadarsini, S.; Meenaxshi, C. In vitro antioxidant and anticancer activity of flavonoid fraction from the aerial parts of Cissus quadrangularis Linn. against human breast carcinoma cell lines. J. Chem. 2013, 2013, 150675. [Google Scholar] [CrossRef] [Green Version]
  46. Shoibe, M.; Chy, M.; Uddin, N.; Alam, M.; Adnan, M.; Islam, M.; Nihar, S.W.; Rahman, N.; Suez, E. In vitro and in vivo biological activities of Cissus adnata (Roxb.). Biomedicines 2017, 5, 63. [Google Scholar] [CrossRef] [Green Version]
  47. Park, E.; Hong, K.; Kwon, B.-M.; Kim, Y.; Kim, J.-H. Jaceosidin ameliorates insulin resistance and kidney dysfunction by enhancing insulin receptor signaling and the antioxidant defense system in type 2 diabetic mice. J. Med. Food 2020, 23, 1083–1092. [Google Scholar] [CrossRef]
  48. Chiou, W.F.; Lee, C.H.; Liao, J.F.; Chen, C.C. 8-Prenylkaempferol accelerates osteoblast maturation through bone morphogenetic protein-2/p38 pathway to activate Runx2 transcription. Life Sci. 2011, 88, 335–342. [Google Scholar] [CrossRef]
  49. Rehman, S.U.; Choe, K.; Yoo, H.H. Review on a traditional herbal medicine, Eurycoma longifolia Jack (Tongkat Ali): Its traditional uses, chemistry, evidence-based pharmacology and toxicology. Molecules 2016, 21, 331. [Google Scholar] [CrossRef] [Green Version]
  50. Li, L.; Lv, Y.; Xu, L.; Zheng, Q. Quantitative efficacy of soy isoflavones on menopausal hot flashes. Br. J. Clin. Pharmacol. 2015, 79, 593–604. [Google Scholar] [CrossRef] [Green Version]
  51. Wang, Q.; Lu, W.; Yin, T.; Lu, L. Calycosin suppresses TGF-β-induced epithelial-to-mesenchymal transition and migration by upregulating BATF2 to target PAI-1 via the Wnt and PI3K/Akt signaling pathways in colorectal cancer cells. J. Exp. Clin. Cancer Res. 2019, 38, 240. [Google Scholar] [CrossRef] [Green Version]
  52. Lee, W.; Choo, S.; Sim, H.; Bae, J.-S. Inhibitory activities of ononin on particulate matter-induced oxidative stress. Biotechnol. Bioprocess. Eng. 2021, 26, 208–215. [Google Scholar] [CrossRef]
  53. Behbahani, M. Anti-human immunodeficiency virus-1 activities of pratensein and pratensein glycoside from Alhaji maurorum and its parasite Cuscuta kotchiana. Chin. J. Integr. Med. 2017, 2017, 1–6. [Google Scholar] [CrossRef]
  54. Tan, L.; Wang, M.; Kang, Y.; Azeem, F.; Zhou, Z.; Tuo, D.; María Preciado Rojo, L.; Khan, I.A.; Pan, Z. Biochemical and functional characterization of anthocyanidin reductase (ANR) from Mangifera indica L. Molecules 2018, 23, 2876. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Jiang, T.; Guo, K.; Liu, L.; Tian, W.; Xie, X.; Wen, S.; Wen, C. Integrated transcriptomic and metabolomic data reveal the flavonoid biosynthesis metabolic pathway in Perilla frutescens (L.) leaves. Sci. Rep. 2020, 10, 16207. [Google Scholar] [CrossRef] [PubMed]
  56. Yuan, Y.; Zhang, J.; Liu, X.; Meng, M.; Wang, J.; Lin, J. Tissue-specific transcriptome for Dendrobium officinale reveals genes involved in flavonoid biosynthesis. Genomics 2020, 112, 1781–1794. [Google Scholar] [CrossRef] [PubMed]
  57. Goto-Yamamoto, N.; Wan, G.H.; Masaki, K.; Kobayashi, S. Structure and transcription of three chalcone synthase genes of grapevine (Vitis vinifera). Plant Sci. 2002, 162, 867–872. [Google Scholar] [CrossRef]
  58. Muir, S.R.; Collins, G.J.; Robinson, S.; Hughes, S.; Bovy, A.; De Vos, C.H.R.; van Tunen, A.J.; Verhoeyen, M.E. Overexpression of petunia chalcone isomerase in tomato results in fruit containing increased levels of flavonols. Nat. Biotechnol. 2001, 19, 470–474. [Google Scholar] [CrossRef]
  59. Jez, J.M.; Bowman, M.E.; Dixon, R.A.; Noel, J.P. Structure and mechanism of the evolutionarily unique plant enzyme chalcone isomerase. Nat. Struct. Biol. 2000, 7, 786–791. [Google Scholar]
  60. Forkmann, G.; Martens, S. Metabolic engineering and applications of flavonoids. Curr. Opin. Biotechnol. 2001, 12, 155–160. [Google Scholar] [CrossRef]
  61. Zhang, S.; Li, H.; Liang, X.; Yan, Y.; Xia, P.; Jia, Y.; Liang, Z. Enhanced production of phenolic acids in Salvia miltiorrhiza hairy root cultures by combing the RNAi-mediated silencing of chalcone synthase gene with salicylic acid treatment. Biochem. Eng. J. 2015, 103, 185–192. [Google Scholar] [CrossRef]
  62. Cho, K.; Cho, K.-S.; Sohn, H.-B.; Ha, I.J.; Hong, S.-Y.; Lee, H.; Kim, Y.-M.; Nam, M.H. Network analysis of the metabolome and transcriptome reveals novel regulation of potato pigmentation. J. Exp. Bot. 2016, 67, 1519–1533. [Google Scholar] [CrossRef] [Green Version]
  63. Zhao, M.; Li, J.; Zhu, L.; Chang, P.; Li, L.; Zhang, L. Identification and characterization of MYB-bHLH-WD40 regulatory complex members controlling anthocyanidin biosynthesis in blueberry fruits development. Genes 2019, 10, 496. [Google Scholar] [CrossRef] [Green Version]
  64. Schaart, J.G.; Dubos, C.; Romero De La Fuente, I.; van Houwelingen, A.M.M.L.; de Vos, R.C.H.; Jonker, H.H.; Bovy, A.; Xu, W.; Routaboul, J.-M.; Lepiniec, L. Identification and characterization of MYB-bHLH-WD 40 regulatory complexes controlling proanthocyanidin biosynthesis in strawberry (F ragaria× ananassa) fruits. New Phytol. 2013, 197, 454–467. [Google Scholar] [CrossRef]
  65. Chen, Y.; Pan, W.; Jin, S.; Lin, S. Combined metabolomic and transcriptomic analysis reveals key candidate genes involved in the regulation of flavonoid accumulation in Anoectochilus roxburghii. Process Biochem. 2020, 91, 339–351. [Google Scholar] [CrossRef]
  66. Liang, W.; Ni, L.; Carballar-Lejarazú, R.; Zou, X.; Sun, W.; Wu, L.; Yuan, X.; Mao, Y.; Huang, W.; Zou, S. Comparative transcriptome among Euscaphis konishii Hayata tissues and analysis of genes involved in flavonoid biosynthesis and accumulation. BMC Genom. 2019, 20, 1–14. [Google Scholar] [CrossRef]
  67. Zou, Y.; Chang, S.K.C.; Gu, Y.; Qian, S.Y. Antioxidant activity and phenolic compositions of lentil (Lens culinaris var. Morton) extract and its fractions. J. Agric. Food Chem. 2011, 59, 2268–2276. [Google Scholar] [CrossRef] [Green Version]
  68. Chomczynski, P.; Sacchi, N. Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem. 1987, 162, 156–159. [Google Scholar] [CrossRef]
  69. Schmittgen, T.D.; Livak, K.J. Analyzing real-time PCR data by the comparative CT method. Nat. Protoc. 2008, 3, 1101–1108. [Google Scholar] [CrossRef]
Metabolites 11 00741 g001 550
Figure 1. Total flavonoids content for C. rotundifolia organs. Flavonoids content was expressed as the rutin equivalent mg/g of the dry weight. The data were expressed as mean ± SD for three replicates. The letters indicate a significant difference at p < 0.05.
Figure 1. Total flavonoids content for C. rotundifolia organs. Flavonoids content was expressed as the rutin equivalent mg/g of the dry weight. The data were expressed as mean ± SD for three replicates. The letters indicate a significant difference at p < 0.05.
Metabolites 11 00741 g001
Metabolites 11 00741 g002 550
Figure 2. Flavonoids and phenolic acids tentatively identified. (A) An integrated view of the identified compounds. The numbers correspond to the compounds in the subclass. (B) Metabolites unique or shared among the C. rotundifolia organs. The numbers represent the compounds in each section.
Figure 2. Flavonoids and phenolic acids tentatively identified. (A) An integrated view of the identified compounds. The numbers correspond to the compounds in the subclass. (B) Metabolites unique or shared among the C. rotundifolia organs. The numbers represent the compounds in each section.
Metabolites 11 00741 g002
Metabolites 11 00741 g003 550
Figure 3. Accumulation profile for tentatively identified compounds. (A) Heatmap showing the hierarchical cluster analysis for metabolites identified in C. rotundifolia tissues. The colors indicate the metabolites relative abundance with red indicating high values and green for metabolites with lower abundance. (B) Stacked bar graph showing the annotation of the identified DAMs between organs. The y-axis represents the number of compounds and the colors indicate the different subclasses.
Figure 3. Accumulation profile for tentatively identified compounds. (A) Heatmap showing the hierarchical cluster analysis for metabolites identified in C. rotundifolia tissues. The colors indicate the metabolites relative abundance with red indicating high values and green for metabolites with lower abundance. (B) Stacked bar graph showing the annotation of the identified DAMs between organs. The y-axis represents the number of compounds and the colors indicate the different subclasses.
Metabolites 11 00741 g003
Metabolites 11 00741 g004 550
Figure 4. A schematic representation of phenylpropanoid and flavonoids biosynthesis in C. rotundifolia. (A) Proposed biosynthetic pathway modified from literature and the KEGG pathway database. The numbers in the brackets after each gene name denotes the number of the respective genes in C. rotundifolia. (B) Heat map representation for flavonoids biosynthetic-related gene expression patterns. PAL: phenylalanine ammonia-lyase; C4H: cinnamic acid 4-hydroxylase; 4CL: 4-Coumaric acid: CoA ligase; CCR: cinnamoyl-CoA; HCT: shikimate-o-hydroxycinnamoyltransferase; CHS: chalcone synthase; CHI: chalcone isomerase; F3H: flavanone 3-hydroxylase; F3′5′H: flavonoid 3′,5′-hydroxylase; FLS: flavonol synthase; DFR: dihydroflavanol 4-reductase; LDOX/ANS: leucoanthocyanidin dioxygenase/anthocyanidin synthase; ANR: anthocyanidin reductase; LAR: leucoanthocyanidin reductase; FNS: flavone synthase.
Figure 4. A schematic representation of phenylpropanoid and flavonoids biosynthesis in C. rotundifolia. (A) Proposed biosynthetic pathway modified from literature and the KEGG pathway database. The numbers in the brackets after each gene name denotes the number of the respective genes in C. rotundifolia. (B) Heat map representation for flavonoids biosynthetic-related gene expression patterns. PAL: phenylalanine ammonia-lyase; C4H: cinnamic acid 4-hydroxylase; 4CL: 4-Coumaric acid: CoA ligase; CCR: cinnamoyl-CoA; HCT: shikimate-o-hydroxycinnamoyltransferase; CHS: chalcone synthase; CHI: chalcone isomerase; F3H: flavanone 3-hydroxylase; F3′5′H: flavonoid 3′,5′-hydroxylase; FLS: flavonol synthase; DFR: dihydroflavanol 4-reductase; LDOX/ANS: leucoanthocyanidin dioxygenase/anthocyanidin synthase; ANR: anthocyanidin reductase; LAR: leucoanthocyanidin reductase; FNS: flavone synthase.
Metabolites 11 00741 g004
Metabolites 11 00741 g005 550
Figure 5. Statistics of differentially expressed genes across three C. rotundifolia tissues. R, S, and L represent root, stem, and leaf respectively.
Figure 5. Statistics of differentially expressed genes across three C. rotundifolia tissues. R, S, and L represent root, stem, and leaf respectively.
Metabolites 11 00741 g005
Table
Table 1. A list of 30 flavonoids tentatively identified in C. rotundifolia tissues.
Table 1. A list of 30 flavonoids tentatively identified in C. rotundifolia tissues.
Metabolite NamePrecursor Ion (Q1) (Da)Product Ion (Q3) (Da)Retention Time (min)Main Fragments (Da)
Diosmetin (5,7,3′-Trihydroxy-4′-methoxyflavone)299.06256.046.1256.04, 284.03, 299.05, 299.13, 299.06
Quercetin-3-o-galactoside (Hyperin)463.1300.034.2300.03, 301.03, 463.09
Catechin gallate441.3289.084.2124.02, 125.03, 169.02, 193.01, 203.07, 245.08, 289.08, 331.05, 441.08
Kaempferol-3-o-rhamnoside (Afzelin)(Kaempferin)431284.045229.05, 227.04, 255.03, 284.04, 285.05, 431.11
Epiafzelechin275.1139.044.2107.05, 111.04, 121.06, 139.04, 149.06, 145.06, 173.06, 191.07, 275.09
Delphinidin303.05137.024.8137.02, 153.02, 165.02, 129.05, 257.04
Isoorientin-7-o-(6″-p-coumaroyl)glucoside757.2637.164.2147.04, 291.08, 309.09, 329.07, 353.07, 431.1, 449.11, 577.13, 637.15, 757.20
Kaempferol-3-o-rutinoside(Nicotiflorin)593.162854.3285.04, 593.15
Kaempferol-3-arabinopyranoside419.1133.055133.05, 287.06
Epicatechin291123.053.8119.05, 123.05, 139.04, 147.04, 165.06, 161.06, 179.07, 207.07, 291.09
Isohemiphloin433.12313.074.2125.02, 152.99, 211.06, 271.06, 331.07, 343.08, 359.15, 433.23, 433.11, 433.2
Calycosin285225.065.5225.09, 242.06, 269.04, 270.05, 285.08
5-Hydroxy-6,7,8,3′,4′-pentamethoxyflavone389.1359.087.2341.08, 359.07, 389.12
Pratensein301.07269.046.4167.04, 181.07, 258.05, 269.08, 286.05, 301.07
Aureusidin287.05153.025.6153.02, 287.06
5-Hydroxyauranetin389.1359.077341.06, 359.07, 389.12
Epigallocatechin305219.073125.02, 137.02, 167.03, 165.02, 179.03, 219.07, 221.05, 305.07
Gallocatechin307163.042.8123.04, 139.04, 163.04, 177.05, 195.06, 233.06
Kaempferol (3,5,7,4′-Tetrahydroxyflavone)285.05229.056.2151.00, 185.06, 211.04, 229.05,
Naringenin-7-o-glucoside (Prunin)4331514.8119.05, 151.00, 177.02, 255.03, 271.06, 284.03, 301.03, 417.08, 433.11
Luteolin-8-c-glucoside (Orientin)449.1329.074287.06, 299.06, 329.07, 353.07, 383.08, 413.09, 431.10, 449.11
Avicularin(Quercetin-3-o-α-L-arabinofuranoside)435.083034.5303.05, 257.04, 229.05
Quercetin-3-o-arabinoside (Guaijaverin)433.08255.034.7151.00, 179.00, 255.03, 271.02, 300.03, 301.04, 433.08
Isoorientin-7-o-glucoside611.1329.073.5299.06, 319.04, 329.07, 383.08, 431.10, 449.11, 465.10, 611.15
Vitexin-7-o-(6″-p-coumaroyl)glucoside741.2415.14.4147.04, 309.11, 313.07, 337.07, 415.10, 433.11, 741.20
Luteolin-7-o-rutinoside595.16287.054.3287.06, 449.11
Phloretin273.08123.046119.05, 123.04, 167.03, 273.08
Kaempferol-3,7-o-dirhamnoside (Kaempferitrin)579.2433.114.3287.05, 433.11
Quercetin303.04229.055.6153.02, 165.02, 229.05, 257.04, 285.04
Kaempferol-3-o-arabinoside (Juglanin)417.1284.034.9227.04, 255.03, 284.04, 285.04, 417.09
Table
Table 2. Summary of annotated central genes involved in flavonoids biosynthesis in C. rotundifolia and D. cambodiana as previously reported [40].
Table 2. Summary of annotated central genes involved in flavonoids biosynthesis in C. rotundifolia and D. cambodiana as previously reported [40].
EnzymeGene DesignationNo. of Annotated Sequences (C. rotundifolia)No. of Annotated Sequences (D. cambodiana)
Phenylalanine ammonia-lyase PAL 76
Cinnamic acid 4-hydroxylase C4H 41
4-Coumaric acid: CoA ligase 4CL 818
Chalcone synthase CHS 110
Chalcone isomerase CHI 76
Cinnamoyl-CoA CCR 2-
Flavanone 3-hydroxylase F3H 17
Flavonoid 3′,5′-hydroxylase F3′5′H 1-
Shikimate-o-hydroxycinnamoyltransferase HCT 6-
Flavonol synthase FLS 510
Dihydroflavanol 4-reductase DFR 116
Leucoanthocyanidin dioxygenase/anthocyanidin synthase LDOX/ANS 2-
Anthocyanidin reductase ANR 1-
Leucoanthocyanidin reductase LAR 61
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gichuki, D.K.; Li, Q.; Hou, Y.; Liu, Y.; Ma, M.; Zhou, H.; Xu, C.; Zhu, Z.; Wang, L.; Musila, F.M.; et al. Characterization of Flavonoids and Transcripts Involved in Their Biosynthesis in Different Organs of Cissus rotundifolia Lam. Metabolites 2021, 11, 741. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo11110741

AMA Style

Gichuki DK, Li Q, Hou Y, Liu Y, Ma M, Zhou H, Xu C, Zhu Z, Wang L, Musila FM, et al. Characterization of Flavonoids and Transcripts Involved in Their Biosynthesis in Different Organs of Cissus rotundifolia Lam. Metabolites. 2021; 11(11):741. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo11110741

Chicago/Turabian Style

Gichuki, Duncan Kiragu, Qingyun Li, Yujun Hou, Yuanshuang Liu, Mengxue Ma, Huimin Zhou, Chen Xu, Zhenfei Zhu, Lina Wang, Fredrick Mutie Musila, and et al. 2021. "Characterization of Flavonoids and Transcripts Involved in Their Biosynthesis in Different Organs of Cissus rotundifolia Lam" Metabolites 11, no. 11: 741. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo11110741

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop