Genome-wide identification of AGO18b-bound miRNAs and phasiRNAs in maize by cRIP-seq

MicroRNAs in Plants Development and Stress Resistance

Plant growth and development are governed by a rigorously timed sequence of ontogenetic programmes. MicroRNAs (miRNAs), a class of short noncoding RNAs, function as master regulators of gene expression by targeting mRNAs for cleavage or direct translational inhibition at the posttranscriptional level in eukaryotes. Numerous miRNA molecules that control significant agronomic properties in plants have been found. On the one hand, miRNAs target transcription factors (TFs) to determine plant structure, such as root development, internode elongation, leaf morphogenesis, sex determination and nutrient transition. On the other hand, miRNAs alter expression levels to adapt to biological and abiotic stresses, including fungi, bacteria, viruses, drought, waterlogging, high temperature, low temperature, salinity, nutrient deficiencies, heavy metals and other abiotic stresses. To fully understand the role of miRNAs in plants, we review the regulatory role of miRNAs in plant development and stress resistance. Beyond that, we propose that the novel miRNA in review can be effectively further studied with artificial miRNA (amiRNA) or short tandem target mimics (STTM) and miRNA delivery in vitro can be used to improve crop yield and agricultural sustainability.

The role of Ancestral MicroRNAs in grass inflorescence development

Plant inflorescences are complex, highly diverse structures whose morphology is determined in meristems that form during reproductive development. Inflorescence structure influences flower formation, and consequently grain number, and yield in crops. Correct inflorescence and flower development require tight control of gene expression via complex interplay between regulatory networks. MicroRNAs (miRNAs) have emerged as fundamental modulators of gene expression at the transcriptional and/or post-transcriptional level in plant inflorescence development. First discovered more than three decades ago, miRNAs have proved to be revolutionary in advancing our mechanistic understanding of gene expression. This review highlights current knowledge of downstream target genes and pathways of some highly conserved miRNAs that regulate the maintenance, identity, and activity of inflorescence and floral meristems in economically and agriculturally important grass species, including rice (Oryza sativa), maize (Zea mays), barley (Hordeum vulgare), and wheat (Triticum aestivum). Furthermore, we summarize emerging regulatory networks of miRNAs and their targets to suggest new avenues and strategies for application of miRNAs as a tool to enhance crop yield and performance.

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.

Genome-wide in-silico analysis of ethylene biosynthesis gene family in Musa acuminata L. and their response under nutrient stress

Ethylene is a gaseous phytohormone involved in plants' growth and developmental processes, including seed germination, root initiation, fruit ripening, flower and leaf senescence, abscission, and stress responses. Ethylene biosynthesis (EB) gene analysis in response to nitrogen (N) and potassium (K) stress has not yet been conducted in Musa acuminata (banana) roots. The genome mining of banana (Musa acuminata L.) revealed 14 putative 1-aminocyclopropane-1-carboxylate synthase (ACS), 10 1-aminocyclopropane-1-carboxylate oxidase (ACO), and 3 Ethylene overproducer 1 (ETO1) genes. ACS, ACO, and ETO1 proteins possessed amino acid residues ranging from 422-684, 636-2670, and 893-969, respectively, with molecular weight (Mw) ranging from 4.93-7.55 kD, 10.1-8.3 kD and 10.1-10.78 kD. The number of introns present in ACS, ACO, and ETO1 gene sequences ranges from 0-14, 1-6, and 0-6, respectively. The cis-regulatory element analysis revealed the presence of light-responsive, abscisic acid, seed regulation, auxin-responsive, gibberellin element, endosperm-specific, anoxic inducibility, low-temperature responsiveness, salicylic acid responsiveness, meristem-specific and stress-responsive elements. Comprehensive phylogenetic analyses ACS, ACO, and ETO1 genes of Banana with Arabidopsis thaliana revealed several orthologs and paralogs assisting in understanding the putative functions of these genes. The expression profile of Musa acuminata genes in root under normal and low levels of nitrogen and potassium shows that MaACS14 and MaACO6 expressed highly at normal nitrogen supply. MaACS1 expression was significantly upregulated at low potassium levels, whereas, MaACO6 gene expression was significantly downregulated. The functional divergence and site-specific selective pressures on specific gene sequences of banana have been investigated. The bioinformatics-based genome-wide assessment of the family of banana attempted in the present study could be a significant step for deciphering novel ACS, ACO, and ETO1 genes based on genome-wide expression profiling.

Beyond Loading: Functions of Plant ARGONAUTE Proteins

ARGONAUTE (AGO) proteins are key components of the RNA-induced silencing complex (RISC) that mediates gene silencing in eukaryotes. Small-RNA (sRNA) cargoes are selectively loaded into different members of the AGO protein family and then target complementary sequences to in-duce transcriptional repression, mRNA cleavage, or translation inhibition. Previous reviews have mainly focused on the traditional roles of AGOs in specific biological processes or on the molecular mechanisms of sRNA sorting. In this review, we summarize the biological significance of canonical sRNA loading, including the balance among distinct sRNA pathways, cross-regulation of different RISC activities during plant development and defense, and, especially, the emerging roles of AGOs in sRNA movement. We also discuss recent advances in novel non-canonical functions of plant AGOs. Perspectives for future functional studies of this evolutionarily conserved eukaryotic protein family will facilitate a more comprehensive understanding of the multi-faceted AGO proteins.

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.

Genome-Wide Association Study and Genomic Prediction on Plant Architecture Traits in Sweet Corn and Waxy Corn

Sweet corn and waxy corn has a better taste and higher accumulated nutritional value than regular maize, and is widely planted and popularly consumed throughout the world. Plant height (PH), ear height (EH), and tassel branch number (TBN) are key plant architecture traits, which play an important role in improving grain yield in maize. In this study, a genome-wide association study (GWAS) and genomic prediction analysis were conducted on plant architecture traits of PH, EH, and TBN in a fresh edible maize population consisting of 190 sweet corn inbred lines and 287 waxy corn inbred lines. Phenotypic data from two locations showed high heritability for all three traits, with significant differences observed between sweet corn and waxy corn for both PH and EH. The differences between the three subgroups of sweet corn were not obvious for all three traits. Population structure and PCA analysis results divided the whole population into three subgroups, i.e., sweet corn, waxy corn, and the subgroup mixed with sweet and waxy corn. Analysis of GWAS was conducted with 278,592 SNPs obtained from resequencing data; 184, 45, and 68 significantly associated SNPs were detected for PH, EH, and TBN, respectively. The phenotypic variance explained (PVE) values of these significant SNPs ranged from 3.50% to 7.0%. The results of this study lay the foundation for further understanding the genetic basis of plant architecture traits in sweet corn and waxy corn. Genomic selection (GS) is a new approach for improving quantitative traits in large plant breeding populations that uses whole-genome molecular markers. The marker number and marker quality are essential for the application of GS in maize breeding. GWAS can choose the most related markers with the traits, so it can be used to improve the predictive accuracy of GS.

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.

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.

Transgenerational conditioned male fertility of HD-ZIP IV transcription factor mutant ocl4: impact on 21-nt phasiRNA accumulation in pre-meiotic maize anthers

Maize Outer cell layer 4 (ocl4) encodes an HD-ZIP IV transcription factor required for robust male fertility and 21-nt phasiRNA biogenesis. ocl4 fertility is favored in warm conditions, and phasiRNAs are partially restored. Environment-sensitive male-sterile plants have been described before and can result from different molecular mechanisms and biological processes, but putative environment-conditioned, transgenerational rescue of their male fertility is a rather new mystery. Here, we report a derivative line of the male-sterile outer cell layer 4 (ocl4) mutant of maize, in which fertility was restored and perpetuated over several generations. Conditioned fertile ocl4 anthers exhibit the anatomical abnormality of a partially duplicated endothecial layer, just like their sterile counterparts. We profiled the dynamics of phased, small interfering RNAs (phasiRNAs) during pre-meiotic development in fully sterile and various grades of semi-fertile ocl4 anthers. The conditioned fertile anthers accumulated significantly higher 21-nt phasiRNAs compared to ocl4 sterile samples, suggesting a partial restoration of phasiRNAs in conditioned fertility. We found that the biogenesis of 21-nt phasiRNAs is largely dependent on Ocl4 at three key steps: (1) production of PHAS precursor transcripts, (2) expression of miR2118 that modulates precursor processing, and (3) accumulation of 21-nt phasiRNAs.

Oxygen, secreted proteins and small RNAs: mobile elements that govern anther development

Correct anther development is essential for male fertility and subsequently agricultural yield. Defects in anther development range from the early stage of stamen formation until the late stage of tapetum degeneration. In particular, the specification of the four distinct somatic layers and the inner sporogenous cells need perfect orchestration relying on precise cell-cell communication. Up to now, several signals, which coordinate the anther´s developmental program, have been identified. Among the known signals are phytohormones, environmental conditions sensed via glutaredoxins, several receptor-like kinases triggered by ligands like MAC1, and small RNAs such as miRNAs and the monocot-prevalent reproductive phasiRNAs. Rather than giving a full review on anther development, here we discuss anther development with an emphasis on mobile elements like ROS/oxygen, secreted proteins and small RNAs (only briefly touching on phytohormones), how they might act and interact, and what the future of this research area might reveal.

Rice-specific Argonaute 17 controls reproductive growth and yield-associated phenotypes

This manuscript describes the functions of an Argonaute protein named AGO17 in rice. AGO17 is required for the development of rice reproductive tissues. Argonaute (AGO) proteins are a well-conserved multigene family of regulators mediating gene silencing across eukaryotes. Monocot plants have additional members of AGO, the functions of which are poorly understood. Among the non-dicot AGO1 clade members in monocots, AGO17 expresses highly in reproductive tissues. Here we show that overexpression of Oryza sativa indica AGO17 in rice resulted in robust growth and increased yield, whereas its silencing resulted in reduced panicle length, less fertility, and poor growth. Small (s)RNA transcriptome analysis revealed misregulation of several miRNAs and other categories of sRNAs in silenced and overexpression lines, in agreement with its likely competition with other AGO1 clade members. Targets of differentially expressed miRNAs included previously unreported target RNAs coding for proteins involved in development, phase transition, and transport. Our results indicate a distinctive role for OsAGO17 in rice reproductive development that could be harnessed to improve yield.

Quantitative, super-resolution localization of small RNAs with sRNA-PAINT

Small RNAs are non-coding RNAs that play important roles in the lives of both animals and plants. They are 21- to 24-nt in length and ∼10 nm in size. Their small size and high diversity have made it challenging to develop detection methods that have sufficient resolution and specificity to multiplex and quantify. We created a method, sRNA-PAINT, for the detection of small RNAs with 20 nm resolution by combining the super-resolution method, DNA-based points accumulation in nanoscale topography (DNA-PAINT), and the specificity of locked nucleic acid (LNA) probes for the in situ detection of multiple small RNAs. The method relies on designing probes to target small RNAs that combine DNA oligonucleotides (oligos) for PAINT with LNA-containing oligos for hybridization; therefore, we developed an online tool called 'Vetting & Analysis of RNA for in situ Hybridization probes' (VARNISH) for probe design. Our method utilizes advances in DNA-PAINT methodologies, including qPAINT for quantification, and Exchange-PAINT for multiplexing. We demonstrated these capabilities of sRNA-PAINT by detecting and quantifying small RNAs in different cell layers of early developmental stage maize anthers that are important for male sexual reproduction.

Loss of function of Oryza sativa Argonaute 18 induces male sterility and reduction in phased small RNAs

In this manuscript, we show that Oryza sativa indica Argonaute protein AGO18 is required for male gametophyte development likely to through a small RNA-mediated mechanism. Monocots have evolved unique gene silencing pathways due to the presence of unique members of Dicer-like and Argonaute (AGO) family members. Among the monocot AGO homologs, AGO18 occupies a unique position. Previous reports have implicated this protein in viral resistance as well as in gametogenesis, likely through its competition with AGO1 clade members for micro(mi)RNAs and other small (s)RNAs. Although expression of rice AGO18 in specific stages of male gametogenesis has been documented, its major functions in plant development remain poorly understood. Here, we show that Oryza sativa indica AGO18 is involved in male gametophyte development. Knockdown (KD) of AGO18 in transgenic rice lines resulted in stunted plants that are male sterile, whereas their carpels were functional. Transcriptome analysis revealed downregulation of several pollen development-associated genes in KD lines. sRNA sequencing in vegetative and reproductive tissues of KD lines indicated reduction of miRNAs and phased secondary sRNAs implicated in male gametophyte development. Our results indicate a distinct role for rice AGO18 in male fertility.

Tapetum-Dependent Male Meiosis Progression in Plants: Increasing Evidence Emerges

In higher plants, male meiosis is a key process during microsporogenesis and is crucial for male fertility and seed set. Meiosis involves a highly dynamic organization of chromosomes and cytoskeleton and specifically takes place within sexual cells. However, studies in multiple plant species have suggested that the normal development of tapetum, the somatic cell layer surrounding the developing male meiocytes, is indispensable for the completion of the male meiotic cell cycle. Disrupted tapetum development causes alterations in the expression of a large range of genes involved in male reproduction. Moreover, recent experiments suggest that small RNAs (sRNAs) present in the anthers, including microRNAs (miRNAs) and phased, secondary, small interfering RNAs (phasiRNAs), play a potential but important role in controlling male meiosis, either by influencing the expression of meiotic genes in the meiocytes or through other unclear mechanisms, supporting the hypothesis that male meiosis is non-cell autonomously regulated. In this mini review, we summarize the recorded meiotic defects that occur in plants with defective tapetum development in both Arabidopsis and crops. Thereafter, we outline the latest understanding on the molecular mechanisms that potentially underpin the tapetum-dependent regulation of male meiosis, and we especially discuss the regulatory role of sRNAs. At the end, we propose several outstanding questions that should be addressed in future studies.