Transcriptomic Analysis for Diurnal Temperature Differences Reveals Gene-Regulation-Network Response to Accumulation of Bioactive Ingredients of Protocorm-like Bodies in Dendrobium officinale

Abstract

Dendrobium officinale Kimura et Migo (D. officinale) is one of the most important traditional Chinese medicinal herbs, celebrated for its abundant bioactive ingredients. This study demonstrated that the diurnal temperature difference (DIF) (T1: 13/13 °C, T2: 25/13 °C, and T3: 25/25 °C) was more favorable for high chlorophyll, increased polysaccharide, and total flavonoid contents compared to constant temperature treatments in D. officinale PLBs. The transcriptome analysis revealed 4251, 4404, and 4536 differentially expressed genes (DEGs) in three different comparisons (A: 25/13 °C vs. 13/13 °C, B: 13/13 °C vs. 25/25 °C, and C: 25/13 °C vs. 25/25 °C, respectively). The corresponding up-/down-regulated DEGs were 1562/2689, 2825/1579, and 2310/2226, respectively. GO and KEGG enrichment analyses of DEGs showed that the pathways of biosynthesis of secondary metabolites, carotenoid biosynthesis, and flavonoid biosynthesis were enriched in the top 20; further analysis of the sugar- and flavonol-metabolism pathways in D. officinale PLBs revealed that the DIF led to a differential gene expression in the enzymes linked to sugar metabolism, as well as to flavonol metabolism. Certain key metabolic genes related to ingredient accumulation were identified, including those involved in polysaccharide metabolism (SUS, SUT, HKL1, HGL, AMY1, and SS3) and flavonol (UGT73C and UGT73D) metabolism. Therefore, these findings indicated that these genes may play an important role in the regulatory network of the DIF in the functional metabolites of D. officinale PLBs. In a MapMan annotation of abiotic stress pathways, the DEGs with significant changes in their expression levels were mainly concentrated in the heat-stress pathways, including heat-shock proteins (HSPs) and heat-shock transcription factors (HSFs). In particular, the expression levels of HSP18.2, HSP70, and HSF1 were significantly increased under DIF treatment, which suggested that HSF1, HSP70 and HSP18.2 may respond to the DIF. In addition, they can be used as candidate genes to study the effect of the DIF on the PLBs of D. officinale. The results of our qPCR analysis are consistent with those of the transcriptome-expression analysis, indicating the reliability of the sequencing. The results of this study revealed the transcriptome mechanism of the DIF on the accumulation of the functional metabolic components of D. officinale. Furthermore, they also provide an important theoretical basis for improving the quality of D. officinale via the DIF in production.