[go: up one dir, main page]

 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (24)

Search Parameters:
Keywords = Dendrobium catenatum

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
17 pages, 5193 KiB  
Article
Whole-Genome Analysis of ZF-HD Genes among Three Dendrobium Species and Expression Patterns in Dendrobium chrysotoxum
by Xin He, Xuewei Zhao, Qinyao Zheng, Meng-Meng Zhang, Ye Huang, Zhong-Jian Liu and Siren Lan
Horticulturae 2024, 10(6), 610; https://doi.org/10.3390/horticulturae10060610 - 8 Jun 2024
Viewed by 395
Abstract
ZF-HD transcription factors, which are unique to land plants, are involved in the regulation of abiotic stress response and related signaling pathways, and play a crucial role in plant growth and development. Dendrobium is one of the largest genera of orchids, with a [...] Read more.
ZF-HD transcription factors, which are unique to land plants, are involved in the regulation of abiotic stress response and related signaling pathways, and play a crucial role in plant growth and development. Dendrobium is one of the largest genera of orchids, with a high ornamental and ecological value. However, the specific functions of the ZF-HDs in Dendrobium remain unknown. In this study, we identified a total of 53 ZF-HDs from D. chrysotoxum (17), D. catenatum (23), and D. huoshanense (13), and analyzed their physicochemical properties, phylogenetic relationships, chromosomal locations, protein structures, conserved motifs, and expression patterns. The phylogenetic relationships revealed that 53 ZF-HDs were classified into six subfamilies (ZHDI–V and MIF), and all ZF-HD proteins contained motif 1 and motif 4 conserved domains, while a minority of these proteins had exons. The analysis of cis-elements in the promoters of ZF-HDs from three Dendrobium species showed that growth- and development-related elements were the most prevalent, followed by hormone response and abiotic stress response elements. Through collinearity analysis, 14 DchZF-HDs were found to be collinear with DhuZF-HDs, and 12 DchZF-HDs were found to be collinear with DcaZF-HDs. Furthermore, RT-qPCR analysis revealed that DchZF-HDs play a regulatory role in the development of lateral organs during the flowering process. The results indicated that DchZHD2 plays a role in the unpigmented bud stage, while DchMIF8 and DchZHD16 play significant roles during the pigmented bud and initial bloom stages. Hence, this study provides a crucial basis for further exploring ZF-HDs functions in regulating the floral organs of orchids. Full article
Show Figures

Figure 1

Figure 1
<p>Phylogenetic tree generated for the 70 <span class="html-italic">ZF-HD</span>s retrieved from <span class="html-italic">D. chrysotoxum</span>, <span class="html-italic">D. catenatum</span>, <span class="html-italic">D. huoshanense</span>, and <span class="html-italic">A. thaliana</span>.</p>
Full article ">Figure 2
<p>The motif and gene structure of <span class="html-italic">ZF-HD</span>s in <span class="html-italic">D. chrysotoxum</span>, <span class="html-italic">D. catenatum</span>, <span class="html-italic">D. huoshanense</span>, and <span class="html-italic">A. thaliana</span>. (<b>A</b>) Phylogenetic tree of 70 <span class="html-italic">ZF-HD</span>s constructed using MEGA7.0. (<b>B</b>) The conserved motif of ZF-HD proteins. (<b>C</b>) The ZF-HDs structure. (<b>D</b>) The sequence information for motif 1 and motif 4 (* indicates highly conserved cysteine residues).</p>
Full article ">Figure 3
<p>The <span class="html-italic">ZF-HD</span>s distribution on the chromosomes of three <span class="html-italic">Dendrobiums</span>. (<b>A</b>) <span class="html-italic">D. chrysotoxum</span>, (<b>B</b>) <span class="html-italic">D. catenatum</span>, and (<b>C</b>) <span class="html-italic">D. huoshanense</span>. Black designates the name of chromosomes, and red is employed to indicate the name of <span class="html-italic">ZF-HD</span>s.</p>
Full article ">Figure 4
<p>Collinearity analysis of <span class="html-italic">ZF-HD</span>s in three <span class="html-italic">Dendrobiums</span>.</p>
Full article ">Figure 5
<p>Analysis of <span class="html-italic">cis</span>-element components of the promoter of ZF-HDs across three <span class="html-italic">Dendrobiums</span>. (<b>A</b>) <span class="html-italic">D. chrysotoxum</span>, (<b>B</b>) <span class="html-italic">D. huoshanense</span>, and (<b>C</b>) <span class="html-italic">D. catenatum</span>. The left segment represents the count of <span class="html-italic">cis</span>-elements in <span class="html-italic">ZF-HDs</span>, while the right segment provides the statistical analysis of different categories of <span class="html-italic">ZF-HDs</span>. The types and numbers of <span class="html-italic">ZF-HDs</span> are listed in <a href="#app1-horticulturae-10-00610" class="html-app">Table S2</a>.</p>
Full article ">Figure 6
<p>The expression levels of 17 <span class="html-italic">ZF-HD</span> genes during various developmental periods (S1: unpigmented bud stage; S2: pigmented bud stage; S3: initial bloom stage) and in distinct components (Pe: petals; Se: sepals; Ov: ovary, Lip: lip; Gy: gynostemium) of <span class="html-italic">D. chrysotoxum</span>. The FPKM values of <span class="html-italic">ZF-HD</span>s observed in <span class="html-italic">D. chrysotoxum</span> are detailed in <a href="#app1-horticulturae-10-00610" class="html-app">Table S5</a>.</p>
Full article ">Figure 7
<p>Real-time reverse transcription quantitative PCR (RT–qPCR) validation of <span class="html-italic">DchZF-HD</span>s. S1 represents the unpigmented bud stage, S2 represents the pigmented bud stage, and S3 represents the initial bloom stage. The Y-axis signifies the relative expression values (2<sup>−ΔΔCT</sup>). Additionally, the red asterisk serves as an indicator of the significance level of the P-value in the respective test (** <span class="html-italic">p</span> &lt; 0.01, **** <span class="html-italic">p</span> &lt; 0.0001). Primers and RT-qPCR analysis of <span class="html-italic">DchZF-HD</span>s are shown in <a href="#app1-horticulturae-10-00610" class="html-app">Table S6</a>.</p>
Full article ">
17 pages, 10310 KiB  
Article
Identification and Analysis of PEPC Gene Family Reveals Functional Diversification in Orchidaceae and the Regulation of Bacterial-Type PEPC
by Ruyi Li, Xuyong Gao, Yuwei Wu, Chunyi Wei, Ming-He Li, Ding-Kun Liu and Zhong-Jian Liu
Int. J. Mol. Sci. 2024, 25(4), 2055; https://doi.org/10.3390/ijms25042055 - 8 Feb 2024
Cited by 1 | Viewed by 1122
Abstract
Phosphoenolpyruvate carboxylase (PEPC) gene family plays a crucial role in both plant growth and response to abiotic stress. Approximately half of the Orchidaceae species are estimated to perform CAM pathway, and the availability of sequenced orchid genomes makes them ideal subjects for investigating [...] Read more.
Phosphoenolpyruvate carboxylase (PEPC) gene family plays a crucial role in both plant growth and response to abiotic stress. Approximately half of the Orchidaceae species are estimated to perform CAM pathway, and the availability of sequenced orchid genomes makes them ideal subjects for investigating the PEPC gene family in CAM plants. In this study, a total of 33 PEPC genes were identified across 15 orchids. Specifically, one PEPC gene was found in Cymbidium goeringii and Platanthera guangdongensis; two in Apostasia shenzhenica, Dendrobium chrysotoxum, D. huoshanense, Gastrodia elata, G. menghaiensis, Phalaenopsis aphrodite, Ph. equestris, and Pl. zijinensis; three in C. ensifolium, C. sinense, D. catenatum, D. nobile, and Vanilla planifolia. These PEPC genes were categorized into four subgroups, namely PEPC-i, PEPC-ii, and PEPC-iii (PTPC), and PEPC-iv (BTPC), supported by the comprehensive analyses of their physicochemical properties, motif, and gene structures. Remarkably, PEPC-iv contained a heretofore unreported orchid PEPC gene, identified as VpPEPC4. Differences in the number of PEPC homolog genes among these species were attributed to segmental duplication, whole-genome duplication (WGD), or gene loss events. Cis-elements identified in promoter regions were predominantly associated with light responsiveness, and circadian-related elements were observed in each PEPC-i and PEPC-ii gene. The expression levels of recruited BTPC, VpPEPC4, exhibited a lower expression level than other VpPEPCs in the tested tissues. The expression analyses and RT-qPCR results revealed diverse expression patterns in orchid PEPC genes. Duplicated genes exhibited distinct expression patterns, suggesting functional divergence. This study offered a comprehensive analysis to unveil the evolution and function of PEPC genes in Orchidaceae. Full article
Show Figures

Figure 1

Figure 1
<p>Phylogenetic trees of PEPC proteins based on 33 orchid PEPC proteins and 4 AtPEPC proteins. The different sizes of the circle on the nodes represent bootstrap percentages; the upper left corner shows the phylogenetic topology structure. PEPC-i, PEPC-ii, PEPC-iii, and PEPC-iv are grouped together as indicated in pale red, pink, blue, and green, respectively. <span class="html-italic">A. shenzhenica</span>, <span class="html-italic">C. ensifolium</span>, <span class="html-italic">C. goeringii</span>, <span class="html-italic">C. sinense, D. catenatum, D. chrysotoxum</span>, <span class="html-italic">D. huoshanense</span>, <span class="html-italic">D. nobile</span>, <span class="html-italic">G. elata</span>, <span class="html-italic">G. menghaiensis</span>, <span class="html-italic">Ph aphrodite</span>, <span class="html-italic">Ph. equestris</span>, <span class="html-italic">Pl guangdongensis</span>, <span class="html-italic">Pl. zijinensis</span>, and <span class="html-italic">V. planifolia</span> are labeled as As, Ce, Cg, Cs, Dc, De, Dh, Dn, Ge, Gm, Pa, Pe, Pg, Pz, and Vp, respectively.</p>
Full article ">Figure 2
<p>Gene structure, conserved motifs, and domains of <span class="html-italic">PEPC</span>s. (<b>A</b>) The NJ tree contains 33 orchid <span class="html-italic">PEPC</span>s. (<b>B</b>) Squares of different colors represent conserved motifs of <span class="html-italic">PEPC</span>s. (<b>C</b>) Squares of different colors represent the gene structures of <span class="html-italic">PEPC</span>s.</p>
Full article ">Figure 3
<p><span class="html-italic">Cis</span>-acting elements in the promoter regions of <span class="html-italic">PEPC</span> genes. (<b>A</b>) Functions of <span class="html-italic">cis</span>-acting elements in different orchid <span class="html-italic">PEPC</span>s. (<b>B</b>) Number of <span class="html-italic">cis</span>-acting elements in different orchid <span class="html-italic">PEPC</span>s.</p>
Full article ">Figure 4
<p>Chromosomal localization of <span class="html-italic">PEPC</span> genes. The scale bars on the left are the length (Mb) of the chromosomes of each orchid. Red lines connect the possible duplicated genes. (<b>A</b>) <span class="html-italic">A. shenzhenica</span>; (<b>B</b>) <span class="html-italic">C. ensifolium</span>; (<b>C</b>) <span class="html-italic">D. huoshanense</span>; (<b>D</b>) <span class="html-italic">C. goeringii</span>; (<b>E</b>) <span class="html-italic">D. chrysotoxum</span>; (<b>F</b>) <span class="html-italic">C. sinense</span>; (<b>G</b>) <span class="html-italic">G. menghaiensis</span>; (<b>H</b>) <span class="html-italic">Pl. zijinensis</span>; (<b>I</b>) <span class="html-italic">G. elata</span>; (<b>J</b>) <span class="html-italic">D. catenatum</span>; (<b>K</b>) <span class="html-italic">Ph. equestris</span>; (<b>L</b>) <span class="html-italic">Pl. guangdongensis</span>; (<b>M</b>) <span class="html-italic">Ph. aphrodite</span>; (<b>N</b>) <span class="html-italic">D. nobile</span>; (<b>O</b>) <span class="html-italic">V. planifolia</span>.</p>
Full article ">Figure 5
<p>Location and orthologs or paralogs of 38 <span class="html-italic">PEPC</span> gene pairs intra- or inter-species of <span class="html-italic">G. menghaiensis</span>, <span class="html-italic">V. planifolia</span>, <span class="html-italic">C. ensifolium</span>, <span class="html-italic">D. chrysotoxum</span>, <span class="html-italic">D. huoshanense</span>, and <span class="html-italic">D. nobile</span> genomes. The chromosomes of <span class="html-italic">G. menghaiensis</span>, <span class="html-italic">V. planifolia</span>, <span class="html-italic">C. ensifolium</span>, <span class="html-italic">D. chrysotoxum</span>, <span class="html-italic">D. huoshanense</span>, and <span class="html-italic">D. nobile</span> were shown with different colors and labeled as Gm, Vp, Ce, De, Dh, and Dn, respectively. The gene pairs among different species are shown with different colors.</p>
Full article ">Figure 6
<p>Temporal expression patterns of <span class="html-italic">PEPC</span> genes in orchids. Gene expression levels were quantified in FPKM. (<b>A</b>) <span class="html-italic">Ph. equestris</span>; (<b>B</b>) <span class="html-italic">V. planifolia</span>.</p>
Full article ">Figure 7
<p>Spatial expression patterns of <span class="html-italic">PEPC</span> genes in orchids. Gene expression levels (FPKM) were normalized by log2 scale. (<b>A</b>) <span class="html-italic">A. shenzhenica</span>; (<b>B</b>) <span class="html-italic">D. catenatum</span>; (<b>C</b>) <span class="html-italic">V. planifolia</span>, mesocarp, placental, and hairs in three months pod; (<b>D</b>) <span class="html-italic">Ph. aphrodite</span>, flower bud: 2.5 cm, stalk: 1.5–3 cm.</p>
Full article ">Figure 8
<p>RT-qPCR analysis of three <span class="html-italic">DcPEPC</span>s genes in <span class="html-italic">D. catenatum</span> in different tissues (stem, leaf, and root). Error bars indicate the SD of three biological replicates.</p>
Full article ">
23 pages, 6599 KiB  
Article
Whole-Transcriptome Sequencing Reveals the Global Molecular Responses and NAC Transcription Factors Involved in Drought Stress in Dendrobium catenatum
by Siqi Zhang, Yuliang Han, Qinzong Zeng, Chenchang Wang, Huizhong Wang, Juncheng Zhang, Maohong Cai, Jiangjie Lu and Tao Chen
Antioxidants 2024, 13(1), 94; https://doi.org/10.3390/antiox13010094 - 12 Jan 2024
Viewed by 1115
Abstract
Dendrobium catenatum is a highly drought-tolerant herb, which usually grows on cliffs or in the branches of trees, yet the underlying molecular mechanisms for its tolerance remain poorly understood. We conducted a comprehensive study utilizing whole-transcriptome sequencing approaches to investigate the molecular response [...] Read more.
Dendrobium catenatum is a highly drought-tolerant herb, which usually grows on cliffs or in the branches of trees, yet the underlying molecular mechanisms for its tolerance remain poorly understood. We conducted a comprehensive study utilizing whole-transcriptome sequencing approaches to investigate the molecular response to extreme drought stress in D. catenatum. A large number of differentially expressed mRNAs, lncRNAs, and circRNAs have been identified, and the NAC transcription factor family was highly enriched. Meanwhile, 46 genes were significantly up-regulated in the ABA-activated signaling pathway. In addition to the 89 NAC family members accurately identified in this study, 32 members were found to have different expressions between the CK and extreme drought treatment. They may regulate drought stress through both ABA-dependent and ABA-independent pathways. Moreover, the 32 analyzed differentially expressed DcNACs were found to be predominantly expressed in the floral organs and roots. The ceRNA regulatory network showed that DcNAC87 is at the core of the ceRNA network and is regulated by miR169, miR393, and four lncRNAs. These investigations provided valuable information on the role of NAC transcription factors in D. catenatum’s response to drought stress. Full article
(This article belongs to the Section ROS, RNS and RSS)
Show Figures

Figure 1

Figure 1
<p>Well-watered plant (<b>A</b>) and drought-treated plant (<b>B</b>) of species <span class="html-italic">D. catenatum</span> (RJ) under 25 °C, 16 h of illumination, and 8 h of darkness conditions.</p>
Full article ">Figure 2
<p>Identification and analysis of differentially expressed mRNAs (DEmRNAs) in response to drought treatment. (<b>A</b>) Sample clustering analysis of three biological replicates of CK and drought-treated samples; (<b>B</b>) volcano plot displaying the identification of DEmRNAs in CK and drought treatment samples with significance plotted against fold change; (<b>C</b>) statistical analysis presenting the number of up- and down-regulated DEmRNAs; (<b>D</b>) heat map showing the expression pattern of DEmRNAs in three replicates of both CK and drought-treated samples; (<b>E</b>) Gene Ontology (GO) classifications of the DEmRNAs.</p>
Full article ">Figure 3
<p>Analysis of DEmRNAs and NAC family in RJ. (<b>A</b>) Column diagram displaying the number of different types of transcription factors in DEmRNAs; (<b>B</b>) phylogenetic tree of the NAC family in <span class="html-italic">D. catenatum</span> was constructed using neighbor-joining method with 1000 bootstrap replications. The three branches were color-coded, and the DE<span class="html-italic">DcNAC</span> genes were labeled in black. (<b>C</b>) Heatmap of the expression patterns of differentially expressed <span class="html-italic">DcNAC</span> (DE<span class="html-italic">DcNAC</span>) gene in the control group and drought group under three biological replicates. The outer circle of the heat map represents the DE<span class="html-italic">DcNAC</span> evolutionary tree. (<b>D</b>) Gene Ontology (GO) classification.</p>
Full article ">Figure 4
<p>The expression patterns of DE<span class="html-italic">DcNAC</span> genes in different tissues of RJ.</p>
Full article ">Figure 5
<p>Expression patterns of 32 DE<span class="html-italic">DcNAC</span> genes under ABA and NaCl treatment. The purple, green, and blue columns represented CK, ABA treatment, and NaCl treatment, respectively. <span class="html-italic">DcACTIN2</span> was used as internal control. Values were presented as means ± SD (<span class="html-italic">n</span> = 3). (* <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> &lt; 0.01, Student’s <span class="html-italic">t</span>-test).</p>
Full article ">Figure 6
<p>ROS levels in RJ under ABA, drought, and NaCl treatments. (<b>A</b>) NBT staining for superoxide; (<b>B</b>) DAB staining for H<sub>2</sub>O<sub>2</sub>; (<b>C</b>) relative NBT staining intensities; (<b>D</b>) relative DAB staining intensities, (Student’s <span class="html-italic">t</span>-test; * <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> &lt; 0.01, <span class="html-italic">n</span> = 10); Bars = 1 cm.</p>
Full article ">Figure 7
<p>Identification and analysis of DElncRNAs under drought stress. (<b>A</b>) Venn diagram of DElncRNAs in CK and drought treatment samples; (<b>B</b>) the quantitative distribution analysis of various types of DElncRNAs; (<b>C</b>) the length distribution of DElncRNAs; (<b>D</b>) the distribution of exon numbers in DElncRNAs; (<b>E</b>) volcano plot of log2 FoldChange and <span class="html-italic">p</span>-values (padj) of DElncRNAs in CK and drought treatment samples; (<b>F</b>) heat map of DElncRNAs in CK and drought treatment samples.</p>
Full article ">Figure 8
<p>Identification and analysis of DEcircRNAs under drought stress. (<b>A</b>) Length distribution of DEcircRNAs; (<b>B</b>) distribution of exon numbers in DEcircRNAs; (<b>C</b>) distribution of predicted circRNAs in each sample; (<b>D</b>) quantitative distribution of different types of DEcircRNAs; (<b>E</b>) Gene Ontology (GO) classification of DEcircRNAs.</p>
Full article ">Figure 9
<p>Identification and analysis of miRNAs under CK condition and drought stress. (<b>A</b>) Sample clustering analysis of three biological replicates of CK and drought-treated samples; (<b>B</b>) length distribution of sRNAs; (<b>C</b>) distribution of predicted miRNAs in each sample; (<b>D</b>) quantitative distribution of different types of sRNAs; (<b>E</b>) top five expressed novel miRNAs in each treatment; (<b>F</b>) heat map of the expression patterns of DEmiRNAs in three replicates of both CK and drought-treated samples.</p>
Full article ">Figure 10
<p>Mini-ceRNA network constructed with 7 known drought-related DEmRNAs, DElncRNAs, DEmiRNAs in RJ (<b>A</b>), and mini-ceRNA network constructed with <span class="html-italic">DcNAC87</span> (<b>B</b>). Colored circles represented the DEmRNAs (purple), DEmiRNAs (blue), and DElncRNAs (green).</p>
Full article ">
16 pages, 9579 KiB  
Article
Metabolic Pathway Engineering Improves Dendrobine Production in Dendrobium catenatum
by Meili Zhao, Yanchang Zhao, Zhenyu Yang, Feng Ming, Jian Li, Demin Kong, Yu Wang, Peng Chen, Meina Wang and Zhicai Wang
Int. J. Mol. Sci. 2024, 25(1), 397; https://doi.org/10.3390/ijms25010397 - 28 Dec 2023
Cited by 1 | Viewed by 1225
Abstract
The sesquiterpene alkaloid dendrobine, widely recognized as the main active compound and a quality control standard of medicinal orchids in the Chinese Pharmacopoeia, demonstrates diverse biological functions. In this study, we engineered Dendrobium catenatum as a chassis plant for the production of dendrobine [...] Read more.
The sesquiterpene alkaloid dendrobine, widely recognized as the main active compound and a quality control standard of medicinal orchids in the Chinese Pharmacopoeia, demonstrates diverse biological functions. In this study, we engineered Dendrobium catenatum as a chassis plant for the production of dendrobine through the screening and pyramiding of key biosynthesis genes. Initially, previously predicted upstream key genes in the methyl-D-erythritol 4-phosphate (MEP) pathway for dendrobine synthesis, including 4-(Cytidine 5′-Diphospho)-2-C-Methyl-d-Erythritol Kinase (CMK), 1-Deoxy-d-Xylulose 5-Phosphate Reductoisomerase (DXR), 2-C-Methyl-d-Erythritol 4-Phosphate Cytidylyltransferase (MCT), and Strictosidine Synthase 1 (STR1), and a few downstream post-modification genes, including Cytochrome P450 94C1 (CYP94C1), Branched-Chain-Amino-Acid Aminotransferase 2 (BCAT2), and Methyltransferase-like Protein 23 (METTL23), were chosen due to their deduced roles in enhancing dendrobine production. The seven genes (SG) were then stacked and transiently expressed in the leaves of D. catenatum, resulting in a dendrobine yield that was two-fold higher compared to that of the empty vector control (EV). Further, RNA-seq analysis identified Copper Methylamine Oxidase (CMEAO) as a strong candidate with predicted functions in the post-modification processes of alkaloid biosynthesis. Overexpression of CMEAO increased dendrobine content by two-fold. Additionally, co-expression analysis of the differentially expressed genes (DEGs) by weighted gene co-expression network analysis (WGCNA) retrieved one regulatory transcription factor gene MYB61. Overexpression of MYB61 increased dendrobine levels by more than two-fold in D. catenatum. In short, this work provides an efficient strategy and prospective candidates for the genetic engineering of D. catenatum to produce dendrobine, thereby improving its medicinal value. Full article
Show Figures

Figure 1

Figure 1
<p>Multigene reconstruction to enhance dendrobine production. (<b>A</b>) Physical map of the multigene constructs used for integration and expression of the synthetic operons. The selected target genes are depicted as light blue boxes. Promotor <span class="html-italic">Prrn</span> is shown in dark yellow and terminator <span class="html-italic">TrbcL</span> in orange. The <span class="html-italic">HPTII</span> selectable marker gene for transformation is represented as a light yellow box. An intercistronic expression element (<span class="html-italic">IEE</span>) was put between <span class="html-italic">STR1</span> and <span class="html-italic">CYP94C1</span> operons to ensure the downstream cistron expression under the same promotor. (<b>B</b>) <span class="html-italic">Not</span> I-digestion analysis of pYLTAC380H-multigene (<span class="html-italic">EV</span>, <span class="html-italic">TG</span>, <span class="html-italic">FG</span>, and <span class="html-italic">SG</span>) constructs. M: DNA ladder marker. (<b>C</b>) Dendrobine content in <span class="html-italic">D. catenatum</span> leaves infiltrated with <span class="html-italic">Agrobacterium tumefaciens</span> carrying multigene constructs. There are three replicates for each sample. (<b>D</b>) The relative expression of each gene in a specific multigene construct. <span class="html-italic">EV</span>: empty vector; <span class="html-italic">TG</span>: two genes; <span class="html-italic">FG</span>: five genes; <span class="html-italic">SG</span>: seven genes. Asterisks indicate significance based on the Student’s <span class="html-italic">t</span>-test. * <span class="html-italic">p</span> ≤ 0.05; ** <span class="html-italic">p</span> ≤ 0.01; *** <span class="html-italic">p</span> ≤ 0.001; **** <span class="html-italic">p</span> ≤ 0.0001.</p>
Full article ">Figure 2
<p>Functional verification of <span class="html-italic">MCT</span> in dendrobine synthesis. (<b>A</b>) <span class="html-italic">MCT</span> expression was checked by qRT-PCR with samples being collected 6 h after infiltration (n = 3). (<b>B</b>) Transiently infiltrated <span class="html-italic">D. catenatum</span> leaves were harvested (5 dpi) for dendrobine measurement (n = 3). (<b>C</b>) CRISPRi was performed to knock down <span class="html-italic">MCT</span> expression. The kinase-dead version of Cas9 (dCas9) was used to block transcription in the promotor region of <span class="html-italic">MCT</span>. Empty vector without dCas9 served as the control. <span class="html-italic">MCT</span> knock-down was verified by qRT-PCR, with samples being collected 6 h after infiltration (n = 3). (<b>D</b>) Leave samples were collected at 5 dpi and subjected to dendrobine measurement (n = 3). ** <span class="html-italic">p</span> ≤ 0.01; *** <span class="html-italic">p</span> ≤ 0.001; **** <span class="html-italic">p</span> ≤ 0.0001.</p>
Full article ">Figure 3
<p>Generation and identification of multigene-transgenic plants. (<b>A</b>) Transgenic <span class="html-italic">D. catenatum</span> plantlets of <span class="html-italic">SG</span>-multigene (11-month-old) and <span class="html-italic">EV</span> control (13-month-old). (<b>B</b>) Analysis of the <span class="html-italic">hptII</span> transgene in <span class="html-italic">SG</span>-multigene and <span class="html-italic">EV</span> control transgenic <span class="html-italic">D. catematum</span> plantlets. DNA isolated from wild-type plants was used as the negative control “−”, while the <span class="html-italic">hptII</span> gene in the <span class="html-italic">EV</span> plasmid was used as the positive control “+”. <span class="html-italic">DcActin</span> was PCR amplified to demonstrate an equal amount of loading. <span class="html-italic">nptII</span> was amplified to avoid bacterial contamination. (<b>C</b>) <span class="html-italic">FPPS</span> expression was checked to demonstrate activation of the dendrobine synthesis pathway. <span class="html-italic">EV</span> transgenic <span class="html-italic">D. catenatum</span> served as the control (Ctrl). (<b>D</b>) <span class="html-italic">SG</span>-transgenic <span class="html-italic">D. catenatum</span> grown in pine-bark pots for 10 months. (<b>E</b>) Molecular characterization of the <span class="html-italic">hptII</span> transgene in <span class="html-italic">SG</span>-transgenic <span class="html-italic">D. catenatum</span> grown in pine-bark pots. (<b>F</b>) Molecular characterization of the <span class="html-italic">hptII</span> transgene in <span class="html-italic">SG</span>-transgenic <span class="html-italic">Arabidopsis</span>. (<b>G</b>) Expression analysis of individual genes in <span class="html-italic">SG</span>-transgenic <span class="html-italic">Arabidopsis</span> by qRT-PCR. Scale bar in (<b>A</b>) represents 1 cm. **** <span class="html-italic">p</span> ≤ 0.0001 represents significance.</p>
Full article ">Figure 4
<p><span class="html-italic">SG</span>-transgenic <span class="html-italic">Arabidopsis</span> tolerant to salinity stress. (<b>A</b>) Plant growth in response to various concentrations of NaCl (0, 100, 120, and 200 mM). (<b>B</b>) Plant growth in terms of fresh weight under 100 mM NaCl for four weeks. (<b>C</b>) Cell damage in terms of MDA release under 100 mM NaCl for four weeks. (<b>D</b>) Representative image showing the transgenic plants growing in earth-pot for one month. (<b>E</b>) Comparison of plant height for the transgenic plants growing in earth-pot for one month. Data are represented as means ± SE from three replicates. *** <span class="html-italic">p</span> ≤ 0.001 and **** <span class="html-italic">p</span> ≤ 0.0001 are of significance compared to <span class="html-italic">EV</span> controls. Scale bar represents 1 cm. The center line in (<b>E</b>) represents median.</p>
Full article ">Figure 5
<p><span class="html-italic">SG</span>-transgenic <span class="html-italic">Arabidopsis</span> was tolerant to drought stress. (<b>A</b>) Plant growth in response to different concentrations of PEG6000 (0, 250, 400, and 550 g/L). (<b>B</b>) Plant growth in terms of fresh weight under varied concentrations of PEG6000 for four weeks. (<b>C</b>) Cell damage in terms of MDA release under varied concentrations of PEG6000 for four weeks. * <span class="html-italic">p</span> ≤ 0.05; ** <span class="html-italic">p</span> ≤ 0.01; **** <span class="html-italic">p</span> ≤ 0.0001 represent significance. Scale bar represents 1 cm.</p>
Full article ">Figure 6
<p>Dendrobine synthesis-related gene screening by Venn diagram and KEGG analysis. (<b>A</b>) Venn distribution of DEGs for transcriptomes of different tissues (stem vs. root; leaf vs. root) from three <span class="html-italic">Dendrobium</span> species (<span class="html-italic">D. houshanense</span>; <span class="html-italic">D. catenatum</span>, and <span class="html-italic">D. moniliforme</span>). (<b>B</b>) KEGG pathway enrichment of 648 DEGs from (<b>A</b>). The x-axis represents the enrichment ratio and the y-axis represents the pathway name. (<b>C</b>) Venn diagram representation of the number of DEGs in samples from protocorm-like bodies (PLBs), samples under <span class="html-italic">MF23</span> treatment, and the 648 DEGs in (<b>A</b>). (<b>D</b>) KEGG pathway enrichment of DEGs from (<b>C</b>).</p>
Full article ">Figure 7
<p>Dendrobine synthesis-related hub gene screening by WGCNA. (<b>A</b>) Hierarchical cluster dendrogram showing six expression modules of co-expressed genes. Each leaf in the tree represents an individual gene, with the branch representing a module of highly connected genes. The designated color rows below correspond to module membership. (<b>B</b>) Scale-free fit index at different threshold values (<span class="html-italic">β</span>). Asterisk indicates the selected soft-thresholding power. (<b>C</b>) Heatmap of connectivity of eigengenes. (<b>D</b>) Module-trait correlations and corresponding <span class="html-italic">p</span>-values (in parenthesis). The color in the box indicates −log(<span class="html-italic">P</span>) and the color scale indicates the <span class="html-italic">p</span>-value from the Fisher exact test. Treatment means <span class="html-italic">MF23</span> infection.</p>
Full article ">Figure 8
<p>Functional verification of downstream genes in dendrobine synthesis. (<b>A</b>) <span class="html-italic">CMEAO</span> overexpression was verified by qRT-PCR. Samples were collected at 24 h post-infiltration. (<b>B</b>) Transiently infiltrated <span class="html-italic">D. catenatum</span> leaves (5 dpi) were harvested and subjected to dendrobine measurement. Empty vectors are used as controls. (<b>C</b>) qRT-PCR verification of <span class="html-italic">MYB61</span> overexpression in transiently infiltrated one-year-old <span class="html-italic">D. catenatum</span> leaves. (<b>D</b>) Dendrobine content in <span class="html-italic">D. catenatum</span> leaves transiently overexpressing <span class="html-italic">MYB61</span>. Statistical significance was demonstrated as following ** <span class="html-italic">p</span> ≤ 0.01, *** <span class="html-italic">p</span> ≤ 0.001, **** <span class="html-italic">p</span> ≤ 0.0001.</p>
Full article ">
14 pages, 8960 KiB  
Article
Genome-Wide Identification of Fatty Acyl-CoA Reductase (FAR) Genes in Dendrobium catenatum and Their Response to Drought Stress
by Yutong Ren, Peng Wang, Tingting Zhang, Wen Liu, Yujuan Wang, Jun Dai and Yang Zhou
Horticulturae 2023, 9(9), 982; https://doi.org/10.3390/horticulturae9090982 - 31 Aug 2023
Cited by 2 | Viewed by 1078
Abstract
Dendrobium catenatum is a high-value medicinal plant that is predominantly found in high mountain areas, thriving amidst cliffs and rock crevices. However, its wild resources face constant threats from adverse environmental conditions, especially drought stress. Fatty acyl-CoA reductase (FAR) is crucial in plant [...] Read more.
Dendrobium catenatum is a high-value medicinal plant that is predominantly found in high mountain areas, thriving amidst cliffs and rock crevices. However, its wild resources face constant threats from adverse environmental conditions, especially drought stress. Fatty acyl-CoA reductase (FAR) is crucial in plant drought resistance, but there is a lack of research on FAR genes in D. catenatum. In this study, the FAR family genes were identified from the D. catenatum genome. Their genomic characteristics were investigated using bioinformatics techniques, and their expression patterns in different tissues and under 20% PEG8000 conditions mimicking drought stress were analyzed using quantitative real-time RT-PCR (RT-qPCR). Seven DcFAR genes were identified from the D. catenatum genome. The encoded amino acids range between 377 and 587 aa, with molecular weights between 43.41 and 66.15 kD and isoelectric points between 5.55 and 9.02. Based on the phylogenetic relationships, the FAR family genes were categorized into three subgroups, each with similar conserved sequences and gene structures. The cis-acting elements of the promoter regions were assessed, and the results reveal that the DcFAR upstream promoter region contains multiple stress-related elements, suggesting its potential involvement in abiotic stress responses. The RT-qPCR results show distinct expression patterns of DcFAR genes in various plant tissues. It was observed that the expression of most DcFAR genes was upregulated under drought stress. Among them, the expression levels of DcFAR2, DcFAR3, DcFAR5, and DcFAR7 genes under drought stress were 544-, 193-, 183-, and 214-fold higher compared to the control, respectively. These results indicate that DcFAR2/3/5/7 might play significant roles in D. catenatum drought tolerance. This research offers insight into the function of DcFAR genes and provides theoretical support for breeding drought-resistant D. catenatum varieties. Full article
Show Figures

Figure 1

Figure 1
<p>Phylogenetic analysis and protein sequence comparison of FAR proteins. (<b>A</b>) Phylogenetic analysis of FAR proteins from <span class="html-italic">D. catenatum</span>, <span class="html-italic">Arabidopsis</span>, rice, and sorghum. The phylogenetic tree was constructed using the Neighbor-Joining (NJ) method in MEGA-X, with default parameters. Red star represents <span class="html-italic">D. catenatum</span> FAR (DcFAR); black triangle represents <span class="html-italic">Arabidopsis</span> FAR (AtFAR); green square represents sorghum FAR; and blue circle represents rice FAR (OsFAR). (<b>B</b>) Protein sequence comparison of FAR proteins using DNAMAN. The three red boxes represent the GXXGXX(G/A) and YXXXK conserved domains.</p>
Full article ">Figure 2
<p>Conserved motif and gene structure analyses of DcFAR family members. (<b>A</b>) Conserved motifs of DcFAR proteins. Rectangular boxes of different colors represent different conserved motifs. (<b>B</b>) Exon/intron structure of <span class="html-italic">DcFAR</span>. UTR(s), exon(s), and intron(s) are represented by green boxes, yellow boxes, and black lines, respectively.</p>
Full article ">Figure 3
<p><span class="html-italic">Cis</span>-acting element analysis in the promoter of <span class="html-italic">DcFAR</span> genes. (<b>A</b>) Locations of <span class="html-italic">cis</span>-acting elements in the promoter of <span class="html-italic">DcFAR</span>. Ellipses of different colors represent different types of <span class="html-italic">cis</span>-acting elements and their positions in each <span class="html-italic">DcFAR</span> gene promoter. (<b>B</b>) Statistics of the number of <span class="html-italic">cis</span>-acting elements in <span class="html-italic">DcFAR</span> promoters. Different colors and numbers represent the number of different <span class="html-italic">cis</span>-acting elements in each <span class="html-italic">DcFAR</span> promoter.</p>
Full article ">Figure 4
<p>Expression analysis of <span class="html-italic">DcFAR</span> genes under drought stress using RT-qPCR. The data are expressed as mean ± standard deviation (<span class="html-italic">n</span> = 3). Vertical bars represent the means of fold change in expression and standard deviations calculated from the replicates. Values of 0, 3, 6, 9, 12, 24, and 48 indicate hours after treatment. Asterisks (* or **) indicate a significant difference at <span class="html-italic">p</span> &lt; 0.05 or 0.01, respectively.</p>
Full article ">Figure 5
<p>Expression analysis of <span class="html-italic">DcFAR</span> genes in different tissues using RT-qPCR. Mean expression values were calculated from three independent biological replicates relative to the value in roots and visualized using TBtools. Green and red indicate low and high levels of expression, respectively. RO: root; ST: stem; LE: leaf; CA: capsule; E: petal; SE: sepal; FS: flower stalk; GY: gynostemia; LI: lip.</p>
Full article ">
16 pages, 10017 KiB  
Article
Genome-Wide Identification and Expression Analysis of the Chalcone Synthase (CHS) Gene Family in Dendrobium catenatum
by Ting Yang, Tingting Zhang, Yuxin Li, Yuqian Kang, Peng Wang, Wen Liu, Yujuan Wang, Libo Tian, Jun Dai and Yang Zhou
Agronomy 2023, 13(6), 1488; https://doi.org/10.3390/agronomy13061488 - 28 May 2023
Cited by 5 | Viewed by 1611
Abstract
Flavonoids are considered to be critical metabolites in regulating plant responses to abiotic stress. Chalcone synthase (CHS) catalyzes the first key step in the flavonoid biosynthesis pathway. However, there is no in-depth information on the sequence and functional analysis of CHS genes in [...] Read more.
Flavonoids are considered to be critical metabolites in regulating plant responses to abiotic stress. Chalcone synthase (CHS) catalyzes the first key step in the flavonoid biosynthesis pathway. However, there is no in-depth information on the sequence and functional analysis of CHS genes in Dendrobium catenatum. In this study, a total of 14 DcCHS genes were identified, encoding proteins of 349–504 amino acids in length, a protein molecular weight ranging between 39.08 and 56.56 kDa, and isoelectric points from 5.64 to 9.63. The DcCHS proteins were then divided into three groups according to their phylogenetic relationships. The members of each group had similar conserved motifs and gene structures. Quantitative real-time RT-PCR (RT-qPCR) analysis revealed that the DcCHS genes exhibited variable expression patterns in the different plant tissues evaluated. Furthermore, six genes were differentially expressed following exposure to abiotic stresses: DcCHS-6, DcCHS-5/-6, DcCHS-13/-14, and DcCHS-6/-8/-9/-13, which were specifically expressed in response to drought, heat, cold, and salt stress, respectively. This is the first genome-wide analysis of the CHS genes in D. catenatum, and our findings can provide essential information for a better understanding of the function of DcCHS genes, thus facilitating further research on D. catenatum stress tolerance. Full article
(This article belongs to the Section Horticultural and Floricultural Crops)
Show Figures

Figure 1

Figure 1
<p>Flavonoid biosynthetic pathway. Enzymatic reactions: CHS: chalcone synthase; CHI: chalcone isomerase; F3′H: flavonoid 3′-hydroxylase; ANR: anthocyanin reductase; FLS: flavonol synthase; DFR: dihydroflavonol 4-reductase; ANS: anthocyanidin synthase.</p>
Full article ">Figure 2
<p>Phylogenetic analyses of CHS proteins. MEGA-X with default parameters was used to construct a phylogenetic tree of CHS proteins. The three groups are shown in different colors. Pink stars indicate <span class="html-italic">Dendrobium catenatum</span> CHSs (DcCHSs); light blue squares indicate <span class="html-italic">Phalaenopsis equestris</span> CHSs (PeCHSs); green triangles represent <span class="html-italic">Apostasia shenzhenica</span> CHSs (AsCHSs); purple circles represent <span class="html-italic">Arabidopsis thaliana</span> CHSs (AtCHSs); blue checkmarks represent <span class="html-italic">Oryza sativa</span> CHSs (OsCHSs).</p>
Full article ">Figure 3
<p>Motifs, gene structures, and phylogenetic relationships of DcCHS family members. (<b>A</b>) Phylogenetic tree of 14 DcCHS proteins generated using Neighbor-Joining (NJ) method. (<b>B</b>) Conserved motifs of DcCHS proteins. (<b>C</b>) Exon/intron structures of <span class="html-italic">DcCHS</span> genes. UTR(s), exon(s), and intron(s) are indicated with green boxes, yellow boxes, and black lines, respectively.</p>
Full article ">Figure 4
<p>Alignments of the conserved domains of DcCHS proteins. (<b>A</b>) Alignments of the Chal_sti_synt_C and Chal_sti_synt_N domains from the MEME results for DcCHS. Motifs 1, 3, and 7 formed the putative Chal_sti_synt_N domain; motifs 2, 5, and 6 formed the putative Chal_sti_synt_C domain. (<b>B</b>) Alignment of the conserved domains of CHS proteins from <span class="html-italic">A. thaliana</span> and <span class="html-italic">D. catenatum</span>. The conserved amino acids are shown in the red box.</p>
Full article ">Figure 5
<p>Analysis of <span class="html-italic">cis</span>-elements in the <span class="html-italic">DcCHS</span> genes’ promoter regions. (<b>A</b>) Phylogenetic tree of 14 DcCHS proteins was generated using the NJ method. (<b>B</b>) Different types of <span class="html-italic">cis</span>-elements and their locations in each <span class="html-italic">DcCHS</span> gene are indicated by different colored blocks. (<b>C</b>) Numbers of different promoter elements in the <span class="html-italic">DcCHS</span> genes indicated by different numbers and colors in the grid. The number of <span class="html-italic">cis</span>-elements varies from few (blue) to many (red).</p>
Full article ">Figure 6
<p>Expression of <span class="html-italic">DcCHS</span> genes in different organs of <span class="html-italic">D. catenatum</span> using RT-qPCR. Mean expression values were obtained using TBtools; green and red indicate low and high levels of expression, respectively. RO: root; ST: stem; LE: leaf; SE: sepal; GY: gynostemia; PE: petal; LI: lip; FS: flower stalk.</p>
Full article ">Figure 7
<p>Expression analysis of the <span class="html-italic">DcCHS</span> genes under drought stress using RT-qPCR. The data are expressed as mean ± standard deviation (<span class="html-italic">n</span> = 3). Vertical bars represent means of fold change in expression and standard deviations calculated from replicates. Values of 0, 3, 6, 9, 12, 24, and 48 indicate hours after treatment. Asterisks (* or **) indicate a significant difference at <span class="html-italic">p</span> &lt; 0.05 or 0.01, respectively.</p>
Full article ">Figure 8
<p>Expression analysis of the <span class="html-italic">DcCHS</span> genes under heat stress using RT-qPCR. The data are expressed as mean ± standard deviation (<span class="html-italic">n</span> = 3). Vertical bars represent means of fold change in expression and standard deviations calculated from replicates. Values of 0, 3, 6, 9, 12, 24, and 48 indicate hours after treatment. Asterisks (* or **) indicate a significant difference at <span class="html-italic">p</span> &lt; 0.05 or 0.01, respectively.</p>
Full article ">Figure 9
<p>Expression analysis of the <span class="html-italic">DcCHS</span> genes under cold stress using RT-qPCR. The data are expressed as mean ± standard deviation (<span class="html-italic">n</span> = 3). Vertical bars represent means of fold change in expression and standard deviations calculated from replicates. Values of 0, 3, 6, 9, 12, 24, and 48 indicate hours after treatment. Asterisks (* or **) indicate a significant difference at <span class="html-italic">p</span> &lt; 0.05 or 0.01, respectively.</p>
Full article ">Figure 10
<p>Expression analysis of the <span class="html-italic">DcCHS</span> genes under salt stress using RT-qPCR. The data are expressed as mean ± standard deviation (<span class="html-italic">n</span> = 3). Vertical bars represent means of fold change in expression and standard deviations calculated from replicates. Values of 0, 3, 6, 9, 12, 24, and 48 indicate hours after treatment. Asterisks (* or **) indicate a significant difference at <span class="html-italic">p</span> &lt; 0.05 or 0.01, respectively.</p>
Full article ">
19 pages, 4435 KiB  
Article
Identification of the bZIP Gene Family and Investigation of Their Response to Drought Stress in Dendrobium catenatum
by Peng Wang, Yuxin Li, Tingting Zhang, Yuqian Kang, Wei Li, Jian Wang, Wengang Yu and Yang Zhou
Agronomy 2023, 13(1), 236; https://doi.org/10.3390/agronomy13010236 - 12 Jan 2023
Cited by 2 | Viewed by 1486
Abstract
Dendrobium catenatum is a valuable Chinese herbal medicine that naturally grows on cliffs and tree trunks and is often threatened by adverse environmental conditions. The bZIP transcription factor is known to play a critical role in the response of plant to stress. However, [...] Read more.
Dendrobium catenatum is a valuable Chinese herbal medicine that naturally grows on cliffs and tree trunks and is often threatened by adverse environmental conditions. The bZIP transcription factor is known to play a critical role in the response of plant to stress. However, the functions of the bZIP gene family in D. catenatum are poorly understood. In this study, 62 bZIP genes were identified from D. catenatum, which encoded proteins with an amino acid number of 130~692, a molecular weight of 15.24 to 74.94 kDa, and an isoelectric point of 5.13 to 11.58. The bZIP family can be divided into 10 subgroups by evolutionary tree analysis, and the conserved motifs of each protein subgroup were similar. The exon number of bZIP genes ranged from 1 to 12 as shown by gene structure analysis. DcbZIP promoter prediction analysis identified 21 cis-acting elements. The expression of DcbZIP genes under drought treatment was analyzed using the public RNA-seq data, and 33 upregulated genes were further screened. A co-expression network analysis revealed that 17 core genes were closely correlated with other genes and their expression was measured using RT-qPCR. The results showed that DcbZIP6, DcbZIP34, DcbZIP42 and DcbZIP47 are the main contributors to drought tolerance in D. catenatum. In summary, we identified candidate bZIP genes in D. catenatum with a apotential contribution to drought stress response, and this study lays the foundation for exploring the functions of bZIP and provides a theoretical basis for improving the drought tolerance of D. catenatum. Full article
(This article belongs to the Section Horticultural and Floricultural Crops)
Show Figures

Figure 1

Figure 1
<p>Phylogenetic analyses of bZIP proteins from <span class="html-italic">D. catenatum</span>, <span class="html-italic">Arabidopsis</span> and rice. A phylogenetic tree of bZIP proteins was constructed using MEGA 6.0 software. The 10 subgroups are indicated by different colors and different colored letters. The blue circles represent <span class="html-italic">A. thaliana</span> bZIPs (AtbZIPs), the red stars represent <span class="html-italic">D. catenatum</span> bZIPs (DcbZIPs) and the black triangles represent <span class="html-italic">O. sativa</span> bZIPs (OsbZIPs).</p>
Full article ">Figure 2
<p>Motifs, structures and phylogenetic relationships of DcbZIP family members. (<b>A</b>) A phylogenetic tree of 62 DcbZIP proteins constructed using the maximum likelihood method. The different subgroups are indicated by different background colors and letters. (<b>B</b>) Conserved motifs of DcbZIP proteins. Different motifs are represented by different colored boxes. (<b>C</b>) Exon/intron structures of <span class="html-italic">DcbZIP</span> genes with UTR(s), exon(s) and intron(s) indicated with green boxes, yellow boxes and black lines, respectively. The phylogenetic tree, conserved motifs and gene structures were predicted with TBtools.</p>
Full article ">Figure 3
<p>Analysis of <span class="html-italic">cis</span>-elements in the <span class="html-italic">DcbZIP</span> genes promoter regions. (<b>A</b>) A phylogenetic tree of 62 DcbZIP proteins was constructed with the ML method. The different subgroups are indicated by different background colors and letters. (<b>B</b>) The different colored blocks represent the different types of <span class="html-italic">cis</span>-elements and their locations in each <span class="html-italic">DcbZIP</span> gene. (<b>C</b>) The different colors and numbers in the grid indicate the numbers of different promoter elements in the <span class="html-italic">DcbZIP</span> genes. The types, numbers and locations of potential elements in the promoter regions 2-kb upstream of the <span class="html-italic">DcbZIP</span> genes were analyzed by PlantCARE.</p>
Full article ">Figure 4
<p>Hierarchical clustering of <span class="html-italic">DcbZIP</span> genes expression profiles across different <span class="html-italic">D. catenatum</span> tissues. Data were normalized relative to each gene’s mean expression value across all tissues and log2−transformed. TPM (transcripts per million) values were used to create heat maps showing the expressions of <span class="html-italic">DcbZIP</span> genes in different tissues. The expression level ranges from low expression (green) to high expression (red).</p>
Full article ">Figure 5
<p>The expression patterns of <span class="html-italic">DcbZIP</span> genes evaluated by RNA-seq. The left heatmap shows the expressions of <span class="html-italic">DcbZIP</span> genes under different volumetric water contents of base material. Data were normalized relative to each gene’s mean expression value across all treatments and log2-transformed. TPM (transcripts per million) values were used to create heat map showing the expression of <span class="html-italic">DcbZIP</span> genes. The expression levels range from low expression (green) to high expression (red). The right heatmap shows the TPM ratios, with high ratios in light orange and low ratios in cyan. T1 was equal to the TPM of 10–15% divided by the TPM of 30–35%, and T2 was equal to the TPM of 0% divided by the TPM of 30–35%. Ranges of 30–35%, 10–15% and 0% represent the volumetric water content of the base material declining to ~30–35%, ~10–15% and ~0%, respectively.</p>
Full article ">Figure 6
<p>Comprehensive analysis of drought-related <span class="html-italic">DcbZIP</span> genes. (<b>A</b>) Correlation analysis of 33 upregulated <span class="html-italic">DcbZIP</span> genes. Red indicates a positive correlation and blue indicates a negative correlation. The circle size indicates the absolute value of the correlation coefficient. (<b>B</b>) Co-expression network analysis. The red lines represent a positive correlation. The size represents the degree calculated by the Cytoscape. The four polygons situated in the center represent core genes with high connectivity.</p>
Full article ">Figure 7
<p>Interaction networks of the DcbZIPs in <span class="html-italic">D. catenatum</span> and their orthologs in <span class="html-italic">Arabidopsis</span>. The network was predicted by STRING. The DcbZIP proteins were shown in the red font below with the <span class="html-italic">Arabidopsis</span> orthologs in the red dashed box. Empty nodes represent proteins with unknown 3D structures. Filled nodes represent known or predicted 3D structures. Light blue lines and rose lines represent known interactions. Green, orange and blue lines represent predicted interactions. Cyan lines represent text mining, black lines represent co-expression and light purple lines represent protein homology.</p>
Full article ">Figure 8
<p>Real-time quantitative PCR analyses of the 17 core <span class="html-italic">DcbZIP</span> genes under drought stress in roots, stems and leaves. The mean expression value was calculated from three replicates. Vertical bars indicate the standard deviation. Values of 0, 3, 6, 9, 12, 24 and 48 indicate hours after treatment. Mean values and standard deviations were calculated according to the data. The unstressed level (0 h) was used as a control. Asterisks (* or **) indicate a significant difference at <span class="html-italic">p</span> &lt; 0.05 or <span class="html-italic">p</span> &lt; 0.01, respectively.</p>
Full article ">
11 pages, 2709 KiB  
Article
Dendrobium Multi-Omics Reveal Lipid Remodeling in Response to Freezing
by Xinqiao Zhan, Yichun Qian and Bizeng Mao
Metabolites 2022, 12(12), 1216; https://doi.org/10.3390/metabo12121216 - 3 Dec 2022
Cited by 3 | Viewed by 1506
Abstract
Freezing damage is a common phenomenon responsible for reduced yields of economic crops. Regulation of lipid metabolism plays an important role in plant growth and adaptation during freezing. We previously carried out transcriptome and untargeted metabolome analyses to determine the regulation of flavonol [...] Read more.
Freezing damage is a common phenomenon responsible for reduced yields of economic crops. Regulation of lipid metabolism plays an important role in plant growth and adaptation during freezing. We previously carried out transcriptome and untargeted metabolome analyses to determine the regulation of flavonol and anthocyanin biosynthesis during freezing treatment (FT) and post-freezing recovery (FR) in Dendrobium catenatum. However, changes in lipid levels are hard to confirm by untargeted metabolomics analysis alone. Regulation of lipid metabolism in response to freezing is largely unknown in Dendrobium. In this study, a multi-omics strategy was used to offer a better means of studying metabolic flow during FT and FR. To this end, 6976 proteins were identified by the 4D_label-free proteome, including 5343 quantified proteins. For each of the two conditions, we enriched differentially accumulated proteins (DAPs) into 15 gene ontology (GO) terms, including primary metabolism, lipid metabolism, and photosynthesis processes. We also identified 7 lipid categories and 3672 lipid species using lipidome assays. We found significant remodeling occurring in the phospholipid category during FT and FR. We also found that most sphingolipids were significantly upregulated. An integrated multi-omics analysis revealed significant changes in the expression levels of 141 mRNAs and encoding proteins under both FT and FR conditions. During FT, phospholipase A (PLA) and phospholipase D (PLD) were associated with phospholipid editing and galactolipid remodeling. These results provide valuable new insights into how the freezing tolerance of D. catenatum might be improved by genetic engineering. Full article
(This article belongs to the Section Plant Metabolism)
Show Figures

Figure 1

Figure 1
<p>Overview of the proteomes during freezing stress. (<b>a</b>) A heatmap of the abundance of proteins in CK, FT, and FR. The bar indicates the significant values. (<b>b</b>) PCA analysis of the proteomes in CK, FT, and FR. (<b>c</b>) Differentially accumulated proteins (DAPs) in three comparisons.</p>
Full article ">Figure 2
<p>Gene ontology (GO) enrichment analysis of DAPs and DEGs in two comparisons. (<b>a</b>) The Go enrichment analysis of DAPs in FT vs. CK. (<b>b</b>) The Go enrichment analysis of DAPs in FR vs. CK. The <span class="html-italic">x</span>-axis represents the z-score of GO enrichment, while the <span class="html-italic">y</span>-axis represents −log<sub>10</sub> adjusted <span class="html-italic">p</span>-value. The dot sizes represent the number of DAPs. GO terms are labeled with a cut-off (adjusted <span class="html-italic">p</span>-value &lt; 0.05). (<b>c</b>) The cellular component enrichment of DEGs from three comparisons. The bar indicates the <span class="html-italic">p</span>-value of KEGG enrichment. (<b>d</b>) Heat maps show the abundance of representative transcripts from cellular component enrichment. The enrichment pathway of DEGs is shown on the right. The bar indicates the significant values of gene expression levels. Gene abbreviations are listed in <a href="#app1-metabolites-12-01216" class="html-app">Table S2</a>.</p>
Full article ">Figure 3
<p>Overview of the lipidomes during freezing stress. (<b>a</b>) A heatmap of the abundance of lipid species in CK, FT, and FR. (<b>b</b>) PCA analysis of lipidomes in CK, FT, and FR. (<b>c</b>) The content of total lipid in CK, FT, and FR. Values are means ± S.D. (<span class="html-italic">n</span> = 3); Student’s <span class="html-italic">t</span>-test. (<b>d</b>) Differentially accumulated lipids (DALs) of two comparisons. Lipid abbreviations are listed in <a href="#app1-metabolites-12-01216" class="html-app">Table S3</a>.</p>
Full article ">Figure 4
<p>Integrated analysis of multi-omics during freezing stress. (<b>a</b>) Scatter plot showing the relationship between changes in protein and mRNA abundances in FT vs. CK. (<b>b</b>) Scatter plot showing the relationship between changes in protein and mRNA abundances in FT vs. CK. The colored points indicate significant upregulation or downregulation of protein and mRNA levels. (<b>c</b>) KEGG enrichment analysis of mRNAs and encoding proteins based on (<b>a</b>,<b>b</b>). (<b>d</b>) Enrichment of mRNAs and encoding proteins in FT vs. CK. (<b>e</b>) Enrichment of mRNAs and encoding proteins in FR vs. CK. (<b>f</b>) Simplified lipid biosynthesis flow based on multi-omics. FCs in lipid, mRNA (underlined on the right), and protein levels (left) are listed and highlighted in different colors according to FC values. The ion characteristic fragments and chromatogram of lipids are shown in <a href="#app1-metabolites-12-01216" class="html-app">Figure S5</a>. Abbreviations: DGDS, digalactoslydiacylglycerol synthase; MGDS, monogalactosyldiacylglycerol synthase; PA, phosphatidic acid; PP, phosphatidate phosphatase; PLA, phospholipase A; PLD, phospholipase D.</p>
Full article ">
15 pages, 4106 KiB  
Article
Population Diversity Analysis Provide Insights into Provenance Identification of Dendrobium catenatum
by Xin-Yi Wu, Ting-Zhang Li, Fang Zheng, Jian-Bing Chen, Yue-Hong Yan and Jiu-Xiang Huang
Genes 2022, 13(11), 2093; https://doi.org/10.3390/genes13112093 - 10 Nov 2022
Viewed by 1631
Abstract
Dendrobium catenatum (Dendrobium officinale) is a valuable genuine herb. The source of this species is difficult to be identified by traditional methods including morphology, spectroscopy, and chromatography. We used the restriction site-associated DNA sequencing (RAD-seq) approach to perform the high-throughput sequencing [...] Read more.
Dendrobium catenatum (Dendrobium officinale) is a valuable genuine herb. The source of this species is difficult to be identified by traditional methods including morphology, spectroscopy, and chromatography. We used the restriction site-associated DNA sequencing (RAD-seq) approach to perform the high-throughput sequencing of 24 D. catenatum provenances. In this study, 371.18 Gb clean data were obtained, and 655,057 high-quality SNPs were selected after their filtration. We used phylogenetic tree, genetic structure, and principal component analyses to examine the genetic diversities and genetic relationships of the 109 accessions. We found that D. catenatum could be divided into two groups, and each group was closely related to the distribution of the sampling sites. At the population level, the average nucleotide diversity (π) of the D. catenatum population mutation parameters was 0.1584 and the expected heterozygosity (HE) was 0.1575. The GXLPTP07 accessions showed the highest genetic diversity in terms of the private allele number, observed heterozygosity, and nucleotide diversity. The Mantel test showed a significant positive correlation between the genetic and geographic distances among the overall distribution. A genetic information database of D. catenatum was established, which confirmed that RAD-seq technology has the potential to be applied in the identification of medicinal Dendrobium of different origins. Full article
(This article belongs to the Special Issue Phylogenetics, Genetics, and Breeding of Medicinal Plants)
Show Figures

Figure 1

Figure 1
<p>Sampling locations of <span class="html-italic">D. catenatum</span> species used in this study.</p>
Full article ">Figure 2
<p>ML phylogenetic tree of the <span class="html-italic">D. catenatum</span> accessions and model-based clustering with K from 2. Numbers near the nodes are bootstrap percentages.</p>
Full article ">Figure 3
<p>Genetic structure of cultivated 109 <span class="html-italic">D. catenatum</span> for K = 2–8 based on the Admixture software (K = 2 with cross validation error is 0.18580, while K = 3 with ross validation error is 0.19141).</p>
Full article ">Figure 4
<p>PCA of 109 samples in <span class="html-italic">D. catenatum</span>. Red shapes indicate Group I, while blue shapes indicate Group II.</p>
Full article ">Figure 5
<p>Linkage disequilibrium decay patterns of different <span class="html-italic">D. catenatum</span> group.</p>
Full article ">Figure 6
<p>Genetic distance (<span class="html-italic">F</span><sub>ST</sub>) for <span class="html-italic">D. catenatum</span>. The samples from Group I were in red font, and the samples from Group II were in blue font.</p>
Full article ">Figure 7
<p>Relationships between the pairwise genetic distance (<span class="html-italic">F</span><sub>ST</sub>) and geographic distance with blue bars and shading represent the 95% highest posterior density.</p>
Full article ">
21 pages, 5093 KiB  
Article
Genome-Wide Identification of AP2/ERF Transcription Factor Family and Functional Analysis of DcAP2/ERF#96 Associated with Abiotic Stress in Dendrobium catenatum
by Yuliang Han, Maohong Cai, Siqi Zhang, Jiawen Chai, Mingzhe Sun, Yingwei Wang, Qinyu Xie, Youheng Chen, Huizhong Wang and Tao Chen
Int. J. Mol. Sci. 2022, 23(21), 13603; https://doi.org/10.3390/ijms232113603 - 6 Nov 2022
Cited by 16 | Viewed by 2143
Abstract
APETALA2/Ethylene Responsive Factor (AP2/ERF) family plays important roles in reproductive development, stress responses and hormone responses in plants. However, AP2/ERF family has not been systematically studied in Dendrobium catenatum. In this study, 120 AP2/ERF family members were identified for the [...] Read more.
APETALA2/Ethylene Responsive Factor (AP2/ERF) family plays important roles in reproductive development, stress responses and hormone responses in plants. However, AP2/ERF family has not been systematically studied in Dendrobium catenatum. In this study, 120 AP2/ERF family members were identified for the first time in D. catenatum, which were divided into four groups (AP2, RAV, ERF and DREB subfamily) according to phylogenetic analysis. Gene structures and conserved motif analysis showed that each DcAP2/ERF family gene contained at least one AP2 domain, and the distribution of motifs varied among subfamilies. Cis-element analysis indicated that DcAP2/ERF genes contained abundant cis-elements related to hormone signaling and stress response. To further identify potential genes involved in drought stress, 12 genes were selected to detect their expression under drought treatment through qRT-PCR analysis and DcAP2/ERF#96, a nuclear localized ethylene-responsive transcription factor, showed a strong response to PEG treatment. Overexpression of DcAP2/ERF#96 in Arabidopsis showed sensitivity to ABA. Molecular, biochemical and genetic assays indicated that DcAP2ERF#96 interacts with DREB2A and directly inhibits the expression of P5CS1 in response to the ABA signal. Taken together, our study provided a molecular basis for the intensive study of DcAP2/ERF genes and revealed the biological function of DcAP2ERF#96 involved in the ABA signal. Full article
Show Figures

Figure 1

Figure 1
<p>Phylogenetic analysis of AP2/ERF in <span class="html-italic">Arabidopsis thaliana</span>, <span class="html-italic">D. catenatum</span> and <span class="html-italic">Phalaenopsis equestris</span>. (<b>A</b>) Phylogenetic trees of three species were constructed by the neighbor-joining method with 1000 bootstrap replications. Each gene cluster was labeled with distinguishable colors. (<b>B</b>) Gene quantity and structure of different subfamilies of <span class="html-italic">AP2/ERF</span> family in three species.</p>
Full article ">Figure 2
<p>Phylogenetic relationships, conserved motifs, protein domain and gene structure in <span class="html-italic">DcAP2/ERF</span> genes. (<b>A</b>) Conserved motif analysis of <span class="html-italic">DcAP2/ERFs</span>. The four subfamilies of genes were colored on the left side, the yellow color represented the ERF subfamily, the green color represented the DREB subfamily, the pink color represented the RAV subfamily, and the blue color represented the AP2 subfamily. (<b>B</b>) Conserved domain analysis of <span class="html-italic">DcAP2/ERFs</span>. (<b>C</b>) Gene structure analysis of <span class="html-italic">DcAP2/ERFs</span>. The black lines in each gene represent the introns.</p>
Full article ">Figure 3
<p>Classification and analysis of cis-acting elements in the promoter regions of <span class="html-italic">DcAP2/ERF</span> genes. The 2 kb region upstream of the genes was analyzed using the PlantCARE website. Graphs in different colors represent different classes of cis-acting elements. The numbers indicate the amount of each cis-acting element.</p>
Full article ">Figure 4
<p>Expression of <span class="html-italic">DcAP2/ERF</span> genes in response to cold stress. (<b>A</b>) Heat map showing the expression pattern <span class="html-italic">DcAP2/ERF</span> genes in leaves under cold stress for 20 h. The red, green, yellow and blue arcs in the outer ring of the heat map represent the range of <span class="html-italic">ERF</span>, <span class="html-italic">DREB</span>, <span class="html-italic">RAV</span> and <span class="html-italic">AP2</span> subfamily, respectively. (<b>B</b>) The volcano plot showing the upregulation and downregulation of genes under low-temperature treatment.</p>
Full article ">Figure 5
<p>Expression of <span class="html-italic">DcAP2/ERF</span> genes in response to MeJA. (<b>A</b>) Heat map showing expression pattern <span class="html-italic">DcAP2/ERF</span> genes under 1 mM MeJA. The red, green, yellow and blue arcs in the outer ring of the heat map represent the range of <span class="html-italic">ERF</span>, <span class="html-italic">DREB</span>, <span class="html-italic">RAV</span> and <span class="html-italic">AP2</span> subfamily, respectively. (<b>B</b>) The volcano plot showing the upregulation and downregulation of genes under 1 mM MeJA treatment.</p>
Full article ">Figure 6
<p>Expression patterns of 12 selected genes in leaves, stems and roots under 20% PEG6000 treatment. The orange color represented leaves, the green color represented stems, and the yellow color represented roots. <span class="html-italic">DcACTIN</span> was used as an internal control. Values are presented as means ± SD (<span class="html-italic">n</span> = 3). (* <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> &lt; 0.01, Student’s <span class="html-italic">t</span>-test).</p>
Full article ">Figure 7
<p>The growth status of <span class="html-italic">D. catenatum</span> stems is strongly affected by drought stress. (<b>A</b>) The growth status of <span class="html-italic">D. catenatum</span> under normal and extreme drought conditions. (<b>B</b>) Statistics on biomass loss under different degrees of drought stress.</p>
Full article ">Figure 8
<p><span class="html-italic">DcAP2ERF#96</span> exhibits significant transcription factor characteristics. (<b>A</b>,<b>B</b>) Subcellular localization of DcAP2ERF#96 protein in <span class="html-italic">D. catenatum</span> protoplasts (<b>A</b>) and in the leaf epidermal cells of <span class="html-italic">N. benthamiana</span> (<b>B</b>). D53-mCherry was used as a nuclear marker. Bars = 20 μm. The fluorescence signals were detected by confocal microscopy. (<b>C</b>) Yeast one-hybrid verification of the binding of DcAP2ERF#96 protein to three tandem repeats of CRT box, DRE box and GCC box. (<b>D</b>) Yeast two-hybrid validation of DcAP2ERF#96 protein self-activation regions. The transformed yeast cells were plated on DDO (SD/-Trp/-Leu) and QDO (SD/-Trp/-Leu/-His/-Ade).</p>
Full article ">Figure 9
<p><span class="html-italic">DcAP2ERF#96</span> overexpression lines showed sensitivity to ABA. (<b>A</b>) 10-d-old seedlings were transferred to vertical plates containing <sup>1</sup>/<sub>2</sub> MS + 10 μM ABA and its control (<sup>1</sup>/<sub>2</sub> MS). Phenotypes were observed and recorded one week later. Bars = 1 cm. (<b>B</b>,<b>C</b>) Leaf area (<b>B</b>) and root length (<b>C</b>) statistics of <span class="html-italic">Col-0</span> and <span class="html-italic">DcAP2ERF#96</span> overexpression lines under 10 μM ABA treatment (<span class="html-italic">n</span> = 9). (<b>D</b>) Relative expression levels of ABA signaling pathway-related genes in <span class="html-italic">Col-0</span> and overexpression lines after 10 μM ABA treatment for one week. <span class="html-italic">AtACTIN2</span> was used as an internal control. Values are presented as means ± SD (<span class="html-italic">n</span> = 3). * <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> &lt; 0.01, Student’s <span class="html-italic">t</span>-test.</p>
Full article ">Figure 10
<p><span class="html-italic">DcAP2ERF#96</span> interacts with <span class="html-italic">DcDREB2A</span> and inhibits <span class="html-italic">P5CS1</span> and <span class="html-italic">RD29A</span> in <span class="html-italic">Arabidopsis thaliana</span>. (<b>A</b>) Distribution of DRE, CRT and GCC motifs in ABA-related gene <span class="html-italic">P5CS1</span>, <span class="html-italic">RD29A</span> and <span class="html-italic">RAB18</span> promoters. (<b>B</b>) Yeast one hybrid verification of the binding of <span class="html-italic">DcAP2ERF#96</span> to <span class="html-italic">P5CS1</span>, <span class="html-italic">RD29A</span> and <span class="html-italic">RAB18</span> promoters. The pB42AD-DcAP2ERF#96 and pLacZi2μ with different promoters were co-transformed into the EGY48 yeast strain. The transformed yeast cells were plated on DDO (SD/-Ura/-Trp) with X-Gal. (<b>C</b>) Yeast two hybrids verified the interaction of DcAP2ERF#96 with DREB2A-1 and DREB2A-2. The transformed yeast cells were plated on DDO (SD/-Trp/-Leu) and QDO (SD/-Trp/-Leu/-His/-Ade). (<b>D</b>) Regulatory model for the regulation of the ABA signal by <span class="html-italic">DcAP2ERF#96</span> in <span class="html-italic">Arabidopsis thaliana</span> and <span class="html-italic">D. catenatum</span>. Overexpression of <span class="html-italic">DcAP2ERF#96</span> inhibits the expression of ABA downstream genes such as <span class="html-italic">RD29A</span> and <span class="html-italic">P5CS1</span> by binding to their promoters, thereby affecting plant sensitivity to ABA signaling. DcAP2ERF#96 can interact with DREB2A protein in <span class="html-italic">D. catenatum</span>, whose homologous gene in <span class="html-italic">Arabidopsis</span> has been reported to positively regulate the expression of <span class="html-italic">RD29A</span>.</p>
Full article ">
23 pages, 5106 KiB  
Article
Comprehensive Genomic Characterization of the NAC Transcription Factors and Their Response to Drought Stress in Dendrobium catenatum
by Yuxin Li, Tingting Zhang, Wenting Xing, Jian Wang, Wengang Yu and Yang Zhou
Agronomy 2022, 12(11), 2753; https://doi.org/10.3390/agronomy12112753 - 5 Nov 2022
Cited by 7 | Viewed by 1592
Abstract
As a large transcription factor family, NAC family proteins (NAM, ATAF1,2, and CUC2) play critical roles in plant growth, development, and response to stresses. Herein, the NAC gene family of Dendrobium catenatum was identified and analyzed by bioinformatics methods. Their expression patterns in [...] Read more.
As a large transcription factor family, NAC family proteins (NAM, ATAF1,2, and CUC2) play critical roles in plant growth, development, and response to stresses. Herein, the NAC gene family of Dendrobium catenatum was identified and analyzed by bioinformatics methods. Their expression patterns in different tissues and under drought stress were analyzed using RNA-seq data and the quantitative real-time reverse transcription-polymerase chain reaction (RT-qPCR) method. A total of 90 NAC genes were identified, encoding amino acids with numbers ranging from 88 to 1065, with protein molecular weight ranging from 10.34 to 119.24 kD, and isoelectric point ranging from 4.5 to 9.99. Phylogenetic analysis showed that the DcNAC proteins could be divided into 17 subgroups, and each subgroup had conserved motif composition and gene structure. Twenty types of cis-elements were identified in the DcNAC promoters. RNA-seq analysis showed that the expression of DcNAC genes had tissue specificity and displayed different expression patterns under drought stress. Co-expression network analysis of the DcNAC genes revealed nine hub genes, and their expression levels were then validated by RT-qPCR. The results showed that DcNAC6, DcNAC18, DcNAC29, DcNAC44, and DcNAC51 (mainly in roots) as well as DcNAC16 and DcNAC64 (mainly in leaves) were considered as the candidate genes for drought tolerance in D. catenatum. Taken together, this study identified candidate NAC genes with potential functions in response to drought stress, which is valuable for development of drought resistance in D. catenatum. Full article
(This article belongs to the Section Horticultural and Floricultural Crops)
Show Figures

Figure 1

Figure 1
<p>Phylogenetic analyses of NAC proteins from <span class="html-italic">D. catenatum</span>, <span class="html-italic">Arabidopsis</span>, and rice. A phylogenetic tree of NAC proteins was constructed using MEGA-X software (version 10.0.5). The NAC domains were classified into two large groups: Groups I and II. Group I was divided into 13 subgroups (TERN, ONAC022, SENU5, NAP, AtNAC3, ATAF, OsNAC3, NAC2, ANAC011, TIP, OsNAC7, NAC1, and NAM). Group II was divided into 4 subgroups (ANAC001, ONAC003, ONAC001, and ANAC063). The 17 subgroups are indicated with different colors. The red pentacles represent <span class="html-italic">A. thaliana</span> NACs (AtNACs), the green boxes represent <span class="html-italic">D. catenatum</span> NACs (DcNACs), and the blue triangles represent <span class="html-italic">O. sativa</span> NACs (OsNACs).</p>
Full article ">Figure 2
<p>Motifs, structures, and phylogenetic relationships of DcNAC family members. (<b>a</b>) A phylogenetic tree of 90 DcNAC proteins was constructed using the maximum likelihood method. The different subgroups are indicated with different background colors and letters. (<b>b</b>) Conserved motifs of DcNAC proteins. Different motifs are represented by different colored boxes. (<b>c</b>) Exon/intron structures of DcNAC genes with UTR(s), exon(s), and intron(s) are indicated with yellow and green boxes, and black lines, respectively. The phylogenetic tree, conserved motifs, and gene structures were predicted with TBtools.</p>
Full article ">Figure 3
<p>NAC domain and TAR in <span class="html-italic">D. catenatum</span> NAC family protein. (<b>a</b>) The predicted domains in the N-terminus and C-terminus. The NAC domain consists of subdomains A to E. The TAR (transcriptional activation region) is the C-terminal activation domain. (<b>b</b>) Alignment of the NAC domain from the MEME results for the DcNACs. The motifs 3, 4, 1, 2, and 5 form the putative DcNAC NAC domain.</p>
Full article ">Figure 4
<p>The number of <span class="html-italic">DcNAC</span> promoters containing various <span class="html-italic">cis</span>-acting elements. Different colors represent different <span class="html-italic">cis</span>-element types. The numbers above the columns represent the number of promoters.</p>
Full article ">Figure 5
<p>Hierarchical clustering of expression profiles of <span class="html-italic">DcNAC</span> gene expression profiles across different <span class="html-italic">D. catenatum</span> tissues. Data were normalized relative to each gene’s mean expression value across all tissues and log2-transformed. TPM (transcripts per million) values were used to create heat maps showing the expression of <span class="html-italic">DcNAC</span> genes in different tissues, ranging from low expression (green) to high expression (red).</p>
Full article ">Figure 6
<p>The expression patterns of <span class="html-italic">DcNAC</span> genes evaluated by RNA-seq. The left heatmap shows the expression of <span class="html-italic">DcNAC</span> genes under different volumetric water contents of base material. Data were normalized relative to each gene’s mean expression value across all treatments and log2-transformed. TPM (transcripts per million) values were used to create a heat map showing the expression of <span class="html-italic">DcNAC</span> genes. The expression level ranged from low expression (green) to high expression (red). The right heatmap showed the TPM ratios, with high ratios in red and low ratios in cyan. T1 was equal to the TPM of 10–15% divided by the TPM of 30–35%, and T2 was equal to the TPM of 0% divided by the TPM of 30–35%. The amounts of 30–35%/10–15%/0% represent the fact that the volumetric water content of the base material declined to ≈30–35%, ≈10–15%, and ≈0%, respectively.</p>
Full article ">Figure 7
<p>Comprehensive analysis of drought-related <span class="html-italic">DcNAC</span> genes. (<b>a</b>) Correlation analysis of upregulated <span class="html-italic">DcNAC</span> genes. Red indicates a positive correlation, and blue indicates a negative correlation. The circle size indicates the absolute value of the correlation coefficient. (<b>b</b>) Co-expression network analysis. The red line represents a positive correlation. The size represents the degree calculated by the Cytoscape.</p>
Full article ">Figure 8
<p>Interaction networks of the selected DcNACs in <span class="html-italic">D. catenatum</span> and related genes in <span class="html-italic">Arabidopsis</span>. The colored, white, empty, and filled nodes represent the first shell of interactors, the second shell of interactors, proteins with unknown 3D structure, and proteins with known or predicted 3D structure, respectively. The light blue and purple lines indicate known interactions. The green, red, and blue lines indicate predicated interactions. The cyan, black, and yellow lines represent protein homology, co-expression, and text mining, respectively.</p>
Full article ">Figure 9
<p>Real-time quantitative PCR analyses of the nine hub <span class="html-italic">DcNAC</span> genes under drought stress in roots and leaves. The mean expression value was calculated from three replicates. Vertical bars indicate the standard deviation. Values of 0, 3, 6, 9, 12, 24, and 48 indicate hours after treatment. Mean values and standard deviations were calculated according to the data. The unstressed level (0 h) was used as a control. Asterisks (* or **) indicate a significant difference at <span class="html-italic">p</span> &lt; 0.05 or <span class="html-italic">p</span> &lt; 0.01, respectively.</p>
Full article ">
20 pages, 2646 KiB  
Review
Glucomannan in Dendrobium catenatum: Bioactivities, Biosynthesis and Perspective
by Luyan Qi, Yan Shi, Cong Li, Jingjing Liu, Sun-Li Chong, Kean-Jin Lim, Jinping Si, Zhigang Han and Donghong Chen
Genes 2022, 13(11), 1957; https://doi.org/10.3390/genes13111957 - 27 Oct 2022
Cited by 11 | Viewed by 2659
Abstract
Dendrobium catenatum is a classical and precious dual-use plant for both medicine and food in China. It was first recorded in Shen Nong’s Herbal Classic, and has the traditional functions of nourishing yin, antipyresis, tonifying the stomach, and promoting fluid production. The [...] Read more.
Dendrobium catenatum is a classical and precious dual-use plant for both medicine and food in China. It was first recorded in Shen Nong’s Herbal Classic, and has the traditional functions of nourishing yin, antipyresis, tonifying the stomach, and promoting fluid production. The stem is its medicinal part and is rich in active polysaccharide glucomannan. As an excellent dietary fiber, glucomannan has been experimentally confirmed to be involved in anti-cancer, enhancing immunity, lowering blood sugar and blood lipids, etc. Here, the status quo of the D. catenatum industry, the structure, bioactivities, biosynthesis pathway and key genes of glucomannan are systematically described to provide a crucial foundation and theoretical basis for understanding the value of D. catenatum and the potential application of glucomannan in crop biofortification. Full article
(This article belongs to the Special Issue Phylogenetics, Genetics, and Breeding of Medicinal Plants)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Epiphytic cultivation model of <span class="html-italic">D. catenatum</span>. (<b>A</b>). Facility-aided cultivation; (<b>B</b>) Rock-dependent eco-cultivation; (<b>C</b>) Trunk-dependent eco-cultivation; (<b>D</b>) Stereo cultivation.</p>
Full article ">Figure 2
<p>Compositions of different glucomannans in <span class="html-italic">D. catenatum</span> ([<a href="#B27-genes-13-01957" class="html-bibr">27</a>,<a href="#B29-genes-13-01957" class="html-bibr">29</a>,<a href="#B30-genes-13-01957" class="html-bibr">30</a>]).</p>
Full article ">Figure 3
<p>The effects of glucomannan from <span class="html-italic">D. catenatum</span> stem. DOP, <span class="html-italic">Dendrobium officinale</span> polysaccharide; TLR4, toll-like receptor4; IL-4, Interleukin-4; IL-6, Interleukin-6; IL-10, Interleukin-10; IL-1β, Interleukin-1β; TNF, tumor necrosis factor; PI3K/AKT, phosphatidylinositol 3 kinase/protein kinase B; CXCR4, chemokine receptor4; CCR7, CC chemokine receptor 7 CC chemokine receptor 7; AMPK, adenosine monophosphate-activated protein kinase; ULK-1, UNC-51 like autophagy activating kinase 1; LC3, light chain 3; SCFAs, short-chain fatty acids; GPRs, g-protein-coupled receptors; MAPK, mitogen-activated protein kinase; PP, Peyer’s patche; MLN, mesenteric lymph node; MNNG, 1-methyl-2-nitro-1-nitroguanidine; GS, glycogen synthase; GP, glycogen phosphorylase; GSK-3β, glycogen synthase kinase 3β; PK, pyruvate kinase; HK, hexokinase; PEPCK, phosphoenolpyruvate carboxykinase; G6Pase, glucose-6-phosphatase; GLP-1, glucagon-like peptide-1.</p>
Full article ">Figure 4
<p>Putative glucomannan biosynthetic pathway in <span class="html-italic">D. catenatum</span> (modified from [<a href="#B41-genes-13-01957" class="html-bibr">41</a>,<a href="#B44-genes-13-01957" class="html-bibr">44</a>,<a href="#B45-genes-13-01957" class="html-bibr">45</a>,<a href="#B88-genes-13-01957" class="html-bibr">88</a>] Glc, glucose; Fru, fructose; Fru-6-P, fructose-6-phosphate; Glc-1-P, glucosophosphate-1-P; Glc-6-P, glucosophosphate-6-P; Man, mannose; Man-6-P, mannose-6-phosphate; Man-1-P, mannose-1-phosphate; AsA, ascorbic acid; SUT, sucrose transporter; STP, sugar transporter proteins; SuS, sucrose synthase; INV, invertase; FRK, fructokinase; HXK, hexokinase; PGI, phosphoglucose isomerase; UGP, UDP-Glc pyrophosphorylase; AGP, ADP-Glc pyrophosphorylase; GGP, GDP-Glc pyrophosphorylase; PGM, phosphoglucomutase; GPT, glucose-6-phosphate transporter; BT1, brittle-1 protein, an ADP-Glc transporter; PMI, phosphate mannose isomerase; PMM, phosphomannomutase; GMP, GDP-mannose pyrophosphorylase; GMT, GDP-mannose transporter; CSLA, cellulose-like synthase A; CSLD, cellulose-like synthase D.</p>
Full article ">Figure 5
<p>Expression profile of glucomannan pathway genes in tissues and organs of <span class="html-italic">D. catenatum</span>. The raw RNA-seq reads of different tissues and organs in <span class="html-italic">D. catenatum</span> are derived from the NCBI database (<a href="http://www.ncbi.nlm.nih.gov" target="_blank">www.ncbi.nlm.nih.gov</a>, accessed on 16 September 2022), including the leaf (SRR4431601), the root (SRR5722140), the green root tip (SRR4431599), the white part of the root (SRR4431598), the stem (SRR4431600), the flower bud (SRR4431603), the sepal (SRR4431597), the labellum (SRR4431602), the pollinia (SRR5722145), and the gynostemium (SRR4431596). A heatmap has been generated via TBtools software [<a href="#B89-genes-13-01957" class="html-bibr">89</a>]. The color scale represents log2 of FPKM expression values; green and red indicate a low and high level of gene expression, respectively.</p>
Full article ">Figure 6
<p>Biofortification strategy of glucomannan based on the in-depth dissection of molecular mechanisms and combined application of new-developed technologies.</p>
Full article ">
9 pages, 1693 KiB  
Communication
Identification of Two GDSL-Type Esterase/Lipase Genes Related to Tissue-Specific Lipolysis in Dendrobium catenatum by Multi-Omics Analysis
by Xinqiao Zhan, Yichun Qian and Bizeng Mao
Life 2022, 12(10), 1563; https://doi.org/10.3390/life12101563 - 9 Oct 2022
Cited by 2 | Viewed by 1517
Abstract
Dendrobium catenatum is an important herb and widely cultivated in China. GDSL-Type Esterase/Lipase proteins (GELPs) are widely distributed in plants and play crucial roles in stress responses, plant growth, and development. However, no identification or functional analysis of GELPs was reported in D. [...] Read more.
Dendrobium catenatum is an important herb and widely cultivated in China. GDSL-Type Esterase/Lipase proteins (GELPs) are widely distributed in plants and play crucial roles in stress responses, plant growth, and development. However, no identification or functional analysis of GELPs was reported in D. catenatum. This study identifies 52 GELPs in D. catenatum genome, which is classified into four groups by phylogenetic analysis. Four conservative blocks (Ser-Gly-Asn-His) are found in most GELP domains. Transcriptome analysis reveals the expression profiles of GELPs in different organs and flowering phases. Co-expression analysis of the transcriptome and lipidome identifies a GELP gene, Dca016600, that positively correlates with 23 lipids. The purified Dca016600 protein shows the optimum pH is active from 8.0 to 8.5, and the optimum temperature is active from 30 °C to 40 °C. The kinetic study provides Vmax (233.43 μmol·min−1·mg−1) and Km (1.49 mM) for substrate p-nitrophenyl palmitate (p-NPP). Integrated analysis of the transcriptome and proteome identifies a GELP gene, Dca005399, which is specially induced by freezing. Interestingly, Dca005399 shows high expression in symbiotic germination seeds and sepals. This study provides new insights into the function of D. catenatum GELPs in plant development and stress tolerance. Full article
(This article belongs to the Special Issue Plant Biotic and Abiotic Stresses)
Show Figures

Figure 1

Figure 1
<p>Comprehensive identification of GELP family. (<b>a</b>) Phylogenetic analysis of GELPs in <span class="html-italic">D. catenatum</span>. A total of 52 GELPs were used to construct the unrooted maximum-likelihood phylogenies. (<b>b</b>) Domain organization of GELPs. Conservative blocks (Block I, II, III, and V) of GELPs were shown from up and down. Conservative amino acid residues Ser-Gly-Asn-His in blocks are marked by black triangles. (<b>c</b>) Expression patterns of GELP genes in different tissues. (<b>d</b>) Expression patterns of GELP genes in three flowering phases. S1, the flower buds were green in the early developmental stage; S2, the flowers had purple pigmentation in the columns and the lips; S3, the sepals and petals had turned yellow and red. (<b>e</b>) Lipidome analysis of four tissues in <span class="html-italic">D. catenatum</span>. Raw data is shown in <a href="#app1-life-12-01563" class="html-app">Table S1</a>. Lipid abbreviations is listed in <a href="#app1-life-12-01563" class="html-app">Table S2</a>. Color scales represented the values of log<sub>2</sub> in gene expression levels or lipid content.</p>
Full article ">Figure 2
<p>Identification and characterization of Dca016600. (<b>a</b>) Correlation analysis of lipidome and GELP expression levels. The dotted and solid lines, respectively, represent positive and negative correlations. The thickness of the line is determined by a Pearson correlation coefficient &gt;0.9 or &lt;−0.9, respectively. The dot sizes and colors represent the correlated number of lipids and genes. (<b>b</b>) Phylogenetic analysis of Dca016600 was intercepted to <a href="#app1-life-12-01563" class="html-app">Supplemental Figure S2</a>. Effects of pH (<b>c</b>), temperature (<b>d</b>), and the effect of substrate concentrations (<b>e</b>) for Dca016600 activity. Data represent the mean ± SD of three independent experiments. Km and Vmax values are determined using Lineweaver-Burk plot in <a href="#app1-life-12-01563" class="html-app">Figure S4</a>.</p>
Full article ">Figure 3
<p>Identification and expression analysis of Dca005399. (<b>a</b>) The relationship between changes in protein and mRNA abundances in FT vs. CK. (<b>b</b>) The relationship between changes in protein and mRNA abundances in FR vs. CK. The colored points indicate significant upregulation or downregulation of protein and mRNA levels. CK, control condition; FT, freezing treatment; FR, post-freezing recovery. (<b>c</b>) Phylogenetic analysis of Dca005399 was intercepted to <a href="#app1-life-12-01563" class="html-app">Supplemental Figure S2</a>. (<b>d</b>) Tissues expression analysis of <span class="html-italic">Dca005399</span>. Data represent the mean ± SD of three independent experiments. MS, mature seed; SGS, symbiotic germination seed; AGS, asymbiotic germination seed.</p>
Full article ">
13 pages, 4464 KiB  
Article
Genome-Wide Analysis of the WOX Transcription Factor Genes in Dendrobium catenatum Lindl.
by Hefan Li, Cheng Li, Yuhua Wang, Xiangshi Qin, Lihua Meng and Xudong Sun
Genes 2022, 13(8), 1481; https://doi.org/10.3390/genes13081481 - 19 Aug 2022
Cited by 6 | Viewed by 1989
Abstract
The WUSCHEL-related homeobox (WOX) proteins are a class of transcription factors exclusive to plants. They can promote cell division or inhibit stem cell differentiation to regulate plant growth and development. However, the WOX transcription factor genes in the monocotyledon Dendrobium catenatum Lindl. remain [...] Read more.
The WUSCHEL-related homeobox (WOX) proteins are a class of transcription factors exclusive to plants. They can promote cell division or inhibit stem cell differentiation to regulate plant growth and development. However, the WOX transcription factor genes in the monocotyledon Dendrobium catenatum Lindl. remain relatively uncharacterized. Specifically, the effects of phytohormones on their expression levels are unclear. In this study, we identified and analyzed 10 candidate DcaWOX transcription factor genes in D. catenatum. The DcaWOX family was divided into the modern/WUS, intermediate, and ancient clades. The subcellular localization analysis detected DcaWOX-GFP fusion proteins in the tobacco epidermal leaf cell nucleus. In DcaWOX, members of the WUS clade with the WUS-box motif can significantly activate the expression of TPL in vivo, while members of the intermediate and ancient clades cannot. The expression of the DcaWOX genes varied among the examined tissues. Moreover, the DcaWOX expression patterns were differentially affected by the phytohormone treatments, with differences detected even between homologs of the same gene. Furthermore, the gene expression patterns were consistent with the predicted cis-acting elements in the promoters. The above results suggest that DcaWOX may have an important role in its growth and development and resistance to stress. The results of this comprehensive investigation of the DcaWOX gene family provide the basis for future studies on the roles of WOX genes in D. catenatum. Full article
(This article belongs to the Section Molecular Genetics and Genomics)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Phylogenetic tree of WOX proteins from <span class="html-italic">D. catenatum</span> and <span class="html-italic">A. thaliana</span>. The tree was constructed according to the neighbor-joining method implemented using MEGA 7.0 software and a bootstrap analysis with 1000 iterations.</p>
Full article ">Figure 2
<p><span class="html-italic">D. catenatum WOX</span> gene structures and the conserved motifs in the encoded proteins. (<b>A</b>) Alignment of the homeodomain. (* Indicates the separation of every ten amino acid residues) (<b>B</b>) Phylogenetic tree of <span class="html-italic">D. catenatum</span> WOX proteins. (<b>C</b>) <span class="html-italic">DcaWOX</span> gene structure. (<b>D</b>) Cis-acting elements in the <span class="html-italic">DcaWOX</span> gene promoters. (<b>E</b>) Homeodomain revealed by SMART. (<b>F</b>) Conserved motifs in <span class="html-italic">D. catenatum</span> WOX proteins identified by MEME. Each motif is indicated by a specific color. (<b>G</b>) Three-dimensional conformation of <span class="html-italic">D. catenatum</span> WOX proteins.</p>
Full article ">Figure 2 Cont.
<p><span class="html-italic">D. catenatum WOX</span> gene structures and the conserved motifs in the encoded proteins. (<b>A</b>) Alignment of the homeodomain. (* Indicates the separation of every ten amino acid residues) (<b>B</b>) Phylogenetic tree of <span class="html-italic">D. catenatum</span> WOX proteins. (<b>C</b>) <span class="html-italic">DcaWOX</span> gene structure. (<b>D</b>) Cis-acting elements in the <span class="html-italic">DcaWOX</span> gene promoters. (<b>E</b>) Homeodomain revealed by SMART. (<b>F</b>) Conserved motifs in <span class="html-italic">D. catenatum</span> WOX proteins identified by MEME. Each motif is indicated by a specific color. (<b>G</b>) Three-dimensional conformation of <span class="html-italic">D. catenatum</span> WOX proteins.</p>
Full article ">Figure 2 Cont.
<p><span class="html-italic">D. catenatum WOX</span> gene structures and the conserved motifs in the encoded proteins. (<b>A</b>) Alignment of the homeodomain. (* Indicates the separation of every ten amino acid residues) (<b>B</b>) Phylogenetic tree of <span class="html-italic">D. catenatum</span> WOX proteins. (<b>C</b>) <span class="html-italic">DcaWOX</span> gene structure. (<b>D</b>) Cis-acting elements in the <span class="html-italic">DcaWOX</span> gene promoters. (<b>E</b>) Homeodomain revealed by SMART. (<b>F</b>) Conserved motifs in <span class="html-italic">D. catenatum</span> WOX proteins identified by MEME. Each motif is indicated by a specific color. (<b>G</b>) Three-dimensional conformation of <span class="html-italic">D. catenatum</span> WOX proteins.</p>
Full article ">Figure 3
<p>Subcellular localization of DcaWOX-GFP in <span class="html-italic">Nicotiana benthamiana</span> leaves. DcaWOX2-GFP, DcaWOX3a-GFP, DcaWOX5-GFP, DcaWOX9-GFP, DcaWOX11b, and DcaWOX13a were localized in the nucleus. Bar = 50 μm.</p>
Full article ">Figure 4
<p>Analysis of the effects of DcaWOX2, DcaWOX3a, DcaWOX5, DcaWOX9, DcaWOX11b, and DcaWOX13a. <span class="html-italic">TPL</span> expression was activated by (<b>A</b>) DcaWOX2, (<b>B</b>) DcaWOX3a, and (<b>C</b>) DcaWOX5. <span class="html-italic">TPL</span> expression was not activated by (<b>D</b>) DcaWOX9, (<b>E</b>) DcaWOX11b, and (<b>F</b>) DcaWOX13a. a, b indicates significant differences among treatments (<span class="html-italic">n</span> = 5).</p>
Full article ">Figure 5
<p>Expression profiles of <span class="html-italic">D. catenatum WOX</span> genes. (<b>A</b>) <span class="html-italic">DcaWOX</span> genes expression in different tissues. (<b>B</b>) <span class="html-italic">DcaWOX</span> genes expression after treatments with ABA, IAA, JA, and SA. ABA: abscisic acid; IAA: indole−3−acetic acid; JA: jasmonic acid; SA: salicylic acid; ck: control check; 3 h and 6 h: time-points after initiating the phytohormone treatments. (<span class="html-italic">n</span> = 3).</p>
Full article ">
15 pages, 4303 KiB  
Article
Metabolic Profiling of Terpene Diversity and the Response of Prenylsynthase-Terpene Synthase Genes during Biotic and Abiotic Stresses in Dendrobium catenatum
by Xinqiao Zhan, Yichun Qian and Bizeng Mao
Int. J. Mol. Sci. 2022, 23(12), 6398; https://doi.org/10.3390/ijms23126398 - 7 Jun 2022
Cited by 10 | Viewed by 2599
Abstract
Dendrobium catenatum is a widely cultivated Chinese orchid herb rich in abundant secondary metabolites, such as terpenes. However, terpene distribution and characterization of terpene biosynthesis-related genes remain unknown in D. catenatum. In this study, metabolic profiling was performed to analyze [...] Read more.
Dendrobium catenatum is a widely cultivated Chinese orchid herb rich in abundant secondary metabolites, such as terpenes. However, terpene distribution and characterization of terpene biosynthesis-related genes remain unknown in D. catenatum. In this study, metabolic profiling was performed to analyze terpene distribution in the root, stem, leaf, and flower of D. catenatum. A total of 74 terpene compounds were identified and classified. Clustering analysis revealed that terpene compounds exhibited a tissue-specific accumulation, including monoterpenes in the flowers, sesquiterpenes in the stems, and triterpenes in the roots. Transcriptome analysis revealed that the ‘terpenoid backbone biosynthesis’ pathway was only significantly enriched in root vs. flower. The expression of terpene biosynthesis-related genes was spatiotemporal in the flowers. Prenylsynthase-terpene synthases (PS-TPSs) are the largest and core enzymes for generating terpene diversity. By systematic sequence analysis of six species, 318 PS-TPSs were classified into 10 groups and 51 DcaPS-TPSs were found in eight of them. Eighteen DcaPS-TPSs were regulated by circadian rhythm under drought stress. Most of the DcaPS-TPSs were influenced by cold stress and fungi infection. The cis-element of the majority of the DcaPS-TPS promoters was related to abiotic stress and plant development. Methyl jasmonate levels were significantly associated with DcaTPSs expression and terpene biosynthesis. These results provide insight into further functional investigation of DcaPS-TPSs and the regulation of terpene biosynthesis in Dendrobium. Full article
Show Figures

Figure 1

Figure 1
<p>The distribution of terpenes in four tissues of <span class="html-italic">D</span>. <span class="html-italic">catenatum</span>. (<b>a</b>) Four tissues were used for terpene metabolic analysis. (<b>b</b>) The total terpene content in four tissues of <span class="html-italic">D</span>. <span class="html-italic">catenatum</span>. Values are means ± S.D. (<span class="html-italic">n</span> = 6). Student’s <span class="html-italic">t</span>-test, ** <span class="html-italic">p</span> &lt; 0.01. (<b>c</b>) PCA analysis of volatile and non-volatile terpenes. (<b>d</b>) Terpene distribution in four tissues of <span class="html-italic">D</span>. <span class="html-italic">catenatum</span>. The color bar represents the normalization for log<sub>2</sub>-metabolite intensity using the Pheatmap software package. (<b>e</b>) Volcano plot showing the DAMs in root vs. flower. (<b>f</b>) Volcano plot showing the DAMs in stem vs. leaf. (<b>g</b>) Top ten DAMs in root vs. flower. (<b>h</b>) Top ten DAMs in stem vs. leaf.</p>
Full article ">Figure 2
<p>Enrichment analysis of terpene biosynthesis pathway in <span class="html-italic">D</span>. <span class="html-italic">catenatum</span> growth. (<b>a</b>) GSEA for root vs. flower comparison. Enrichment scores reveal the expression levels of enriched genes. <span class="html-italic">p</span> value &lt; 0.05 is considered significant. (<b>b</b>) Overview of the two terpene biosynthetic pathways in plants. PS-TPSs are highlighted in black bold. (<b>c</b>) Expression profiles of terpene biosynthesis-related genes in four tissues. (<b>d</b>) Expression profiles of terpene biosynthesis-related genes during the three flowering phases. The color bar represents the normalization for log<sub>2</sub>-FPKM using the Pheatmap software package. Intermediates abbreviations are: DMAPP, Dimethylallyl pyrophosphate; DOXP, 1-deoxy-D-xylulose 5-phosphate; FPP, farnesyl pyrophosphate; GPP, geranyl pyrophosphate; GGPP, geranylgeranyl pyrophosphate; IPP, isopentenyl pyrophosphate; SPP, solanesyl pyrophosphate. Enzymes abbreviations are: AACT, acetoacetyl-CoA thiolase; DXR, DOXP reductoisomerase; DXS, DOXP synthase; FPS, FPP synthase; FK, farnesol kinase; GcpE, (E)-4-hydroxy-3-methylbut-2-enyl-pyrophosphate synthase; GPS, GPP synthase; GGPS, GGPP synthase; HMGR, HMG-CoA reductase; HMGS, HMG-CoA synthase; IDI, isopentenyl pyrophosphate isomerase; IspD, 2-C-methyl-D-erythritol 4-phosphate cytidylyltransferase; IspE, 4-diphosphocytidyl-2-C-methyl-D-erythritol kinase; IspF, 2-C-methyl-D-erythritol 2,4-cyclopyrophosphate synthase; IspH, 4-hydroxy-3-methylbut-2-en-1-yl pyrophosphate reductase; IspS, isoprene synthase; RCE1, prenyl protein peptidase; SPS, SPP synthase; TPS, terpene synthase.</p>
Full article ">Figure 3
<p>Phylogenetic analysis of PS-TPS proteins in <span class="html-italic">Arabidopsis thaliana</span>, <span class="html-italic">Dendrobium catenatum</span>, <span class="html-italic">Oryza sativa</span>, <span class="html-italic">Phalaenopsis equestris</span>, <span class="html-italic">Selaginella moellendorffi</span>, and <span class="html-italic">Chlamydomonas reinhardtii</span>. A total of 71 PSs and 247 TPSs were used to construct the unrooted maximum likelihood phylogenies.</p>
Full article ">Figure 4
<p>Domain organization of DcaPS-TPS. (<b>a</b>) The conserved motifs of DcaPSs. (<b>b</b>) The conserved motifs of DcaTPSs. Asp-rich motifs (DDxxD) are found in DcaPS-TPS, where ‘x’ is any amino acid.</p>
Full article ">Figure 5
<p>Spatiotemporal expression patterns of <span class="html-italic">DcaPS-TPS</span> genes. (<b>a</b>) Expression patterns of <span class="html-italic">DcaPS-TPS</span> genes in different tissues. (<b>b</b>) Expression patterns of <span class="html-italic">DcaPS-TPS</span> genes in three flowering phases. In the early developmental stage (F1), the column was short and lip crest were green; in the medium stage of flower bud development (F2), the column and lip had purple pigmentation; when the flower matured and opened (F3), the sepals and petals turned yellow and red pigmentation. Color scale represents the value of the relative expression level log<sub>2</sub>-FPKM.</p>
Full article ">Figure 6
<p>The response of <span class="html-italic">DcaPS-TPS</span> genes under abiotic and biotic stresses. (<b>a</b>) Expression of <span class="html-italic">DcaPS-TPS</span> genes under drought stress. The seedlings were watered on the 1st day, dried from the 2nd to the 7th day, and then rewatered on the 8th day. Leaves were collected at different times; DR5/DR8, DR6/DR10, and DR7/DR15 indicate sampling at 06:30 and 18:30 on the 2nd, 7th, and 9th days, respectively, and DR11 indicates sampling at 18:30 on the 8th day. (<b>b</b>) Expression of <span class="html-italic">DcaPS-TPS</span> genes under cold stress. Two-year-old cultivated plants were placed at −6 ℃ for freezing treatment (FT) and 8 °C for post freezing-recovery (FR). (<b>c</b>) Expression of <span class="html-italic">DcaPS-TPS</span> genes during <span class="html-italic">Colletotrichum gloeosporioides</span> infection. Two-year-old cultivated plants were infected with <span class="html-italic">C. gloeosporioides</span> for 15 days. Color scale represents the value of the relative expression level log<sub>2</sub>-FPKM.</p>
Full article ">Figure 7
<p>Predicted cis-elements in <span class="html-italic">DcaPS-TPS</span> promoters and the response of <span class="html-italic">DcaPS-TPS</span> under MeJA treatment. (<b>a</b>) Promoter sequences (−2000 bp) of <span class="html-italic">DcaPS-TPS</span> genes were analyzed by PlantCARE. The classification and number of cis elements were shown on top and left, respectively. (<b>b</b>) The endogenous JA content of in Pr and control plants. Values are means ± S.D. (<span class="html-italic">n</span> = 6). Student’s <span class="html-italic">t</span>-test, * <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> &lt; 0.01. (<b>c</b>) The total terpene content of in Pr and control plants. Values are means ± S.D. (<span class="html-italic">n</span> = 6). Student’s <span class="html-italic">t</span>-test, * <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> &lt; 0.01. (<b>d</b>) The expression levels of six <span class="html-italic">DcaPS-TPS</span> genes under MeJA treatment. Values are means ± S.D. (<span class="html-italic">n</span> = 3). n, not detectable.</p>
Full article ">
Back to TopTop