Conserved and non-conserved triggers of 24-nucleotide reproductive phasiRNAs in eudicots

The floral development of the allotetraploid Coffea arabica L. correlates with a small RNA dynamic reprogramming

Noncoding and coding RNAs are key regulators of plant growth, development, and stress responses. To investigate the types of transcripts accumulated during the vegetative to reproductive transition and floral development in the Coffea arabica L., we sequenced small RNA libraries from eight developmental stages, up to anthesis. We combined these data with messenger RNA and PARE sequencing of two important development stages that marks the transition of an apparent latent to a rapid growth stage. In addition, we took advantage of multiple in silico tools to characterize genomic loci producing small RNAs such as phasiRNAs, miRNAs, and tRFs. Our differential and co-expression analysis showed that some types of small RNAs such as tRNAs, snoRNAs, snRNAs, and phasiRNAs preferentially accumulate in a stage-specific manner. Members of the miR482/miR2118 superfamily and their 21-nucleotide phasiRNAs originating from resistance genes show a robust co-expression pattern that is maintained across all the evaluated developmental stages. Finally, the majority of miRNAs accumulate in a family stage-specific manner, related to modulated hormonal responses and transcription factor expression.

Phased secondary small interfering RNAs in Camellia sinensis var. assamica

Phased secondary small interfering RNAs (phasiRNAs) in plants play important roles in regulating genome stability, plant development and stress adaption. Camellia sinensis var. assamica has immense economic, medicinal and cultural significance. However, there are still no studies of phasiRNAs and their putative functions in this valuable plant. We identified 476 and 43 PHAS loci which generated 4290 twenty one nucleotide (nt) and 264 twenty four nt phasiRNAs, respectively. Moreover, the analysis of degradome revealed more than 35000 potential targets for these phasiRNAs. We identified several conserved 21 nt phasiRNA generation pathways in tea plant, including miR390 → TAS3, miR482/miR2118 → NB-LRR, miR393 → F-box, miR828 → MYB/TAS4, and miR7122 → PPR in this study. Furthermore, we found that some transposase and plant mobile domain genes could generate phasiRNAs. Our results show that phasiRNAs target genes in the same family in cis- or trans-manners, and different members of the same gene family may generate the same phasiRNAs. The phasiRNAs, generated by transposase and plant mobile domain genes, and their targets, suggest that phasiRNAs may be involved in the inhibition of transposable elements in tea plant. To summarize, these results provide a comprehensive view of phasiRNAs in Camellia sinensis var. assamica.

Non-coding RNAs-mediated environmental surveillance determines male fertility in plants

Plants are continuously exposed to environmental stresses leading to significant yield losses. With the changing climatic conditions, the intensity and duration of these stresses are expected to increase, posing a severe threat to crop productivity worldwide. Male gametogenesis is one of the most sensitive developmental stages. Exposure to environmental stresses during this stage leads to male sterility and yield loss. Elucidating the underlying molecular mechanism of environment-affected male sterility is essential to address this challenge. High-throughput RNA sequencing studies, loss-of-function phenotypes of sRNA biogenesis genes and functional genomics studies with non-coding RNAs have started to unveil the roles of small RNAs, long non-coding RNAs and the complex regulatory interactions between them in regulating male fertility under different growth regimes. Here, we discuss the current understanding of the non-coding RNA-mediated environmental stress surveillance and regulation of male fertility in plants. The candidate ncRNAs emerging from these studies can be leveraged to generate environment-sensitive male sterile lines for hybrid breeding or mitigate the impact of climate change on male fertility, as the situation demands.

CRISPR/Cas9-targeted mutagenesis of TaDCL4, TaDCL5 and TaRDR6 induces male sterility in common wheat

Phased, small interfering RNAs (phasiRNAs) are important for plant anther development, especially for male sterility. PhasiRNA biogenesis is dependent on genes like RNA polymerase 6 (RDR6), DICER-LIKE 4 (DCL4), or DCL5 to produce 21- or 24 nucleotide (nt) double-strand small RNAs. Here, we generated mutants of DCL4, DCL5 and RDR6 using CRISPR/Cas9 system and studied their effects on plant reproductive development and phasiRNA production in wheat. We found that RDR6 mutation caused sever consequence throughout plant development starting from seed germination and the dcl4 mutants grew weaker with thorough male sterility, while dcl5 plants developed normally but exhibited male sterility. Correspondingly, DCL4 and DCL5, respectively, specified 21- and 24-nt phasiRNA biogenesis, while RDR6 contributed to both. Also, the three key genes evolved differently in wheat, with TaDCL5-A/B becoming non-functioning and TaRDR6-A being lost after polyploidization. Furthermore, we found that PHAS genes (phasiRNA precursors) identified via phasiRNAs diverged rapidly among sub-genomes of polyploid wheat. Despite no similarity being found among phasiRNAs of grasses, their targets were enriched for similar biological functions. In light of the important roles of phasiRNA pathways in gametophyte development, genetic dissection of the function of key genes may help generate male sterile lines suitable for hybrid wheat breeding.

Genetic architecture underlying variation in floral meristem termination in Aquilegia

Floral organs are produced by floral meristems (FMs), which harbor stem cells in their centers. Since each flower only has a finite number of organs, the stem cell activity of an FM will always terminate at a specific time point, a process termed floral meristem termination (FMT). Variation in the timing of FMT can give rise to floral morphological diversity, but how this process is fine-tuned at a developmental and evolutionary level is poorly understood. Flowers from the genus Aquilegia share identical floral organ arrangement except for stamen whorl number (SWN), making Aquilegia a well-suited system for investigation of this process: differences in SWN between species represent differences in the timing of FMT. By crossing A. canadensis and A. brevistyla, quantitative trait locus (QTL) mapping has revealed a complex genetic architecture with seven QTL. We explored potential candidate genes under each QTL and characterized novel expression patterns of select loci of interest using in situ hybridization. To our knowledge, this is the first attempt to dissect the genetic basis of how natural variation in the timing of FMT is regulated, and our results provide insight into how floral morphological diversity can be generated at the meristematic level.

24-nt phasiRNAs move from tapetal to meiotic cells in maize anthers

In maize, 24-nt phased, secondary small interfering RNAs (phasiRNAs) are abundant in meiotic stage anthers, but their distribution and functions are not precisely known. Using laser capture microdissection, we analyzed tapetal cells, meiocytes and other somatic cells at several stages of anther development to establish the timing of 24-PHAS precursor transcripts and the 24-nt phasiRNA products. By integrating RNA and small RNA profiling plus single-molecule and small RNA FISH (smFISH or sRNA-FISH) spatial detection, we demonstrate that the tapetum is the primary site of 24-PHAS precursor and Dcl5 transcripts and the resulting 24-nt phasiRNAs. Interestingly, 24-nt phasiRNAs accumulate in all cell types, with the highest levels in meiocytes, followed by tapetum. Our data support the conclusion that 24-nt phasiRNAs are mobile from tapetum to meiocytes and to other somatic cells. We discuss possible roles for 24-nt phasiRNAs in anther cell types.

Dynamically expressed small RNAs, substantially driven by genomic structural variants, contribute to transcriptomic changes during tomato domestication

Tomato has undergone extensive selections during domestication. Recent progress has shown that genomic structural variants (SVs) have contributed to gene expression dynamics during tomato domestication, resulting in changes of important traits. Here, we performed comprehensive analyses of small RNAs (sRNAs) from nine representative tomato accessions. We demonstrate that SVs substantially contribute to the dynamic expression of the three major classes of plant sRNAs: microRNAs (miRNAs), phased secondary short interfering RNAs (phasiRNAs), and 24-nucleotide heterochromatic siRNAs (hc-siRNAs). Changes in the abundance of phasiRNAs and 24-nucleotide hc-siRNAs likely contribute to the alteration of mRNA gene expression in cis during tomato domestication, particularly for genes associated with biotic and abiotic stress tolerance. We also observe that miRNA expression dynamics are associated with imprecise processing, alternative miRNA-miRNA* selections, and SVs. SVs mainly affect the expression of less-conserved miRNAs that do not have established regulatory functions or low abundant members in highly expressed miRNA families. Our data highlight different selection pressures on miRNAs compared to phasiRNAs and 24-nucleotide hc-siRNAs. Our findings provide insights into plant sRNA evolution as well as SV-based gene regulation during crop domestication. Furthermore, our dataset provides a rich resource for mining the sRNA regulatory network in tomato.

Heat-responsive microRNAs and phased small interfering RNAs in reproductive development of flax

Plants will face increased heat stress due to rising global temperatures. Heat stress affects plant reproductive development and decreases productivity; however, the underlying molecular mechanisms of these processes are poorly characterized. Plant small RNAs (sRNAs) have important regulatory roles in plant reproductive development following abiotic stress responses. We generated sRNA transcriptomes of reproductive bud stages at three different time points to identify sRNA-mediated pathways responsive to heat stress in flax (Linum usitatissimum). With added sRNA transcriptomes of vegetative tissues, we comprehensively annotated miRNA and phasiRNA-encoding genes (PHAS) in flax. We identified 173 miRNA genes, of which 42 are newly annotated. Our analysis revealed that 141 miRNA genes were differentially accumulated between tissue types, while 18 miRNA genes were differentially accumulated in reproductive tissues following heat stress, including members of miR482/2118 and miR2275 families, known triggers of reproductive phasiRNAs. Furthermore, we identified 68 21-PHAS flax loci from protein-coding and noncoding regions, four 24-PHAS loci triggered by miR2275, and 658 24-PHAS-like loci with unknown triggers, derived mostly from noncoding regions. The reproductive phasiRNAs are mostly downregulated in response to heat stress. Overall, we found that several previously unreported miRNAs and phasiRNAs are responsive to heat stress in flax reproductive tissues.

Spatiotemporal regulation and roles of reproductive phasiRNAs in plants

Since co-suppression was discovered as a pioneer silencing phenomenon of RNA interference (RNAi) in petunia in 1990, many types of small RNAs have been identified in the RNAi pathway among various eukaryotes. In plants, a large number of 21- or 24-nucleotide (nt) phased small interfering RNAs (phasiRNAs) are produced via processing of long RNA precursors by Dicer-like proteins. However, the roles of phasiRNAs remain largely unknown. The development of imaging technology and RNA profiling has clarified the spatiotemporal regulation of phasiRNAs, and subsequently the different functions of 21-nt trans-acting phasiRNAs and 24-nt cis-regulatory phasiRNAs during male organ development. This review focuses on the biogenesis, diversification, spatiotemporal expression pattern and function of phasiRNAs in plants.