Phosphorylation of 399S at CsHsp70 of Cymbidium sinense is essential to maintain chlorophyll stability.

The Chinese orchids symbolise nobility and gentility in China, and the variation of leaf color makes Cymbidium sinense more diversified and valuable. However, its color variations especially at the protein level still remain largely unexplored. In this study, the proteomics and phosphoproteomics of Cymbidium sinense leaf color variation mutants were studied. A total of 1059 differentially abundant proteins (DAPs) and 1127 differentially abundant phosphorylation sites belonging to 644 phosphoproteins (DAPPs) were identified in the yellow section of leaf variegation mutant of Cymbidium sinense (MY) compared with the green section (MG). Moreover, 349 co-expressing proteins were found in both omics' datasets, while only 26 proteins showed the same expression patterns in the two omics. The interaction network analysis of kinases and phosphatases showed that DAPs and DAPPs in photosynthesis, response to hormones, pigment metabolic process, phosphorylation, glucose metabolic process, and dephosphorylation might contribute to leaf color variation. The abundance of 28 Hsps and 28 phosphorylation sites belonging to 10 Hsps showed significant differences between MG and MY. CsHsp70 was selected to explore the function in Cymbidium sinense leaf variegation. The results showed CsHsp70 is essential for maintaining photosynthetic pigment content and the 399S phosphorylation site is crucial to the function of CsHsp70. Collectively, our findings construct a comprehensive coverage of protein and protein phosphorylation in leaf variegation of C. sinense, providing valuable insights into its formation mechanisms.

Genome-Wide Identification and Drought Stress Response Pattern of the NF-Y Gene Family in Cymbidium sinense.

Cymbidium sinense, a type of orchid plant, is more drought-resistant and ornamental than other terrestrial orchids. Research has shown that many members of the NUCLEAR FACTOR Y (NF-Y) transcription factor family are responsive to plant growth, development, and abiotic stress. However, the mechanism of the NF-Y gene family's response to abiotic stress in orchids has not yet been reported. In this study, phylogenetic analysis allowed for 27 CsNF-Y genes to be identified (5 CsNF-YAs, 9 CsNF-YBs, and 13 CsNF-YC subunits), and the CsNF-Ys were homologous to those in Arabidopsis and Oryza. Protein structure analysis revealed that different subfamilies contained different motifs, but all of them contained Motif 2. Secondary and tertiary protein structure analysis indicated that the CsNF-YB and CsNF-YC subfamilies had a high content of alpha helix structures. Cis-element analysis showed that elements related to drought stress were mainly concentrated in the CsNF-YB and CsNF-YC subfamilies, with CsNF-YB3 and CsNF-YC12 having the highest content. The results of a transcriptome analysis showed that there was a trend of downregulation of almost all CsNF-Ys in leaves under drought stress, while in roots, most members of the CsNF-YB subfamily showed a trend of upregulation. Additionally, seven genes were selected for real-time reverse transcription quantitative PCR (qRT-PCR) experiments. The results were generally consistent with those of the transcriptome analysis. The regulatory roles of CsNF-YB 1, 2, and 4 were particularly evident in the roots. The findings of our study may make a great contribution to the understanding of the role of CsNF-Ys in stress-related metabolic processes.

Genome-Wide Identification, Expression, and Molecular Characterization of the CONSTANS-like Gene Family in Seven Orchid Species.

The orchid is one of the most distinctive and highly valued flowering plants. Nevertheless, the CONSTANS-like (COL) gene family plays significant roles in the control of flowering, and its functions in Orchidaceae have been minimally explored. This research identified 68 potential COL genes within seven orchids' complete genome, divided into three groups (groups I, II, and III) via a phylogenetic tree. The modeled three-dimensional structure and the conserved domains exhibited a high degree of similarity among the orchid COL proteins. The selection pressure analysis showed that all orchid COLs suffered a strong purifying selection. Furthermore, the orchid COL genes exhibited functional and structural heterogeneity in terms of collinearity, gene structure, cis-acting elements within their promoters, and expression patterns. Moreover, we identified 50 genes in orchids with a homology to those involved in the COL transcriptional regulatory network in Arabidopsis. Additionally, the first overexpression of CsiCOL05 and CsiCOL09 in Cymbidium sinense protoplasts suggests that they may antagonize the regulation of flowering time and gynostemium development. Our study will undoubtedly provide new resources, ideas, and values for the modern breeding of orchids and other plants.

Identification of the CONSTANS-like family in Cymbidium sinense, and their functional characterization.

BACKGROUND: Cymbidium sinense is an orchid that is typically used as a potted plant, given its high-grade ornamental characteristics, and is most frequently distributed in China and SE Asia. The inability to strictly regulate flowering in this economically important potted and cut-flower orchid is a bottleneck that limits its industrial development. Studies on C. sinense flowering time genes would help to elucidate the mechanism regulating flowering. There are very few studies on the genetic regulation of flowering pathways in C. sinense. Photoperiod significantly affects the flowering of C. sinense, but it was unknown how the CONSTANS gene family is involved in regulating flowering. RESULTS: In this study, eight CONSTANS-like genes were identified and cloned. They were divided into three groups based on a phylogenetic analysis. Five representative CsCOL genes (CsCOL3/4/6/8/9) were selected from the three groups to perform expression characterization and functional study. CsCOL3/4/6/8/9 are nucleus-localized proteins, and all five CsCOL genes were expressed in all organs, mainly in leaves followed by sepals. The expression levels of CsCOL3/4 (group I) were higher in all organs than other CsCOL genes. Developmental stage specific expression revealed that the expression of CsCOL3/4/9 peaked at the initial flowering stage. In contrast, the transcript level of CsCOL6/8 was highest at the pedicel development stage. Photoperiodic experiments demonstrated that the transcripts of the five CsCOL genes exhibited distinct diurnal rhythms. Under LD conditions, the overexpression of CsCOL3/4 promoted early flowering, and CsCOL6 had little effect on flowering time, whereas CsCOL8 delayed flowering of Arabidopsis thaliana. However, under SD conditions, overexpression of CsCOL4/6/8 promoted early flowering and the rosette leaves growth, and CsCOL3 induced flower bud formation in transgenic Arabidopsis. CONCLUSION: The phylogenetic analysis, temporal and spatial expression patterns, photoperiodic rhythms and functional study indicate that CsCOL family members in C. sinense were involved in growth, development and flowering regulation through different photoperiodic pathway. The results will be useful for future research on mechanisms pertaining to photoperiod-dependent flowering, and will also facilitate genetic engineering-based research that uses Cymbidium flowering time genes.

Genome-Wide Identification Analysis of the R2R3-MYB Transcription Factor Family in Cymbidium sinense for Insights into Drought Stress Responses.

Cymbidium sinense represents a distinctive Orchidaceae plant that is more tolerant than other terrestrial orchids. Studies have shown that many members of the MYB transcription factor (TF) family, especially the R2R3-MYB subfamily, are responsive to drought stress. This study identified 103 CsMYBs; phylogenetic analysis classified these genes into 22 subgroups with Arabidopsis thaliana. Structural analysis showed that most CsMYB genes contained the same motifs, three exons and two introns, and showed a helix-turn-helix 3D structure in each R repeat. However, the members of subgroup 22 contained only one exon and no intron. Collinear analysis revealed that C. sinense had more orthologous R2R3-MYB genes with wheat than A. thaliana and rice. Ka/Ks ratios indicated that most CsMYB genes were under purifying negative selection pressure. Cis-acting elements analysis revealed that drought-related elements were mainly focused on subgroups 4, 8, 18, 20, 21, and 22, and Mol015419 (S20) contained the most. The transcriptome analysis results showed that expression patterns of most CsMYB genes were upregulated in leaves in response to slight drought stress and downregulated in roots. Among them, members in S8 and S20 significantly responded to drought stress in C. sinense. In addition, S14 and S17 also participated in these responses, and nine genes were selected for the real-time reverse transcription quantitative PCR (RT-qPCR) experiment. The results were roughly consistent with the transcriptome. Our results, thus, provide an important contribution to understanding the role of CsMYBs in stress-related metabolic processes.

Genome-wide identification of Cymbidium sinense WRKY gene family and the importance of its Group III members in response to abiotic stress.

Transcription factors (TFs) of the WRKY family play pivotal roles in defense responses and secondary metabolism of plants. Although WRKY TFs are well documented in numerous plant species, no study has performed a genome-wide investigation of the WRKY gene family in Cymbidium sinense. In the present work, we found 64 C. sinense WRKY (CsWRKY) TFs, and they were further divided into eight subgroups. Chromosomal distribution of CsWRKYs revealed that the majority of these genes were localized on 16 chromosomes, especially on Chromosome 2. Syntenic analysis implied that 13 (20.31%) genes were derived from segmental duplication events, and 17 orthologous gene pairs were identified between Arabidopsis thaliana WRKY (AtWRKY) and CsWRKY genes. Moreover, 55 of the 64 CsWRKYs were detectable in different plant tissues in response to exposure to plant hormones. Among them, Group III members were strongly induced in response to various hormone treatments, indicating their potential essential roles in hormone signaling. We subsequently analyzed the function of CsWRKY18 in Group III. The CsWRKY18 was localized in the nucleus. The constitutive expression of CsWRKY18 in Arabidopsis led to enhanced sensitivity to ABA-mediated seed germination and root growth and elevated plant tolerance to abiotic stress within the ABA-dependent pathway. Overall, our study represented the first genome-wide characterization and functional analysis of WRKY TFs in C. sinense, which could provide useful clues about the evolution and functional description of CsWRKY genes.

First Report of stem rot in Cymbidium sinense Caused by Fusarium oxysporum (FOSC) in China.

Cymbidium sinense (Jackson ex Andr.) Willd is a perennial terrestrial plant in the orchid family mainly distributed in China, Japan, India and Southeast Asia that occupies a strong position in the flower market due to its bright green leaves and fragrant flowers (Zhang et al. 2013). Cymbidium sinense is not only valued by people for its ornamental and economic value, but its roots have antiasthmatic medicinal properties (Ke et al. 2004). In August 2020, about 15% stem rot on two-year old C. sinense with varying severity was observed in five nursery gardens located in Enshi city (N 30° 16', E 109° 29'), Hubei province, China. Typical symptoms of C. sinense included roots and inner part of the pseudobulbs changing from white to brown and rotting. Leaves became brown and withered from bottom to top, and there was an obvious blight yellow halo at the junction of diseased and healthy tissue, which eventually caused the whole plant to wilt and die (Fig. 1d). To isolate the pathogen, a total of 15 leaf tissues from the disease-health junction (3 × 3 mm) from 5 individual plants (3 leaves/plant) with symptoms were surface sterilized with 75% ethanol for 30 s and 2% sodium hypochlorite (NaOCl) for 3 min. The sterilized tissue was rinsed three times with sterilized water, and then placed on potato dextrose agar (PDA) for incubation at 28°C in the dark for 5 days. Isolated colonies were subcultured by a hyphal tip protocol. Thirteen fungal isolates were obtained. Through preliminary pathogenicity tests, we found that ten isolates induced leaf blight. These ten isolates with pathogenicity showed similar morphological characteristics, with initial white-flocculent aerial mycelium that secreted a lavender pigment and produced colonies with an irregular edge after 3 days on PDA. The ten strains were cultured on PDA plates at 28℃ for 5 and 15 days to observe colony and conidial characteristics. The ten strains were identified as Fusarium based on morphological characteristics (Leslie and Summerell 2006). Strain ML0303 was selected for further identification. Macroconidia were falciform, hyaline, slightly pointed at both ends with two to four septa, 24.0 ± 5.6 µm × 4.7 ± 0.8 µm (n = 50). Microconidia were hyaline, oval, globose, with zero to one septum, 5.5 ± 1.3 µm × 2.2 ± 0.5 µm (n = 50) (Fig. 1c). Total genomic DNA of strain ML0303 was extracted with a CTAB protocol (Stenglein and Balatti 2006). The translation elongation factor (EF-1α), RNA polymerase II second largest subunit (RPB2) and ß-tubulin (Tub2) genes were amplified respectively using primer pairs EF1/EF2, RPB2-5F2/RPB2-7cR and T1/T22 respectively (O'Donnell. et al. 2010, O'Donnell. et al. 1997). The EF-1α, RPB2 and Tub2 (accession numbers-MW719874, OL614838, OL689398, respectively) gene sequences were submitted to GenBank. EF-1α, RPB2 and Tub2 sequences of ML0303 showed 99.5% - 100% identity respectively with Fusarium oxysporum in the Genbank and FUSARIUM-ID databases. The multilocus sequence data was used to infer a phylogenetic tree via a Neighbor-joining (NJ), Maximum-likelihood (ML) and Maximum-Parsimony(MP) together with reference sequences from GenBank. The topology of the three trees was similar; only the NJ tree is presented here. Strain ML0303 and F. oxysporum formed a clade supported with high values (NJ/ML/MP: 96,95,97). The results indicated that the fungus was F. oxysporum based on the phylogenetic analysis and BLASTn queries. For pathogenicity tests, conidia of strain ML0303 were collected by rinsing PDA plates. Two-year-old C. sinense grown in plastic pots filled with sterilized autoclaved sandy loam soil were used for the tests. Three pots (two plants/pot) were included in each treatment. Spore suspensions (106spores/ml) of strain ML0303 were used to irrigate the stem-zone of the plants, and sterile water was used as control. The two treatments were placed in a greenhouse and incubated at 28±2℃ with a 14-hour light/10-hour dark cycle. The experiment was repeated twice. After three weeks, stem rot symptoms were observed on C. sinense inoculated with ML0303, that were the as same as observed in the nursery (Fig. 1e-h). No symptoms were observed on the negative control. Fusarium oxysporum was re-isolated from the infected plants to fulfill Koch's postulates. Partial EF-1α and RPB2 gene sequences were used for molecular identification. Members of the FOSC are notorious for causing many diseases, which includes stem rot of Sulcorebutia heliosa and root rot of Torreya grandis (Garibaldi et al. 2020; Zhang et al. 2016). To our knowledge, this is the first report of stem rot by F. oxysporum on C. sinense in China. The finding of this pathogen provides a clear target for stem rot control.

Selection and validation of reference genes for quantitative real-time PCR in Cymbidium sinense.

Selection of reference genes (RGs) is important for the accurate analysis of real-time quantitative PCR (RT-qPCR) results. This study screened RGs of Cymbidium sinense for more accurate quantification of target genes. The two most stable RGs for all tissues were ACT and EF1α, those for vegetative organs were UBQ3 and ACT, while those for reproductive organs as well as organs in the full flowering stage were EF1α and ACT. The AGAMOUS (CsAG1) expression level was verified using EF1α, ACT, GAPDH, UBQ2 and UBQ3 as RG. The expression profile of CsAG1 was consistent when normalized with EF1α, ACT and UBQ3. The results have practical value for the expression of key genes during the development of C. sinense.

The genome of Cymbidium sinense revealed the evolution of orchid traits.

The Orchidaceae is of economic and ecological importance and constitutes ˜10% of all seed plant species. Here, we report a genome physical map for Cymbidium sinense, a well-known species belonging to genus Cymbidium that has thousands of natural variation varieties of flower organs, flower and leaf colours and also referred as the King of Fragrance, which make it arose into a unique cultural symbol in China. The high-quality chromosome-scale genome assembly was 3.52 Gb in size, 29 638 protein-coding genes were predicted, and evidence for whole-genome duplication shared with other orchids was provided. Marked amplification of cytochrome- and photosystem-related genes was observed, which was consistent with the shade tolerance and dark green leaves of C. sinense. Extensive duplication of MADS-box genes, and the resulting subfunctional and expressional differentiation, was associated with regulation of species-specific flower traits, including wild-type and mutant-type floral patterning, seasonal flowering and ecological adaption. CsSEP4 was originally found to positively regulate gynostemium development. The CsSVP genes and their interaction proteins CsAP1 and CsSOC1 were significantly expanded and involved in the regulation of low-temperature-dependent flowering. Important genetic clues to the colourful leaf traits, purple-black flowers and volatile trait in C. sinense were also found. The results provide new insights into the molecular mechanisms of important phenotypic traits of Cymbidium and its evolution and serve as a powerful platform for future evolutionary studies and molecular breeding of orchids.

Comparative Transcriptomic and Metabolic Analyses Reveal the Molecular Mechanism of Ovule Development in the Orchid, Cymbidium sinense.

Ovule development is pivotal to plant reproduction and seed development. Cymbidium sinense (Orchidaceae) has high ornamental value due to its pleasant aroma and elegant floral morphology. The regulatory mechanism underlying ovule development in orchids, especially C. sinense, is largely unknown and information on the C. sinense genome is very scarce. In this study, a combined analysis was performed on the transcriptome and non-targeted metabolomes of 18 C. sinense 'Qi Jian Hei Mo' ovule samples. Transcriptome analysis assembled gene-related information related to six growth stages of C. sinense ovules (S1-S6, equivalent to 30, 35, 42, 46, 53, and 60 days after pollination). Illumina sequencing technology was used to obtain the complete set of transcriptome sequences of the 18 samples. A total of 81,585 unigene sequences were obtained after assembly, 24,860 (30.47%) of which were functionally annotated. Using transcriptome sequencing technology, a total of 9845 differentially expressed unigenes (DEUs) were identified in C. sinense ovules that were assigned to specific metabolic pathways according to the Kyoto Encyclopedia of Genes and Genomes (KEGG). DEUs associated with transcription factors (TFs) and phytohormones were identified and analyzed. The TFs homeobox and MADS-box were associated with C. sinense ovule development. In particular, the phytohormones associated with DEUs such as indole-3-acetic acid (IAA), cytokinin (CK), gibberellin (GA), abscisic acid (ABA), brassinosteroid (BR), and jasmonate (JA), may have important regulatory effects on C. sinense ovule development. Metabolomic analysis showed an inconsistent number of KEGG annotations of differential metabolites across comparisons (S2_vs_S4, S2_vs_S5, and S4_vs_S5 contained 23, 26, and 3 annotations, respectively) in C. sinense ovules. This study provides a valuable foundation for further understanding the regulation of orchid ovule development and formation, and establishes a theoretical background for future practical applications during orchid cultivation.

Complete chloroplast genome sequence of an orchid hybrid Cymbidium sinense (♀) × C. goeringii (♂).

The complete chloroplast genome sequence of the Cymbidium hybrid, C. sinense (♀) × C. goeringii (♂) was assembled in this study. The circular genome was 150,149 bp in length with an overall GC content of 37.1% and consisted of a pair of 25,691 bp inverted repeats, and two single-copy regions that were 84,987 bp and 13,780 bp, respectively. Gene annotation analysis identified 109 genes including 75 protein-coding genes, 30 transfer RNA, and 4 ribosomal RNA genes. Phylogenetic analysis showed its closest relationship to Cymbidium sinense, reflecting a maternal inheritance of chloroplasts.

Comparative Metabolomic Analysis Reveals Distinct Flavonoid Biosynthesis Regulation for Leaf Color Development of Cymbidium sinense 'Red Sun'.

The colorful leaf is an important ornamental character of Cymbidium sinense (C. sinense), especially the red leaf, which has always been attracted by breeders and consumers. However, little is documented on the formation mechanism of the red leaf of C. sinense. In this study, the changing patterns of flavonoid-related metabolites, corresponding enzyme activities and genes expression in the leaves of C. sinense 'Red Sun' from red to yellow and finally to green was investigated. A total of 196 flavonoid-related metabolites including 11 anthocyanins metabolites were identified using UPLC-MS/MS-based approach. In the process of leaf color change, 42 metabolites were identified as having significantly different contents and the content of 28 differential metabolites turned to zero. In anthocyanin biosynthetic pathway, content of all 15 identified metabolites showed downregulation trend in the process of leaf color change. Among the 15 metabolites, the contents of Naringenin chalcone, Pelargonidin O-acetylhexoside and Anthocyanin 3-O-beta-d-glucoside decreased to zero in the green leaf stage. The changing pattern of enzyme activity of 10 enzymes involved in the anthocyanin biosynthetic pathway showed different trends from red leaves that have turned yellow and finally green, while the expression of genes encoding these enzymes was all down-regulated in the process of leaf color change. The results of this study revealed the types of flavonoid-related metabolites and the comprehensive analysis of metabolites content, enzyme activities and genes expression providing a new reference for breeders to improve the leaf color of C. sinense 'Red Sun'.