Search

In monocots other than maize (Zea mays) and rice (Oryza sativa), the repertoire and diversity of microRNAs (miRNAs) and the populations of phased, secondary, small interfering RNAs (phasiRNAs) are poorly characterized. To remedy this, we sequenced small RNAs (sRNA) from vegetative and dissected inflorescence tissue in 28 phylogenetically diverse monocots and from several early-diverging angiosperm lineages, as well as publicly available data from 10 additional monocot species. We annotated miRNAs, small interfering RNAs (siRNAs) and phasiRNAs across the monocot phylogeny, identifying miRNAs apparently lost or gained in the grasses relative to other monocot families, as well as a number of transfer RNA fragments misannotated as miRNAs. Using our miRNA database cleaned of these misannotations, we identified conservation at the 8th, 9th, 19th, and 3′-end positions that we hypothesize are signatures of selection for processing, targeting, or Argonaute sorting. We show that 21-nucleotide (nt) reproductive phasiRNAs are far more numerous in grass genomes than other monocots. Based on sequenced monocot genomes and transcriptomes, DICER-LIKE5, important to 24-nt phasiRNA biogenesis, likely originated via gene duplication before the diversification of the grasses. This curated database of phylogenetically diverse monocot miRNAs, siRNAs, and phasiRNAs represents a large collection of data that should facilitate continued exploration of sRNA diversification in flowering plants.

Isolated or syndromic congenital cataracts are heterogeneous developmental defects, making the identification of the associated genes challenging. In the past, mouse lens expression microarrays have been successfully applied in bioinformatics tools (e.g., iSyTE) to facilitate human cataract-associated gene discovery. To develop a new resource for geneticists, we report high-throughput RNA sequencing (RNA-seq) profiles of mouse lens at key embryonic stages (E)10.5 (lens pit), E12.5 (primary fiber cell differentiation), E14.5 and E16.5 (secondary fiber cell differentiation). These stages capture important events as the lens develops from an invaginating placode into a transparent tissue. Previously, in silico whole-embryo body (WB)-subtraction-based “lens-enriched” expression has been effective in prioritizing cataract-linked genes. To apply an analogous approach, we generated new mouse WB RNA-seq datasets and show that in silico WB subtraction of lens RNA-seq datasets successfully identifies key genes based on lens-enriched expression. At ≥2 counts-per-million expression, ≥1.5 log2 fold-enrichment (p

Summary and keywordsLittle is known about the characteristics and function of reproductive phased, secondary, small interfering RNAs (phasiRNAs) in the Poaceae, despite the availability of significant genomic resources, experimental data, and a growing number of computational tools. We utilized machine-learning methods to identify sequence-based and structural features that distinguish phasiRNAs in rice and maize from other small RNAs (sRNAs).We developed Random Forest classifiers that can distinguish reproductive phasiRNAs from other sRNAs in complex sets of sequencing data, utilizing sequence-based (k-mers) and features describing position-specific sequence biases.The classification performance attained is >80% in accuracy, sensitivity, specificity, and positive predicted value. Feature selection identified important features in both ends of phasiRNAs. We demonstrated that phasiRNAs have strand specificity and position-specific nucleotide biases potentially influencing AGO sorting; we also predicted targets to infer functions of phasiRNAs, and computationally-assessed their sequence characteristics relative to other sRNAs.Our results demonstrate that machine-learning methods effectively identify phasiRNAs despite the lack of characteristic features typically present in precursor loci of other small RNAs, such as sequence conservation or structural motifs. The 5’-end features we identified provide insights into AGO-phasiRNA interactions; we describe a hypothetical model of competition for AGO loading between phasiRNAs of different nucleotide compositions.

In grasses, two pathways that generate diverse and numerous 21-nt (premeiotic) and 24-nt (meiotic) phased siRNAs are highly enriched in anthers, the male reproductive organs. These “phasiRNAs” are analogous to mammalian piRNAs, yet their functions and evolutionary origins remain largely unknown. The 24-nt meiotic phasiRNAs have only been described in grasses, wherein their biogenesis is dependent on a specialized Dicer (DCL5). To assess how evolution gave rise to this pathway, we examined reproductive phasiRNA pathways in nongrass monocots: garden asparagus, daylily, and lily. The common ancestors of these species diverged approximately 115–117 million years ago (MYA). We found that premeiotic 21-nt and meiotic 24-nt phasiRNAs were abundant in all three species and displayed spatial localization and temporal dynamics similar to grasses. The miR2275-triggered pathway was also present, yielding 24-nt reproductive phasiRNAs, and thus originated more than 117 MYA. In asparagus, unlike in grasses, these siRNAs are largely derived from inverted repeats (IRs); analyses in lily identified thousands of precursor loci, and many were also predicted to form foldback substrates for Dicer processing. Additionally, reproductive phasiRNAs were present in female reproductive organs and thus may function in both male and female germinal development. These data describe several distinct mechanisms of production for 24-nt meiotic phasiRNAs and provide new insights into the evolution of reproductive phasiRNA pathways in monocots.

Although successful in identifying new cataract-linked genes, the previous version of the database iSyTE (integrated Systems Tool for Eye gene discovery) was based on expression information on just three mouse lens stages and was functionally limited to visualization by only UCSC-Genome Browser tracks. To increase its efficacy, here we provide an enhanced iSyTE version 2.0 (URL: http://research.bioinformatics.udel.edu/iSyTE) based on well-curated, comprehensive genome-level lens expression data as a one-stop portal for the effective visualization and analysis of candidate genes in lens development and disease. iSyTE 2.0 includes all publicly available lens Affymetrix and Illumina microarray datasets representing a broad range of embryonic and postnatal stages from wild-type and specific gene-perturbation mouse mutants with eye defects. Further, we developed a new user-friendly web interface for direct access and cogent visualization of the curated expression data, which supports convenient searches and a range of downstream analyses. The utility of these new iSyTE 2.0 features is illustrated through examples of established genes associated with lens development and pathobiology, which serve as tutorials for its application by the end-user. iSyTE 2.0 will facilitate the prioritization of eye development and disease-linked candidate genes in studies involving transcriptomics or next-generation sequencing data, linkage analysis and GWAS approaches.

Root-knot nematodes (RKNs) are one of the most important biotic factors limiting crop productivity in many crop plants. The major RKN control strategies include development of resistant cultivars, application of nematicides and crop rotation, but each has its own limitations. In recent years, RNA interference (RNAi) has become a powerful approach for developing nematode resistance. The two housekeeping genes, splicing factor and integrase, of Meloidogyne incognita were targeted for engineering nematode resistance using a host-delivered RNAi (HD-RNAi) approach. Splicing factor and integrase genes are essential for nematode development as they are involved in RNA metabolism. Stable homozygous transgenic Arabidopsis lines expressing dsRNA for both genes were generated. In RNAi lines of splicing factor gene, the number of galls, females and egg masses was reduced by 71.4, 74.5 and 86.6%, respectively, as compared with the empty vector controls. Similarly, in RNAi lines of the integrase gene, the number of galls, females and egg masses was reduced up to 59.5, 66.8 and 63.4%, respectively, compared with the empty vector controls. Expression analysis revealed a reduction in mRNA abundance of both targeted genes in female nematodes feeding on transgenic plants expressing dsRNA constructs. The silencing of housekeeping genes in the nematodes through HD-RNAi significantly reduced root-knot nematode infectivity and suggests that they will be useful in developing RKN resistance in crop plants.

Plant small RNAs are ∼20 to 24 nucleotide noncoding RNAs that typically have repressive regulatory roles in gene expression, functioning at the transcriptional or post-transcriptional level. This influence on regulation of developmental and physiological processes has direct effects on phenotype. High-throughput sequencing technologies have enabled the sequencing of millions of small RNAs. Along with decreased sequencing costs, recent improvements in small RNA library construction have facilitated the ability to use minimal amounts of input RNA for analysis. This unit describes steps to isolate total RNA from limited amounts of plant tissue to construct small RNA libraries and perform small RNA data processing. © 2017 by John Wiley & Sons, Inc.

The Lemnaceae family comprises aquatic plants of angiosperms gaining attention due to their utility in wastewater treatment, and rapid production of biomass that can be used as feed, fuel, or food. Moreover, it can serve as a model species for neotenous growth and environmental adaptation. The latter properties are subject to post-transcriptional regulation of gene expression, meriting investigation of how miRNAs in Spirodela polyrhiza, the most basal and most thoroughly sequenced member of the family, are expressed under different growth conditions. To further scientific understanding of its capacity to adapt to environmental cues, we measured miRNA expression and processing of their target sequences under different temperatures, and in the presence of abscisic acid, copper, kinetin, nitrate, and sucrose. Using two small RNA sequencing experiments and one degradome sequencing experiment, we provide evidence for 108 miRNAs. Sequencing cleaved mRNAs validated 42 conserved miRNAs with 83 targets and 24 novel miRNAs regulating 66 targets and created a list of 575 predicted and verified targets. These analyses revealed condition-induced changes in miRNA expression and cleavage activity, and resulted in the addition of stringently reviewed miRNAs to miRBase. This combination of small RNA and degradome sequencing provided not only high confidence predictions of conserved and novel miRNAs and targets, but also a view of the post-transcriptional regulation of adaptations. A unique aspect is the role of miR156 and miR172 expression and activity in its clonal propagation and neoteny. Additionally, low levels of 24 nt sRNAs were observed, despite the lack of recent retrotransposition.

Opacification of the ocular lens, termed cataract, is a common cause of blindness. To become transparent, lens fiber cells undergo degradation of their organelles, including their nuclei, presenting a fundamental question: does signaling/transcription sufficiently explain differentiation of cells progressing toward compromised transcriptional potential? We report that a conserved RNA-binding protein Celf1 post-transcriptionally controls key genes to regulate lens fiber cell differentiation. Celf1-targeted knockout mice and celf1-knockdown zebrafish and Xenopus morphants have severe eye defects/cataract. Celf1 spatiotemporally down-regulates the cyclin-dependent kinase (Cdk) inhibitor p27Kip1 by interacting with its 5’ UTR and mediating translation inhibition. Celf1 deficiency causes ectopic up-regulation of p21Cip1. Further, Celf1 directly binds to the mRNA of the nuclease Dnase2b to maintain its high levels. Together these events are necessary for Cdk1-mediated lamin A/C phosphorylation to initiate nuclear envelope breakdown and DNA degradation in fiber cells. Moreover, Celf1 controls alternative splicing of the membrane-organization factor beta-spectrin and regulates F-actin-crosslinking factor Actn2 mRNA levels, thereby controlling fiber cell morphology. Thus, we illustrate new Celf1-regulated molecular mechanisms in lens development, suggesting that post-transcriptional regulatory RNA-binding proteins have evolved conserved functions to control vertebrate oculogenesis., Author summary Besides associated with aging, the eye disease cataract can occur early in life because of defects in lens development. Lens fiber cells degrade their nuclei to achieve lens transparency, which poses a fundamental question: how is differentiation regulated in a cell progressing toward compromised transcriptional potential? We demonstrate that this is achieved by distinct post-transcriptional gene expression control mechanisms mediated by a conserved RNA-binding protein Celf1. Celf1 regulates key factors required for normal nuclear degradation and cell morphology of fiber cells by controlling abundance of target mRNAs and/or their translation into protein. These data identify new Celf1-regulated molecular mechanisms necessary for lens transparency and suggest that conserved post-transcriptional regulatory RNA-binding proteins have evolved to control eye development in vertebrates.

The root-knot nematode (RKN), Meloidogyne incognita, is an obligate, sedentary endoparasite that infects a large number of crops and severely affects productivity. The commonly used nematode control strategies have their own limitations. Of late, RNA interference (RNAi) has become a popular approach for the development of nematode resistance in plants. Transgenic crops capable of expressing dsRNAs, specifically in roots for disrupting the parasitic process, offer an effective and efficient means of producing resistant crops. We identified nematode-responsive and root-specific (NRRS) promoters by using microarray data from the public domain and known conserved cis-elements. A set of 51 NRRS genes was identified which was narrowed down further on the basis of presence of cis-elements combined with minimal expression in the absence of nematode infection. The comparative analysis of promoters from the enriched NRRS set, along with earlier reported nematode-responsive genes, led to the identification of specific cis-elements. The promoters of two candidate genes were used to generate transgenic plants harboring promoter GUS constructs and tested in planta against nematodes. Both promoters showed preferential expression upon nematode infection, exclusively in the root in one and galls in the other. One of these NRRS promoters was used to drive the expression of splicing factor, a nematode-specific gene, for generating host-delivered RNAi-mediated nematode-resistant plants. Transgenic lines expressing dsRNA of splicing factor under the NRRS promoter exhibited upto a 32% reduction in number of galls compared to control plants.

Although majority of the genes linked to early-onset cataract exhibit lens fiber cell-enriched expression, our understanding of gene regulation in these cells is limited to function of just eight transcription factors and largely in the context of crystallins. We report on small Maf transcription factors Mafg and Mafk as regulators of several non-crystallin human cataract-associated genes in fiber cells and establish their significance to this disease. We applied a bioinformatics tool for cataract gene discovery iSyTE to identify Mafg and its co-regulators in the lens, and generated various null-allelic combinations of Mafg:Mafk mouse mutants for phenotypic and molecular analysis. By age 4 months, Mafg−/−:Mafk+/− mutants exhibit lens defects that progressively develop into cataract. High-resolution phenotypic characterization of Mafg−/−:Mafk+/− mouse lens reveals severely disorganized fiber cells, while microarray-based expression profiling identifies 97 differentially regulated genes (DRGs). Integrative analysis of Mafg−/−:Mafk+/− lens-DRGs with (1) binding motifs and genomic targets of small Mafs and their regulatory partners, (2) iSyTE lens expression data, and (3) interactions between DRGs in the String database, unravel a detailed small Maf regulatory network in the lens, several nodes of which are linked to cataract. This approach identifies 36 high-priority candidates from the original 97 DRGs. Significantly, 8/36 (22 %) DRGs are associated with cataracts in human (GSTO1, MGST1, SC4MOL, UCHL1) or mouse (Aldh3a1, Crygf, Hspb1, Pcbd1), suggesting a multifactorial etiology that includes oxidative stress and misregulation of sterol synthesis. These data identify Mafg and Mafk as new cataract-associated candidates and define their function in regulating largely non-crystallin genes linked to human cataract.

Small RNAs are ubiquitous, versatile repressors and include (1) microRNAs (miRNAs), processed from mRNA forming stem-loops; and (2) small interfering RNAs (siRNAs), the latter derived in plants by a process typically requiring an RNA-dependent RNA polymerase. We constructed and analyzed an expression atlas of soybean (Glycine max) small RNAs, identifying over 500 loci generating 21-nucleotide phased siRNAs (phasiRNAs; from PHAS loci), of which 483 overlapped annotated protein-coding genes. Via the integration of miRNAs with parallel analysis of RNA end (PARE) data, 20 miRNA triggers of 127 PHAS loci were detected. The primary class of PHAS loci (208 or 41% of the total) corresponded to NB-LRR genes; some of these small RNAs preferentially accumulate in nodules. Among the PHAS loci, novel representatives of TAS3 and noncanonical phasing patterns were also observed. A noncoding PHAS locus, triggered by miR4392, accumulated preferentially in anthers; the phasiRNAs are predicted to target transposable elements, with their peak abundance during soybean reproductive development. Thus, phasiRNAs show tremendous diversity in dicots. We identified novel miRNAs and assessed the veracity of soybean miRNAs registered in miRBase, substantially improving the soybean miRNA annotation, facilitating an improvement of miRBase annotations and identifying at high stringency novel miRNAs and their targets.

Phased, secondary siRNAs (phasiRNAs) are found widely in plants, from protein-coding transcripts and long, non-coding RNAs; animal piRNAs are also phased. Integrated methods characterizing “PHAS” loci are unavailable, and existing methods are quite limited and inefficient in handling large volumes of sequencing data. The PHASIS suite described here provides complete tools for the computational characterization of PHAS loci, with an emphasis on plants, in which these loci are numerous. Benchmarked comparisons demonstrate that PHASIS is sensitive, highly scalable and fast. Importantly, PHASIS eliminates the requirement of a sequenced genome and PARE/degradome data for discovery of phasiRNAs and their miRNA triggers.

Summary: We describe miTRATA, the first web-based tool for microRNA Truncation and Tailing Analysis—the analysis of 3′ modifications of microRNAs including the loss or gain of nucleotides relative to the canonical sequence. miTRATA is implemented in Python (version 3) and employs parallel processing modules to enhance its scalability when analyzing multiple small RNA (sRNA) sequencing datasets. It utilizes miRBase, currently version 21, as a source of known microRNAs for analysis. miTRATA notifies user(s) via email to download as well as visualize the results online. miTRATA’s strengths lie in (i) its biologist-focused web interface, (ii) improved scalability via parallel processing and (iii) its uniqueness as a webtool to perform microRNA truncation and tailing analysis. Availability and implementation: miTRATA is developed in Python and PHP. It is available as a web-based application from https://wasabi.dbi.udel.edu/∼apps/ta/. Contact: meyers@dbi.udel.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Sex chromosomes evolved from autosomes many times across the eukaryote phylogeny. Several models have been proposed to explain this transition, some involving male and female sterility mutations linked in a region of suppressed recombination between X and Y (or Z/W, U/V) chromosomes. Comparative and experimental analysis of a reference genome assembly for a double haploid YY male garden asparagus (Asparagus officinalis L.) individual implicates separate but linked genes as responsible for sex determination. Dioecy has evolved recently within Asparagus and sex chromosomes are cytogenetically identical with the Y, harboring a megabase segment that is missing from the X. We show that deletion of this entire region results in a male-to-female conversion, whereas loss of a single suppressor of female development drives male-to-hermaphrodite conversion. A single copy anther-specific gene with a male sterile Arabidopsis knockout phenotype is also in the Y-specific region, supporting a two-gene model for sex chromosome evolution.

We study the molecular mechanisms that are essential for development and maintenance of transparency in the eye lens. Toward this goal, we have applied the bioinformatics-based gene discovery tool ...

Plant architecture is determined by meristems that initiate leaves during vegetative development and flowers during reproductive development. Maize (Zea mays) inflorescences are patterned by a series of branching events, culminating in floral meristems that produce sexual organs. The maize fuzzy tassel (fzt) mutant has striking inflorescence defects with indeterminate meristems, fasciation, and alterations in sex determination. fzt plants have dramatically reduced plant height and shorter, narrower leaves with leaf polarity and phase change defects. We positionally cloned fzt and discovered that it contains a mutation in a dicer-like1 homolog, a key enzyme required for microRNA (miRNA) biogenesis. miRNAs are small noncoding RNAs that reduce target mRNA levels and are key regulators of plant development and physiology. Small RNA sequencing analysis showed that most miRNAs are moderately reduced in fzt plants and a few miRNAs are dramatically reduced. Some aspects of the fzt phenotype can be explained by reduced levels of known miRNAs, including miRNAs that influence meristem determinacy, phase change, and leaf polarity. miRNAs responsible for other aspects of the fzt phenotype are unknown and likely to be those miRNAs most severely reduced in fzt mutants. The fzt mutation provides a tool to link specific miRNAs and targets to discrete phenotypes and developmental roles.

MicroRNAs are ubiquitous in plant genomes but vary greatly in their abundance within and conservation among plant lineages. To gain insight into the evolutionary birth/death dynamics of microRNA families, we sequenced small RNA and 5′-end PARE libraries generated from two closely related species of Gossypium. Here, we demonstrate that 33 microRNA families, with similar copy numbers and average evolutionary rates, are conserved in the two congeneric cottons. Analysis of the presence/absence of these microRNA families in other land plants sheds light on their depth of phylogenetic origin and lineage-specific loss/gain. Conserved microRNA families in Gossypium exhibit a striking interspecific asymmetry in expression, potentially connected to relative proximity to neighboring transposable elements. A complex correlated expression pattern of microRNA target genes with their controlling microRNAs indicates that possible functional divergence of conserved microRNA families can also exist even within a single plant genus.

Parallel analysis of RNA ends (PARE) is a technique utilizing high-throughput sequencing to profile uncapped, mRNA cleavage or decay products on a genome-wide basis. Tools currently available to validate miRNA targets using PARE data employ only annotated genes, whereas important targets may be found in unannotated genomic regions. To handle such cases and to scale to the growing availability of PARE data and genomes, we developed a new tool, 'sPARTA' (small RNA-PARE target analyzer) that utilizes a built-in, plant-focused target prediction module (aka 'miRferno'). sPARTA not only exhibits an unprecedented gain in speed but also it shows greater predictive power by validating more targets, compared to a popular alternative. In addition, the novel 'seed-free' mode, optimized to find targets irrespective of complementarity in the seed-region, identifies novel intergenic targets. To fully capitalize on the novelty and strengths of sPARTA, we developed a web resource, 'comPARE', for plant miRNA target analysis; this facilitates the systematic identification and analysis of miRNA-target interactions across multiple species, integrated with visualization tools. This collation of high-throughput small RNA and PARE datasets from different genomes further facilitates re-evaluation of existing miRNA annotations, resulting in a 'cleaner' set of microRNAs.