Organ-Specific Gene Expression Reveals the Role of the Cymbidium ensifolium-miR396/Growth-Regulating Factors Module in Flower Development of the Orchid Plant Cymbidium ensifolium

Identification of microRNAs and their target genes associated with chasmogamous and cleistogamous flower development in Viola prionantha

MAIN CONCLUSION

Differentially expressed microRNAs were found associated with the development of chasmogamous and cleistogamous flowers in Viola prionantha, revealing potential roles of microRNAs in the developmental evolution of dimorphic flowers. In Viola prionantha, chasmogamous (CH) flowers are induced by short daylight, while cleistogamous (CL) flowers are triggered by long daylight. How environmental factors and microRNAs (miRNAs) affect dimorphic flower formation remains unknown. In this study, small RNA sequencing was performed on CH and CL floral buds at different developmental stages in V. prionantha, differentially expressed miRNAs (DEmiRNAs) were identified, and their target genes were predicted. In CL flowers, Viola prionantha miR393 (vpr-miR393a/b) and vpr-miRN3366 were highly expressed, while in CH flowers, vpr-miRN2005, vpr-miR172e-2, vpr-miR166m-3, vpr-miR396f-2, and vpr-miR482d-2 were highly expressed. In the auxin-activated signaling pathway, vpr-miR393a/b and vpr-miRN2005 could target Vpr-TIR1/AFB and Vpr-ARF2, respectively, and other DEmiRNAs could target genes involved in the regulation of transcription, e.g., Vpr-AP2-7. Moreover, Vpr-UFO and Vpr-YAB5, the main regulators in petal and stamen development, were co-expressed with Vpr-TIR1/AFB and Vpr-ARF2 and showed lower expression in CL flowers than in CH flowers. Some V. prionantha genes relating to the stress/defense responses were co-expressed with Vpr-TIR1/AFB, Vpr-ARF2, and Vpr-AP2-7 and highly expressed in CL flowers. Therefore, in V. prionantha, CH-CL flower development may be regulated by the identified DEmiRNAs and their target genes, thus providing the first insight into the formation of dimorphic flowers in Viola.

Screening and validation of optimal miRNA reference genes in different developing stages and tissues of Lilium henryi Baker

Dynamic miRNA detection using the qRT-PCR technique requires appropriate reference genes to ensure data reliability. Previous studies have screened internal reference genes in plants during embryonic development and various stress treatment, involving relatively few tissues and organs. There is no relevant miRNA study in Lilium henryi Baker and limited research on the optimal miRNA reference genes in lilies, such as 5S, 18S, U6 and Actin. Twelve genes were selected as candidate reference genes whose expression stability was analyzed in petals at different developmental stages and other tissues using various algorithms, such as geNorm, NormFinder, BestKeeper, and Delta CT. The results revealed that the optimal combination of reference genes for Lilium henryi Baker petals at different developmental stages was osa-miR166m and osa-miR166a-3p, while that for different tissues of Lilium henryi Baker was osa-miR166g-3p and osa-miR166a-3p.Four important genes related to growth and development regulation, namely, osa-miR156a, osa-miR395b, osa-miR396a-3p, and osa-miR396a-5p, were selected for validation. The findings of the present study could contribute to future investigations onmiRNA expression and the related functions in Lilium henryi Baker while providing important references for the normalization of the miRNA expression in other varieties of lily.

PeGRF6-PeGIF1 complex regulates cell proliferation in the leaf of Phalaenopsis equestris

Phalaenopsis equestris is an ornamental plant with very large leaves. In this study, we identified genes related to the regulation of leaf development in Phalaenopsis and explored their mechanism of action. Sequence alignment and phylogenetic analyses revealed that PeGRF6 in the PeGRF family of P. equestris has similarities with the Arabidopsis genes AtGRF1 and AtGRF2, which are known to be involved in the regulation of leaf development. Among the PeGRFs, PeGRF6 was continuously and stably expressed at various stages of leaf development. The functions of PeGRF6 and of its complex formed with PeGIF1 in leaf development were verified by virus-induced gene silencing (VIGS) technology. The results show that the PeGRF6-PeGIF1 complex forms in the nucleus and positively regulates leaf cell proliferation via influencing cell size. Interestingly, VIGS suppression of PeGRF6 resulted in anthocyanin accumulation in Phalaenopsis leaves. Analyses of the regulatory mechanism of the miR396-PeGRF6 model based on the P. equestris small RNA library constructed here suggested that PeGRF6 transcripts are cleaved by Peq-miR396. These results show that, compared with PeGRF6 or PeGIF1 alone, the PeGRF6-PeGIF1 complex plays a more important role in the leaf development of Phalaenopsis, possibly by regulating the expression of cell cycle-related genes.

Identification and analysis of miRNAs differentially expressed in male and female Trichosanthes kirilowii maxim

Trichosanthes kirilowii Maxim. (TK) is a dioecious plant in the Cucurbitaceae family of which different sexes have separate medicinal uses. We used Illumina high-throughput sequencing technology to sequence miRNAs from male and female flower buds of TK. We performed bioinformatics analysis, miRNA identification, and target gene prediction on the data obtained from sequencing, and association analysis was performed in combination with the results of a previous transcriptome sequencing study. As a result, there were 80 differentially expressed miRNAs (DESs) between the female and male plants (48 upregulated and 32 downregulated in female plants). Moreover, 27 novel miRNAs in DESs were predicted to have 282 target genes, and 51 known miRNAs were predicted to have 3418 target genes. By establishing a regulatory network between miRNAs and target genes, 12 core genes were screened, including 7 miRNAs and 5 target genes. Among them, tkmiR157a-5p, tkmiR156c, tkmiR156_2, and tkmiR156k_2 jointly target the regulation of tkSPL18 and tkSPL13B. These two target genes are specifically expressed in male and female plants, respectively, and are involved in the biosynthesis process of BR, which is closely related to the sex differentiation process of TK. The identification of these miRNAs will provide a reference for the analysis of the sex differentiation mechanism of TK.

The Integrated mRNA and miRNA Approach Reveals Potential Regulators of Flowering Time in Arundina graminifolia

Orchids are among the most precious flowers in the world. Regulation of flowering time is one of the most important targets to enhance their ornamental value. The beauty of Arundina graminifolia is its year-round flowering, although the molecular mechanism of this flowering ability remains masked. Therefore, we performed a comprehensive assessment to integrate transcriptome and miRNA sequencing to disentangle the genetic regulation of flowering in this valuable species. Clustering analyses provided a set of molecular regulators of floral transition and floral morphogenesis. We mined candidate floral homeotic genes, including FCA, FPA, GI, FT, FLC, AP2, SOC1, SVP, GI, TCP, and CO, which were targeted by a variety of miRNAs. MiR11091 targeted the highest number of genes, including candidate regulators of phase transition and hormonal control. The conserved miR156-miR172 pathway of floral time regulation was evident in our data, and we found important targets of these miRNAs in the transcriptome. Moreover, endogenous hormone levels were determined to decipher the hormonal control of floral buds in A. graminifolia. The qRT-PCR analysis of floral and hormonal integrators validated the transcriptome expression. Therefore, miRNA-mediated mining of candidate genes with hormonal regulation forms the basis for comprehending the complex regulatory network of perpetual flowering in precious orchids. The findings of this study can do a great deal to broaden the breeding programs for flowering time manipulation of orchids.

Transcriptome mining of hormonal and floral integrators in the leafless flowers of three cymbidium orchids

Flowering is the most studied ornamental trait in orchids where long vegetative phase may span up to three years. Cymbidium orchids produce beautiful flowers with astonishing shapes and pleasant scent. However, an unusually long vegetative phase is a major drawback to their ornamental value. We observed that under certain culture conditions, three cymbidium species (Cymbidium ensifolium, C. goeringii and C. sinense) skipped vegetative growth phase and directly flowered within six months, that could be a breakthrough for future orchids with limited vegetative growth. Hormonal and floral regulators could be the key factors arresting vegetative phase. Therefore, transcriptomic analyses were performed for leafless flowers and normal vegetative leaves to ascertain differentially expressed genes (DEGs) related to hormones (auxin, cytokinin, gibberellin, abscisic acid and ethylene), floral integrators and MADS-box genes. A significant difference of cytokinin and floral regulators was observed among three species as compared to other hormones. The MADS-box genes were significantly expressed in the leafless flowers of C. sinense as compared to other species. Among the key floral regulators, CONSTANS and AGAMOUS-like genes showed the most differential expression in the leafless flowers as compared to leaves where the expression was negligible. However, CONSTANS also showed downregulation. Auxin efflux carriers were mainly downregulated in the leafless flowers of C. ensifolium and C. sinense, while they were upregulated in C. goeringii. Moreover, gibberellin and cytokinin genes were also downregulated in C. ensifolium and C. sinense flowers, while they were upregulated in C. goeringii, suggesting that species may vary in their responses. The data mining thus, outsources the valuable information to direct future research on orchids at industrial levels.

Identification of long non-coding RNAs and microRNAs involved in anther development in the tropical Camellia oleifera

BACKGROUND

Explored the molecular science of anther development is important for improving productivity and overall yield of crops. Although the role of regulatory RNAs, including long non-coding RNAs (lncRNAs) and microRNAs (miRNAs), in regulating anther development has been established, their identities and functions in Camellia oleifera, an important industrial crop, have yet not been clearly explored. Here, we report the identification and characterization of genes, lncRNAs and miRNAs during three stages of the tropical C. oleifera anther development by single-molecule real-time sequencing, RNA sequencing and small RNA sequencing, respectively.

RESULTS

These stages, viz. the pollen mother cells stage, tetrad stage and uninucleate pollen stage, were identified by analyzing paraffin sections of floral buds during rapid expansion periods. A total of 18,393 transcripts, 414 putative lncRNAs and 372 miRNAs were identified, of which 5,324 genes, 115 lncRNAs, and 44 miRNAs were differentially accumulated across three developmental stages. Of these, 44 and 92 genes were predicted be regulated by 37 and 30 differentially accumulated lncRNAs and miRNAs, respectively. Additionally, 42 differentially accumulated lncRNAs were predicted as targets of 27 miRNAs. Gene ontology enrichment indicated that potential target genes of lncRNAs were enriched in photosystem II, regulation of autophagy and carbohydrate phosphatase activity, which are essential for anther development. Functional annotation of genes targeted by miRNAs indicated that they are relevant to transcription and metabolic processes that play important roles in microspore development. An interaction network was built with 2 lncRNAs, 6 miRNAs and 10 mRNAs. Among these, miR396 and miR156 family were up-regulated, while their targets, genes (GROWTH REGULATING FACTORS and SQUAMOSA PROMOTER BINDING PROTEIN-LIKE genes) and lncRNAs, were down-regulated. Further, the trans-regulated targets of these lncRNAs, like wall-associated kinase2 and phosphomannose isomerase1, are involved in pollen wall formation during anther development.

CONCLUSIONS

This study unravels lncRNAs, miRNAs and miRNA-lncRNA-mRNA networks involved in development of anthers of the tropical C. oleifera lays a theoretical foundation for further elucidation of regulatory roles of lncRNAs and miRNAs in anther development.

Transcriptional Proposition for Uniquely Developed Protocorm Flowering in Three Orchid Species: Resources for Innovative Breeding

During orchid seed culture, seeds germinate as protocorms, and protocorms normally develop into plant with leaves and roots. Orchids require many years of vegetative development for flowering. However, under a certain combination of growth cultures, we observed that protocorms can directly flower without leaves and roots. Therefore, we performed comparative transcriptome analysis to identify the different transcriptional regulators of two types of protocorms of Cymbidium ensifolium, Cymbidium sinense, and Cymbidium goeringii. Zinc finger, MYB, AP2, and bHLH were the most abundant transcription factor (TF) families in the transcriptome. Weighted gene coexpression network analysis (WGCNA) was performed to identify hub genes related to leaf and flower development. The key hubs included SPL6, SVP, SEP2, KNOX1, AP2, OFP1, COL12, MYB13, MYB36, MYB59, bHLH086, and ARF7. The hub genes were further validated through statistical tools to propose the roles of key TFs. Therefore, this study initiates to answer that why there is no leaf initiation and root development and how can protocorm bypass the vegetative phase to flower? The outcomes can direct future research on short-span flowering in orchids through protocorms.