The novel developed microsatellite markers revealed potential hybridization among Cymbidium species and the interspecies sub-division of C. goeringii and C. ensifolium

BACKGROUND

Orchids (Cymbidium spp.) exhibit significant variations in floral morphology, pollinator relations, and ecological habitats. Due to their exceptional economic and ornamental value, Cymbidium spp. have been commercially cultivated for centuries. SSR markers are extensively used genetic tools for biology identification and population genetics analysis.

RESULT

In this study, nine polymorphic EST-SSR loci were isolated from Cymbidium goeringii using RNA-Seq technology. All nine SSR loci showed transferability in seven other congeneric species, including 51 cultivars. The novel SSR markers detected inter-species gene flow among the Cymbidium species and intra-species sub-division of C. goeringii and C. ensifolium, as revealed by neighborhood-joining and Structure clustering analyses.

CONCLUSION

In this study, we developed nine microsatellites using RNA-Seq technology. These SSR markers aided in detecting potential gene flow among Cymbidium species and identified the intra-species sub-division of C. goeringii and C. ensifolium.

Genome-based identification of the CYP75 gene family in Orchidaceae and its expression patterns in Cymbidium goeringii

With a great diversity of species, Orchidaceae stands out as an essential component of plant biodiversity, making it a primary resource for studying angiosperms evolution and genomics. This study focuses on 13 published orchid genomes to identify and analyze the CYP75 gene family belonging to the cytochrome P450 superfamily, which is closely related to flavonoid biosynthetic enzymes and pigment regulation. We found 72 CYP75s in the 13 orchid genomes and further classified them into two classes: CYP75A and CYP75B subfamily, the former synthesizes blue anthocyanins, while the latter is involved in the production of red anthocyanins. Furthermore, the amount of CYP75Bs (53/72) greatly exceeds the amount of CYP75As (19/72) in orchids. Our findings suggest that CYP75B genes have a more important evolutionary role, as red plants are more common in nature than blue plants. We also discovered unique conserved motifs in each subfamily that serve as specific recognition features (motif 19 belong to CYP75A; motif 17 belong to CYP75B). Two diverse-colored varieties of C. goeringii were selected for qRT-PCR experiments. The expression of CgCYP75B1 was significantly higher in the purple-red variant compared to the yellow-green variant, while CgCYP75A1 showed no significant difference. Based on transcriptomic expression analysis, CYP75Bs are more highly expressed than CYP75As in floral organs, especially in colorful petals and lips. These results provide valuable information for future studies on CYP75s in orchids and other angiosperms.

Microsatellite Dataset for Cultivar Discrimination in Spring Orchid (Cymbidium goeringii)

Cymbidium goeringii Reichb. fil., locally known as the spring orchid in the Republic of Korea, is one of the most important and popular horticultural species in the family Orchidaceae. C. goeringii cultivars originated from plants with rare phenotypes in wild mountains where pine trees commonly grow. This study aimed to determine the cultivar-specific combined genotypes (CGs) of short sequence repeats (SSRs) by analyzing multiple samples per cultivar of C. goeringii. In this study, we collected more than 4000 samples from 67 cultivars and determined the genotypes of 12 SSRs. Based on the most frequent combined genotypes (CG1s), the average observed allele number and combined matching probability were 11.8 per marker and 3.118 × 10-11, respectively. Frequencies of the CG1 in 50 cultivars (n ≥ 10) ranged from 40.9% to 100.0%, with an average of 70.1%. Assuming that individuals with the CG1 are genuine in the corresponding cultivars, approximately 30% of C. goeringii on the farms and markets may be not genuine. The dendrogram of the phylogenetic tree and principal coordinate analysis largely divided the cultivars into three groups according to their countries of origin; however, the genetic distances were not great among the cultivars. In conclusion, this dataset of C. goeringii cultivar-specific SSR profiles could be used for ecogenetic studies and forensic authentication. This study suggests that genetic authentication should be introduced for the sale of expensive C. goeringii cultivars. We believe that this study will help establish a genetic method for the forensic authentication of C. goeringii cultivars.

Genome-wide analysis of the TCP gene family and their expression pattern in Cymbidium goeringii

TCP gene family are specific transcription factors for plant, and considered to play an important role in development and growth. However, few related studies investigated the TCP gene trait and how it plays a role in growth and development of Orchidaceae. In this study, we obtained 14 TCP genes (CgTCPs) from the Spring Orchid Cymbidium goeringii genome. The classification results showed that 14 CgTCPs were mainly divided into two clades as follows: four PCF genes (Class I), nine CIN genes and one CYC gene (Class II). The sequence analysis showed that the TCP proteins of C. goeringii contain four conserved regions (basic Helix-Loop-Helix) in the TCP domain. The exon-intron structure varied in the clade according to a comparative investigation of the gene structure, and some genes had no introns. There are fewer CgTCP homologous gene pairs compared with Dendrobium catenatum and Phalaenopsis equestris, suggesting that the TCP genes in C. goeringii suffered more loss events. The majority of the cis-elements revealed to be enriched in the function of light responsiveness, followed by MeJA and ABA responsiveness, demonstrating their functions in regulating by light and phytohormones. The collinearity study revealed that the TCPs in D. catenatum, P. equestris and C. goeringii almost 1:1. The transcriptomic data and real-time reverse transcription-quantitative PCR (RT-qPCR) expression profiles showed that the flower-specific expression of the TCP class II genes (CgCIN2, CgCIN5 and CgCIN6) may be related to the regulation of florescence. Altogether, this study provides a comprehensive analysis uncovering the underlying function of TCP genes in Orchidaceae.

Genome-wide identification and analysis of anthocyanin synthesis-related R2R3-MYB genes in Cymbidium goeringii

The MYB gene family plays a vital regulatory role in plant metabolism, stress response, and floral color. The R2R3-MYB gene family of C. goeringii was identified, and its expression was analyzed using bioinformatics in this article. The R2R3-MYB genes of Arabidopsis thaliana were used as a reference to determine 104 CgMYB genes and categorize them into 22 subfamilies. Exon/intron organizations and conserved motif analysis revealed that the majority of CgMYB genes were highly conserved, and chromosome localization and collinearity analysis provided evidence of tandem duplication and segmental duplication events, indicating the phenomenon of gene family expansion and contraction. The function of CgMYB genes was analyzed by cis-acting element and gene ontology (GO) enrichment. In addition, we selected CgMYB91 and CgMYB32 for RT-qPCR, suggesting that CgMYB91 and CgMYB32 are associated with anthocyanin formation. In short, this study provides a comprehensive and specific function of the R2R3-MYB transcription factors (TFs) in orchids.

miR390 family of Cymbidium goeringii is involved in the development of reproductive organs in transgenic Arabidopsis

BACKGROUND

miR390s is an ancient family with a high level of conservation among plant miRNAs. Through the auxin signal transduction pathway, miR390 participates in diverse biological processes of plant growth and development. As an important Chinese traditional orchid, Cymbidium goeringii has unique flower shape and elegant fragrance. But its development has been greatly restricted because of the low flower bud differentiation and the difficult reproduction. This study aims to provide guidance for the role of cgo-miR390 in reproductive organ development to enhance the ornamental and economic value of Cymbidium.

RESULTS

MIR390a, MIR390b and MIR390c of C. goeringii were cloned, and their length ranged from 130 to 150 nt. Each precursor sequence of cgo-miR390 contains 2 to 3 mature miRNAs. Three kinds of cgo-miR390s displayed distinct temporal and spatial expression patterns during floral development in C. goeringii. The overexpression of MIR390s alters morphology and function of stamens and pistils in Arabidopsis, such as enlargement of anther aspect ratio and separation of stylar and stigmas, which affects the development of fruits and seeds. In particular, the pollen amount decreased and the seed abortion rate increased in cgo-MIR390c-overexpressed plants.

CONCLUSIONS

cgo-miR390 family affected the development of reproductive organs in transgenic Arabidopsis. The study provides references for the genetic improvement for orchid with potentially great economic benefit.

Functional analysis of ARF1 from Cymbidium goeringii in IAA response during leaf development

Background

Cymbidium is an economically important genus of flowering orchids cultivated in China because of showing graceful leaf shapes and elegant flower coloration. However, the deterioration of the ecological environment and the difficulty of conservation management have become increasing challenges for maintaining its germplasm resources. ARFs are critical transcription factors in the auxin signaling pathway and have been found to play pivotal roles in leaf growth and development in previous studies. However, their functions and mechanisms in Cymbidium goeringii remain to be clarified.

Methods

The sequence of the CgARF1 gene was analyzed by bioinformatics. The expression of this gene in different tissues and under IAA treatment was detected by quantitative real-time PCR analysis. The CgARF1 gene was overexpressed in wild-type Arabidopsis and Nicotiana benthamiana via the Agrobacterium infection method. Acetone-ethanol solvent extraction was applied for the determination of chlorophyll contents, and the contents of endogenous hormones were determined using the enzyme-linked immunosorbent assay technique.

Results

CgARF1 cloned from C. goeringii 'Songmei' was 2,049 bp in length and encoded 682 amino acids containing three typical domains: a B3 DNA binding domain, an Aux_resp domain and an AUX/IXX family domain. The expression pattern of CgARF1 in different tissues of C. goeringii showed that its expression was highest in the leaves and changed greatly under IAA treatment. Subcellular localization studies showed that the protein was mainly localized in the cell nucleus. CgARF1-overexpressing lines exhibited leaf senescence and a decreased chlorophyll content. Under IAA treatment, CgARF1 regulates the rooting length, rooting number and rooting rate from detached leaves. The levels of endogenous hormones in transgenic leaves were also significantly changed.

Conclusion

These results indicated that CgARF1 overexpression is responsive to IAA treatment during leaf development. This study provides a foundation for future research on the function of the ARF gene family in C. goeringii.

The Leafless Orchid Cymbidium macrorhizon Performs Photosynthesis in the Pericarp during the Fruiting Season

Photosynthesis with highly photoreactive chlorophyll (Chl) provides energy for plant growth but with simultaneous risk of photooxidative damage and photoprotection costs. Although the leafless orchid Cymbidium macrorhizon mostly depends on mycorrhizal fungi for carbon, it accumulates Chl particularly during fruiting and may not be fully mycoheterotrophic. In fact, stable isotopic analysis suggested that the fruiting C. macrorhizon specimens obtain a significant proportion of its carbon demands through photosynthesis. However, actual photosynthetic characteristics of this leafless orchid are unknown. To reveal the functionality of photosynthetic electron transport in C. macrorhizon, we compared its photosynthetic properties with those of its relative mixotrophic orchid Cymbidium goeringii and the model plant Arabidopsis thaliana. Compared with C. goeringii and A. thaliana, maximum photochemical efficiency of PSII was substantially low in C. macrorhizon. Chl fluorescence induction kinetics revealed that the electron transport capacity of PSII was limited in C. macrorhizon. Chl fluorescence analysis at 77 K suggested partial energetic disconnection of the light-harvesting antenna from the PSII reaction center in C. macrorhizon. Despite its low PSII photochemical efficiency, C. macrorhizon showed photosynthetic electron transport activity both in the field and under laboratory conditions. Cymbidium macrorhizon developed strong nonphotochemical quenching in response to increased light intensity as did C. goeringii, suggesting the functionality of photoprotective systems in this orchid. Moreover, C. macrorhizon fruit developed stomata on the pericarp and showed net O2-evolving activity. Our data demonstrate that C. macrorhizon can perform photosynthetic electron transport in the pericarp, although its contribution to net carbon acquisition may be limited.

Functional analysis of CgWRKY57 from Cymbidium goeringii in ABA response

Background

The orchid is one of the top ten Chinese flowers and has high ornamental value and elegant color. However, orchids are vulnerable to abiotic stresses during their growth and development, and the molecular mechanism of the abiotic stress response in orchids is unclear. WRKY proteins belong to a transcription factor family that plays important roles in biotic stress, abiotic stress, growth and development in plants, but little is known about the WRKY family in Cymbidium goeringii.

Methods

The specific fragment of the CgWRKY57 gene of C. goeringii was analyzed by bioinformatics. The expression of the CgWRKY57 gene of C. goeringii under 4 °C, 42 °C water and ABA stress as well as different tissues was detected by real-time fluorescence quantitative PCR. CgWRKY57 gene was overexpressed in wild type Arabidopsis thaliana by inflorescence infection method, and the function of transgenic lines under ABA stress was analyzed.

Results

CgWRKY57 was cloned from C. goeringii and found to encode 303 amino acids. The CgWRKY57 protein is an acidic, nonsecreted hydrophilic protein without a signal peptide or transmembrane domain. The CgWRKY57 protein is located to the nucleus and may function intracellularly according to its predicted subcellular localization. A domain analysis and homology comparison showed that the CgWRKY57 protein has a "WRKYGQK" domain and belongs to Group III of the WRKY family, and a phylogenetic analysis demonstrated that CgWRKY57 is closely related to OsWRKY47. CgWRKY57 was expressed in the roots, stems, leaves and floral organs of C. goeringii, and its expression level was highest in the roots according to real-time qPCR analysis. There were significant differences in CgWRKY57 expression under 4 °C, 42 °C ABA and water stress treatments, and its expression changed greatly under ABA stress. The expression of CgWRKY57 in transgenic plants was significantly higher than that in wild type plants under ABA stress, and the root length and germination rate were reduced in transgenic plants compared to wild type plants.

Conclusions

These results indicate that CgWRKY57 overexpression is responsive to ABA stress, and they provide a foundation for future analyses of the biological functions of the WRKY family in C. goeringii.

nov., isolated from Cymbidium goeringii roots

A Gram-stain-variable, aerobic, rod-shaped, endospore-forming strain R196^T) was isolated from internal tissues of roots of Cymbidium goeringii. Cells were motile with peritrichous flagella. The colonies were light pink on tryptone soya agar medium. Phylogenetic analysis based on 16S rRNA gene sequences showed that strain R196^T fell into a phylogenetic cluster belonging to the genus Paenibacillus. Strain R196^T was closely related to Paenibacillus cavernae C4-5^T and Paenibacillus contaminans CKOBP-6^T with 93.6 and 93.3% sequence similarities, respectively. The major cellular polar lipids were diphosphatidylglycerol, phosphatidylethanolamine, phosphatidylglycerol, three unidentified phospholipids, two unidentified aminophospholipids and an unidentified aminolipid. The dominant respiratory quinone was MK-7. The main cellular fatty acids were anteiso-C_{15 : 0} (53.01%), C_{16 : 0} (13.04%) and iso-C_{16 : 0} (10.80%). The genome size of R196^T was 9.45 Mb, containing 7617 predicted protein-coding genes, with a DNA G+C content of 57.7 mol%. Based on the results of phenotypic, chemotaxonomic and whole-genome analyses, strain R196^T represents a novel species of the genus Paenibacillus, for which the name Paenibacillus cymbidii sp. nov. is proposed. The type strain is R196^T (=ACCC 61713^T=KCTC 33718^T).

RNA sequencing analysis of Cymbidium goeringii identifies floral scent biosynthesis related genes

BACKGROUND

Cymbidium goeringii belongs to the Orchidaceae, which is one of the most abundant angiosperm families. Cymbidium goeringii consist with high economic value and characteristics include fragrance and multiple flower colors. Floral scent is one of the important strategies for ensuring fertilization. However, limited genetic data is available in this non-model plant, and little known about the molecular mechanism responsible for floral scent in this orchid. Transcriptome and expression profiling data are needed to identify genes and better understand the biological mechanisms of floral scents in this species. Present transcriptomic data provides basic information on the genes and enzymes related to and pathways involved in flower secondary metabolism in this plant.

RESULTS

In this study, RNA sequencing analyses were performed to identify changes in gene expression and biological pathways related scent metabolism. Three cDNA libraries were obtained from three developmental floral stages: closed bud, half flowering stage and full flowering stage. Using Illumina technique 159,616,374 clean reads were obtained and were assembled into 85,868 final unigenes (average length 1194 nt), 33.85% of which were annotated in the NCBI non redundant protein database. Among this unigenes 36,082 were assigned to gene ontology and 23,164 were combined with COG groups. Total 33,417 unigenes were assigned in 127 pathways according to the Kyoto Encyclopedia of Genes and Genomes pathway database. According these transcriptomic data we identified number of candidates genes which differentially expressed in different developmental stages of flower related to fragrance biosynthesis. In q-RT-PCR most of the fragrance related genes highly expressed in half flowering stage.

CONCLUSIONS

RNA-seq and DEG data provided comprehensive gene expression information at the transcriptional level that could be facilitate the molecular mechanisms of floral biosynthesis pathways in three developmental phase's flowers in Cymbidium goeringii, moreover providing useful information for further analysis on C. goeringii, and other plants of genus Cymbidium.

Low-temperature-induced changes in the transcriptome reveal a major role of CgSVP genes in regulating flowering of Cymbidium goeringii

BACKGROUND

Cymbidium goeringii is one of the most horticulturally important and popular ornamental plants in the orchid family (Orchidaceae). It blooms in winter during January-March and a period of low temperature is necessary for its normal flowering, otherwise there is flower bud abortion, which seriously affects the economic benefits. However, the molecular mechanism underlying winter-blooming behavior in C. goeringii is unclear.

RESULTS

In this research, we firstly study the flowering physiology of C. goeringii by cytobiology observations and physiological experiments. Using comparative transcriptome analysis, we identified 582 differentially expressed unigenes responding to cold treatment that were involved in metabolic process, flowering time, hormone signaling, stress response, and cell cycle, implying their potential roles in regulating winter-blooming of C. goeringii. Twelve MADS-box genes among them were investigated by full-length cDNA sequence analysis and expression validation, which indicated that three genes within the SHORT VEGETATIVE PHASE (SVP) sub-group had the most significant repressed expression after cold treatment. Further analysis revealed that the SVP genes showed population variation in expression that correlated with cold-regulated flowering and responded to low temperature earlier than the flowering pathway integrators CgAP1, CgSOC1, and CgLFY, suggesting a potential role of CgSVP genes in the early stage of low-temperature-induced blooming of C. goeringii. Moreover, a yeast two-hybrid experiment confirmed that CgSVP proteins interacted with CgAP1 and CgSOC1, suggesting that they may synergistically control the process of C. goeringii flowering in winter.

CONCLUSIONS

This study represents the first exploration of flowering physiology of C. goeringii and provides gene expression information that could facilitate our understanding of molecular regulation of orchid plant winter-flowering, which could provide new insights and practical guidance for improving their flowering regulation and molecular breeding.

Integrated mRNA and microRNA transcriptome variations in the multi-tepal mutant provide insights into the floral patterning of the orchid Cymbidium goeringii

Background

Cymbidium goeringii is a very famous traditional orchid plant in China, which is well known for its spectacular and diverse flower morphology. In particular, the multi-tepal mutants have considerable ecological and cultural value. However, the current understanding of the molecular mechanisms of floral patterning and multi-tepal development is limited. In this study, we performed expression profiling of both microRNA (miRNA) and mRNA from wild-type and typical multi-tepal-mutant flowers of C. goeringii for the first time, to identify the genes and pathways regulating floral morphogenesis in C. goeringii.

Results

Total clean reads of 98,988,774 and 100,188,534 bp were obtained from the wild-type and mutant library, respectively, and de novo assembled into 98,446 unigenes, with an average length of 989 bp. Among them, 18,489 were identified as differentially expressed genes between the two libraries according to comparative transcript profiling. The majority of the gene ontology terms and Kyoto Encyclopedia of Genes and Genomes pathway enrichment responses were for membrane-building and ploidy-related processes, consistent with the excessive floral organs and altered cell size observed in the mutant. There were 29 MADS-box genes, as well as a large number of floral-related regulators and hormone-responsive genes, considered as candidates regulating floral patterning of C. goeringii. Small RNA sequencing revealed 132 conserved miRNA families expressed in flowers of C. goeringii, and 11 miRNAs corresponding to 455 putative target genes were considered to be responsible for multi-tepal development. Importantly, integrated analysis of mRNA and miRNA sequencing data showed two transcription factor/microRNA-based genetic pathways contributing to the multi-tepal trait: well-known floral-related miR156/SPL and miR167/ARF regulatory modes involved in reproductive organ development; and the miR319/TCP4–miR396/GRF regulatory cascade probably regulating cell proliferation of the multi-tepal development.

Conclusions

Integrated mRNA and miRNA profiling data provided comprehensive gene expression information on the wild-type and multi-tepal mutant at the transcriptional level that could facilitate our understanding of the molecular mechanisms of floral patterning of C. goeringii. These data could also be used as an important resource for investigating the genetics of floral morphogenesis and various biological mechanisms of orchid plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3756-9) contains supplementary material, which is available to authorized users.