Analysis of microRNA reveals cleistogamous and chasmogamous floret divergence in dimorphic plant

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.

The roles of non-coding RNAs in male reproductive development and abiotic stress responses during this unique process in flowering plants

Successful male reproductive development is the guarantee for sexual reproduction of flowering plants. Male reproductive development is a complicated and multi-stage process that integrates physiological processes and adaptation and tolerance to a myriad of environmental stresses. This well-coordinated process is governed by genetic and epigenetic machineries. Non-coding RNAs (ncRNAs) play pleiotropic roles in the plant growth and development. The identification, characterization and functional analysis of ncRNAs and their target genes have opened a new avenue for comprehensively revealing the regulatory network of male reproductive development and its response to environmental stresses in plants. This review briefly addresses the types, origin, biogenesis and mechanisms of ncRNAs in plants, highlights important updates on the roles of ncRNAs in regulating male reproductive development and emphasizes the contribution of ncRNAs, especially miRNAs and lncRNAs, in responses to abiotic stresses during this unique process in flowering plants.

Identification of candidate genes responsible for chasmogamy in wheat

BACKGROUND

The flowering biology of wheat plants favours self-pollination which causes obstacles in wheat hybrid breeding. Wheat flowers can be divided into two groups, the first one is characterized by flowering and pollination within closed flowers (cleistogamy), while the second one possesses the ability to open flowers during processes mentioned above (chasmogamy). The swelling of lodicules is involved in the flowering of cereals and among others their morphology, calcium and potassium content differentiate between cleistogamic and non-cleistogamous flowers. A better understanding of the chasmogamy mechanism can lead to the development of tools for selection of plants with the desired outcrossing rate. To learn more, the sequencing of transcriptomes (RNA-Seq) and Representational Difference Analysis products (RDA-Seq) were performed to investigate the global transcriptomes of wheat lodicules in two highly chasmogamous (HCH, Piko and Poezja) and two low chasmogamous (LCH, Euforia and KWS Dacanto) varieties at two developmental stages-pre-flowering and early flowering.

RESULTS

The differentially expressed genes were enriched in five, main pathways: "metabolism", "organismal systems", "genetic information processing", "cellular processes" and "environmental information processing", respectively. Important genes with opposite patterns of regulation between the HCH and LCH lines have been associated with the lodicule development i.e. expression levels of MADS16 and MADS58 genes may be responsible for quantitative differences in chasmogamy level in wheat.

CONCLUSIONS

We conclude that the results provide a new insight into lodicules involvement in the wheat flowering process. This study generated important genomic information to support the exploitation of the chasmogamy in wheat hybrid breeding programs.

Transcriptomic heterochrony and completely cleistogamous flower development in the mycoheterotrophic orchid Gastrodia

Cleistogamy, in which plants can reproduce via self-fertilization within permanently closed flowers, has evolved in > 30 angiosperm lineages; however, consistent with Darwin's doubts about its existence, complete cleistogamy - the production of only cleistogamous flowers - has rarely been recognized. Thus far, the achlorophyllous orchid genus, Gastrodia, is the only known genus with several plausible completely cleistogamous species. Here, we analyzed the floral developmental transcriptomes of two recently evolved, completely cleistogamous Gastrodia species and their chasmogamous sister species to elucidate the possible changes involved in producing common cleistogamous traits. The ABBA-BABA test did not support introgression and protein sequence convergence as evolutionary mechanisms leading to cleistogamy, leaving convergence in gene expression as a plausible mechanism. Regarding transcriptomic differentiation, the two cleistogamous species had common modifications in the expression of developmental regulators, exhibiting a gene family-wide signature of convergent expression changes in MADS-box genes. Our transcriptomic pseudotime analysis revealed a prolonged juvenile state and eventual maturation, a heterochronic pattern consistent with partial neoteny, in cleistogamous flower development. These findings indicate that transcriptomic partial neoteny, arising from changes in the expression of conserved developmental regulators, might have contributed to the rapid and repeated evolution of cleistogamous flowers in Gastrodia.

The emerging role of small RNAs in ovule development, a kind of magic

In plants, small RNAs have been recognized as key genetic and epigenetic regulators of development. Small RNAs are usually 20 to 30 nucleotides in length and they control, in a sequence specific manner, the transcriptional or post-transcriptional expression of genes. In this review, we present a comprehensive overview of the most recent findings about the function of small RNAs in ovule development, including megasporogenesis and megagametogenesis, both in sexual and apomictic plants. We discuss recent studies on the role of miRNAs, siRNAs and trans-acting RNAs (ta-siRNAs) in early female germline differentiation. The mechanistic complexity and unique regulatory features are reviewed, and possible directions for future research are provided.

The genome of Cleistogenes songorica provides a blueprint for functional dissection of dimorphic flower differentiation and drought adaptability

Cleistogenes songorica (2n = 4x = 40) is a desert grass with a unique dimorphic flowering mechanism and an ability to survive extreme drought. Little is known about the genetics underlying drought tolerance and its reproductive adaptability. Here, we sequenced and assembled a high-quality chromosome-level C. songorica genome (contig N50 = 21.28 Mb). Complete assemblies of all telomeres, and of ten chromosomes were derived. C. songorica underwent a recent tetraploidization (~19 million years ago) and four major chromosomal rearrangements. Expanded genes were significantly enriched in fatty acid elongation, phenylpropanoid biosynthesis, starch and sucrose metabolism, and circadian rhythm pathways. By comparative transcriptomic analysis we found that conserved drought tolerance related genes were expanded. Transcription of CsMYB genes was associated with differential development of chasmogamous and cleistogamous flowers, as well as drought tolerance. Furthermore, we found that regulation modules encompassing miRNA, transcription factors and target genes are involved in dimorphic flower development, validated by overexpression of CsAP2_9 and its targeted miR172 in rice. Our findings enable further understanding of the mechanisms of drought tolerance and flowering in C. songorica, and provide new insights into the adaptability of native grass species in evolution, along with potential resources for trait improvement in agronomically important species.

Global identification and characterization of miRNA family members responsive to potassium deprivation in wheat (Triticum aestivum L.)

Potassium (K) is essential for plant growth and stress responses. MicroRNAs (miRNAs) are involved in adaptation to nutrient deprivation through modulating gene expression. Here, we identified the miRNAs responsive to K deficiency in Triticum aestivum based on high-throughput small RNA sequencing analyses. Eighty-nine miRNAs, including 68 previously reported ones and 21 novel ones, displayed differential expression under K deficiency. In Gene Ontology and Kyoto Encyclopedia and Genome analyses, the putative target genes of the differentially expressed miRNAs were categorized into functional groups associated with ADP-binding activity, secondary metabolic pathways, and biosynthesis and metabolism. Functional characterization of tae-miR408, an miRNA significantly down-regulated under K deficiency, revealed its important role in mediating low-K tolerance. Compared with wild type, transgenic tobacco lines overexpressing tae-miR408 showed significantly improved K uptake, biomass, photosynthesis, and reactive oxygen species scavenging under K deficiency. These results show that distinct miRNAs function in the plant response to K deficiency through regulating target genes involved in energy metabolism and various secondary metabolic pathways. Our findings shed light on the plant response to K deficiency mediated by miRNAs in T. aestivum. Distinct miRNAs, such as tae-miR408, are valuable targets for generating crop varieties with improved K-use efficiency.

Genome-Wide Analysis of the Role of NAC Family in Flower Development and Abiotic Stress Responses in Cleistogenes songorica

Plant-specific NAC (NAM, ATAF, CUC) transcription factor (TF) family plays important roles in biological processes such as plant growth and response to stress. Nevertheless, no information is known about NAC TFs in Cleistogenes songorica, a prominent xerophyte desert grass in northwestern China. In this study, 162 NAC genes were found from the Cleistogenes songorica genome, among which 156 C. songoricaNAC (CsNAC) genes (96.3%) were mapped onto 20 chromosomes. The phylogenetic tree constructed by CsNAC and rice NAC TFs can be separated into 14 subfamilies. Syntenic and Ka/Ks analyses showed that CsNACs were primarily expanded by genomewide replication events, and purifying selection was the primary force driving the evolution of CsNAC family genes. The CsNAC gene expression profiles showed that 36 CsNAC genes showed differential expression between cleistogamous (CL) and chasmogamous (CH) flowers. One hundred and two CsNAC genes showed differential expression under heat, cold, drought, salt and ABA treatment. Twenty-three CsNAC genes were commonly differentially expressed both under stress responses and during dimorphic floret development. Gene Ontology (GO) annotation, coexpression network and qRT-PCR tests revealed that these CsNAC genes may simultaneously regulate dimorphic floret development and the response to stress. Our results may help to characterize the NAC transcription factors in C. songorica and provide new insights into the functional research and application of the NAC family in crop improvement, especially in dimorphic floret plants.

Comprehensive analysis of bZIP transcription factors uncovers their roles during dimorphic floret differentiation and stress response in Cleistogenes songorica

BACKGROUND

Transcription factors act as important regulators of transcription networks. Basic leucine zipper (bZIP) transcription factors have been shown to be involved in multiple biological processes in plants. However, no information is available for the bZIP family in Cleistogenes songorica, which is an important xerophytic and allotetraploid grass in desert grasslands.

RESULTS

In this study, 86 CsbZIPs were identified in the allotetraploid C. songorica genome. For location analysis, CsbZIPs were distributed evenly across two subgenomes of C. songorica. Phylogenetic tree analysis among three species indicated that CsbZIPs were evolutionarily more closely related to OsbZIPs than AtbZIPs. Syntenic and phylogenetic analyses confirmed that the CsbZIPs were mainly expanded by whole-genome duplication events. Furthermore, it was determined that rice and C. songorica might have undergone purified selection during their long evolutionary history by calculating the Ks values and Ka/Ks ratios of orthologous gene pairs. By analysing the expression patterns of CsbZIPs in different tissues and under abiotic stresses, 21 CsbZIP genes were differentially expressed between chasmogamous (CH) and cleistogamous (CL) flowers, including two FLOWERING LOCUS D (FD) genes. In shoots and roots, 79.1 and 87.2% of the CsbZIP genes, respectively, displayed transcript changes under at least one stress treatment, such as heat, cold, drought and salt. Strikingly, 17 common CsbZIP genes showed differential expression under stress response and during CL flowering. Co-expression network, GO annotation and real-time quantitative reverse transcription PCR (qRT-PCR) analyses revealed a close relationship between CL flowering-associated genes and abiotic stress-related genes.

CONCLUSIONS

BZIP TFs were comprehensively analysed and identified in allotetraploid C. songorica. Our results provide insights into the evolutionary history of the bZIP family in C. songorica and provide abiotic stress-responsive and CL-associated candidate CsbZIP genes for potential applications in the genetic improvement of plants.

Differential co-expression networks of long non-coding RNAs and mRNAs in Cleistogenes songorica under water stress and during recovery

BACKGROUND

Water stress seriously constrains plant growth and yield. Long non-coding RNAs (lncRNAs) serve as versatile regulators in various biological regulatory processes. To date, the systematic screening and potential functions of lncRNA have not yet been characterized in Cleistogenes songorica, especially under water stress conditions.

RESULTS

In this study, we obtained the root and shoot transcriptomes of young C. songorica plants subjected to different degrees of water stress and recovery treatments by Illumina-based RNA-seq. A total of 3397 lncRNAs were identified through bioinformatics analysis. LncRNA differential expression analysis indicated that the higher response of roots compared to shoots during water stress and recovery. We further identified the 1644 transcription factors, 189 of which were corresponded to 163 lncRNAs in C. songorica. Though comparative analyses with major Poaceae species based on blast, 81 water stress-related orthologues regulated to lncRNAs were identified as a core of evolutionary conserved genes important to regulate water stress responses in the family. Among these target genes, two genes were found to be involved in the abscisic acid (ABA) signalling pathway, and four genes were enriched for starch and sucrose metabolism. Additionally, the 52 lncRNAs were predicted as target mimics for microRNAs (miRNAs) in C. songorica. RT-qPCR results suggested that MSTRG.43964.1 and MSTRG.4400.2 may regulate the expression of miRNA397 and miRNA166, respectively, as target mimics under water stress and during recovery. Finally, a co-expression network was constructed based on the lncRNAs, miRNAs, protein-coding genes (PCgenes) and transcription factors under water stress and during recovery in C. songorica.

CONCLUSIONS

In C. songorica, lncRNAs, miRNAs, PCgenes and transcription factors constitute a complex transcriptional regulatory network which lncRNAs can regulate PCgenes and miRNAs under water stress and recovery. This study provides fundamental resources to deepen our knowledge on lncRNAs during ubiquitous water stress.