ShortStack: comprehensive annotation and quantification of small RNA genes

A low-bias and sensitive small RNA library preparation method using randomized splint ligation

Small RNAs are important regulators of gene expression and are involved in human development and disease. Next generation sequencing (NGS) allows for scalable, genome-wide studies of small RNA; however, current methods are challenged by low sensitivity and high bias, limiting their ability to capture an accurate representation of the cellular small RNA population. Several studies have shown that this bias primarily arises during the ligation of single-strand adapters during library preparation, and that this ligation bias is magnified by 2′-O-methyl modifications (2′OMe) on the 3′ terminal nucleotide. In this study, we developed a novel library preparation process using randomized splint ligation with a cleavable adapter, a design which resolves previous challenges associated with this ligation strategy. We show that a randomized splint ligation based workflow can reduce bias and increase the sensitivity of small RNA sequencing for a wide variety of small RNAs, including microRNA (miRNA) and tRNA fragments as well as 2′OMe modified RNA, including Piwi-interacting RNA and plant miRNA. Finally, we demonstrate that this workflow detects more differentially expressed miRNA between tumorous and matched normal tissues. Overall, this library preparation process allows for highly accurate small RNA sequencing and will enable studies of 2′OMe modified RNA with new levels of detail.

Systematic analysis of long intergenic non-coding RNAs in C. elegans germline uncovers roles in somatic growth

Long intergenic non-coding RNAs (lincRNAs) are transcripts longer than 200 nucleotides that are transcribed from non-coding loci yet undergo biosynthesis similar to coding mRNAs. The disproportional number of lincRNAs expressed in testes suggests that lincRNAs are important during gametogenesis, but experimental evidence has implicated very few lincRNAs in this process. We took advantage of the relatively limited number of lincRNAs in the genome of the nematode Caenorhabditis elegans to systematically analyse the functions of lincRNAs during meiosis. We deleted six lincRNA genes that are highly and dynamically expressed in the C. elegans gonad and tested the effects on central meiotic processes. Surprisingly, whereas the lincRNA deletions did not strongly impact fertility, germline apoptosis, crossovers, or synapsis, linc-4 was required for somatic growth. Slower growth was observed in linc-4-deletion mutants and in worms depleted of linc-4 using RNAi, indicating that linc-4 transcripts are required for this post-embryonic process. Unexpectedly, analysis of worms depleted of linc-4 in soma versus germline showed that the somatic role stems from linc-4 expression in germline cells. This unique feature suggests that some lincRNAs, like some small non-coding RNAs, are required for germ-soma interactions.

Identification of small non-coding RNAs responsive to GUN1 and GUN5 related retrograde signals in Arabidopsis thaliana

Chloroplast perturbations activate retrograde signalling pathways, causing dynamic changes of gene expression. Besides transcriptional control of gene expression, different classes of small non-coding RNAs (sRNAs) act in gene expression control, but comprehensive analyses regarding their role in retrograde signalling are lacking. We performed sRNA profiling in response to norflurazon (NF), which provokes retrograde signals, in Arabidopsis thaliana wild type (WT) and the two retrograde signalling mutants gun1 and gun5. The RNA samples were also used for mRNA and long non-coding RNA profiling to link altered sRNA levels to changes in the expression of their cognate target RNAs. We identified 122 sRNAs from all known sRNA classes that were responsive to NF in the WT. Strikingly, 142 and 213 sRNAs were found to be differentially regulated in both mutants, indicating a retrograde control of these sRNAs. Concomitant with the changes in sRNA expression, we detected about 1500 differentially expressed mRNAs in the NF-treated WT and around 900 and 1400 mRNAs that were differentially regulated in the gun1 and gun5 mutants, with a high proportion (~30%) of genes encoding plastid proteins. Furthermore, around 20% of predicted miRNA targets code for plastid-localised proteins. Among the sRNA-target pairs, we identified pairs with an anticorrelated expression as well pairs showing other expressional relations, pointing to a role of sRNAs in balancing transcriptional changes upon retrograde signals. Based on the comprehensive changes in sRNA expression, we assume a considerable impact of sRNAs in retrograde-dependent transcriptional changes to adjust plastidic and nuclear gene expression.

Long-Read cDNA Sequencing Enables a "Gene-Like" Transcript Annotation of Transposable Elements

Transcript-based annotations of genes facilitate both genome-wide analyses and detailed single-locus research. In contrast, transposable element (TE) annotations are rudimentary, consisting of information only on TE location and type. The repetitiveness and limited annotation of TEs prevent the ability to distinguish between potentially functional expressed elements and degraded copies. To improve genome-wide TE bioinformatics, we performed long-read sequencing of cDNAs from Arabidopsis (Arabidopsis thaliana) lines deficient in multiple layers of TE repression. These uniquely mapping transcripts were used to identify the set of TEs able to generate polyadenylated RNAs and create a new transcript-based annotation of TEs that we have layered upon the existing high-quality community standard annotation. We used this annotation to reduce the bioinformatic complexity associated with multimapping reads from short-read RNA sequencing experiments, and we show that this improvement is expanded in a TE-rich genome such as maize (Zea mays). Our TE annotation also enables the testing of specific standing hypotheses in the TE field. We demonstrate that inaccurate TE splicing does not trigger small RNA production, and the cell more strongly targets DNA methylation to TEs that have the potential to make mRNAs. This work provides a transcript-based TE annotation for Arabidopsis and maize, which serves as a blueprint to reduce the bioinformatic complexity associated with repetitive TEs in any organism.

Spatiotemporal Resolved Leaf Angle Establishment Improves Rice Grain Yield via Controlling Population Density

Leaf angle is mainly determined by the lamina joint (LJ) and contributes to ideal crop architecture for high yield. Here, we dissected five successive stages with distinct cytological features of LJs spanning organogenesis to leaf angle formation and obtained the underlying stage-specific mRNAs and small RNAs, which well explained the cytological dynamics during LJ organogenesis and leaf angle plasticity. Combining the gene coexpression correlation with high-throughput promoter analysis, we identified a set of transcription factors (TFs) determining the stage- and/or cytological structure-specific profiles. The functional studies of these TFs demonstrated that cytological dynamics determined leaf angle and that the knockout rice of these TFs with erect leaves significantly enhanced yield by maintaining the proper tiller number under dense planting. This work revealed the high-resolution mechanisms of how the cytological dynamics of LJ determined leaf erectness and served as a valuable resource to remodel rice architecture for high yield by controlling population density.

Extracellular vesicles from Heligmosomoides bakeri and Trichuris muris contain distinct microRNA families and small RNAs that could underpin different functions in the host

Extracellular vesicles (EVs) have emerged as a ubiquitous component of helminth excretory-secretory products that can deliver parasite molecules to host cells to elicit immunomodulatory effects. RNAs are one type of cargo molecule that can underpin EV functions, hence there is extensive interest in characterising the RNAs that are present in EVs from different helminth species. Here we outline methods for identifying all of the small RNAs (sRNA) in helminth EVs and address how different methodologies may influence the sRNAs detected. We show that different EV purification methods introduce relatively little variation in the sRNAs that are detected, and that different RNA library preparation methods yielded larger differences. We compared the EV sRNAs in the gastrointestinal nematode Heligmosomoides bakeri with those in EVs from the distantly related gastrointestinal nematode Trichuris muris, and found that many of the sRNAs in both organisms derive from repetitive elements or intergenic regions. However, only in H. bakeri do these RNAs contain a 5' triphosphate, and Guanine (G) starting nucleotide, consistent with their biogenesis by RNA-dependent RNA polymerases (RdRPs). Distinct microRNA (miRNA) families are carried in EVs from each parasite, with H. bakeri EVs specific for miR-71, miR-49, miR-63, miR-259 and miR-240 gene families, and T. muris EVs specific for miR-1, miR-1822 and miR-252, and enriched for miR-59, miR-72 and miR-44 families, with the miR-9, miR-10, miR-80 and let-7 families abundant in both. We found a larger proportion of miRNA reads derive from the mouse host in T. muris EVs, compared with H. bakeri EVs. Our report underscores potential biases in the sRNAs sequenced based on library preparation methods, suggests specific nematode lineages have evolved distinct sRNA synthesis/export pathways, and highlights specific differences in EV miRNAs from H. bakeri and T. muris that may underpin functional adaptation to their host niches.

Pan-Genome of Wild and Cultivated Soybeans

Soybean is one of the most important vegetable oil and protein feed crops. To capture the entire genomic diversity, it is needed to construct a complete high-quality pan-genome from diverse soybean accessions. In this study, we performed individual de novo genome assemblies for 26 representative soybeans that were selected from 2,898 deeply sequenced accessions. Using these assembled genomes together with three previously reported genomes, we constructed a graph-based genome and performed pan-genome analysis, which identified numerous genetic variations that cannot be detected by direct mapping of short sequence reads onto a single reference genome. The structural variations from the 2,898 accessions that were genotyped based on the graph-based genome and the RNA sequencing (RNA-seq) data from the representative 26 accessions helped to link genetic variations to candidate genes that are responsible for important traits. This pan-genome resource will promote evolutionary and functional genomics studies in soybean.

Identification of small RNAs during cold acclimation in Arabidopsis thaliana

BACKGROUND

Cold stress causes dynamic changes in gene expression that are partially caused by small non-coding RNAs since they regulate protein coding transcripts and act in epigenetic gene silencing pathways. Thus, a detailed analysis of transcriptional changes of small RNAs (sRNAs) belonging to all known sRNA classes such as microRNAs (miRNA) and small interfering RNA (siRNAs) in response to cold contributes to an understanding of cold-related transcriptome changes.

RESULT

We subjected A. thaliana plants to cold acclimation conditions (4 °C) and analyzed the sRNA transcriptomes after 3 h, 6 h and 2 d. We found 93 cold responsive differentially expressed miRNAs and only 14 of these were previously shown to be cold responsive. We performed miRNA target prediction for all differentially expressed miRNAs and a GO analysis revealed the overrepresentation of miRNA-targeted transcripts that code for proteins acting in transcriptional regulation. We also identified a large number of differentially expressed cis- and trans-nat-siRNAs, as well as sRNAs that are derived from long non-coding RNAs. By combining the results of sRNA and mRNA profiling with miRNA target predictions and publicly available information on transcription factors, we reconstructed a cold-specific, miRNA and transcription factor dependent gene regulatory network. We verified the validity of links in the network by testing its ability to predict target gene expression under cold acclimation.

CONCLUSION

In A. thaliana, miRNAs and sRNAs derived from cis- and trans-NAT gene pairs and sRNAs derived from lncRNAs play an important role in regulating gene expression in cold acclimation conditions. This study provides a fundamental database to deepen our knowledge and understanding of regulatory networks in cold acclimation.

Abundant expression of maternal siRNAs is a conserved feature of seed development

Small RNAs are abundant in plant reproductive tissues, especially 24-nucleotide (nt) small interfering RNAs (siRNAs). Most 24-nt siRNAs are dependent on RNA Pol IV and RNA-DEPENDENT RNA POLYMERASE 2 (RDR2) and establish DNA methylation at thousands of genomic loci in a process called RNA-directed DNA methylation (RdDM). In Brassica rapa, RdDM is required in the maternal sporophyte for successful seed development. Here, we demonstrate that a small number of siRNA loci account for over 90% of siRNA expression during B. rapa seed development. These loci exhibit unique characteristics with regard to their copy number and association with genomic features, but they resemble canonical 24-nt siRNA loci in their dependence on RNA Pol IV/RDR2 and role in RdDM. These loci are expressed in ovules before fertilization and in the seed coat, embryo, and endosperm following fertilization. We observed a similar pattern of 24-nt siRNA expression in diverse angiosperms despite rapid sequence evolution at siren loci. In the endosperm, siren siRNAs show a marked maternal bias, and siren expression in maternal sporophytic tissues is required for siren siRNA accumulation. Together, these results demonstrate that seed development occurs under the influence of abundant maternal siRNAs that might be transported to, and function in, filial tissues.

Corrigendum: Novel Meiotic miRNAs and Indications for a Role of PhasiRNAs in Meiosis

[This corrects the article DOI: 10.3389/fpls.2016.00762.].

sRNAbench and sRNAtoolbox 2019: intuitive fast small RNA profiling and differential expression

Since the original publication of sRNAtoolbox in 2015, small RNA research experienced notable advances in different directions. New protocols for small RNA sequencing have become available to address important issues such as adapter ligation bias, PCR amplification artefacts or to include internal controls such as spike-in sequences. New microRNA reference databases were developed with different foci, either prioritizing accuracy (low number of false positives) or completeness (low number of false negatives). Additionally, other small RNA molecules as well as microRNA sequence and length variants (isomiRs) have continued to gain importance. Finally, the number of microRNA sequencing studies deposited in GEO nearly triplicated from 2014 (280) to 2018 (764). These developments imply that fast and easy-to-use tools for expression profiling and subsequent downstream analysis of miRNA-seq data are essential to many researchers. Key features in this sRNAtoolbox release include addition of all major RNA library preparation protocols to sRNAbench and improvements in sRNAde, a tool that summarizes several aspects of small RNA sequencing studies including the detection of consensus differential expression. A special emphasis was put on the user-friendliness of the tools, for instance sRNAbench now supports parallel launching of several jobs to improve reproducibility and user time efficiency.

LuluDB-The Database Created Based on Small RNA, Transcriptome, and Degradome Sequencing Shows the Wide Landscape of Non-coding and Coding RNA in Yellow Lupine (Lupinus luteus L.) Flowers and Pods

Yellow lupine (Lupinus luteus L.) belongs to a legume family that benefits from symbiosis with nitrogen-fixing bacteria. Its seeds are rich in protein, which makes it a valuable food source for animals and humans. Yellow lupine is also the model plant for basic research on nodulation or abscission of organs. Nevertheless, the knowledge about the molecular regulatory mechanisms of its generative development is still incomplete. The RNA-Seq technique is becoming more prominent in high-throughput identification and expression profiling of both coding and non-coding RNA sequences. However, the huge amount of data generated with this method may discourage other scientific groups from making full use of them. To overcome this inconvenience, we have created a database containing analysis-ready information about non-coding and coding L. luteus RNA sequences (LuluDB). LuluDB was created on the basis of RNA-Seq analysis of small RNA, transcriptome, and degradome libraries obtained from yellow lupine cv. Taper flowers, pod walls, and seeds in various stages of development, flower pedicels, and pods undergoing abscission or maintained on the plant. It contains sequences of miRNAs and phased siRNAs identified in L. luteus, information about their expression in individual samples, and their target sequences. LuluDB also contains identified lncRNAs and protein-coding RNA sequences with their organ expression and annotations to widely used databases like GO, KEGG, NCBI, Rfam, Pfam, etc. The database also provides sequence homology search by BLAST using, e.g., an unknown sequence as a query. To present the full capabilities offered by our database, we performed a case study concerning transcripts annotated as DCL 1–4 (DICER LIKE 1–4) homologs involved in small non-coding RNA biogenesis and identified miRNAs that most likely regulate DCL1 and DCL2 expression in yellow lupine. LuluDB is available at http://luluseqdb.umk.pl/basic/web/index.php.

Plant 22-nt siRNAs mediate translational repression and stress adaptation

Small interfering RNAs (siRNAs) are essential for proper development and immunity in eukaryotes¹. Plants produce siRNAs with lengths of 21, 22 or 24 nucleotides. The 21- and 24-nucleotide species mediate cleavage of messenger RNAs and DNA methylation^2,3, respectively, but the biological functions of the 22-nucleotide siRNAs remain unknown. Here we report the identification and characterization of a group of endogenous 22-nucleotide siRNAs that are generated by the DICER-LIKE 2 (DCL2) protein in plants. When cytoplasmic RNA decay and DCL4 are deficient, the resulting massive accumulation of 22-nucleotide siRNAs causes pleiotropic growth disorders, including severe dwarfism, meristem defects and pigmentation. Notably, two genes that encode nitrate reductases-NIA1 and NIA2-produce nearly half of the 22-nucleotide siRNAs. Production of 22-nucleotide siRNAs triggers the amplification of gene silencing and induces translational repression both gene specifically and globally. Moreover, these 22-nucleotide siRNAs preferentially accumulate upon environmental stress, especially those siRNAs derived from NIA1/2, which act to restrain translation, inhibit plant growth and enhance stress responses. Thus, our research uncovers the unique properties of 22-nucleotide siRNAs, and reveals their importance in plant adaptation to environmental stresses.

Computational tools for plant small RNA detection and categorization

Small RNAs (sRNAs) are important short-length molecules with regulatory functions essential for plant development and plasticity. High-throughput sequencing of total sRNA populations has revealed that the largest share of sRNA remains uncategorized. To better understand the role of sRNA-mediated cellular regulation, it is necessary to create accurate and comprehensive catalogues of sRNA and their sequence features, a task that currently relies on nontrivial bioinformatic approaches. Although a large number of computational tools have been developed to predict features of sRNA sequences, these tools are mostly dedicated to microRNAs and none integrates the functionalities necessary to describe units from all sRNA pathways thus far discovered in plants. Here, we review the different classes of sRNA found in plants and describe available bioinformatics tools that can help in their detection and categorization.

Arabidopsis retrotransposon virus-like particles and their regulation by epigenetically activated small RNA

In Arabidopsis, LTR retrotransposons are activated by mutations in the chromatin gene DECREASE in DNA METHYLATION 1 (DDM1), giving rise to 21- to 22-nt epigenetically activated siRNA (easiRNA) that depend on RNA DEPENDENT RNA POLYMERASE 6 (RDR6). We purified virus-like particles (VLPs) from ddm1 and ddm1rdr6 mutants in which genomic RNA is reverse transcribed into complementary DNA. High-throughput short-read and long-read sequencing of VLP DNA (VLP DNA-seq) revealed a comprehensive catalog of active LTR retrotransposons without the need for mapping transposition, as well as independent of genomic copy number. Linear replication intermediates of the functionally intact COPIA element EVADE revealed multiple central polypurine tracts (cPPTs), a feature shared with HIV in which cPPTs promote nuclear localization. For one member of the ATCOPIA52 subfamily (SISYPHUS), cPPT intermediates were not observed, but abundant circular DNA indicated transposon "suicide" by auto-integration within the VLP. easiRNA targeted EVADE genomic RNA, polysome association of GYPSY (ATHILA) subgenomic RNA, and transcription via histone H3 lysine-9 dimethylation. VLP DNA-seq provides a comprehensive landscape of LTR retrotransposons and their control at transcriptional, post-transcriptional, and reverse transcriptional levels.

Genome wide identification, characterization and validation of novel miRNA-based SSR markers in pomegranate (Punica granatum L.)

A total of 17,439 mature miRNAs (~ 21 nt) earlier generated through RNA seq in the pomegranate were used for in silico analysis. After complexity reduction, a total of 1922 representative mature miRNAs were selected and used as query sequences against pomegranate genome to retrieve 2540 homologous contigs with flanking regions (~ 800). By using pre-miRNA prediction web server, a total of 1028 true contigs harbouring pri-miRNAs encoding 1162 pre-miRNAs were identified. Survey of these sequences for SSRs yielded a total of 1358 and 238 SSRs specific to pri-miRNA and pre-miRNAs, respectively. Of these, primer pairs were designed for 897 pri-miRNA and 168 pre-miRNA SSRs. In pri-miRNA sequences, hexa-nucleotides repeats were found to be most abundant (44.18%) followed by mono- (18.41%) and di-nucleotide (17.01%), which is also observed in pre-miRNA sequences. Further, a set of 51 randomly selected pre-miRNA-SSRs was examined for marker polymorphism. The experimental validation of these markers on eight pomegranate genotypes demonstrated 92.15% polymorphism. Utility of these functional markers was confirmed via examination of genetic diversity of 18 pomegranate genotypes using 15 miRNA-SSRs. Further, potential application of miRNA-SSRs for discovery of trait specific candidate genes was showed by validating 51 mature miRNA against publically available 2047 EST sequences of pomegranate by target and network analysis. In summary, the current study offers novel functional molecular markers for pomegranate genetic improvement.

Mobile genomics: tools and techniques for tackling transposons

Next-generation sequencing approaches have fundamentally changed the types of questions that can be asked about gene function and regulation. With the goal of approaching truly genome-wide quantifications of all the interaction partners and downstream effects of particular genes, these quantitative assays have allowed for an unprecedented level of detail in exploring biological interactions. However, many challenges remain in our ability to accurately describe and quantify the interactions that take place in those hard to reach and extremely repetitive regions of our genome comprised mostly of transposable elements (TEs). Tools dedicated to TE-derived sequences have lagged behind, making the inclusion of these sequences in genome-wide analyses difficult. Recent improvements, both computational and experimental, allow for the better inclusion of TE sequences in genomic assays and a renewed appreciation for the importance of TE biology. This review will discuss the recent improvements that have been made in the computational analysis of TE-derived sequences as well as the areas where such analysis still proves difficult. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

Nanopore sequencing and Hi-C scaffolding provide insight into the evolutionary dynamics of transposable elements and piRNA production in wild strains of Drosophila melanogaster

Illumina sequencing has allowed for population-level surveys of transposable element (TE) polymorphism via split alignment approaches, which has provided important insight into the population dynamics of TEs. However, such approaches are not able to identify insertions of uncharacterized TEs, nor can they assemble the full sequence of inserted elements. Here, we use nanopore sequencing and Hi-C scaffolding to produce de novo genome assemblies for two wild strains of Drosophila melanogaster from the Drosophila Genetic Reference Panel (DGRP). Ovarian piRNA populations and Illumina split-read TE insertion profiles have been previously produced for both strains. We find that nanopore sequencing with Hi-C scaffolding produces highly contiguous, chromosome-length scaffolds, and we identify hundreds of TE insertions that were missed by Illumina-based methods, including a novel micropia-like element that has recently invaded the DGRP population. We also find hundreds of piRNA-producing loci that are specific to each strain. Some of these loci are created by strain-specific TE insertions, while others appear to be epigenetically controlled. Our results suggest that Illumina approaches reveal only a portion of the repetitive sequence landscape of eukaryotic genomes and that population-level resequencing using long reads is likely to provide novel insight into the evolutionary dynamics of repetitive elements.

Integrated annotations and analyses of small RNA-producing loci from 47 diverse plants

Plant endogenous small RNAs (sRNAs) are important regulators of gene expression. There are two broad categories of plant sRNAs: microRNAs (miRNAs) and endogenous short interfering RNAs (siRNAs). MicroRNA loci are relatively well-annotated but compose only a small minority of the total sRNA pool; siRNA locus annotations have lagged far behind. Here, we used a large data set of published and newly generated sRNA sequencing data (1333 sRNA-seq libraries containing more than 20 billion reads) and a uniform bioinformatic pipeline to produce comprehensive sRNA locus annotations of 47 diverse plants, yielding more than 2.7 million sRNA loci. The two most numerous classes of siRNA loci produced mainly 24- and 21-nucleotide (nt) siRNAs, respectively. Most often, 24-nt-dominated siRNA loci occurred in intergenic regions, especially at the 5′-flanking regions of protein-coding genes. In contrast, 21-nt-dominated siRNA loci were most often derived from double-stranded RNA precursors copied from spliced mRNAs. Genic 21-nt-dominated loci were especially common from disease resistance genes, including from a large number of monocots. Individual siRNA sequences of all types showed very little conservation across species, whereas mature miRNAs were more likely to be conserved. We developed a web server where our data and several search and analysis tools are freely accessible.

Evolution of PHAS loci in the young spike of Allohexaploid wheat

BACKGROUND

PhasiRNAs (phased secondary siRNAs) play important regulatory roles in the development processes and biotic or abiotic stresses in plants. Some of phasiRNAs involve in the reproductive development in grasses, which include two categories, 21-nt (nucleotide) and 24-nt phasiRNAs. They are triggered by miR2118 and miR2275 respectively, in premeiotic and meiotic anthers of rice, maize and other grass species. Wheat (Triticum aestivum) with three closely related subgenomes (subA, subB and subD), is a model of allopolyploid in plants. Knowledge about the role of phasiRNAs in the inflorescence development of wheat is absent until now, and the evolution of PHAS loci in polyploid plants is also unavailable.

RESULTS

Using 261 small RNA expression datasets from various tissues, a batch of PHAS (phasiRNA precursors) loci were identified in the young spike of wheat, most of which were regulated by miR2118 and miR2275 in their target site regions. Dissection of PHAS and their trigger miRNAs among the diploid (AA and DD), tetraploid (AABB) and hexaploid (AABBDD) genomes of Triticum indicated that distribution of PHAS loci were dominant randomly in local chromosomes, while miR2118 was dominant only in the subB genome. The diversity of PHAS loci in the three subgenomes of wheat and their progenitor genomes (AA, DD and AABB) suggested that they originated or diverged at least before the occurrence of the tetraploid AABB genome. The positive correlation between the PHAS loci or the trigger miRNAs and the ploidy of genome indicated the expansion of genome was the major drive force for the increase of PHAS loci and their trigger miRNAs in Triticum. In addition, the expression profiles of the PHAS transcripts suggested they responded to abiotic stresses such as cold stress in wheat.

CONCLUSIONS

Altogether, non-coding phasiRNAs are conserved transcriptional regulators that display quick plasticity in Triticum genome. They may be involved in reproductive development and abiotic stress in wheat. It could be referred to molecular research on male reproductive development in Triticum.

Translational Landscape of Protein-Coding and Non-Protein-Coding RNAs upon Light Exposure in Arabidopsis

Light is one of the most essential environmental clues for plant growth and morphogenesis. Exposure to blue monochromatic light from darkness is a turning point for plant biological activity, and as a result dramatic changes in gene expression occur. To understand the translational impacts of blue light, we have performed ribosome profiling analysis and called translated open reading frames (ORFs) de novo within not only mRNAs but also non-coding RNAs (ncRNAs). Translation efficiency of 3,823 protein-coding ORFs, such as nuclear chloroplast-related genes, was up-regulated by blue light exposure. Moreover, the translational activation of the microRNA biogenesis-related genes, DCL1 and HYL1, was induced by blue light. Considering the 3-nucleotide codon periodicity of ribosome footprints, a few hundred short ORFs lying on ncRNAs and upstream ORFs (uORFs) on mRNAs were found that had differential translation status between blue light and dark. uORFs are known to have a negative effect on the expression of the main ORFs (mORFs) on the same mRNAs. Our analysis suggests that the translation of uORFs is likely to be more stimulated than that of the corresponding mORFs, and uORF-mediated translational repression of the mORFs in five genes was alleviated by blue light exposure. With data-based annotation of the ORFs, our analysis provides insights into the translatome in response to environmental changes, such as those involving light.

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.

Arabidopsis RNA Polymerase IV generates 21-22 nucleotide small RNAs that can participate in RNA-directed DNA methylation and may regulate genes

The plant-specific RNA Polymerase IV (Pol IV) transcribes heterochromatic regions, including many transposable elements (TEs), with the well-described role of generating 24 nucleotide (nt) small interfering RNAs (siRNAs). These siRNAs target DNA methylation back to TEs to reinforce the boundary between heterochromatin and euchromatin. In the male gametophytic phase of the plant life cycle, pollen, Pol IV switches to generating primarily 21-22 nt siRNAs, but the biogenesis and function of these siRNAs have been enigmatic. In contrast to being pollen-specific, we identified that Pol IV generates these 21-22 nt siRNAs in sporophytic tissues, likely from the same transcripts that are processed into the more abundant 24 nt siRNAs. The 21-22 nt forms are specifically generated by the combined activities of DICER proteins DCL2/DCL4 and can participate in RNA-directed DNA methylation. These 21-22 nt siRNAs are also loaded into ARGONAUTE1 (AGO1), which is known to function in post-transcriptional gene regulation. Like other plant siRNAs and microRNAs incorporated into AGO1, we find a signature of genic mRNA cleavage at the predicted target site of these siRNAs, suggesting that Pol IV-generated 21-22 nt siRNAs may function to regulate gene transcript abundance. Our data provide support for the existing model that in pollen Pol IV functions in gene regulation. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

AUSPP: A universal short-read pre-processing package

There are many short-read aligners that can map short reads to a reference genome/sequence, and most of them can directly accept a FASTQ file as the input query file. However, the raw data usually need to be pre-processed. Few software programs specialize in pre-processing raw data generated by a variety of next-generation sequencing (NGS) technologies. Here, we present AUSPP, a Perl script-based pipeline for pre-processing and automatic mapping of NGS short reads. This pipeline encompasses quality control, adaptor trimming, collapsing of reads, structural RNA removal, length selection, read mapping, and normalized wiggle file creation. It facilitates the processing from raw data to genome mapping and is therefore a powerful tool for the steps before meta-analysis. Most importantly, since AUSPP has default processing pipeline settings for many types of NGS data, most of the time, users will simply need to provide the raw data and genome. AUSPP is portable and easy to install, and the source codes are freely available at https://github.com/highlei/AUSPP.

Deep-sequencing of Solanum commersonii small RNA libraries reveals riboregulators involved in cold stress response

Among wild species used in potato breeding, Solanum commersonii displays the highest tolerance to low temperatures under both acclimated (ACC) and non-acclimated (NACC) conditions. It is also the first wild potato relative with a known whole genome sequence. Recent studies have shown that abiotic stresses induce changes in the expression of many small non-coding RNA (sncRNA). We determined the small non-coding RNA (sncRNAome) of two clones of S. commersonii contrasting in their cold response phenotypes via smRNAseq. Differential analysis provided evidence that expression of several miRNAs changed in response to cold stress conditions. Conserved miR408a and miR408b changed their expression under NACC conditions, whereas miR156 and miR169 were differentially expressed only under ACC conditions. We also report changes in tasiRNA and secondary siRNA expression under both stress conditions. Our results reveal possible roles of sncRNA in the regulatory networks associated with tolerance to low temperatures and provide useful information for a more strategic use of genomic resources in potato breeding.

Integrated Genome-Scale Analysis and Northern Blot Detection of Retrotransposon siRNAs Across Plant Species

Cells have sophisticated RNA-directed mechanisms to regulate genes, destroy viruses, or silence transposable elements (TEs). In terrestrial plants, a specialized non-coding RNA machinery involving RNA polymerase IV (Pol IV) and small interfering RNAs (siRNAs) targets DNA methylation and silencing to TEs. Here, we present a bioinformatics protocol for annotating and quantifying siRNAs that derive from long terminal repeat (LTR) retrotransposons. The approach was validated using small RNA northern blot analyses, comparing the species Arabidopsis thaliana and Brachypodium distachyon. To assist hybridization probe design, we configured a genome browser to show small RNA-seq mappings in distinct colors and shades according to their nucleotide lengths and abundances, respectively. Samples from wild-type and pol IV mutant plants, cross-species negative controls, and a conserved microRNA control validated the detected siRNA signals, confirming their origin from specific TEs and their Pol IV-dependent biogenesis. Moreover, an optimized labeling method yielded probes that could detect low-abundance siRNAs from B. distachyon TEs. The integration of de novo TE annotation, small RNA-seq profiling, and northern blotting, as outlined here, will facilitate the comparative genomic analysis of RNA silencing in crop plants and non-model species.

Small noncoding RNA discovery and profiling with sRNAtools based on high-throughput sequencing

Small noncoding RNAs (sRNA/sncRNAs) are generated from different genomic loci and play important roles in biological processes, such as cell proliferation and the regulation of gene expression. Next-generation sequencing (NGS) has provided an unprecedented opportunity to discover and quantify diverse kinds of sncRNA, such as tRFs (tRNA-derived small RNA fragments), phasiRNAs (phased, secondary, small-interfering RNAs), Piwi-interacting RNA (piRNAs) and plant-specific 24-nt short interfering RNAs (siRNAs). However, currently available web-based tools do not provide approaches to comprehensively analyze all of these diverse sncRNAs. This study presents a novel integrated platform, sRNAtools (https://bioinformatics.caf.ac.cn/sRNAtools), that can be used in conjunction with high-throughput sequencing to identify and functionally annotate sncRNAs, including profiling microRNAss, piRNAs, tRNAs, small nuclear RNAs, small nucleolar RNAs and rRNAs and discovering isomiRs, tRFs, phasiRNAs and plant-specific 24-nt siRNAs for up to 21 model organisms. Different modules, including single case, batch case, group case and target case, are developed to provide users with flexible ways of studying sncRNA. In addition, sRNAtools supports different ways of uploading small RNA sequencing data in a very interactive queue system, while local versions based on the program package/Docker/virtureBox are also available. We believe that sRNAtools will greatly benefit the scientific community as an integrated tool for studying sncRNAs.

Tissue-specific transposon-associated small RNAs in the gymnosperm tree, Norway spruce

BACKGROUND

Small RNAs (sRNAs) are regulatory molecules impacting on gene expression and transposon activity. MicroRNAs (miRNAs) are responsible for tissue-specific and environmentally-induced gene repression. Short interfering RNAs (siRNA) are constitutively involved in transposon silencing across different type of tissues. The male gametophyte in angiosperms has a unique set of sRNAs compared to vegetative tissues, including phased siRNAs from intergenic or genic regions, or epigenetically activated siRNAs. This is contrasted by a lack of knowledge about the sRNA profile of the male gametophyte of gymnosperms.

RESULTS

Here, we isolated mature pollen from male cones of Norway spruce and investigated its sRNA profiles. While 21-nt sRNAs is the major size class of sRNAs in needles, in pollen 21-nt and 24-nt sRNAs are the most abundant size classes. Although the 24-nt sRNAs were exclusively derived from TEs in pollen, both 21-nt and 24-nt sRNAs were associated with TEs. We also investigated sRNAs from somatic embryonic callus, which has been reported to contain 24-nt sRNAs. Our data show that the 24-nt sRNA profiles are tissue-specific and differ between pollen and cell culture.

CONCLUSION

Our data reveal that gymnosperm pollen, like angiosperm pollen, has a unique sRNA profile, differing from vegetative leaf tissue. Thus, our results reveal that angiosperm and gymnosperm pollen produce new size classes not present in vegetative tissues; while in angiosperm pollen 21-nt sRNAs are generated, in the gymnosperm Norway spruce 24-nt sRNAs are generated. The tissue-specific production of distinct TE-derived sRNAs in angiosperms and gymnosperms provides insights into the diversification process of sRNAs in TE silencing pathways between the two groups of seed plants.

Integrated Analysis of Small RNA, Transcriptome and Degradome Sequencing Provides New Insights into Floral Development and Abscission in Yellow Lupine (Lupinus luteus L.)

The floral development in an important legume crop yellow lupine (Lupinus luteus L., Taper cv.) is often affected by the abscission of flowers leading to significant economic losses. Small non-coding RNAs (sncRNAs), which have a proven effect on almost all developmental processes in other plants, might be of key players in a complex net of molecular interactions regulating flower development and abscission. This study represents the first comprehensive sncRNA identification and analysis of small RNA, transcriptome and degradome sequencing data in lupine flowers to elucidate their role in the regulation of lupine generative development. As shedding in lupine primarily concerns flowers formed at the upper part of the inflorescence, we analyzed samples from extreme parts of raceme separately and conducted an additional analysis of pedicels from abscising and non-abscising flowers where abscission zone forms. A total of 394 known and 28 novel miRNAs and 316 phased siRNAs were identified. In flowers at different stages of development 59 miRNAs displayed differential expression (DE) and 46 DE miRNAs were found while comparing the upper and lower flowers. Identified tasiR-ARFs were DE in developing flowers and were strongly expressed in flower pedicels. The DEmiR-targeted genes were preferentially enriched in the functional categories related to carbohydrate metabolism and plant hormone transduction pathways. This study not only contributes to the current understanding of how lupine flowers develop or undergo abscission but also holds potential for research aimed at crop improvement.

Secretion of an Argonaute protein by a parasitic nematode and the evolution of its siRNA guides

Extracellular RNA has been proposed to mediate communication between cells and organisms however relatively little is understood regarding how specific sequences are selected for export. Here, we describe a specific Argonaute protein (exWAGO) that is secreted in extracellular vesicles (EVs) released by the gastrointestinal nematode Heligmosomoides bakeri, at multiple copies per EV. Phylogenetic and gene expression analyses demonstrate exWAGO orthologues are highly conserved and abundantly expressed in related parasites but highly diverged in free-living genus Caenorhabditis. We show that the most abundant small RNAs released from the nematode parasite are not microRNAs as previously thought, but rather secondary small interfering RNAs (siRNAs) that are produced by RNA-dependent RNA Polymerases. The siRNAs that are released in EVs have distinct evolutionary properties compared to those resident in free-living or parasitic nematodes. Immunoprecipitation of exWAGO demonstrates that it specifically associates with siRNAs from transposons and newly evolved repetitive elements that are packaged in EVs and released into the host environment. Together this work demonstrates molecular and evolutionary selectivity in the small RNA sequences that are released in EVs into the host environment and identifies a novel Argonaute protein as the mediator of this.

Isolating Male Meiocytes from Maize and Wheat for "-Omics" Analyses

Genome-wide gene expression studies have become a routine approach due to the advances in sequencing technologies, their ease of use, and increasing affordability. Simultaneous investigation of small RNA expression adds further valuable information but is not adopted as widely yet. Both RNA-seq and small RNA-seq benefit from the use of specific cell types. Here, we describe a protocol for the isolation of male meiotic cells from maize or wheat plants, along with the application of downstream RNA sequencing, extendable to other -omics approaches.

Germ Granules Govern Small RNA Inheritance

In C. elegans nematodes, components of liquid-like germ granules were shown to be required for transgenerational small RNA inheritance. Surprisingly, we show here that mutants with defective germ granules can nevertheless inherit potent small RNA-based silencing responses, but some of the mutants lose this ability after many generations of homozygosity. Animals mutated in pptr-1, which is required for stabilization of P granules in the early embryo, display extraordinarily strong heritable RNAi responses, lasting for tens of generations. Intriguingly, the RNAi capacity of descendants derived from mutants defective in the core germ granule proteins MEG-3 and MEG-4 is determined by the genotype of the ancestors and changes transgenerationally. Further, whether the meg-3/4 mutant alleles were present in the paternal or maternal lineages leads to different transgenerational consequences. Small RNA inheritance, rather than maternal contribution of the germ granules themselves, mediates the transgenerational defects in RNAi of meg-3/4 mutants and their progeny. Accordingly, germ granule defects lead to heritable genome-wide mis-expression of endogenous small RNAs. Upon disruption of germ granules, hrde-1 mutants can inherit RNAi, although HRDE-1 was previously thought to be absolutely required for RNAi inheritance. We propose that germ granules sort and shape the RNA pool, and that small RNA inheritance maintains this activity for multiple generations.

SPAR: small RNA-seq portal for analysis of sequencing experiments

The introduction of new high-throughput small RNA sequencing protocols that generate large-scale genomics datasets along with increasing evidence of the significant regulatory roles of small non-coding RNAs (sncRNAs) have highlighted the urgent need for tools to analyze and interpret large amounts of small RNA sequencing data. However, it remains challenging to systematically and comprehensively discover and characterize sncRNA genes and specifically-processed sncRNA products from these datasets. To fill this gap, we present Small RNA-seq Portal for Analysis of sequencing expeRiments (SPAR), a user-friendly web server for interactive processing, analysis, annotation and visualization of small RNA sequencing data. SPAR supports sequencing data generated from various experimental protocols, including smRNA-seq, short total RNA sequencing, microRNA-seq, and single-cell small RNA-seq. Additionally, SPAR includes publicly available reference sncRNA datasets from our DASHR database and from ENCODE across 185 human tissues and cell types to produce highly informative small RNA annotations across all major small RNA types and other features such as co-localization with various genomic features, precursor transcript cleavage patterns, and conservation. SPAR allows the user to compare the input experiment against reference ENCODE/DASHR datasets. SPAR currently supports analyses of human (hg19, hg38) and mouse (mm10) sequencing data. SPAR is freely available at https://www.lisanwanglab.org/SPAR.

Roles of Argonautes and Dicers on Sclerotinia sclerotiorum Antiviral RNA Silencing

RNA silencing or RNA interference (RNAi) is an essential mechanism in animals, fungi, and plants that functions in gene regulation and defense against foreign nucleic acids. In fungi, RNA silencing has been shown to function primarily in defense against invasive nucleic acids. We previously determined that mycoviruses are triggers and targets of RNA silencing in Sclerotinia sclerotiorum. However, recent progresses in RNAi or dsRNA-based pest control requires more detailed characterization of the RNA silencing pathways in S. sclerotiorum to investigate the utility of dsRNA-based strategy for white mold control. This study elucidates the roles of argonaute enzymes, agl-2 and agl-4, in small RNA metabolism in S. sclerotiorum. Gene disruption mutants of agl-2 and agl-4 were compared for changes in phenotype, virulence, viral susceptibility, and small RNA profiles. The Δagl-2 mutant but not the Δagl-4 mutant had significantly slower growth and virulence prior to virus infection. Similarly, the Δagl-2 mutant but not the Δagl-4 mutant, showed greater debilitation under virus infection compared to uninfected strains. The responses were confirmed in complementation studies and revealed the antiviral role of agl-2. Gene disruption mutants of agl-2, agl-4, Dicer-like (dcl)-1, and dcl-2 did not change the stability of the most abundant endogenous small RNAs, which suggests the existence of alternative enzymes/pathways for small RNA biogenesis in S. sclerotiorum. Furthermore, in vitro synthesized dsRNA targeting agl-2 showed a significantly reduced average lesion diameter (P < 0.05) on canola leaves with agl-2 down-regulated compared to controls. This is the first report describing the effectiveness of RNA pesticides targeting S. sclerotiorum RNA silencing pathway for the control of the economically important pathogen.

Small RNA discovery in the interaction between barley and the powdery mildew pathogen

Background

Plants encounter pathogenic and non-pathogenic microorganisms on a nearly constant basis. Small RNAs such as siRNAs and miRNAs/milRNAs influence pathogen virulence and host defense responses. We exploited the biotrophic interaction between the powdery mildew fungus, Blumeria graminis f. sp. hordei (Bgh), and its diploid host plant, barley (Hordeum vulgare) to explore fungal and plant sRNAs expressed during Bgh infection of barley leaf epidermal cells.

Results

RNA was isolated from four fast-neutron immune-signaling mutants and their progenitor over a time course representing key stages of Bgh infection, including appressorium formation, penetration of epidermal cells, and development of haustorial feeding structures. The Cereal Introduction (CI) 16151 progenitor carries the resistance allele Mla6, while Bgh isolate 5874 harbors the AVR_a6 avirulence effector, resulting in an incompatible interaction. Parallel Analysis of RNA Ends (PARE) was used to verify sRNAs with likely transcript targets in both barley and Bgh. Bgh sRNAs are predicted to regulate effectors, metabolic genes, and translation-related genes. Barley sRNAs are predicted to influence the accumulation of transcripts that encode auxin response factors, NAC transcription factors, homeodomain transcription factors, and several splicing factors. We also identified phasing small interfering RNAs (phasiRNAs) in barley that overlap transcripts that encode receptor-like kinases (RLKs) and nucleotide-binding, leucine-rich domain proteins (NLRs).

Conclusions

These data suggest that Bgh sRNAs regulate gene expression in metabolism, translation-related, and pathogen effectors. PARE-validated targets of predicted Bgh milRNAs include both EKA (effectors homologous to AVR_k1 and AVR_a10) and CSEP (candidate secreted effector protein) families. We also identified barley phasiRNAs and miRNAs in response to Bgh infection. These include phasiRNA loci that overlap with a significant proportion of receptor-like kinases, suggesting an additional sRNA control mechanism may be active in barley leaves as opposed to predominant R-gene phasiRNA overlap in many eudicots. In addition, we identified conserved miRNAs, novel miRNA candidates, and barley genome mapped sRNAs that have PARE validated transcript targets in barley. The miRNA target transcripts are enriched in transcription factors, signaling-related proteins, and photosynthesis-related proteins. Together these results suggest both barley and Bgh control metabolism and infection-related responses via the specific accumulation and targeting of genes via sRNAs.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5947-z) contains supplementary material, which is available to authorized users.

Small RNA-seq reveals novel regulatory components for apomixis in Paspalum notatum

Background

Apomixis is considered an evolutionary deviation of the sexual reproductive pathway leading to the generation of clonal maternal progenies by seeds. Recent evidence from model and non-model species suggested that this trait could be modulated by epigenetic mechanisms involving small RNAs (sRNAs). Here we profiled floral sRNAs originated from apomictic and sexual Paspalum notatum genotypes in order to identify molecular pathways under epigenetic control that might be involved in the transition from sexuality to agamospermy.

Results

The mining of genes participating in sRNA-directed pathways from floral Paspalum transcriptomic resources showed these routes are functional during reproductive development, with several members differentially expressed in apomictic and sexual plants. Triplicate floral sRNA libraries derived from apomictic and a sexual genotypes were characterized by using high-throughput sequencing technology. EdgeR was apply to compare the number of sRNA reads between sexual and apomictic libraries that map over all Paspalum floral transcripts. A total of 1525 transcripts showed differential sRNA representation, including genes related to meiosis, plant hormone signaling, biomolecules transport, transcription control and cell cycle. Survey for miRNA precursors on transcriptome and genome references allowed the discovery of 124 entities, including 40 conserved and 8 novel ones. Fifty-six clusters were differentially represented in apomictic and sexual plants. All differentially expressed miRNAs were up-regulated in apomictic libraries but miR2275, which showed different family members with opposed representation. Examination of predicted miRNAs targets detected 374 potential candidates. Considering sRNA, miRNAs and target surveys together, 14 genes previously described as related with auxin metabolism, transport and signaling were detected, including AMINO ACID/AUXIN PERMEASE 15, IAA-AMIDO SYNTHETASE GH3–8, IAA30, miR160, miR167, miR164, miR319, ARF2, ARF8, ARF10, ARF12, AFB2, PROLIFERATING CELL FACTOR 6 and NITRATE TRANSPORTER 1.1.

Conclusions

This work provides a comprehensive survey of the sRNA differential representation in flowers of sexual and apomictic Paspalum notatum plants. An integration of the small RNA profiling data presented here and previous transcriptomic information suggests that sRNA-mediated regulation of auxin pathways is pivotal in promoting apomixis. These results will underlie future functional characterization of the molecular components mediating the switch from sexuality to apomixis.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5881-0) contains supplementary material, which is available to authorized users.

Neuronal Small RNAs Control Behavior Transgenerationally

It is unknown whether the activity of the nervous system can be inherited. In Caenorhabditis elegans nematodes, parental responses can transmit heritable small RNAs that regulate gene expression transgenerationally. In this study, we show that a neuronal process can impact the next generations. Neurons-specific synthesis of RDE-4-dependent small RNAs regulates germline amplified endogenous small interfering RNAs (siRNAs) and germline gene expression for multiple generations. Further, the production of small RNAs in neurons controls the chemotaxis behavior of the progeny for at least three generations via the germline Argonaute HRDE-1. Among the targets of these small RNAs, we identified the conserved gene saeg-2, which is transgenerationally downregulated in the germline. Silencing of saeg-2 following neuronal small RNA biogenesis is required for chemotaxis under stress. Thus, we propose a small-RNA-based mechanism for communication of neuronal processes transgenerationally.

•

C. elegans neuronal small RNAs are characterized by RNA sequencing

•

RDE-4-dependent neuronal endogenous small RNAs communicate with the germline

•

Germline HRDE-1 mediates transgenerational regulation by neuronal small RNAs

•

Neuronal small RNAs regulate germline genes to control behavior transgenerationally

Heat shock protein DNAJA1 stabilizes PIWI proteins to support regeneration and homeostasis of planarian Schmidtea mediterranea

PIWI proteins are key regulators of germline and somatic stem cells throughout different evolutionary lineages. However, how PIWI proteins themselves are regulated remains largely unknown. To identify candidate proteins that interact with PIWI proteins and regulate their stability, here we established a yeast two-hybrid (Y2H) assay in the planarian species Schmidtea mediterranea We show that DNAJA1, a heat shock protein 40 family member, interacts with the PIWI protein SMEDWI-2, as validated by the Y2H screen and co-immunoprecipitation assays. We found that DNAJA1 is enriched in planarian adult stem cells, the nervous system, and intestinal tissues. DNAJA1-knockdown abolished planarian regeneration and homeostasis, compromised stem cell maintenance and PIWI-interacting RNA (piRNA) biogenesis, and deregulated SMEDWI-1/2 target genes. Mechanistically, we observed that DNAJA1 is required for the stability of SMEDWI-1 and SMEDWI-2 proteins. Furthermore, we noted that human DNAJA1 binds to Piwi-like RNA-mediated gene silencing 1 (PIWIL1) and is required for PIWIL1 stability in human gastric cancer cells. In summary, our results reveal not only an evolutionarily conserved functional link between PIWI and DNAJA1 that is essential for PIWI protein stability and piRNA biogenesis, but also an important role of DNAJA1 in the control of proteins involved in stem cell regulation.

Small RNA-Seq Analysis Reveals miRNA Expression Dynamics Across Tissues in the Malaria Vector, Anopheles gambiae

Malaria continues to be a major global health problem, where disease transmission is deeply linked to the repeated blood feeding nature of the anautogenous mosquito. Given the tight link between blood feeding and disease transmission, understanding basic biology behind mosquito physiology is a requirement for developing effective vector-borne disease control strategies. In the mosquito, numerous loss of function studies with notable phenotypes demonstrate microRNAs (miRNAs) play significant roles in mosquito physiology. While the field appreciates the importance of a handful of miRNAs, we still need global mosquito tissue miRNA transcriptome studies. To address this need, our goal was to determine the miRNA transcriptome for multiple tissues of the pre-vitellogenic mosquito. To this end, by using small RNA-Seq analysis, we determined miRNA transcriptomes in tissues critical for mosquito reproduction and immunity including (i) fat body-abdominal wall enriched tissues, (ii) midguts, (iii) ovaries, and (iv) remaining tissues comprised of the head and thorax. We found numerous examples of miRNAs exhibiting pan-tissue high- or low- expression, tissue exclusion, and tissue enrichment. We also updated and consolidated the miRNA catalog and provided a detailed genome architecture map for the malaria vector, Anopheles gambiae. This study aims to build a foundation for future research on how miRNAs and potentially other small RNAs regulate mosquito physiology as it relates to vector-borne disease transmission.

Automated analysis of small RNA datasets with RAPID

Understanding the role of short-interfering RNA (siRNA) in diverse biological processes is of current interest and often approached through small RNA sequencing. However, analysis of these datasets is difficult due to the complexity of biological RNA processing pathways, which differ between species. Several properties like strand specificity, length distribution, and distribution of soft-clipped bases are few parameters known to guide researchers in understanding the role of siRNAs. We present RAPID, a generic eukaryotic siRNA analysis pipeline, which captures information inherent in the datasets and automatically produces numerous visualizations as user-friendly HTML reports, covering multiple categories required for siRNA analysis. RAPID also facilitates an automated comparison of multiple datasets, with one of the normalization techniques dedicated for siRNA knockdown analysis, and integrates differential expression analysis using DESeq2.

Too Many False Targets for MicroRNAs: Challenges and Pitfalls in Prediction of miRNA Targets and Their Gene Ontology in Model and Non-model Organisms

Short ("seed") or extended base pairing between microRNAs (miRNAs) and their target RNAs enables post-transcriptional silencing in many organisms. These interactions allow the computational prediction of potential targets. In model organisms, predicted targets are frequently validated experimentally; hence meaningful miRNA-regulated processes are reported. However, in non-models, these reports mostly rely on computational prediction alone. Many times, further bioinformatic analyses such as Gene Ontology (GO) enrichment are based on these in silico projections. Here such approaches are reviewed, their caveats are highlighted and the ease of picking false targets from predicted lists is demonstrated. Discoveries that shed new light on how miRNAs evolved to regulate targets in various phyletic groups are discussed, in addition to the pitfalls of target identification in non-model organisms. The goal is to prevent the misuse of bioinformatic tools, as they cannot bypass the biological understanding of miRNA-target regulation.

H3K9me3 is required for inheritance of small RNAs that target a unique subset of newly evolved genes

In Caenorhabditis elegans, RNA interference (RNAi) responses can transmit across generations via small RNAs. RNAi inheritance is associated with Histone-3-Lysine-9 tri-methylation (H3K9me3) of the targeted genes. In other organisms, maintenance of silencing requires a feed-forward loop between H3K9me3 and small RNAs. Here, we show that in C. elegans not only is H3K9me3 unnecessary for inheritance, the modification's function depends on the identity of the RNAi-targeted gene. We found an asymmetry in the requirement for H3K9me3 and the main worm H3K9me3 methyltransferases, SET-25 and SET-32. Both methyltransferases promote heritable silencing of the foreign gene gfp, but are dispensable for silencing of the endogenous gene oma-1. Genome-wide examination of heritable endogenous small interfering RNAs (endo-siRNAs) revealed that endo-siRNAs that depend on SET-25 and SET-32 target newly acquired and highly H3K9me3 marked genes. Thus, 'repressive' chromatin marks could be important specifically for heritable silencing of genes which are flagged as 'foreign', such as gfp.

Editorial note

This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Temporal Small RNA Expression Profiling under Drought Reveals a Potential Regulatory Role of Small Nucleolar RNAs in the Drought Responses of Maize

Small RNAs (sRNAs) are short noncoding RNAs that play roles in many biological processes, including drought responses in plants. However, how the expression of sRNAs dynamically changes with the gradual imposition of drought stress in plants is largely unknown. We generated time-series sRNA sequence data from maize ( L.) seedlings under drought stress (DS) and under well-watered (WW) conditions at the same time points. Analyses of length, functional annotation, and abundance of 736,372 nonredundant sRNAs from both DS and WW data, as well as genome copy numbers at the corresponding genomic regions, revealed distinct patterns of abundance and genome organization for different sRNA classes. The analysis identified 6646 sRNAs whose regulation was altered in response to drought stress. Among drought-responsive sRNAs, 1325 showed transient downregulation by the seventh day, coinciding with visible symptoms of drought stress. The profiles revealed drought-responsive microRNAs, as well as other sRNAs that originated from ribosomal RNAs (rRNAs), splicing small nuclear RNAs, and small nucleolar RNAs (snoRNA). Expression profiles of their sRNA derivers indicated that snoRNAs might play a regulatory role through regulating the stability of rRNAs and splicing small nuclear RNAs under drought condition.

Identification and profiling of narrow-leafed lupin (Lupinus angustifolius) microRNAs during seed development

BACKGROUND

Whilst information regarding small RNAs within agricultural crops is increasing, the miRNA composition of the nutritionally valuable pulse narrow-leafed lupin (Lupinus angustifolius) remains unknown.

RESULTS

By conducting a genome- and transcriptome-wide survey we identified 7 Dicer-like and 16 Argonaute narrow-leafed lupin genes, which were highly homologous to their legume counterparts. We identified 43 conserved miRNAs belonging to 16 families, and 13 novel narrow-leafed lupin-specific miRNAs using high-throughput sequencing of small RNAs from foliar and root and five seed development stages. We observed up-regulation of members of the miRNA families miR167, miR399, miR156, miR319 and miR164 in narrow-leafed lupin seeds, and confirmed expression of miR156, miR166, miR164, miR1507 and miR396 using quantitative RT-PCR during five narrow-leafed lupin seed development stages. We identified potential targets for the conserved and novel miRNAs and were able to validate targets of miR399 and miR159 using 5' RLM-RACE. The conserved miRNAs are predicted to predominately target transcription factors and 93% of the conserved miRNAs originate from intergenic regions. In contrast, only 43% of the novel miRNAs originate from intergenic regions and their predicted targets were more functionally diverse.

CONCLUSION

This study provides important insights into the miRNA gene regulatory networks during narrow-leafed lupin seed development.

Meiocyte-Specific and AtSPO11-1-Dependent Small RNAs and Their Association with Meiotic Gene Expression and Recombination

Meiotic recombination ensures accurate chromosome segregation and results in genetic diversity in sexually reproducing eukaryotes. Over the last few decades, the genetic regulation of meiotic recombination has been extensively studied in many organisms. However, the role of endogenous meiocyte-specific small RNAs (ms-sRNAs; 21-24 nucleotide [nt]) and their involvement in meiotic recombination are unclear. Here, we sequenced the total small RNA (sRNA) and messenger RNA populations from meiocytes and leaves of wild type Arabidopsis (Arabidopsis thaliana) and meiocytes of spo11-1, a mutant defective in double-strand break formation, and we discovered 2,409 ms-sRNA clusters, 1,660 of which areSPORULATION 11-1 (AtSPO11-1)-dependent. Unlike mitotic small interfering RNAs that are enriched in intergenic regions and associated with gene silencing, ms-sRNAs are significantly enriched in genic regions and exhibit a positive correlation with genes that are preferentially expressed in meiocytes (i.e. Arabidopsis SKP1-LIKE1 and RAD51), in a fashion unrelated to DNA methylation. We also found that AtSPO11-1-dependent sRNAs have distinct characteristics compared with ms-sRNAs and tend to be associated with two known types of meiotic recombination hotspot motifs (i.e. CTT-repeat and A-rich motifs). These results reveal different meiotic and mitotic sRNA landscapes and provide new insights into how sRNAs relate to gene expression in meiocytes and meiotic recombination.

Plant Regulomics Portal (PRP): a comprehensive integrated regulatory information and analysis portal for plant genomes

Gene regulation is a highly complex and networked phenomenon where multiple tiers of control determine the cell state in a spatio-temporal manner. Among these, the transcription factors, DNA and histone modifications, and post-transcriptional control by small RNAs like miRNAs serve as major regulators. An understanding of the integrative and spatio-temporal impact of these regulatory factors can provide better insights into the state of a ‘cell system’. Yet, there are limited resources available to this effect. Therefore, we hereby report an integrative information portal (Plant Regulomics Portal; PRP) for plants for the first time. The portal has been developed by integrating a huge amount of curated data from published sources, RNA-, methylome- and sRNA/miRNA sequencing, histone modifications and repeats, gene ontology, digital gene expression and characterized pathways. The key features of the portal include a regulatory search engine for fetching numerous analytical outputs and tracks of the abovementioned regulators and also a genome browser for integrated visualization of the search results. It also has numerous analytical features for analyses of transcription factors (TFs) and sRNA/miRNA, spot-specific methylation, gene expression and interactions and details of pathways for any given genomic element. It can also provide information on potential RdDM regulation, while facilitating enrichment analysis, generation of visually rich plots and downloading of data in a selective manner. Visualization of intricate biological networks is an important feature which utilizes the Neo4j Graph database making analysis of relationships and long-range system viewing possible. Till date, PRP hosts 571-GB processed data for four plant species namely Arabidopsis thaliana, Oryza sativa subsp. japonica, Zea mays and Glycine max.

Database URL: https://scbb.ihbt.res.in/PRP

Characterization of siRNAs clusters in Arabidopsis thaliana galls induced by the root-knot nematode Meloidogyne incognita

Background

Root-knot nematodes (RKN), genus Meloidogyne, are plant parasitic worms that have the ability to transform root vascular cylinder cells into hypertrophied, multinucleate and metabolically over-active feeding cells. Redifferentiation into feeding cells is the result of a massive transcriptional reprogramming of root cells targeted by RKN. Since RKN are able to induce similar feeding cells in roots of thousands of plant species, these worms are thought to manipulate essential and conserved plant molecular pathways.

Results

Small non-coding RNAs of uninfected roots and infected root galls induced by M. incognita from Arabidopsis thaliana were sequenced by high throughput sequencing. SiRNA populations were analysed by using the Shortstack algorithm. We identified siRNA clusters that are differentially expressed in infected roots and evidenced an over-representation of the 23–24 nt siRNAs in infected tissue. This size corresponds to heterochromatic siRNAs (hc-siRNAs) which are known to regulate expression of transposons and genes at the transcriptional level, mainly by inducing DNA methylation.

Conclusions

Correlation of siRNA clusters expression profile with transcriptomic data identified several protein coding genes that are candidates to be regulated by siRNAs at the transcriptional level by RNA directed DNA methylation (RdDM) pathway either directly or indirectly via silencing of neighbouring transposable elements.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5296-3) contains supplementary material, which is available to authorized users.

Identification and characterisation of microRNAs and their target genes in phosphate-starved Nicotiana benthamiana by small RNA deep sequencing and 5'RACE analysis

BACKGROUND

Phosphorus is an important macronutrient that is severely lacking in soils. In plants, specific microRNAs (miRNAs) essential for nutrient management and the regulation of stress responses are responsible for the control of many phosphate starvation responses. Further understanding of conserved and species-specific microRNA species has potential implications for the development of crops tolerant to soils with low phosphate.

RESULTS

This study identified and characterised phosphate starvation-responsive miRNAs in the native Australian tobacco Nicotiana benthamiana. Small RNA libraries were constructed and sequenced from phosphate-starved plant leaves, stems and roots. Twenty-four conserved miRNA families and 36 species-specific miRNAs were identified. The majority of highly phosphate starvation-responsive miRNAs were highly conserved, comprising of members from the miR399, miR827, and miR2111 families. In addition, two miRNA-star species were identified to be phosphate starvation-responsive. A total of seven miRNA targets were confirmed using RLM-5'RACE to be cleaved by five miRNA families, including two confirmed cleavage targets for Nbe-miR399 species, one for Nbe-miR2111, and two for Nbe-miR398. A number of N. benthamiana-specific features for conserved miRNAs were identified, including species-specific miRNA targets predicted or confirmed for miR399, miR827, and miR398.

CONCLUSIONS

Our results give an insight into the phosphate starvation-responsive miRNAs of Nicotiana benthamiana, and indicate that the phosphate starvation response pathways in N. benthamiana contain both highly conserved and species-specific components.

Transcriptional and Small RNA Responses of the White Mold Fungus Sclerotinia sclerotiorum to Infection by a Virulence-Attenuating Hypovirus

Mycoviruses belonging to the family Hypoviridae cause persistent infection of many different host fungi. We previously determined that the white mold fungus, Sclerotiniasclerotiorum, infected with Sclerotinia sclerotiorum hypovirus 2-L (SsHV2-L) exhibits reduced virulence, delayed/reduced sclerotial formation, and enhanced production of aerial mycelia. To gain better insight into the cellular basis for these changes, we characterized changes in mRNA and small RNA (sRNA) accumulation in S.sclerotiorum to infection by SsHV2-L. A total of 958 mRNAs and 835 sRNA-producing loci were altered after infection by SsHV2-L, among which >100 mRNAs were predicted to encode proteins involved in the metabolism and trafficking of carbohydrates and lipids. Both S. sclerotiorum endogenous and virus-derived sRNAs were predominantly 22 nt in length suggesting one dicer-like enzyme cleaves both. Novel classes of endogenous small RNAs were predicted, including phasiRNAs and tRNA-derived small RNAs. Moreover, S. sclerotiorum phasiRNAs, which were derived from noncoding RNAs and have the potential to regulate mRNA abundance in trans, showed differential accumulation due to virus infection. tRNA fragments did not accumulate differentially after hypovirus infection. Hence, in-depth analysis showed that infection of S. sclerotiorum by a hypovirulence-inducing hypovirus produced selective, large-scale reprogramming of mRNA and sRNA production.

Downregulation of RdDM during strawberry fruit ripening

Background

Recently, DNA methylation was proposed to regulate fleshy fruit ripening. Fleshy fruits can be distinguished by their ripening process as climacteric fruits, such as tomatoes, or non-climacteric fruits, such as strawberries. Tomatoes undergo a global decrease in DNA methylation during ripening, due to increased expression of a DNA demethylase gene. The dynamics and biological relevance of DNA methylation during the ripening of non-climacteric fruits are unknown.

Results

Here, we generate single-base resolution maps of the DNA methylome in immature and ripe strawberry. We observe an overall loss of DNA methylation during strawberry fruit ripening. Thus, ripening-induced DNA hypomethylation occurs not only in climacteric fruit, but also in non-climacteric fruit. Application of a DNA methylation inhibitor causes an early ripening phenotype, suggesting that DNA hypomethylation is important for strawberry fruit ripening. The mechanisms underlying DNA hypomethylation during the ripening of tomato and strawberry are distinct. Unlike in tomatoes, DNA demethylase genes are not upregulated during the ripening of strawberries. Instead, genes involved in RNA-directed DNA methylation are downregulated during strawberry ripening. Further, ripening-induced DNA hypomethylation is associated with decreased siRNA levels, consistent with reduced RdDM activity. Therefore, we propose that a downregulation of RdDM contributes to DNA hypomethylation during strawberry ripening.

Conclusions

Our findings provide new insight into the DNA methylation dynamics during the ripening of non-climacteric fruit and suggest a novel function of RdDM in regulating an important process in plant development.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1587-x) contains supplementary material, which is available to authorized users.