Elsevier

Process Biochemistry

Volume 91, April 2020, Pages 339-351
Process Biochemistry

Combined metabolomic and transcriptomic analysis reveals key candidate genes involved in the regulation of flavonoid accumulation in Anoectochilus roxburghii

https://doi.org/10.1016/j.procbio.2020.01.004Get rights and content

Highlights

  • A. roxburghii was first characterized using high-throughput sequencing technologies and metabolomics.

  • Correlation analysis between metabolites and regulatory genes revealed flavonoid accumulation-related regulatory genes.

  • This study screened and obtained key genes potentially involved in the accumulation of flavonoids in A. roxburghii.

  • This study provides an in-depth understanding of the molecular mechanisms underlying flavonoid accumulation in A. roxburghii.

Abstract

Anoectochilus roxburghii, a species used in Chinese herbal medicine, has unique characteristics. This plant has important medicinal and ornamental value and is distributed primarily in China, including Fujian, Zhejiang, Jiangxi, and Guizhou provinces, and in Taiwan. Flavonoids are involved in leaf pigment formation and are major pharmacodynamic substances. However, the molecular mechanisms that regulate accumulation of flavonoids remain unclear, which has significantly limited their application. To elucidate the molecular mechanisms underlying the regulation of flavonoid accumulation in A. roxburghii, root, stem and leaf samples were collected for constructing transcriptomic and metabolomic datasets using RNA sequencing and liquid chromatography-mass spectrometry (LC–MS) techniques. Sequencing of the transcriptomes of the different organs generated 60.2 Gb of data, which were assembled into 186,865 unigenes. Metabolomic analysis resulted in the identification of 10,690 metabolites. Based on the ESI+ mode, 4,010, 4,008 and 4,013 metabolites were annotated in roots, stems, and leaves, respectively; based on the ESI– mode, 1,530, 1,530, and 1,531 metabolites were annotated, respectively. Differential analyses of the transcriptome and flavonoid metabolism of the different organs revealed the greatest significant differences between roots and leaves, followed by differences between stems and leaves; differences between roots and stems were the smallest. According to analysis of differentially expressed genes (DEGs), the secondary metabolism-related DEGs between roots and stems, between roots and leaves, and between leaves and stems were classified into 16, 21, and 19 secondary metabolic Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, respectively. Among these pathways, flavonoid biosynthesis, flavone and flavonol biosynthesis, and isoflavonoid biosynthesis significantly differed in terms of both gene expression and secondary metabolism. Moreover, the key genes involved in flavonoid accumulation were comprehensively analyzed using metabolic and transcriptomic data, and ten transcription factor-encoding genes and fourteen flavonoid biosynthesis genes were identified. Specifically, four of the ten transcription factor-encoding genes appear to activate CHS, CHI, CYP75A and ANR gene expression. The combined approach used in this study provides an in-depth understanding of the molecular regulatory mechanisms underlying flavonoid accumulation in A. roxburghii and is a powerful tool that can uncover valuable information for plant breeders.

Introduction

Anoectochilus roxburghii (Wall.) Lindl is a perennial evergreen plant species of the Orchidaceae family. A. roxburghii has unique characteristics and value for use in Chinese herbal medicine; indeed, it is often used as an important medicinal and health care product in China and other Asian countries [1]. This plant is known as “the king of medicine” because of its multiple pharmacological effects and is applied for treating liver diseases, cardiovascular diseases, diabetes, cancer, and nephritis [[1], [2], [3], [4]]. Because it is a small, exquisite plant with gracefully shaped leaves and networks of golden yellow veins, A. roxburghii also has high ornamental value and is prized for its indoor foliage [2]. Owing to its important medicinal and ornamental value, A. roxburghii is receiving increased attention.

The secondary metabolites of A. roxburghii, especially flavonoids, are involved in leaf pigment formation and are a major source of the pharmacological activities of this species [5,6]. Therefore, it is important to study the molecular regulatory mechanisms underlying flavonoid accumulation in A. roxburghii. Furthermore, the contents of flavonols, such as quercetin, kaempferol and isorhamnetin, are often used as indicators to measure the quality of A. roxburghii [1,6]. Additionally, A. roxburghii is rich in flavonoid compounds, including dihydroquercetin, quercetin, kaempferol, and myricetin [2,6], and differences in the contents of quercetin, kaempferol, myricetin, and total flavonoids can affect pigment formation [6,7].

In recent years, long-term harvesting of A. roxburghii has led to a scarcity of wild resources of this species, pushing it to the verge of extinction. Although wild A. roxburghii was once the only source, the success of artificial cultivation techniques has revolutionized the procurement of the plant [8,9]. Recent studies regarding the medicinal value of A. roxburghii have focused primarily on resource collection of varieties, quality evaluation, active ingredient analysis, and pharmacological effects. Conversely, few studies have investigated the molecular genetics of A. roxburghii, including molecular identification and genetic diversity analyses. Genomic information concerning A. roxburghii is also not available, which has been the primary bottleneck for molecular genetics research in this species [2,8,9]. However, there has been marked progress in plant functional genomics in recent years. In general, transcriptomics is widely used to elucidate the biosynthesis pathways and regulatory mechanisms of key drug-related metabolites in different medicinal plant species [10,11]. Additionally, metabolomic analyses of medicinal plants have greatly facilitated the identification of the metabolic pathways of active pharmaceutical compounds [12]. Hence, integration of transcriptomics and metabolomics has been widely applied to identify the accumulation mechanisms underlying key metabolic pathways, especially the biosynthetic mechanisms of nonmodel plant species [11,13].

In this study, a transcriptomic database of A. roxburghii organs was constructed using a sequencing platform and bioinformatic analysis. Analysis of the expression profiles of different organs utilizing this database will aid in the identification of functional genes in A. roxburghii, laying the foundation for molecular genetics investigations in this species. To date, the molecular regulatory mechanisms of flavonoid accumulation in A. roxburghii affecting leaf pigment production and drug efficacy remain unclear, yet the metabolic pathways and metabolites in different organs, including roots, stems, and leaves, can be analyzed via metabolomics to identify differences in flavonoid accumulation. The data obtained in this study provide important technical support for future research on the accumulation and metabolic regulation of flavonoids in different organs of A. roxburghii and will assist in the identification of other active drug compounds.

Section snippets

Sample collection

Sample materials were collected from Nanjing, Zhangzhou, Fujian Province, China, and identified as A. roxburghii. In March 2016, A. roxburghii was cultivated in Nanjing according to the underwood planting method. The planting site is located at 24′ 26″ north, 117′ 36″ east, at approximately 300 m above sea level, and is characterized as a southern Fujian hilly landform. The annual average temperature is 20.4 °C, and the annual rainfall is 2,001.3 mm. The forest is dominated by evergreen

Construction and functional analysis of the transcriptome libraries of A. roxburghii roots, stems and leaves

The transcriptomes of the root, stem, and leaf libraries of A. roxburghii were sequenced, resulting in 60.2 Gb of data. After raw data were filtered and removed, the average sizes of clean reads of the transcriptomic data were 44.60, 44.58, and 44.61 Mb for roots, stems, and leaves, accounting for 78.87 %, 78.68 %, and 77.24 % of the total reads, respectively. The Q20 (%) clean reads accounted for 97.99 %, 97.70 %, and 97.69 % of the clean reads from the root, stem, and leaf libraries,

Discussion

Transcription factors are involved in many biological processes, and their gene expression is likely related to different traits of different organs; however, flavonoid accumulation is also involved in many biological processes. Therefore, genes related to the different traits of different organs might also be linked to flavonoid accumulation processes [24,[1], [2], [3], [4], [5], [6], [7], [8], [9], [10], [11], [12], [13], [14], [15], [16], [17], [18], [19], [20], [21], [22], [23], [24], [25],

Conclusion

In this study, A. roxburghii was for the first time characterized by high-throughput sequencing and metabolomic techniques. In total, 186,865 unigenes and 10,690 metabolites were identified by transcriptomic and metabolomic approaches. By comprehensively analyzing metabolomic and transcriptomic data, we investigated the regulatory network underlying flavonoid accumulation in different A. roxburghii organs. Correlation analysis between the metabolites and genes revealed regulatory genes related

Author contributions statement

SL obtained the funding and designed and supervised all the experiments. YC performed most of the experiments, analyzed the results, and wrote the article. JS performed some of the experiments. WP critically reviewed the manuscript. All authors read and approved the final manuscript.

Funding

This work was financially supported by grants from the Department of Technology and Science of Fujian Provincial Government (Grant No. 2017J01614), the Fujian Department of Finance Scientific Research Project (Grant No. k81139231) and the Education Department Project of Fujian Province (Grant No. JAT160169).

Declaration of Competing Interest

The authors declare that they have no competing interests.

Acknowledgments

The authors would like to acknowledge the staff of the State Forestry Administration Engineering Research Center of Chinese Fir for their excellent technical support.

References (57)

  • W. Huang et al.

    Isolation and molecular characterisation of flavonoid 3’-hydroxylase and flavonoid 3’, 5’-hydroxylase genes from a traditional Chinese medicinal plant, Epimedium sagittatum

    Gene

    (2012)
  • G. Forkmann et al.

    Metabolic engineering and applications of flavonoids

    Curr. Opin. Biotechnol.

    (2001)
  • K. Matsui et al.

    Buckwheat R2R3 MYB transcription factor FeMYBF1 regulates flavonol biosynthesis

    Plant Sci.

    (2018)
  • X. Wang et al.

    Simultaneous Structural Identification of Natural Products in Fractions of Crude Extract of the Rare Endangered Plant Anoectochilus roxburghii Using 1H NMR/RRLC-MS Parallel Dynamic Spectroscopy

    Int. J. Mol. Sci.

    (2011)
  • W. Wang et al.

    Effects of supplemental lighting with different light qualities on growth and secondary metabolite content of Anoectochilus roxburghii

    PeerJ

    (2018)
  • S. Ye et al.

    Effects of Light Quality on Morphology, Enzyme Activities, and Bioactive Compound Contents in Anoectochilus roxburghii

    Front. Plant Sci.

    (2017)
  • C.W. Yu et al.

    The complete chloroplast genome sequence of Anoectochilus roxburghii

    Mitochondrial DNA A DNA Mapp. Seq. Anal.

    (2016)
  • T. Lv et al.

    DNA barcodes for the identification of Anoectochilus roxburghii and its adulterants

    Planta

    (2015)
  • A. Rai et al.

    High-throughput sequencing and de novo transcriptome assembly of Swertia japonica to identify genes involved in the biosynthesis of therapeutic metabolites

    Plant Cell Rep.

    (2016)
  • W. Gao et al.

    Combining metabolomics and transcriptomics to characterize tanshinone biosynthesis in Salvia miltiorrhiza

    BMC Genomics

    (2014)
  • M. Zhu et al.

    Current advances in the metabolomics study on Lotus seeds

    Front. Plant Sci.

    (2016)
  • K. Cho et al.

    Network analysis of the metabolome and transcriptome reveals novel regulation of potato pigmentation

    J. Exp. Bot.

    (2016)
  • S. Ka et al.

    HLAscan: genotyping of the HLA region using next-generation sequencing data

    BMC Bioinformatics

    (2017)
  • S.C. Niu et al.

    De novo transcriptome assembly databases for the butterfly orchid Phalaenopsis equestris

    Sci. Data

    (2016)
  • Z. Chu et al.

    De novo assembly and comparative analysis of the transcriptome of embryogenic callus formation in bread wheat (Triticum aestivum L.)

    BMC Plant Biol.

    (2017)
  • L. Wang et al.

    DEGseq: an R package for identifying differentially expressed genes from RNA-seq data

    Bioinformatics

    (2010)
  • H. Li et al.

    De novo assembly and discovery of genes that are involved in drought tolerance in Tibetan Sophora moorcroftiana

    PLoS One

    (2015)
  • H. Ghosson et al.

    Simultaneous untargeted and targeted metabolomics profiling of underivatized primary metabolites in sulfur-deficient barley by ultra-high performance liquid chromatography-quadrupole/time-of-flight mass spectrometry

    Plant Methods

    (2018)
  • Cited by (11)

    • Combined analysis of transcriptomics and metabolomics revealed complex metabolic genes for diterpenoids biosynthesis in different organs of Anoectochilus roxburghii

      2023, Chinese Herbal Medicines
      Citation Excerpt :

      Therefore, the research on the diterpenoids biosynthesis of A. roxburghii is beneficial to further explore its formation mechanism of medicinal value. In recent years, with the development of biological technology in the application of medicinal plants, the combination of transcriptomics and metabolomics analysis method has been widely used in the research of biosynthesis and regulatory mechanisms of active ingredients in medicinal plants (Chen, Yao, et al., 2020; Chen, Pan, Jin, & Lin, 2020; Gong et al., 2020; Lin et al., 2019; Nett, Lau, & Sattely, 2020; Zhang, Li, Chen, Meng, & Guo, 2020). Nett et al. established a metabolic route of tropolone-containing colchicine alkaloids by transcriptomics, metabolomics analysis and pathway reconstruction, which almost completely elucidated the biosynthetic pathway of colchicine (Nett, Lau, & Sattely, 2020).

    View all citing articles on Scopus
    View full text