Spatiotemporally dynamic, cell-type-dependent premeiotic and meiotic phasiRNAs in maize anthers

Pseudouridine guides germline small RNA transport and epigenetic inheritance

Developmental epigenetic modifications in plants and animals are mostly reset during gamete formation but some are inherited from the germline. Small RNAs guide these epigenetic modifications but how inherited small RNAs are distinguished in plants and animals is unknown. Pseudouridine (Ψ) is the most abundant RNA modification but has not been explored in small RNAs. Here, we develop assays to detect Ψ in short RNA sequences, demonstrating its presence in mouse and Arabidopsis microRNAs. Germline small RNAs, namely epigenetically activated small interfering RNAs (easiRNAs) in Arabidopsis pollen and Piwi-interacting RNAs in mouse testes, are enriched for Ψ. In pollen, pseudouridylated easiRNAs are transported to sperm cells from the vegetative nucleus, and PAUSED/HEN5 (PSD), the plant homolog of Exportin-t, interacts genetically with Ψ and is required for this transport. We further show that Exportin-t is required for the triploid block: small RNA dosage-dependent seed lethality that is epigenetically inherited from pollen. Thus, Ψ has a conserved role in marking inherited small RNAs in the germline.

Epigenetic gene regulation in plants and its potential applications in crop improvement

DNA methylation, also known as 5-methylcytosine, is an epigenetic modification that has crucial functions in plant growth, development and adaptation. The cellular DNA methylation level is tightly regulated by the combined action of DNA methyltransferases and demethylases. Protein complexes involved in the targeting and interpretation of DNA methylation have been identified, revealing intriguing roles of methyl-DNA binding proteins and molecular chaperones. Structural studies and in vitro reconstituted enzymatic systems have provided mechanistic insights into RNA-directed DNA methylation, the main pathway catalysing de novo methylation in plants. A better understanding of the regulatory mechanisms will enable locus-specific manipulation of the DNA methylation status. CRISPR-dCas9-based epigenome editing tools are being developed for this goal. Given that DNA methylation patterns can be stably transmitted through meiosis, and that large phenotypic variations can be contributed by epimutations, epigenome editing holds great promise in crop breeding by creating additional phenotypic variability on the same genetic material.

The role of microRNAs in NBS-LRR gene expression and its implications for plant immunity and crop development

Plants evolved, over millions of years, complex defense systems against pathogens. Once infected, the interaction between pathogen effector molecules and host receptors triggers plant immune responses, which include apoptosis, systemic immune response, among others. An important protein family responsible for pathogen effector recognition is the nucleotide binding site-leucine repeat rich (NBS-LRR) proteins. The NBS-LRR gene family is the largest disease resistance gene class in plants. These proteins are widely distributed in vascular plants and have a complex multigenic cluster distribution in plant genomes. To counteract the genetic load of such a large gene family on fitness cost, plants evolved a mechanism using post transcriptional gene silencing induced by small RNAs, particularly microRNAs. For the NBS-LRR gene family, the small RNAs involved in this silencing mechanism are mainly the microRNA482/2118 superfamily. This suppression mechanism is relieved upon pathogen infection, thus allowing increased NBS-LRR expression and triggering plant immunity. In this review, we will discuss the biogenesis of microRNAs and secondary RNAs involved in this silencing mechanism, biochemical and structural features of NBS-LRR proteins in response to pathogen effectors and the evolution of microRNA-based silencing mechanism with a focus on the miR482/2118 family. Furthermore, the biotechnological manipulation of microRNA expression, using both transgenic or genome editing approaches to improve cultivated plants will be discussed, with a focus on the miR482/2118 family in soybean.

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.

Premeiotic 24-nt phasiRNAs are present in the Zea genus and unique in biogenesis mechanism and molecular function

Reproductive phasiRNAs (phased, small interfering RNAs) are broadly present in angiosperms and play crucial roles in sustaining male fertility. While the premeiotic 21-nt (nucleotides) phasiRNAs and meiotic 24-nt phasiRNA pathways have been extensively studied in maize (Zea mays) and rice (Oryza sativa), a third putative category of reproductive phasiRNAs-named premeiotic 24-nt phasiRNAs-have recently been reported in barley (Hordeum vulgare) and wheat (Triticum aestivum). To determine whether premeiotic 24-nt phasiRNAs are also present in maize and related species and begin to characterize their biogenesis and function, we performed a comparative transcriptome and degradome analysis of premeiotic and meiotic anthers from five maize inbred lines and three teosinte species/subspecies. Our data indicate that a substantial subset of the 24-nt phasiRNA loci in maize and teosinte are already highly expressed at the premeiotic phase. The premeiotic 24-nt phasiRNAs are similar to meiotic 24-nt phasiRNAs in genomic origin and dependence on DCL5 (Dicer-like 5) for biogenesis, however, premeiotic 24-nt phasiRNAs are unique in that they are likely i) not triggered by microRNAs, ii) not loaded by AGO18 proteins, and iii) not capable of mediating PHAS precursor cleavage. In addition, we also observed a group of premeiotic 24-nt phasiRNAs in rice using previously published data. Together, our results indicate that the premeiotic 24-nt phasiRNAs constitute a unique class of reproductive phasiRNAs and are present more broadly in the grass family (Poaceae) than previously known.

Premeiotic 24-nt phasiRNAs are present in the Zea genus and unique in biogenesis mechanism and molecular function

Reproductive phasiRNAs are broadly present in angiosperms and play crucial roles in sustaining male fertility. While the premeiotic 21-nt phasiRNAs and meiotic 24-nt phasiRNA pathways have been extensively studied in maize (Zea mays) and rice (Oryza sativa), a third putative category of reproductive phasiRNAs-named premeiotic 24-nt phasiRNAs-have recently been reported in barley (Hordeum vulgare) and wheat (Triticum aestivum). To determine whether premeiotic 24-nt phasiRNAs are also present in maize and related species and begin to characterize their biogenesis and function, we performed a comparative transcriptome and degradome analysis of premeiotic and meiotic anthers from five maize inbred lines and three teosinte species/subspecies. Our data indicate that a substantial subset of the 24-nt phasiRNA loci in maize and teosinte are already highly expressed at premeiotic phase. The premeiotic 24-nt phasiRNAs are similar to meiotic 24-nt phasiRNAs in genomic origin and dependence on DCL5 for biogenesis, however, premeiotic 24-nt phasiRNAs are unique in that they are likely (i) not triggered by microRNAs, (ii) not loaded by AGO18 proteins, and (iii) not capable of mediating cis-cleavage. In addition, we also observed a group of premeiotic 24-nt phasiRNAs in rice using previously published data. Together, our results indicate that the premeiotic 24-nt phasiRNAs constitute a unique class of reproductive phasiRNAs and are present more broadly in the grass family (Poaceae) than previously known.

Review: Plant microRNAs in pathogen defense: A panacea or a piece of the puzzle?

Plant microRNAs (miRNAs) control key agronomic traits that are associated with their conserved role(s) in development. However, despite a multitude of studies, the utility of miRNAs in plant-pathogen resistance remains less certain. Reviewing the literature identifies three general classes of miRNAs regarding plant pathogen defense. Firstly, a number of evolutionary dynamic 22 nucleotide miRNA families that repress large numbers of plant immunity genes, either directly, or through triggering the biogenesis of secondary siRNAs. However, understanding of their role in defense and of their manipulation to enhance pathogen resistance are still lacking. Secondly, highly conserved miRNAs that indirectly impact disease resistance through their targets that are primarily regulating development or hormone signaling. Any alteration of these miRNAs usually results in pleiotropic impacts, which may alter disease resistance in some plant species, and against some pathogens. Thirdly, are the comparatively diverse and evolutionary dynamic set of non-conserved miRNAs, some of which contribute to pathogen resistance, but whose narrow evolutionary presence will likely restrict their utility. Therefore, reflecting the diverse and evolving nature of plant-pathogen interactions, a complex interplay of plant miRNAs with pathogen responses exists. Any miRNA-based solution for pathogen resistance will likely be highly specific, rather than a general panacea.

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.

sRNAminer: A multifunctional toolkit for next-generation sequencing small RNA data mining in plants

Small RNAs (sRNAs), found extensively in plants, play an essential role in plant growth and development. Although various sRNA analysis tools have been developed for plants, the use of most of them depends on programming and command-line environments, which is a challenge for many wet-lab biologists. Furthermore, current sRNA analysis tools mostly focus on the analysis of certain type of sRNAs and are resource-intensive, normally demanding an immense amount of time and effort to learn the use of numerous tools or scripts and assemble them into a workable pipeline to get the final results. Here, we present sRNAminer, a powerful stand-alone toolkit with a user-friendly interface that integrates all common functions for the analysis of three major types of plant sRNAs: microRNAs (miRNAs), phased small interfering RNAs (phasiRNAs), and heterochromatic siRNAs (hc-siRNAs). We constructed a curated or "golden" set of MIRNA and PHAS loci, which was used to assess the performance of sRNAminer in comparison to other existing tools. The results showed that sRNAminer outperformed these tools in multiple aspects, highlighting its functionality. In addition, to enable an efficient evaluation of sRNA annotation results, we developed Integrative Genomics Viewer (IGV)-sRNA, a modified genome browser optimized from IGV and we incorporated it as a functional module in sRNAminer. IGV-sRNA can display a wealth of sRNA-specific features, enabling a more comprehensive understanding of sRNA data. sRNAminer and IGV-sRNA are both platform-independent software that can be run under all operating systems. They are now freely available at https://github.com/kli28/sRNAminer and https://gitee.com/CJchen/IGV-sRNA.

mRNA cleavage by 21-nucleotide phasiRNAs determines temperature-sensitive male sterility in rice

Temperature-sensitive male sterility is one of the core components for hybrid rice (Oryza sativa) breeding based on the 2-line system. We previously found that knockout of ARGONAUTE 1d (AGO1d) causes temperature-sensitive male sterility in rice by influencing phased small interfering RNA (phasiRNA) biogenesis and function. However, the specific phasiRNAs and their targets underlying the temperature-sensitive male sterility in the ago1d mutant remain unknown. Here, we demonstrate that the ago1d mutant displays normal female fertility but complete male sterility at low temperature. Through a multiomics analysis of small RNA (sRNA), degradome, and transcriptome, we found that 21-nt phasiRNAs account for the greatest proportion of the 21-nt sRNA species in rice anthers and are sensitive to low temperature and markedly downregulated in the ago1d mutant. Moreover, we found that 21-nt phasiRNAs are essential for the mRNA cleavage of a set of fertility- and cold tolerance-associated genes, such as Earlier Degraded Tapetum 1 (EDT1), Tapetum Degeneration Retardation (TDR), OsPCF5, and OsTCP21, directly or indirectly determined by AGO1d-mediated gene silencing. The loss of function of 21-nt phasiRNAs can result in upregulation of their targets and causes varying degrees of defects in male fertility and grain setting. Our results highlight the essential functions of 21-nt phasiRNAs in temperature-sensitive male sterility in rice and suggest their promising application in 2-line hybrid rice breeding in the future.

Receptor-like kinases and their signaling cascades for plant male fertility: loyal messengers

Receptor-like kinases (RLKs) are evolved for plant cell-cell communications. The typical RLK protein contains an extracellular and hypervariable N-terminus to perceive various signals, a transmembrane domain to anchor into plasma membrane, and a cytoplasmic, highly conserved kinase domain to phosphorylate target proteins. To date, RLKs have manifested their significance in a myriad of biological processes during plant reproductive growth, especially in male fertility. This review first summarizes a recent update on RLKs and their interacting protein partners controlling anther and pollen development, pollen release from dehisced anther, and pollen function during pollination and fertilization. Then, regulatory networks of RLK signaling pathways are proposed. In addition, we predict RLKs in maize and rice genome, obtain homologs of well-studied RLKs from phylogeny of three subfamilies and then analyze their expression patterns in developing anthers of maize and rice to excavate potential RLKs regulating male fertility in crops. Finally, current challenges and future prospects regarding RLKs are discussed. This review will contribute to a better understanding of plant male fertility control by RLKs, creating potential male sterile lines, and inspiring innovative crop breeding methods.

The plant siRNA landscape

Whereas micro (mi)RNAs are considered the clean, noble side of the small RNA world, small interfering (si)RNAs are often seen as a noisy set of molecules whose barbarian acronyms reflect a large diversity of often elusive origins and functions. Twenty-five years after their discovery in plants, however, new classes of siRNAs are still being identified, sometimes in discrete tissues or at particular developmental stages, making the plant siRNA world substantially more complex and subtle than originally anticipated. Focusing primarily on the model Arabidopsis, we review here the plant siRNA landscape, including transposable elements (TE)-derived siRNAs, a vast array of non-TE-derived endogenous siRNAs, as well as exogenous siRNAs produced in response to invading nucleic acids such as viruses or transgenes. We primarily emphasize the extraordinary sophistication and diversity of their biogenesis and, secondarily, the variety of their known or presumed functions, including via non-cell autonomous activities, in the sporophyte, gametophyte, and shortly after fertilization.

Genome-Wide Identification and Expression Analyses of the FAR1/FHY3 Gene Family Provide Insight into Inflorescence Development in Maize

As transcription factors derived from transposase, FAR-RED IMPAIRED RESPONSE1 (FAR1) and its homolog FHY3 play crucial roles in the regulation of light signaling and various stress responses by coordinating the expression of downstream target genes. Despite the extensive investigation of the FAR1/FHY3 family in Arabidopsis thaliana and other species, a comprehensive examination of these genes in maize has not been conducted thus far. In this study, we employed a genomic mining approach to identify 16 ZmFAR1 genes in the maize inbred line B73, which were further classified into five subgroups based on their phylogenetic relationships. The present study characterized the predicted polypeptide sequences, molecular weights, isoelectric points, chromosomal distribution, gene structure, conserved motifs, subcellular localizations, phylogenetic relationships, and cis-regulatory elements of all members belonging to the ZmFAR1 family. Furthermore, the tissue-specific expression of the 16 ZmFAR1 genes was analyzed using RNA-seq, and their expression patterns under far-red light conditions were validated in the ear and tassel through qRT-qPCR. The observed highly temporal and spatial expression patterns of these ZmFAR1 genes were likely associated with their specific functional capabilities under different light conditions. Further analysis revealed that six ZmFAR1 genes (ZmFAR1-1, ZmFAR1-10, ZmFAR1-11, ZmFAR1-12, ZmFAR1-14, and ZmFAR1-15) exhibited a response to simulated shading treatment and actively contributed to the development of maize ears. Through the integration of expression quantitative trait loci (eQTL) analyses and population genetics, we identified the presence of potential causal variations in ZmFAR1-14 and ZmFAR1-9, which play a crucial role in regulating the kernel row number and kernel volume weight, respectively. In summary, this study represents the initial identification and characterization of ZmFAR1 family members in maize, uncovering the functional variation in candidate regulatory genes associated with the improvement of significant agronomic traits during modern maize breeding.

Toward Transgene-Free Transposon-Mediated Biological Mutagenesis for Plant Breeding

Genetic diversity is a key factor for plant breeding. The birth of novel genic and genomic variants is also crucial for plant adaptation in nature. Therefore, the genomes of almost all living organisms possess natural mutagenic mechanisms. Transposable elements (TEs) are a major mutagenic force driving genetic diversity in wild plants and modern crops. The relatively rare TE transposition activity during the thousand-year crop domestication process has led to the phenotypic diversity of many cultivated species. The utilization of TE mutagenesis by artificial and transient acceleration of their activity in a controlled mode is an attractive foundation for a novel type of mutagenesis called TE-mediated biological mutagenesis. Here, I focus on TEs as mutagenic sources for plant breeding and discuss existing and emerging transgene-free approaches for TE activation in plants. Furthermore, I also review the non-randomness of TE insertions in a plant genome and the molecular and epigenetic factors involved in shaping TE insertion preferences. Additionally, I discuss the molecular mechanisms that prevent TE transpositions in germline plant cells (e.g., meiocytes, pollen, egg and embryo cells, and shoot apical meristem), thereby reducing the chances of TE insertion inheritance. Knowledge of these mechanisms can expand the TE activation toolbox using novel gene targeting approaches. Finally, the challenges and future perspectives of plant populations with induced novel TE insertions (iTE plant collections) are discussed.

Heat-responsive microRNAs participate in regulating the pollen fertility stability of CMS-D2 restorer line under high-temperature stress

Anther development and pollen fertility of cytoplasmic male sterility (CMS) conditioned by Gossypium harknessii cytoplasm (CMS-D2) restorer lines are susceptible to continuous high-temperature (HT) stress in summer, which seriously hinders the large-scale application of "three-line" hybrids in production. Here, integrated small RNA, transcriptome, degradome, and hormone profiling was performed to explore the roles of microRNAs (miRNAs) in regulating fertility stability in mature pollens of isonuclear alloplasmic near-isogenic restorer lines NH and SH under HT stress at two environments. A total of 211 known and 248 novel miRNAs were identified, of which 159 were differentially expressed miRNAs (DEMs). Additionally, 45 DEMs in 39 miRNA clusters (PmCs) were also identified, and most highly expressed miRNAs were significantly induced in SH under extreme HT, especially four MIR482 and six MIR6300 family miRNAs. PmC28 was located in the fine-mapped interval of the Rf1 gene and contained two DEMs, gra-miR482_L-2R + 2 and gma-miR2118a-3p_R + 1_1ss18TG. Transcriptome sequencing identified 6281 differentially expressed genes, of which heat shock protein (HSP)-related genes, such as HSP70, HSP22, HSP18.5-C, HSP18.2 and HSP17.3-B, presented significantly reduced expression levels in SH under HT stress. Through integrating multi-omics data, we constructed a comprehensive molecular network of miRNA-mRNA-gene-KEGG containing 35 pairs of miRNA/target genes involved in regulating the pollen development in response to HT, among which the mtr-miR167a_R + 1, tcc-miR167c and ghr-miR390a, tcc-miR396c_L-1 and ghr-MIR169b-p3_1ss6AG regulated the pollen fertility by influencing ARF8 responsible for the auxin signal transduction, ascorbate and aldarate metabolism, and the sugar and lipid metabolism and transport pathways, respectively. Further combination with hormone analysis revealed that HT-induced jasmonic acid signaling could activate the expression of downstream auxin synthesis-related genes and cause excessive auxin accumulation, followed by a cascade of auxin signal transduction, ultimately resulting in pollen abortion. The results provide a new understanding of how heat-responsive miRNAs regulate the stability of fertility restoration for CMS-D2 cotton under heat stress.

Pseudouridine guides germline small RNA transport and epigenetic inheritance

Epigenetic modifications that arise during plant and animal development, such as DNA and histone modification, are mostly reset during gamete formation, but some are inherited from the germline including those marking imprinted genes1. Small RNAs guide these epigenetic modifications, and some are also inherited by the next generation2,3. In C. elegans, these inherited small RNAs have poly (UG) tails4, but how inherited small RNAs are distinguished in other animals and plants is unknown. Pseudouridine (Ψ) is the most abundant RNA modification but has not been explored in small RNAs. Here, we develop novel assays to detect Ψ in short RNA sequences, demonstrating its presence in mouse and Arabidopsis microRNAs and their precursors. We also detect substantial enrichment in germline small RNAs, namely epigenetically activated siRNAs (easiRNAs) in Arabidopsis pollen, and piwi-interacting piRNAs in mouse testis. In pollen, pseudouridylated easiRNAs are localized to sperm cells, and we found that PAUSED/HEN5 (PSD), the plant homolog of Exportin-t, interacts genetically with Ψ and is required for transport of easiRNAs into sperm cells from the vegetative nucleus. We further show that Exportin-t is required for the triploid block: chromosome dosage-dependent seed lethality that is epigenetically inherited from pollen. Thus, Ψ has a conserved role in marking inherited small RNAs in the germline.

PhasiHunter: a robust phased siRNA regulatory cascade mining tool based on multiple reference sequences

SUMMARY

In recent years, phased small interfering RNA has been found to play crucial roles in many biological processes in plants. However, efficiently predicting phasiRNA regulatory cascades with computational methods is still challenging. Here, we introduce PhasiHunter, a phasiRNA regulatory network prediction tool that has several distinctive features compared to existing tools: (i) PhasiHunter employs two major phasiRNA prediction algorithms, namely phase score and hypergeometric distribution-based methods, to ensure the integrity and accuracy of prediction; (ii) PhasiHunter can identify phasiRNAs and their regulatory networks based on multiple reference sequences and the predicted results can be automatically integrated; (iii) PhasiHunter can efficiently identify the phasiRNAs generated through alternative splicing events; and (iv) the excellent data structure and parallel computing architecture allow PhasiHunter to predict phasiRNAs and their regulatory pathways with high efficiency.

AVAILABILITY AND IMPLEMENTATION

PhasiHunter is an open-source tool that is available at https://github.com/HuangLab-CBI/PhasiHunter.

Small RNAs: Promising Molecules to Tackle Climate Change Impacts in Coffee Production

Over the centuries, human society has evolved based on the ability to select and use more adapted species for food supply, which means making plant species tastier and more productive in particular environmental conditions. However, nowadays, this scenario is highly threatened by climate change, especially by the changes in temperature and greenhouse gasses that directly affect photosynthesis, which highlights the need for strategic studies aiming at crop breeding and guaranteeing food security. This is especially worrying for crops with complex phenology, genomes with low variability, and the ones that support a large production chain, such as Coffea sp. L. In this context, recent advances shed some light on the genome function and transcriptional control, revealing small RNAs (sRNAs) that are responsible for environmental cues and could provide variability through gene expression regulation. Basically, sRNAs are responsive to environmental changes and act on the transcriptional and post-transcriptional gene silencing pathways that regulate gene expression and, consequently, biological processes. Here, we first discuss the predicted impact of climate changes on coffee plants and coffee chain production and then the role of sRNAs in response to environmental changes, especially temperature, in different species, together with their potential as tools for genetic improvement. Very few studies in coffee explored the relationship between sRNAs and environmental cues; thus, this review contributes to understanding coffee development in the face of climate change and towards new strategies of crop breeding.

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.

Epigenetic control of transposons during plant reproduction: From meiosis to hybrid seeds

The regulation of transposable elements (TEs) requires overlapping epigenetic modifications that must be reinforced every cell division and generation. In plants, this is achieved by multiple pathways including small RNAs, DNA methylation, and repressive histone marks that act together to control TE expression and activity throughout the entire life cycle. However, transient TE activation is observed during reproductive transitions as a result of epigenome reprogramming, thus providing windows of opportunity for TE proliferation and epigenetic novelty. Ultimately, these events may originate complex TE-driven transcriptional networks or cell-to-cell communication strategies via mobile small RNAs. In this review, we discuss recent findings and current understanding of TE regulation during sexual plant reproduction, and its implications for fertility, early seed development, and epigenetic inheritance.

Noncoding RNAs and their roles in regulating the agronomic traits of crops

Molecular breeding is one of the most effective methods for improving the performance of crops. Understanding the genome features of crops, especially the physiological functions of individual genes, is of great importance to molecular breeding. Evidence has shown that genomes of both animals and plants transcribe numerous non-coding RNAs, which are involved in almost every aspect of development. In crops, an increasing number of studies have proven that non-coding RNAs are new genetic resources for regulating crop traits. In this review, we summarize the current knowledge of non-coding RNAs, which are potential crop trait regulators, and focus on the functions of long non-coding RNAs (lncRNAs) in determining crop grain yield, phased small-interfering RNAs (phasiRNAs) in regulating fertility, small interfering RNAs (siRNAs) and microRNAs (miRNAs) in facilitating plant immune response and disease resistance, and miRNAs mediating nutrient and metal stress. Finally, we also discuss the next-generation method for ncRNA application in crop domestication and breeding.

ZmAGO18b negatively regulates maize resistance against southern leaf blight

KEY MESSAGE

Here, we report that ZmAGO18b encoding an argonaute protein is a negative regulator of maize resistance against southern leaf blight. Southern leaf blight caused by fungal pathogen Cochliobolus heterostrophus is a destructive disease on maize throughout the world. Argonaute (AGO) proteins, key regulators in small RNA pathway, play important roles in plant defense. But whether they have function in maize resistance against C. heterostrophus is unknown. Association analysis between the nucleic variation of 18 ZmAGO loci with disease phenotype against C. heterostrophus was performed, and the ZmAGO18b locus was identified to be associated with resistance against C. heterostrophus. Overexpression of ZmAGO18b gene suppresses maize resistance against C. heterostrophus, and mutation of ZmAGO18b enhances maize resistance against C. heterostrophus. Further, we identified the resistant haplotype of ZmAGO18b by association analysis of natural variation in ZmAGO18b genomic DNA sequences with seedling resistance phenotypes against C. heterostrophus and confirmed the resistant haplotype is co-segregated with resistance phenotypes against C. heterostrophus in two F₂ populations. In sum, this study reports that ZmAGO18b negatively regulates maize resistance against C. heterostrophus.

Spatial distribution of three ARGONAUTEs regulates the anther phasiRNA pathway

Argonaute protein (AGO) in association with small RNAs is the core machinery of RNA silencing, an essential mechanism for precise development and defense against pathogens in many organisms. Here, we identified two AGOs in rice anthers, AGO1b and AGO1d, that interact with phased small interfering RNAs (phasiRNAs) derived from numerous long non-coding RNAs. Moreover, 3D-immunoimaging and mutant analysis indicated that rice AGO1b and AGO1d cell type-specifically regulate anther development by acting as mobile carriers of these phasiRNAs from the somatic cell layers to the germ cells in anthers. Our study also highlights a new mode of reproductive RNA silencing via the specific nuclear and cytoplasmic localization of three AGOs, AGO1b, AGO1d, and MEL1, in rice pollen mother cells.

Phylogenetic analyses of seven protein families refine the evolution of small RNA pathways in green plants

Several protein families participate in the biogenesis and function of small RNAs (sRNAs) in plants. Those with primary roles include Dicer-like (DCL), RNA-dependent RNA polymerase (RDR), and Argonaute (AGO) proteins. Protein families such as double-stranded RNA-binding (DRB), SERRATE (SE), and SUPPRESSION OF SILENCING 3 (SGS3) act as partners of DCL or RDR proteins. Here, we present curated annotations and phylogenetic analyses of seven sRNA pathway protein families performed on 196 species in the Viridiplantae (aka green plants) lineage. Our results suggest that the RDR3 proteins emerged earlier than RDR1/2/6. RDR6 is found in filamentous green algae and all land plants, suggesting that the evolution of RDR6 proteins coincides with the evolution of phased small interfering RNAs (siRNAs). We traced the origin of the 24-nt reproductive phased siRNA-associated DCL5 protein back to the American sweet flag (Acorus americanus), the earliest diverged, extant monocot species. Our analyses of AGOs identified multiple duplication events of AGO genes that were lost, retained, or further duplicated in subgroups, indicating that the evolution of AGOs is complex in monocots. The results also refine the evolution of several clades of AGO proteins, such as AGO4, AGO6, AGO17, and AGO18. Analyses of nuclear localization signal sequences and catalytic triads of AGO proteins shed light on the regulatory roles of diverse AGOs. Collectively, this work generates a curated and evolutionarily coherent annotation for gene families involved in plant sRNA biogenesis/function and provides insights into the evolution of major sRNA pathways.

Small RNA-mediated DNA methylation during plant reproduction

Reproductive tissues are a rich source of small RNAs, including several classes of short interfering (si)RNAs that are restricted to this stage of development. In addition to RNA polymerase IV-dependent 24-nt siRNAs that trigger canonical RNA-directed DNA methylation, abundant reproductive-specific siRNAs are produced from companion cells adjacent to the developing germ line or zygote and may move intercellularly before inducing methylation. In some cases, these siRNAs are produced via non-canonical biosynthesis mechanisms or from sequences with little similarity to transposons. While the precise role of these siRNAs and the methylation they trigger is unclear, they have been implicated in specifying a single megaspore mother cell, silencing transposons in the male germ line, mediating parental dosage conflict to ensure proper endosperm development, hypermethylation of mature embryos, and trans-chromosomal methylation in hybrids. In this review, we summarize the current knowledge of reproductive siRNAs, including their biosynthesis, transport, and function.

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.

IQ67 DOMAIN protein 21 is critical for indentation formation in pavement cell morphogenesis

In plants, cortical microtubules anchor to the plasma membrane in arrays and play important roles in cell shape. However, the molecular mechanism of microtubule binding proteins, which connect the plasma membrane and cortical microtubules in cell morphology remains largely unknown. Here, we report that a plasma membrane and microtubule dual-localized IQ67 domain protein, IQD21, is critical for cotyledon pavement cell (PC) morphogenesis in Arabidopsis. iqd21 mutation caused increased indentation width, decreased lobe length, and similar lobe number of PCs, whereas IQD21 overexpression had a different effect on cotyledon PC shape. Weak overexpression led to increased lobe number, decreased indentation width, and similar lobe length, while moderate or great overexpression resulted in decreased lobe number, indentation width, and lobe length of PCs. Live-cell observations revealed that IQD21 accumulation at indentation regions correlates with lobe initiation and outgrowth during PC development. Cell biological and genetic approaches revealed that IQD21 promotes transfacial microtubules anchoring to the plasma membrane via its polybasic sites and bundling at the indentation regions in both periclinal and anticlinal walls. IQD21 controls cortical microtubule organization mainly through promoting Katanin 1-mediated microtubule severing during PC interdigitation. These findings provide the genetic evidence that transfacial microtubule arrays play a determinant role in lobe formation, and the insight into the molecular mechanism of IQD21 in transfacial microtubule organization at indentations and puzzle-shaped PC development.

Plant non-coding RNAs function in pollen development and male sterility

Male sterility is classified as either cytoplasmic male sterility (CMS) or genic male sterility (GMS). Generally, CMS involves mitochondrial genomes interacting with the nuclear genome, while GMS is caused by nuclear genes alone. Male sterility is regulated by multilevel mechanisms in which non-coding RNAs (ncRNAs), including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and phased small interfering RNAs (phasiRNAs), which have been proven to be critical elements. The development of high-throughput sequencing technology offers new opportunities to evaluate the genetic mechanism of ncRNAs in plant male sterility. In this review, we summarize the critical ncRNAs that regulate gene expression in ways dependent on or independent of hormones, which involve the differentiation of the stamen primordia, degradation of the tapetum, formation of microspores, and the release of pollen. In addition, the key mechanisms of the miRNA-lncRNA-mRNA interaction networks mediating male sterility in plants are elaborated. We present a different perspective on exploring the ncRNA-mediated regulatory pathways that control CMS in plants and create male-sterile lines through hormones or genome editing. A refined understanding of the ncRNA regulatory mechanisms in plant male sterility for the development of new sterile lines would be conducive to improve hybridization breeding.

Mobile ARGONAUTE 1d binds 22-nt miRNAs to generate phasiRNAs important for low-temperature male fertility in rice

Phased small interfering RNAs (phasiRNAs) are abundantly expressed in anthers and linked to environment-related male fertility in grasses, yet how they function under different environmental conditions remains unclear. Here, we identified a rice (Oryza sativa) low temperature-induced Argonaute (AGO) protein, OsAGO1d, that is responsible for generating phasiRNAs and preserving male fertility at low temperature. Loss of OsAGO1d function causes low-temperature male sterility associated with delayed programmed cell death of tapetal cells during anther development. OsAGO1d binds miR2118 and miR2275 family members and triggers phasiRNA biogenesis; it also binds 21-nt phasiRNAs with a 5' terminal U. In total, phasiRNAs from 972 loci are OsAGO1d-dependent. OsAGO1d protein moves from anther wall cells into meiocytes, where it loads miR2275 to produce 24-nt phasiRNAs. Together, our results show that OsAGO1d acts as a mobile signal to fine-tune phasiRNA production and this function is important for male fertility at low temperature.

Rice GLUCAN SYNTHASE-LIKE5 promotes anther callose deposition to maintain meiosis initiation and progression

Callose is a plant cell wall polysaccharide whose deposition is spatiotemporally regulated in various developmental processes and environmental stress responses. The appearance of callose in premeiotic anthers is a prominent histological hallmark for the onset of meiosis in flowering plants; however, the biological role of callose in meiosis remains unknown. Here, we show that rice (Oryza sativa) GLUCAN SYNTHASE LIKE5 (OsGSL5), a callose synthase, localizes on the plasma membrane of pollen mother cells (PMCs) and is responsible for biogenesis of callose in anther locules through premeiotic and meiotic stages. In Osgsl5 mutant anthers mostly lacking callose deposition, aberrant PMCs accompanied by aggregated, unpaired, or multivalent chromosomes were frequently observed and, furthermore, a considerable number of mutant PMCs had untimely progress into meiosis compared to that of wild-type PMCs. Immunostaining of meiosis-specific protein HOMOLOGOUS PAIRING ABERRATION IN RICE MEIOSIS2 in premeiotic PMCs revealed precocious meiosis entry in Osgsl5 anthers. These findings provide insights into the function of callose in controlling the timing of male meiosis initiation and progression, in addition to roles in microsporogenesis, in flowering plants.

Coexpression network and trans-activation analyses of maize reproductive phasiRNA loci

The anther-enriched phased, small interfering RNAs (phasiRNAs) play vital roles in sustaining male fertility in grass species. Their long non-coding precursors are synthesized by RNA polymerase II and are likely regulated by transcription factors (TFs). A few putative transcriptional regulators of the 21- or 24-nucleotide phasiRNA loci (referred to as 21- or 24-PHAS loci) have been identified in maize (Zea mays), but whether any of the individual TFs or TF combinations suffice to activate any PHAS locus is unclear. Here, we identified the temporal gene coexpression networks (modules) associated with maize anther development, including two modules highly enriched for the 21- or 24-PHAS loci. Comparisons of these coexpression modules and gene sets dysregulated in several reported male sterile TF mutants provided insights into TF timing with regard to phasiRNA biogenesis, including antagonistic roles for OUTER CELL LAYER4 and MALE STERILE23. Trans-activation assays in maize protoplasts of individual TFs using bulk-protoplast RNA-sequencing showed that two of the TFs coexpressed with 21-PHAS loci could activate several 21-nucleotide phasiRNA pathway genes but not transcription of 21-PHAS loci. Screens for combinatorial activities of these TFs and, separately, the recently reported putative transcriptional regulators of 24-PHAS loci using single-cell (protoplast) RNA-sequencing, did not detect reproducible activation of either 21-PHAS or 24-PHAS loci. Collectively, our results suggest that the endogenous transcriptional machineries and/or chromatin states in the anthers are necessary to activate reproductive PHAS loci.

The ZmMYB84-ZmPKSB regulatory module controls male fertility through modulating anther cuticle-pollen exine trade-off in maize anthers

Anther cuticle and pollen exine are two crucial lipid layers that ensure normal pollen development and pollen-stigma interaction for successful fertilization and seed production in plants. Their formation processes share certain common pathways of lipid biosynthesis and transport across four anther wall layers. However, molecular mechanism underlying a trade-off of lipid-metabolic products to promote the proper formation of the two lipid layers remains elusive. Here, we identified and characterized a maize male-sterility mutant pksb, which displayed denser anther cuticle but thinner pollen exine as well as delayed tapetal degeneration compared with its wild type. Based on map-based cloning and CRISPR/Cas9 mutagenesis, we found that the causal gene (ZmPKSB) of pksb mutant encoded an endoplasmic reticulum (ER)-localized polyketide synthase (PKS) with catalytic activities to malonyl-CoA and midchain-fatty acyl-CoA to generate triketide and tetraketide α-pyrone. A conserved catalytic triad (C171, H320 and N353) was essential for its enzymatic activity. ZmPKSB was specifically expressed in maize anthers from stages S8b to S9-10 with its peak at S9 and was directly activated by a transcription factor ZmMYB84. Moreover, loss function of ZmMYB84 resulted in denser anther cuticle but thinner pollen exine similar to the pksb mutant. The ZmMYB84-ZmPKSB regulatory module controlled a trade-off between anther cuticle and pollen exine formation by altering expression of a series of genes related to biosynthesis and transport of sporopollenin, cutin and wax. These findings provide new insights into the fine-tuning regulation of lipid-metabolic balance to precisely promote anther cuticle and pollen exine formation in plants.

DNA methylation dynamics during germline development

DNA methylation plays essential homeostatic functions in eukaryotic genomes. In animals, DNA methylation is also developmentally regulated and, in turn, regulates development. In the past two decades, huge research effort has endorsed the understanding that DNA methylation plays a similar role in plant development, especially during sexual reproduction. The power of whole-genome sequencing and cell isolation techniques, as well as bioinformatics tools, have enabled recent studies to reveal dynamic changes in DNA methylation during germline development. Furthermore, the combination of these technological advances with genetics, developmental biology and cell biology tools has revealed functional methylation reprogramming events that control gene and transposon activities in flowering plant germlines. In this review, we discuss the major advances in our knowledge of DNA methylation dynamics during male and female germline development in flowering plants.

Anther development-The long road to making pollen

Anthers express the most genes of any plant organ, and their development involves sequential redifferentiation of many cell types to perform distinctive roles from inception through pollen dispersal. Agricultural yield and plant breeding depend on understanding and consequently manipulating anthers, a compelling motivation for basic plant biology research to contribute. After stamen initiation, two theca form at the tip, and each forms an adaxial and abaxial lobe composed of pluripotent Layer 1-derived and Layer 2-derived cells. After signal perception or self-organization, germinal cells are specified from Layer 2-derived cells, and these secrete a protein ligand that triggers somatic differentiation of their neighbors. Historically, recovery of male-sterile mutants has been the starting point for studying anther biology. Many genes and some genetic pathways have well-defined functions in orchestrating subsequent cell fate and differentiation events. Today, new tools are providing more detailed information; for example, the developmental trajectory of germinal cells illustrates the power of single cell RNA-seq to dissect the complex journey of one cell type. We highlight ambiguities and gaps in available data to encourage attention on important unresolved issues.

Maize sterility gene DRP1 encodes a desiccation-related protein that is critical for Ubisch bodies and pollen exine development

Desiccation tolerance is a remarkable feature of pollen, seeds, and resurrection-type plants. Exposure to desiccation stress can cause sporophytic defects, resulting in male sterility. Here, we report the novel maize sterility gene DRP1 (Desiccation-Related Protein 1), which was identified by bulked-segregant analysis sequencing and encodes a desiccation-related protein. Loss of function of DRP1 results in abnormal Ubisch bodies, defective tectum of the pollen exine, and complete male sterility. Our results suggest that DRP1 may facilitate anther dehydration to maintain appropriate water status. DRP1 is a secretory protein that is specifically expressed in the tapetum and microspore from the tetrad to the uninucleate microspore stage. Differentially expressed genes in drp1 are enriched in Gene Ontology terms for pollen exine formation, polysaccharide catabolic process, extracellular region, and response to heat. In addition, DRP1 is a target of selection that appears to have played an important role in the spread of maize from tropical/subtropical to temperate regions. Taken together, our results suggest that DRP1 encodes a desiccation-related protein whose loss of function causes male sterility. Our findings provide a potential genetic resource that may be used to design crops for heterosis utilization.

Temperature-sensitive male sterility in rice determined by the roles of AGO1d in reproductive phasiRNA biogenesis and function

Phased secondary siRNAs (phasiRNAs) are broadly present in the reproductive tissues of flowering plants, with spatial-temporal specificity. However, the ARGONAUTE (AGO) proteins associated with phasiRNAs and their miRNA triggers remain elusive. Here, through histological and high-throughput sequencing analyses, we show that rice AGO1d, which is specifically expressed in anther wall cells before and during meiosis, associates with both miR2118 and miR2275 to mediate phasiRNA biogenesis. AGO1d preferentially binds to miR2118-triggered 21-nucleotide (nt) phasiRNAs with a 5'-terminal uridine, suggesting a dual role in phasiRNA biogenesis and function. Depletion of AGO1d causes a reduction of 21- and 24-nt phasiRNAs and temperature-sensitive male sterility. At lower temperatures, anthers of the ago1d mutant predominantly show excessive tapetal cells with little starch accumulation during pollen formation, possibly caused by the dysregulation of cell metabolism. These results uncover an essential role of AGO1d in rice anther development at lower temperatures and demonstrate coordinative roles of AGO proteins during reproductive phasiRNA biogenesis and function.

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.

Expansion and contraction of small RNA and methylation machinery throughout plant evolution

The revolution in sequencing has created a wealth of plant genomes that can be mined to understand the evolution of biological complexity. Complexity is often driven by gene duplication, which allows paralogs to specialize in an activity of the ancestral gene or acquire novel functions. Angiosperms encode a variety of gene silencing pathways that share related machinery for small RNA biosynthesis and function. Recent phylogenetic analysis of these gene families plots the expansion, specialization, and occasional contraction of this core machinery. This analysis reveals the ancient origin of RNA-directed DNA Methylation in early land plants, or possibly their algal ancestors, as well as ongoing duplications that evolve novel small RNA pathways.

3D multiple immunoimaging using whole male organs in rice

Spatiotemporal regulation of proteins and RNAs is essential for the precise development of reproductive tissues in many organisms. The anther, a prominent part of the male reproductive organ in plants, contains several somatic cell layers named the anther wall and, within it, the germ cells. Here, we successfully developed a simple 3D organ-immunoimaging technique for rice anthers, which distinguishes each individual cell from the four somatic cell layers and germ cells without the need for transformation, embedding, sectioning, or clearing. The 3D immunostaining method is also applicable to the intracellular localization of meiosis-specific proteins in meiocytes, as exemplified by MEL1, a germ cell-specific ARGONAUTE in the cytoplasm, and ZEP1, a pachytene marker on meiotic chromosomes. Our 3D multiple immunostaining method with single-cell and intracellular resolution will contribute to a comprehensive organ-level elucidation of molecular mechanisms and cellular connectivity.

Epigenetic Dynamics and Regulation of Plant Male Reproduction

Flowering plant male germlines develop within anthers and undergo epigenetic reprogramming with dynamic changes in DNA methylation, chromatin modifications, and small RNAs. Profiling the epigenetic status using different technologies has substantially accumulated information on specific types of cells at different stages of male reproduction. Many epigenetically related genes involved in plant gametophyte development have been identified, and the mutation of these genes often leads to male sterility. Here, we review the recent progress on dynamic epigenetic changes during pollen mother cell differentiation, microsporogenesis, microgametogenesis, and tapetal cell development. The reported epigenetic variations between male fertile and sterile lines are summarized. We also summarize the epigenetic regulation-associated male sterility genes and discuss how epigenetic mechanisms in plant male reproduction can be further revealed.

Transcriptome-wide analyses of RNA m6A methylation in hexaploid wheat reveal its roles in mRNA translation regulation

N6-methyladenosine (m6A) is the most abundant RNA modification in eukaryotic messenger RNAs. m6A was discovered in wheat about 40 years ago; however, its potential roles in wheat remain unknown. In this study, we profiled m6As in spikelets transcriptome at the flowering stage of hexaploid wheat and found that m6As are evenly distributed across the A, B, and D subgenomes but their extents and locations vary across homeologous genes. m6As are enriched in homeologous genes with close expression levels and the m6A methylated genes are more conserved. The extent of m6A methylation is negatively correlated with mRNA expression levels and its presence on mRNAs has profound impacts on mRNA translation in a location-dependent manner. Specifically, m6As within coding sequences and 3'UTRs repress the translation of mRNAs while the m6As within 5'UTRs and start codons could promote it. The m6A-containing mRNAs are significantly enriched in processes and pathways of "translation" and "RNA transport," suggesting the potential role of m6As in regulating the translation of genes involved in translation regulation. Our data also show a stronger translation inhibition by small RNAs (miRNA and phasiRNA) than by m6A methylation, and no synergistical effect between the two was observed. We propose a secondary amplification machinery of translation regulation triggered by the changes in m6A methylation status. Taken together, our results suggest translation regulation as a key role played by m6As in hexaploid wheat.

Spotting the Targets of the Apospory Controller TGS1 in Paspalum notatum

Sexuality and apomixis are interconnected plant reproductive routes possibly behaving as polyphenic traits under the influence of the environment. In the subtropical grass Paspalum notatum, one of the controllers of apospory, a main component of gametophytic apomixis reproduction, is TRIMETHYLGUANOSINE SYNTHASE 1 (TGS1), a multifunctional gene previously associated with RNA cleavage regulation (including mRNA splicing as well as rRNA and miRNA processing), transcriptional modulation and the establishment of heterochromatin. In particular, the downregulation of TGS1 induces a sexuality decline and the emergence of aposporous-like embryo sacs. The present work was aimed at identifying TGS1 target RNAs expressed during reproductive development of Paspalum notatum. First, we mined available RNA databases originated from spikelets of sexual and apomictic plants, which naturally display a contrasting TGS1 representation, to identify differentially expressed mRNA splice variants and miRNAs. Then, the role of TGS1 in the generation of these particular molecules was investigated in antisense tgs1 sexual lines. We found that CHLOROPHYLL A-B BINDING PROTEIN 1B-21 (LHC Ib-21, a component of the chloroplast light harvesting complex), QUI-GON JINN (QGJ, encoding a MAP3K previously associated with apomixis) and miR2275 (a meiotic 24-nt phasi-RNAs producer) are directly or indirectly targeted by TGS1. Our results point to a coordinated control exercised by signal transduction and siRNA machineries to induce the transition from sexuality to apomixis.

3D Imaging and In Situ Hybridization for Uncovering the Functions of MicroRNA in Rice Anther

Small RNAs specifically expressed in reproductive tissues are key regulators of germline development in eukaryotes. Rice microRNA2118 (miR2118), which is enriched during reproduction in grasses, is a trigger to produce phased small interfering RNAs (phasiRNAs). These phasiRNAs demonstrate the temporal regulation with premeiotic phasiRNAs and meiotic phasiRNAs in anther development. Furthermore, the site-specific regulation via miR2118 and phasiRNAs is of importance in soma and germ development in anthers. Accordingly, histological imaging methods are essential tools for understanding spatiotemporal regulation during reproduction and elucidating the reproductive roles of miRNAs and phasiRNAs. We successfully developed a method to visualize the three-dimensional (3D) structure of entire rice anthers, which can also be used for distinguishing the internal structure of the anthers in other plants. Here, we describe the detailed methods of in situ hybridization for miR2118 localization and the visualization of the 3D structure of entire anthers of rice.

Exploring key developmental phases and phase-specific genes across the entirety of anther development in maize

Anther development from stamen primordium to pollen dispersal is complex and essential to sexual reproduction. How this highly dynamic and complex developmental process is controlled genetically is not well understood, especially for genes involved in specific key developmental phases. Here we generated RNA sequencing libraries spanning 10 key stages across the entirety of anther development in maize (Zea mays). Global transcriptome analyses revealed distinct phases of cell division and expansion, meiosis, pollen maturation, and mature pollen, for which we detected 50, 245, 42, and 414 phase-specific marker genes, respectively. Phase-specific transcription factor genes were significantly enriched in the phase of meiosis. The phase-specific expression of these marker genes was highly conserved among the maize lines Chang7-2 and W23, indicating they might have important roles in anther development. We explored a desiccation-related protein gene, ZmDRP1, which was exclusively expressed in the tapetum from the tetrad to the uninucleate microspore stage, by generating knockout mutants. Notably, mutants in ZmDRP1 were completely male-sterile, with abnormal Ubisch bodies and defective pollen exine. Our work provides a glimpse into the gene expression dynamics and a valuable resource for exploring the roles of key phase-specific genes that regulate anther development.

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.

Accumulation dynamics of ARGONAUTE proteins during meiosis in Arabidopsis

Meiosis is a specialized cell division that is key for reproduction and genetic diversity in sexually reproducing plants. Recently, different RNA silencing pathways have been proposed to carry a specific activity during meiosis, but the pathways involved during this process remain unclear. Here, we explored the subcellular localization of different ARGONAUTE (AGO) proteins, the main effectors of RNA silencing, during male meiosis in Arabidopsis thaliana using immunolocalizations with commercially available antibodies. We detected the presence of AGO proteins associated with posttranscriptional gene silencing (AGO1, 2, and 5) in the cytoplasm and the nucleus, while AGOs associated with transcriptional gene silencing (AGO4 and 9) localized exclusively in the nucleus. These results indicate that the localization of different AGOs correlates with their predicted roles at the transcriptional and posttranscriptional levels and provide an overview of their timing and potential role during meiosis.

Unravelling the Role of Epigenetic Modifications in Development and Reproduction of Angiosperms: A Critical Appraisal

Epigenetics are the heritable changes in gene expression patterns which occur without altering DNA sequence. These changes are reversible and do not change the sequence of the DNA but can alter the way in which the DNA sequences are read. Epigenetic modifications are induced by DNA methylation, histone modification, and RNA-mediated mechanisms which alter the gene expression, primarily at the transcriptional level. Such alterations do control genome activity through transcriptional silencing of transposable elements thereby contributing toward genome stability. Plants being sessile in nature are highly susceptible to the extremes of changing environmental conditions. This increases the likelihood of epigenetic modifications within the composite network of genes that affect the developmental changes of a plant species. Genetic and epigenetic reprogramming enhances the growth and development, imparts phenotypic plasticity, and also ensures flowering under stress conditions without changing the genotype for several generations. Epigenetic modifications hold an immense significance during the development of male and female gametophytes, fertilization, embryogenesis, fruit formation, and seed germination. In this review, we focus on the mechanism of epigenetic modifications and their dynamic role in maintaining the genomic integrity during plant development and reproduction.

Small regulatory RNAs in rice epigenetic regulation

Plant small RNAs (sRNAs) are short non-coding RNAs that are implicated in various regulatory processes involving post-transcriptional gene silencing and epigenetic gene regulation. In epigenetic regulation, sRNAs are primarily involved in RNA-directed DNA methylation (RdDM) pathways. sRNAs in the RdDM pathways play a role not only in the suppression of transposable element (TE) activity but also in gene expression regulation. Although the major components of the RdDM pathways have been well studied in Arabidopsis, recent studies have revealed that the RdDM pathways in rice have important biological functions in stress response and developmental processes. In this review, we summarize and discuss recent literature on sRNA-mediated epigenetic regulation in rice. First, we describe the RdDM mechanisms in plants. We then introduce recent discoveries on the biological roles of rice genes involved in the RdDM pathway and TE-derived sRNAs working at specific genomic loci for epigenetic control in rice.