Reducing ligation bias of small RNAs in libraries for next generation sequencing

Identification and Characterization of Cyprinid Herpesvirus-3 (CyHV-3) Encoded MicroRNAs

MicroRNAs (miRNAs) are a class of small non-coding RNAs involved in post-transcriptional gene regulation. Some viruses encode their own miRNAs and these are increasingly being recognized as important modulators of viral and host gene expression. Cyprinid herpesvirus 3 (CyHV-3) is a highly pathogenic agent that causes acute mass mortalities in carp (Cyprinus carpio carpio) and koi (Cyprinus carpio koi) worldwide. Here, bioinformatic analyses of the CyHV-3 genome suggested the presence of non-conserved precursor miRNA (pre-miRNA) genes. Deep sequencing of small RNA fractions prepared from in vitro CyHV-3 infections led to the identification of potential miRNAs and miRNA–offset RNAs (moRNAs) derived from some bioinformatically predicted pre-miRNAs. DNA microarray hybridization analysis, Northern blotting and stem-loop RT-qPCR were then used to definitively confirm that CyHV-3 expresses two pre-miRNAs during infection in vitro. The evidence also suggested the presence of an additional four high-probability and two putative viral pre-miRNAs. MiRNAs from the two confirmed pre-miRNAs were also detected in gill tissue from CyHV-3-infected carp. We also present evidence that one confirmed miRNA can regulate the expression of a putative CyHV-3-encoded dUTPase. Candidate homologues of some CyHV-3 pre-miRNAs were identified in CyHV-1 and CyHV-2. This is the first report of miRNA and moRNA genes encoded by members of the Alloherpesviridae family, a group distantly related to the Herpesviridae family. The discovery of these novel CyHV-3 genes may help further our understanding of the biology of this economically important virus and their encoded miRNAs may have potential as biomarkers for the diagnosis of latent CyHV-3.

The genomes of two key bumblebee species with primitive eusocial organization

Background

The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats.

Results

We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits.

Conclusions

These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0623-3) contains supplementary material, which is available to authorized users.

A non-canonical RNA silencing pathway promotes mRNA degradation in basal Fungi

The increasing knowledge on the functional relevance of endogenous small RNAs (esRNAs) as riboregulators has stimulated the identification and characterization of these molecules in numerous eukaryotes. In the basal fungus Mucor circinelloides, an emerging opportunistic human pathogen, esRNAs that regulate the expression of many protein coding genes have been described. These esRNAs share common machinery for their biogenesis consisting of an RNase III endonuclease Dicer, a single Argonaute protein and two RNA-dependent RNA polymerases. We show in this study that, besides participating in this canonical dicer-dependent RNA interference (RNAi) pathway, the rdrp genes are involved in a novel dicer-independent degradation process of endogenous mRNAs. The analysis of esRNAs accumulated in wild type and silencing mutants demonstrates that this new rdrp-dependent dicer-independent regulatory pathway, which does not produce sRNA molecules of discrete sizes, controls the expression of target genes promoting the specific degradation of mRNAs by a previously unknown RNase. This pathway mainly regulates conserved genes involved in metabolism and cellular processes and signaling, such as those required for heme biosynthesis, and controls responses to specific environmental signals. Searching the Mucor genome for candidate RNases to participate in this pathway, and functional analysis of the corresponding knockout mutants, identified a new protein, R3B2. This RNase III-like protein presents unique domain architecture, it is specifically found in basal fungi and, besides its relevant role in the rdrp-dependent dicer-independent pathway, it is also involved in the canonical dicer-dependent RNAi pathway, highlighting its crucial role in the biogenesis and function of regulatory esRNAs. The involvement of RdRPs in RNA degradation could represent the first evolutionary step towards the development of an RNAi mechanism and constitutes a genetic link between mRNA degradation and post-transcriptional gene silencing.

Most eukaryotic organisms produce different classes of endogenous small RNA (esRNA) molecules that suppress gene expression through RNA interference (RNAi) pathways. These pathways, which may differ among organisms, are normally involved in genome defense, heterochromatin formation and regulation of genes involved in multiple cellular functions. In the basal fungus Mucor circinelloides, an opportunistic human pathogen, we previously demonstrated that biogenesis of a large group of esRNA molecules requires a basic RNAi machinery consisting of a Dicer-like protein, an Argonaute nuclease and two RNA-dependent RNA polymerases. This canonical dicer-dependent pathway regulates different cellular processes, such as vegetative sporulation. Besides those esRNAs generated by this canonical RNAi pathway, we have identified a new rdrp-dependent dicer-independent esRNA class. These esRNAs are produced by a degradation pathway in which the RdRP proteins signal specific transcripts that will be degraded by a newly identified RNase. This RNase, named R3B2, presents unique domain architecture, can only be found in basal fungi and it is also involved in the canonical dicer-dependent RNAi pathway. Our results expand the role of RdRPs in gene silencing and reveal the involvement of these proteins in a new RNA degradation process that could represent the first step in the evolution of RNAi.

An improved protocol for small RNA library construction using High Definition adapters

Next generation sequencing of small RNA (sRNA) libraries is widely used for studying sRNAs in various biological systems. However, cDNA libraries of sRNAs are biased for molecules that are ligated to adapters more or less efficiently than other molecules. One approach to reduce this ligation bias is to use a pool of adapters instead of a single adapter sequence, which allows many sRNAs to be ligated efficiently. We previously developed High Definition (HD) adapters for the Illumina sequencing platform, which contain degenerate nucleotides at the ligating ends of the adapters. However, the current commercial kits produced a large amount of 5’ adapter – 3’ adapter ligation product without the cDNA insert when HD adapters were used to replace the kit adapters. Here, we report a protocol to generate sRNA libraries using HD adapters with greatly reduced proportion of adapter-adapter products due to the degradation of nonligated 3’ adapters. The libraries can be completed within two days and can be used for various biological and clinical samples. As examples for using this protocol, we constructed sRNA libraries using total RNA extracted from cultured mammalian cells and plant leaf tissue.

Regulatory role of cellular and viral microRNAs in insect-virus interactions

The discovery of microRNAs (miRNAs) and their functions has led to a paradigm shift in our understanding of the regulation of gene expression, adding an extra layer of complexity for the mechanisms of gene expression. Both cellular and virus encoded miRNAs play important roles in virus-host interactions that may affect virus replication and the outcome of infection. Recent developments in RNA-seq platforms and bioinformatics tools have accelerated the discovery of miRNAs, their targets, and a myriad of associated research in various species. Here, recent findings and developments in miRNA research pertinent to insect host-virus interactions are reviewed and analyzed.

Deep sequencing of virus-derived small interfering RNAs and RNA from viral particles shows highly similar mutational landscapes of a plant virus population

RNA viruses exist within a host as a population of mutant sequences, often referred to as quasispecies. Within a host, sequences of RNA viruses constitute several distinct but interconnected pools, such as RNA packed in viral particles, double-stranded RNA, and virus-derived small interfering RNAs. We aimed to test if the same representation of within-host viral population structure could be obtained by sequencing different viral sequence pools. Using ultradeep Illumina sequencing, the diversity of two coexisting Potato virus Y sequence pools present within a plant was investigated: RNA isolated from viral particles and virus-derived small interfering RNAs (the derivatives of a plant RNA silencing mechanism). The mutational landscape of the within-host virus population was highly similar between both pools, with no notable hotspots across the viral genome. Notably, all of the single-nucleotide polymorphisms with a frequency of higher than 1.6% were found in both pools. Some unique single-nucleotide polymorphisms (SNPs) with very low frequencies were found in each of the pools, with more of them occurring in the small RNA (sRNA) pool, possibly arising through genetic drift in localized virus populations within a plant and the errors introduced during the amplification of silencing signal. Sequencing of the viral particle pool enhanced the efficiency of consensus viral genome sequence reconstruction. Nonhomologous recombinations were commonly detected in the viral particle pool, with a hot spot in the 3' untranslated and coat protein regions of the genome. We stress that they present an important but often overlooked aspect of virus population diversity.

IMPORTANCE

This study is the most comprehensive whole-genome characterization of a within-plant virus population to date and the first study comparing diversity of different pools of viral sequences within a host. We show that both virus-derived small RNAs and RNA from viral particles could be used for diversity assessment of within-plant virus population, since they show a highly congruent portrayal of the virus mutational landscape within a plant. The study is an important baseline for future studies of virus population dynamics, for example, during the adaptation to a new host. The comparison of the two virus sequence enrichment techniques, sequencing of virus-derived small interfering RNAs and RNA from purified viral particles, shows the strength of the latter for the detection of recombinant viral genomes and reconstruction of complete consensus viral genome sequence.

Bioinformatics approaches for viral metagenomics in plants using short RNAs: model case of study and application to a Cicer arietinum population

Over the past years deep sequencing experiments have opened novel doors to reconstruct viral populations in a high-throughput and cost-effective manner. Currently a substantial number of studies have been performed which employ next generation sequencing techniques to either analyze known viruses by means of a reference-guided approach or to discover novel viruses using a de novo-based strategy. Taking advantage of the well-known Cymbidium ringspot virus we have carried out a comparison of different bioinformatics tools to reconstruct the viral genome based on 21–27 nt short (s)RNA sequencing with the aim to identify the most efficient pipeline. The same approach was applied to a population of plants constituting an ancient variety of Cicer arietinum with red seeds. Among the discovered viruses, we describe the presence of a Tobamovirus referring to the Tomato mottle mosaic virus (NC_022230), which was not yet observed on C. arietinum nor revealed in Europe and a viroid referring to Hop stunt viroid (NC_001351.1) never reported in chickpea. Notably, a reference sequence guided approach appeared the most efficient in such kind of investigation. Instead, the de novo assembly reached a non-appreciable coverage although the most prominent viral species could still be identified. Advantages and limitations of viral metagenomics analysis using sRNAs are discussed.

Computational Prediction of miRNA Genes from Small RNA Sequencing Data

Next-generation sequencing now for the first time allows researchers to gage the depth and variation of entire transcriptomes. However, now as rare transcripts can be detected that are present in cells at single copies, more advanced computational tools are needed to accurately annotate and profile them. microRNAs (miRNAs) are 22 nucleotide small RNAs (sRNAs) that post-transcriptionally reduce the output of protein coding genes. They have established roles in numerous biological processes, including cancers and other diseases. During miRNA biogenesis, the sRNAs are sequentially cleaved from precursor molecules that have a characteristic hairpin RNA structure. The vast majority of new miRNA genes that are discovered are mined from small RNA sequencing (sRNA-seq), which can detect more than a billion RNAs in a single run. However, given that many of the detected RNAs are degradation products from all types of transcripts, the accurate identification of miRNAs remain a non-trivial computational problem. Here, we review the tools available to predict animal miRNAs from sRNA sequencing data. We present tools for generalist and specialist use cases, including prediction from massively pooled data or in species without reference genome. We also present wet-lab methods used to validate predicted miRNAs, and approaches to computationally benchmark prediction accuracy. For each tool, we reference validation experiments and benchmarking efforts. Last, we discuss the future of the field.

Copy number variable microRNAs in schizophrenia and their neurodevelopmental gene targets

BACKGROUND

MicroRNAs (miRNAs) are key regulators of gene expression in the human genome and may contribute to risk for neuropsychiatric disorders. miRNAs play an acknowledged role in the strongest of genetic risk factors for schizophrenia, 22q11.2 deletions. We hypothesized that in schizophrenia there would be an enrichment of other rare copy number variants (CNVs) that overlap miRNAs.

METHODS

Using high-resolution genome-wide microarrays and rigorous methods, we compared the miRNA content of rare CNVs in well-characterized cohorts of schizophrenia cases (n = 420) and comparison subjects, excluding 22q11.2 CNVs. We also performed a gene-set enrichment analysis of the predicted miRNA target genes.

RESULTS

The schizophrenia group was enriched for the proportion of individuals with a rare CNV overlapping a miRNA (3.29-fold increase over comparison subjects, p < .0001). The presence of a rare CNV overlapping a miRNA remained a significant predictor of schizophrenia case status (p = .0072) in a multivariate logistic regression model correcting for total CNV size. In contrast, comparable analyses correcting for CNV size showed no enrichment of rare CNVs overlapping protein-coding genes. A gene-set enrichment analysis indicated that predicted target genes of recurrent CNV-overlapped miRNAs in schizophrenia may be functionally enriched for neurodevelopmental processes, including axonogenesis and neuron projection development. Predicted gene targets driving these results included CAPRIN1, NEDD4, NTRK2, PAK2, RHOA, and SYNGAP1.

CONCLUSIONS

These data are the first to demonstrate a genome-wide role for CNVs overlapping miRNAs in the genetic risk for schizophrenia. The results provide support for an expanded multihit model of causation, with potential implications for miRNA-based therapeutics.

Drawing siRNAs of viral origin out from plant siRNAs libraries

Viruses are obligate intracellular entities that infect all forms of life. In plants, invading viral nucleic acids trigger RNA silencing machinery and it results in the accumulation of viral short interfering RNAs (v-siRNAs). The study of v-siRNAs population in biological samples has become a major part of many research projects aiming to identify viruses infecting them, including unknown viruses, even at extremely low titer. Currently, siRNA populations are investigated by high-throughput sequencing approaches, which generate very large data sets. The major difficulty in these studies is to properly analyze such huge amount of data. In this regard, easy-to-use bioinformatics tools to groom and decipher siRNA libraries and to draw out v-siRNAs are needed. Here we describe a workflow, which permit users with little experience in bioinformatics to draw out v-siRNAs from raw data sequences obtained by Illumina technology. Such pipeline has been released in the context of Galaxy, an open source Web-based platform for bioinformatics analyses.

An optimized kit-free method for making strand-specific deep sequencing libraries from RNA fragments

Deep sequencing of strand-specific cDNA libraries is now a ubiquitous tool for identifying and quantifying RNAs in diverse sample types. The accuracy of conclusions drawn from these analyses depends on precise and quantitative conversion of the RNA sample into a DNA library suitable for sequencing. Here, we describe an optimized method of preparing strand-specific RNA deep sequencing libraries from small RNAs and variably sized RNA fragments obtained from ribonucleoprotein particle footprinting experiments or fragmentation of long RNAs. Our approach works across a wide range of input amounts (400 pg to 200 ng), is easy to follow and produces a library in 2–3 days at relatively low reagent cost, all while giving the user complete control over every step. Because all enzymatic reactions were optimized and driven to apparent completion, sequence diversity and species abundance in the input sample are well preserved.

Interactions of Rice tungro bacilliform pararetrovirus and its protein P4 with plant RNA-silencing machinery

Small interfering RNA (siRNA)-directed gene silencing plays a major role in antiviral defense. Virus-derived siRNAs inhibit viral replication in infected cells and potentially move to neighboring cells, immunizing them from incoming virus. Viruses have evolved various ways to evade and suppress siRNA production or action. Here, we show that 21-, 22-, and 24-nucleotide (nt) viral siRNAs together constitute up to 19% of total small RNA population of Oryza sativa plants infected with Rice tungro bacilliform virus (RTBV) and cover both strands of the RTBV DNA genome. However, viral siRNA hotspots are restricted to a short noncoding region between transcription and reverse-transcription start sites. This region generates double-stranded RNA (dsRNA) precursors of siRNAs and, in pregenomic RNA, forms a stable secondary structure likely inaccessible to siRNA-directed cleavage. In transient assays, RTBV protein P4 suppressed cell-to-cell spread of silencing but enhanced cell-autonomous silencing, which correlated with reduced 21-nt siRNA levels and increased 22-nt siRNA levels. Our findings imply that RTBV generates decoy dsRNA that restricts siRNA production to the structured noncoding region and thereby protects other regions of the viral genome from repressive action of siRNAs, while the viral protein P4 interferes with cell-to-cell spread of antiviral silencing.

MicroRNA Dysregulation, Gene Networks, and Risk for Schizophrenia in 22q11.2 Deletion Syndrome

The role of microRNAs (miRNAs) in the etiology of schizophrenia is increasingly recognized. Microdeletions at chromosome 22q11.2 are recurrent structural variants that impart a high risk for schizophrenia and are found in up to 1% of all patients with schizophrenia. The 22q11.2 deletion region overlaps gene DGCR8, encoding a subunit of the miRNA microprocessor complex. We identified miRNAs overlapped by the 22q11.2 microdeletion and for the first time investigated their predicted target genes, and those implicated by DGCR8, to identify targets that may be involved in the risk for schizophrenia. The 22q11.2 region encompasses seven validated or putative miRNA genes. Employing two standard prediction tools, we generated sets of predicted target genes. Functional enrichment profiles of the 22q11.2 region miRNA target genes suggested a role in neuronal processes and broader developmental pathways. We then constructed a protein interaction network of schizophrenia candidate genes and interaction partners relevant to brain function, independent of the 22q11.2 region miRNA mechanisms. We found that the predicted gene targets of the 22q11.2 deletion miRNAs, and targets of the genome-wide miRNAs predicted to be dysregulated by DGCR8 hemizygosity, were significantly represented in this schizophrenia network. The findings provide new insights into the pathway from 22q11.2 deletion to expression of schizophrenia, and suggest that hemizygosity of the 22q11.2 region may have downstream effects implicating genes elsewhere in the genome that are relevant to the general schizophrenia population. These data also provide further support for the notion that robust genetic findings in schizophrenia may converge on a reasonable number of final pathways.

Small RNAs from plants, bacteria and fungi within the order Hypocreales are ubiquitous in human plasma

BACKGROUND

The human microbiome plays a significant role in maintaining normal physiology. Changes in its composition have been associated with bowel disease, metabolic disorders and atherosclerosis. Sequences of microbial origin have been observed within small RNA sequencing data obtained from blood samples. The aim of this study was to characterise the microbiome from which these sequences are derived.

RESULTS

Abundant non-human small RNA sequences were identified in plasma and plasma exosomal samples. Assembly of these short sequences into longer contigs was the pivotal novel step in ascertaining their origin by BLAST searches. Most reads mapped to rRNA sequences. The taxonomic profiles of the microbes detected were very consistent between individuals but distinct from microbiomes reported at other sites. The majority of bacterial reads were from the phylum Proteobacteria, whilst for 5 of 6 individuals over 90% of the more abundant fungal reads were from the phylum Ascomycota; of these over 90% were from the order Hypocreales. Many contigs were from plants, presumably of dietary origin. In addition, extremely abundant small RNAs derived from human Y RNAs were detected.

CONCLUSIONS

A characteristic profile of a subset of the human microbiome can be obtained by sequencing small RNAs present in the blood. The source and functions of these molecules remain to be determined, but the specific profiles are likely to reflect health status. The potential to provide biomarkers of diet and for the diagnosis and prognosis of human disease is immense.

Evolutionary dynamics of coding and non-coding transcriptomes

Gene expression changes may underlie much of phenotypic evolution. The development of high-throughput RNA sequencing protocols has opened the door to unprecedented large-scale and cross-species transcriptome comparisons by allowing accurate and sensitive assessments of transcript sequences and expression levels. Here, we review the initial wave of the new generation of comparative transcriptomic studies in mammals and vertebrate outgroup species in the context of earlier work. Together with various large-scale genomic and epigenomic data, these studies have unveiled commonalities and differences in the dynamics of gene expression evolution for various types of coding and non-coding genes across mammalian lineages, organs, developmental stages, chromosomes and sexes. They have also provided intriguing new clues to the regulatory basis and phenotypic implications of evolutionary gene expression changes.

MicroRNAs influence reproductive responses by females to male sex peptide in Drosophila melanogaster

Across taxa, female behavior and physiology change significantly following the receipt of ejaculate molecules during mating. For example, receipt of sex peptide (SP) in female Drosophila melanogaster significantly alters female receptivity, egg production, lifespan, hormone levels, immunity, sleep, and feeding patterns. These changes are underpinned by distinct tissue- and time-specific changes in diverse sets of mRNAs. However, little is yet known about the regulation of these gene expression changes, and hence the potential role of microRNAs (miRNAs), in female postmating responses. A preliminary screen of genomic responses in females to receipt of SP suggested that there were changes in the expression of several miRNAs. Here we tested directly whether females lacking four of the candidate miRNAs highlighted (miR-279, miR-317, miR-278, and miR-184) showed altered fecundity, receptivity, and lifespan responses to receipt of SP, when mated once or continually to SP null or control males. The results showed that miRNA-lacking females mated to SP null males exhibited altered receptivity, but not reproductive output, in comparison to controls. However, these effects interacted significantly with the genetic background of the miRNA-lacking females. No significant survival effects were observed in miRNA-lacking females housed continually with SP null or control males. However, continual exposure to control males that transferred SP resulted in significantly higher variation in miRNA-lacking female lifespan than did continual exposure to SP null males. The results provide the first insight into the effects and importance of miRNAs in regulating postmating responses in females.

High-efficiency RNA cloning enables accurate quantification of miRNA expression by deep sequencing

Small RNA cloning and sequencing is uniquely positioned as a genome-wide approach to quantify miRNAs with single-nucleotide resolution. However, significant biases introduced by RNA ligation in current protocols lead to inaccurate miRNA quantification by 1000-fold. Here we report an RNA cloning method that achieves over 95% efficiency for both 5′ and 3′ ligations. It achieves accurate quantification of synthetic miRNAs with less than two-fold deviation from the anticipated value and over a dynamic range of four orders of magnitude. Taken together, this high-efficiency RNA cloning method permits accurate genome-wide miRNA profiling from total RNAs.

Induction and suppression of tick cell antiviral RNAi responses by tick-borne flaviviruses

Arboviruses are transmitted by distantly related arthropod vectors such as mosquitoes (class Insecta) and ticks (class Arachnida). RNA interference (RNAi) is the major antiviral mechanism in arthropods against arboviruses. Unlike in mosquitoes, tick antiviral RNAi is not understood, although this information is important to compare arbovirus/host interactions in different classes of arbovirus vectos. Using an Ixodes scapularis-derived cell line, key Argonaute proteins involved in RNAi and the response against tick-borne Langat virus (Flaviviridae) replication were identified and phylogenetic relationships characterized. Analysis of small RNAs in infected cells showed the production of virus-derived small interfering RNAs (viRNAs), which are key molecules of the antiviral RNAi response. Importantly, viRNAs were longer (22 nucleotides) than those from other arbovirus vectors and mapped at highest frequency to the termini of the viral genome, as opposed to mosquito-borne flaviviruses. Moreover, tick-borne flaviviruses expressed subgenomic flavivirus RNAs that interfere with tick RNAi. Our results characterize the antiviral RNAi response in tick cells including phylogenetic analysis of genes encoding antiviral proteins, and viral interference with this pathway. This shows important differences in antiviral RNAi between the two major classes of arbovirus vectors, and our data broadens our understanding of arthropod antiviral RNAi.

Small RNA profile in moso bamboo root and leaf obtained by high definition adapters

Moso bamboo (Phyllostachy heterocycla cv. pubescens L.) is an economically important fast-growing tree. In order to gain better understanding of gene expression regulation in this important species we used next generation sequencing to profile small RNAs in leaf and roots of young seedlings. Since standard kits to produce cDNA of small RNAs are biased for certain small RNAs, we used High Definition adapters that reduce ligation bias. We identified and experimentally validated five new microRNAs and a few other small non-coding RNAs that were not microRNAs. The biological implication of microRNA expression levels and targets of microRNAs are discussed.

An efficient and sensitive method for preparing cDNA libraries from scarce biological samples

The preparation and high-throughput sequencing of cDNA libraries from samples of small RNA is a powerful tool to quantify known small RNAs (such as microRNAs) and to discover novel RNA species. Interest in identifying the small RNA repertoire present in tissues and in biofluids has grown substantially with the findings that small RNAs can serve as indicators of biological conditions and disease states. Here we describe a novel and straightforward method to clone cDNA libraries from small quantities of input RNA. This method permits the generation of cDNA libraries from sub-picogram quantities of RNA robustly, efficiently and reproducibly. We demonstrate that the method provides a significant improvement in sensitivity compared to previous cloning methods while maintaining reproducible identification of diverse small RNA species. This method should have widespread applications in a variety of contexts, including biomarker discovery from scarce samples of human tissue or body fluids.

Evaluating bias-reducing protocols for RNA sequencing library preparation

Background

Next-generation sequencing does not yield fully unbiased estimates for read abundance, which may impact on the conclusions that can be drawn from sequencing data. The ligation step in RNA sequencing library generation is a known source of bias, motivating developments in enzyme technology and library construction protocols. We present the first comparison of the standard duplex adaptor protocol supplied by Life Technologies for use on the Ion Torrent PGM with an alternate single adaptor approach involving CircLigase (CircLig protocol).

A correlation between over-representation in sequenced libraries and degree of secondary structure has been reported previously, therefore we also investigated whether bias could be reduced by ligation with an enzyme that functions at a temperature not permissive for such structure.

Results

A pool of small RNA fragments of known composition was converted into a sequencing library using one of three protocols and sequenced on an Ion Torrent PGM. The CircLig protocol resulted in less over-representation of specific sequences than the standard protocol. Over-represented sequences are more likely to be predicted to have secondary structure and to co-fold with adaptor sequences. However, use of the thermostable ligase Methanobacterium thermoautotrophicum RNA ligase K97A (Mth K97A) was not sufficient to reduce bias.

Conclusions

The single adaptor CircLigase-based approach significantly reduces, but does not eliminate, bias in Ion Torrent data. Ligases that function at temperatures to remove the possible influence of secondary structure on library generation may be of value, although Mth K97A is not effective in this case.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-569) contains supplementary material, which is available to authorized users.

Conserved microRNA editing in mammalian evolution, development and disease

Background

Mammalian microRNAs (miRNAs) are sometimes subject to adenosine-to-inosine RNA editing, which can lead to dramatic changes in miRNA target specificity or expression levels. However, although a few miRNAs are known to be edited at identical positions in human and mouse, the evolution of miRNA editing has not been investigated in detail. In this study, we identify conserved miRNA editing events in a range of mammalian and non-mammalian species.

Results

We demonstrate deep conservation of several site-specific miRNA editing events, including two that date back to the common ancestor of mammals and bony fishes some 450 million years ago. We also find evidence of a recent expansion of an edited miRNA family in placental mammals and show that editing of these miRNAs is associated with changes in target mRNA expression during primate development and aging. While global patterns of miRNA editing tend to be conserved across species, we observe substantial variation in editing frequencies depending on tissue, age and disease state: editing is more frequent in neural tissues compared to heart, kidney and testis; in older compared to younger individuals; and in samples from healthy tissues compared to tumors, which together suggests that miRNA editing might be associated with a reduced rate of cell proliferation.

Conclusions

Our results show that site-specific miRNA editing is an evolutionarily conserved mechanism, which increases the functional diversity of mammalian miRNA transcriptomes. Furthermore, we find that although miRNA editing is rare compared to editing of long RNAs, miRNAs are greatly overrepresented among conserved editing targets.

Identification of tubular injury microRNA biomarkers in urine: comparison of next-generation sequencing and qPCR-based profiling platforms

Background

MicroRNAs (miRNAs) are small, non-coding RNAs that regulate protein levels post-transcriptionally. miRNAs play important regulatory roles in many cellular processes and have been implicated in several diseases. Recent studies have reported significant levels of miRNAs in a variety of body fluids, raising the possibility that miRNAs could serve as useful biomarkers. Next-generation sequencing (NGS) is increasingly employed in biomedical investigations. Although concordance between this platform and qRT-PCR based assays has been reported in high quality specimens, information is lacking on comparisons in biofluids especially urine. Here we describe the changes in miRNA expression patterns in a rodent model of renal tubular injury (gentamicin). Our aim is to compare RNA sequencing and qPCR based miRNA profiling in urine specimen from control and rats with confirmed tubular injury.

Results

Our preliminary examination of the concordance between miRNA-seq and qRT-PCR in urine specimen suggests minimal agreement between platforms probably due to the differences in sensitivity. Our results suggest that although miRNA-seq has superior specificity, it may not detect low abundant miRNAs in urine samples. Specifically, miRNA-seq did not detect some sequences which were identified by qRT-PCR. On the other hand, the qRT-PCR analysis was not able to detect the miRNA isoforms, which made up the majority of miRNA changes detected by NGS.

Conclusions

To our knowledge, this is the first time that miRNA profiling platforms including NGS have been compared in urine specimen. miRNAs identified by both platforms, let-7d, miR-203, and miR-320, may potentially serve as promising novel urinary biomarkers for drug induced renal tubular epithelial injury.

Importance of the RNA secondary structure for the relative accumulation of clustered viral microRNAs

Micro (mi)RNAs are small non-coding RNAs with key regulatory functions. Recent advances in the field allowed researchers to identify their targets. However, much less is known regarding the regulation of miRNAs themselves. The accumulation of these tiny regulators can be modulated at various levels during their biogenesis from the transcription of the primary transcript (pri-miRNA) to the stability of the mature miRNA. Here, we studied the importance of the pri-miRNA secondary structure for the regulation of mature miRNA accumulation. To this end, we used the Kaposi's sarcoma herpesvirus, which encodes a cluster of 12 pre-miRNAs. Using small RNA profiling and quantitative northern blot analysis, we measured the absolute amount of each mature miRNAs in different cellular context. We found that the difference in expression between the least and most expressed viral miRNAs could be as high as 60-fold. Using high-throughput selective 2′-hydroxyl acylation analyzed by primer extension, we then determined the secondary structure of the long primary transcript. We found that highly expressed miRNAs derived from optimally structured regions within the pri-miRNA. Finally, we confirmed the importance of the local structure by swapping stem-loops or by targeted mutagenesis of selected miRNAs, which resulted in a perturbed accumulation of the mature miRNA.

Small RNA cloning and sequencing strategy affects host and viral microRNA expression signatures

The establishment of the microRNA (miRNA) expression signatures is the basic element to investigate the role played by these regulatory molecules in the biology of an organism. Marek's disease virus 1 (MDV-1) is an avian herpesvirus that naturally infects chicken and induces T cells lymphomas. During latency, MDV-1, like other herpesviruses, expresses a limited subset of transcripts. These include three miRNA clusters. Several studies identified the expression of virus and host encoded miRNAs from MDV-1 infected cell cultures and chickens. But a high discrepancy was observed when miRNA cloning frequencies obtained from different cloning and sequencing protocols were compared. Thus, we analyzed the effect of small RNA library preparation and sequencing on the miRNA frequencies obtained from the same RNA samples collected during MDV-1 infection of chicken at different steps of the oncoviral pathogenesis. Qualitative and quantitative variations were found in the data, depending on the strategy used. One of the mature miRNA derived from the latency-associated-transcript (LAT), mdv1-miR-M7-5p, showed the highest variation. Its cloning frequency was 50% of the viral miRNA counts when a small scale sequencing approach was used. Its frequency was 100 times less abundant when determined through the deep sequencing approach. Northern blot analysis showed a better correlation with the miRNA frequencies found by the small scale sequencing approach. By analyzing the cellular miRNA repertoire, we also found a gap between the two sequencing approaches. Collectively, our study indicates that next-generation sequencing data considered alone are limited for assessing the absolute copy number of transcripts. Thus, the quantification of small RNA should be addressed by compiling data obtained by using different techniques such as microarrays, qRT-PCR and NB analysis in support of high throughput sequencing data. These observations should be considered when miRNA variations are studied prior addressing functional studies.

Elimination of ligation dependent artifacts in T4 RNA ligase to achieve high efficiency and low bias microRNA capture

Adapter ligation is a critical first step in many microRNA analysis methods including microarray, qPCR, and sequencing. Previous studies have shown that ligation bias can have dramatic effects on both the fidelity of expression profiles and reproducibility across samples. We have developed a method for high efficiency and low bias microRNA capture by 3′ adapter ligation using T4 RNA ligase that does not require pooled adapters. Using a panel of 20 microRNA, we investigated the effects of ligase type, PEG concentration, ligase amount, adapter concentration, incubation time, incubation temperature, and adapter design on capture efficiency and bias. Of these factors, high PEG% was found to be critical in suppressing ligation bias. We obtained high average capture efficiency and low CV across the 20 microRNA panel, both in idealized buffer conditions (86%±10%) and total RNA spiking conditions (64%±17%). We demonstrate that this method is reliable across microRNA species that previous studies have had difficulty capturing and that our adapter design performs significantly better than the common adapter designs. Further, we demonstrate that the optimization methodology must be specifically designed for minimizing bias in order to obtain the ideal reaction parameters.

A meta-analysis revealed insights into the sources, conservation and impact of microRNA 5'-isoforms in four model species

MicroRNA (miRNA) 5'-isoforms, or 5'-isomiRs, are small-RNA species that originate from the same genomic loci as the major miRNAs with their 5' ends shifted from the 5' ends of the miRNAs by a few nucleotides. Although 5'-isomiRs have been reported, their origins, properties and potential functions remain to be examined. We systematically studied 5'-isomiRs in human, mouse, fruitfly and worm by analysing a large collection of small non-coding RNA and mRNA profiling data. The results revealed a broad existence of 5'-isomiRs in the four species, many of which were conserved and could arise from genomic loci of canonical and non-canonical miRNAs. The well-conserved 5'-isomiRs have several features, including a preference of the 3p over the 5p arms of hairpins of conserved mammalian miRNAs, altered 5'-isomiRs across species and across tissues, and association with structural variations of miRNA hairpins. Importantly, 5'-isomiRs and their major miRNAs may have different mRNA targets and thus potentially play distinct roles of gene regulation, as shown by an integrative analysis combining miRNA and mRNA profiling data from psoriatic and normal human skin and from murine miRNA knockout assays. Indeed, 18 5'-isomiRs had aberrant expression in psoriatic human skin, suggesting their potential function in psoriasis pathogenesis. The results of the current study deepened our understanding of the diversity and conservation of miRNAs, their plasticity in gene regulation and potential broad function in complex diseases.

Library construction for next-generation sequencing: overviews and challenges

High-throughput sequencing, also known as next-generation sequencing (NGS), has revolutionized genomic research. In recent years, NGS technology has steadily improved, with costs dropping and the number and range of sequencing applications increasing exponentially. Here, we examine the critical role of sequencing library quality and consider important challenges when preparing NGS libraries from DNA and RNA sources. Factors such as the quantity and physical characteristics of the RNA or DNA source material as well as the desired application (i.e., genome sequencing, targeted sequencing, RNA-seq, ChIP-seq, RIP-seq, and methylation) are addressed in the context of preparing high quality sequencing libraries. In addition, the current methods for preparing NGS libraries from single cells are also discussed.

Library preparation methods for next-generation sequencing: tone down the bias

Next-generation sequencing (NGS) has caused a revolution in biology. NGS requires the preparation of libraries in which (fragments of) DNA or RNA molecules are fused with adapters followed by PCR amplification and sequencing. It is evident that robust library preparation methods that produce a representative, non-biased source of nucleic acid material from the genome under investigation are of crucial importance. Nevertheless, it has become clear that NGS libraries for all types of applications contain biases that compromise the quality of NGS datasets and can lead to their erroneous interpretation. A detailed knowledge of the nature of these biases will be essential for a careful interpretation of NGS data on the one hand and will help to find ways to improve library quality or to develop bioinformatics tools to compensate for the bias on the other hand. In this review we discuss the literature on bias in the most common NGS library preparation protocols, both for DNA sequencing (DNA-seq) as well as for RNA sequencing (RNA-seq). Strikingly, almost all steps of the various protocols have been reported to introduce bias, especially in the case of RNA-seq, which is technically more challenging than DNA-seq. For each type of bias we discuss methods for improvement with a view to providing some useful advice to the researcher who wishes to convert any kind of raw nucleic acid into an NGS library.

Identification and characterization of miRNAome in root, stem, leaf and tuber developmental stages of potato (Solanum tuberosum L.) by high-throughput sequencing

BACKGROUND

MicroRNAs (miRNAs) are ubiquitous components of endogenous plant transcriptome. miRNAs are small, single-stranded and ~21 nt long RNAs which regulate gene expression at the post-transcriptional level and are known to play essential roles in various aspects of plant development and growth. Previously, a number of miRNAs have been identified in potato through in silico analysis and deep sequencing approach. However, identification of miRNAs through deep sequencing approach was limited to a few tissue types and developmental stages. This study reports the identification and characterization of potato miRNAs in three different vegetative tissues and four stages of tuber development by high throughput sequencing.

RESULTS

Small RNA libraries were constructed from leaf, stem, root and four early developmental stages of tuberization and subjected to deep sequencing, followed by bioinformatics analysis. A total of 89 conserved miRNAs (belonging to 33 families), 147 potato-specific miRNAs (with star sequence) and 112 candidate potato-specific miRNAs (without star sequence) were identified. The digital expression profiling based on TPM (Transcripts Per Million) and qRT-PCR analysis of conserved and potato-specific miRNAs revealed that some of the miRNAs showed tissue specific expression (leaf, stem and root) while a few demonstrated tuberization stage-specific expressions. Targets were predicted for identified conserved and potato-specific miRNAs, and predicted targets of four conserved miRNAs, miR160, miR164, miR172 and miR171, which are ARF16 (Auxin Response Factor 16), NAM (NO APICAL MERISTEM), RAP1 (Relative to APETALA2 1) and HAM (HAIRY MERISTEM) respectively, were experimentally validated using 5' RLM-RACE (RNA ligase mediated rapid amplification of cDNA ends). Gene ontology (GO) analysis for potato-specific miRNAs was also performed to predict their potential biological functions.

CONCLUSIONS

We report a comprehensive study of potato miRNAs at genome-wide level by high-throughput sequencing and demonstrate that these miRNAs have tissue and/or developmental stage-specific expression profile. Also, predicted targets of conserved miRNAs were experimentally confirmed for the first time in potato. Our findings indicate the existence of extensive and complex small RNA population in this crop and suggest their important role in pathways involved in diverse biological processes, including tuber development.

RNA-binding proteins in regulation of alternative cleavage and polyadenylation

Almost all eukaryotic pre-mRNAs are processed at the 3' end by the cleavage and polyadenylation (C/P) reaction, which preludes termination of transcription and gives rise to the poly(A) tail of mature mRNA. Genomic studies in recent years have indicated that most eukaryotic mRNA genes have multiple cleavage and polyadenylation sites (pAs), leading to alternative cleavage and polyadenylation (APA) products. APA isoforms generally differ in their 3' untranslated regions (3' UTRs), but can also have different coding sequences (CDSs). APA expands the repertoire of transcripts expressed from the genome, and is highly regulated under various physiological and pathological conditions. Growing lines of evidence have shown that RNA-binding proteins (RBPs) play important roles in regulation of APA. Some RBPs are part of the machinery for C/P; others influence pA choice through binding to adjacent regions. In this chapter, we review cis elements and trans factors involved in C/P, the significance of APA, and increasingly elucidated roles of RBPs in APA regulation. We also discuss analysis of APA using transcriptome-wide techniques as well as molecular biology approaches.

Reverse transcription through a bulky triazole linkage in RNA: implications for RNA sequencing

Reverse-transcriptase reads through the bulky artificial triazole linkage with omission of one nucleotide with implications for RNA isolation and analysis.

Discovery of novel small RNAs in the quest to unravel genome complexity

The different RNA populations within the cell, or the 'RNAome', form a complex and integral layer of gene regulation. In particular, small non-coding RNAs such as microRNAs have been shown to regulate as much as half of all genes expressed in mammals. The recent discovery of novel small RNAs derived from tRNAs, snoRNAs and Y RNAs, as well as the development of techniques that can reduce sequencing bias of these molecules, is slowly helping us to unveil the secrets of the genome.

Identification of RNase L-dependent, 3'-end-modified, viral small RNAs in Sindbis virus-infected mammalian cells

Small RNAs play a critical role in host-pathogen interaction. Indeed, small RNA-mediated silencing or RNA interference (RNAi) is one of the earliest forms of antiviral immunity. Although it represents the main defense system against viruses in many organisms, the antiviral role of RNAi has not been clearly proven in higher vertebrates. However, it is well established that their response to viral infection relies on the recognition of viral RNAs by host pattern recognition receptors (PRRs) to trigger activation of the interferon pathway. In the present work, we report the existence of a novel small noncoding RNA population produced in mammalian cells upon RNA virus infection. Using Sindbis virus (SINV) as a prototypic arbovirus model, we profiled the small RNA population of infected cells in both human and African green monkey cell lines. Here, we provide evidence for the presence of discrete small RNAs of viral origin that are not associated with the RNA-induced silencing complex (RISC), that are highly expressed and detected by Northern blot analysis, and that accumulate as 21- to 28-nucleotide (nt) species during infection. We report that the cellular antiviral endoribonuclease RNase L cleaves the viral genome, producing in turn the small RNAs. Surprisingly, we uncovered the presence of a modification on the 3′-end nucleotide of SINV-derived viral small RNAs (SvsRNAs) that might be at the origin of their stability. Altogether, our findings show that stable modified small viral RNAs could represent a novel way to modulate host-virus interaction upon SINV infection.

In a continuous arms race, viruses have to deal with host antiviral responses in order to successfully establish an infection. In mammalian cells, the host defense mechanism relies on the recognition of viral RNAs, resulting in the activation of type I interferons (IFNs). In turn, the expression of many interferon-stimulated genes (ISGs) is induced to inhibit viral replication. Here we report that the cytoplasmic, interferon-induced, cellular endoribonuclease RNase L is involved in the accumulation of a novel small RNA population of viral origin. These small RNAs are produced upon SINV infection of mammalian cells and are stabilized by a 3′-end modification. Altogether, our findings indicate that in our system RNA silencing is not active against Sindbis virus (SINV) and might open the way to a better understanding of the antiviral response mediated by a novel class of small RNAs.

Biases in small RNA deep sequencing data

High-throughput RNA sequencing (RNA-seq) is considered a powerful tool for novel gene discovery and fine-tuned transcriptional profiling. The digital nature of RNA-seq is also believed to simplify meta-analysis and to reduce background noise associated with hybridization-based approaches. The development of multiplex sequencing enables efficient and economic parallel analysis of gene expression. In addition, RNA-seq is of particular value when low RNA expression or modest changes between samples are monitored. However, recent data uncovered severe bias in the sequencing of small non-protein coding RNA (small RNA-seq or sRNA-seq), such that the expression levels of some RNAs appeared to be artificially enhanced and others diminished or even undetectable. The use of different adapters and barcodes during ligation as well as complex RNA structures and modifications drastically influence cDNA synthesis efficacies and exemplify sources of bias in deep sequencing. In addition, variable specific RNA G/C-content is associated with unequal polymerase chain reaction amplification efficiencies. Given the central importance of RNA-seq to molecular biology and personalized medicine, we review recent findings that challenge small non-protein coding RNA-seq data and suggest approaches and precautions to overcome or minimize bias.

A study of small RNAs from cerebral neocortex of pathology-verified Alzheimer's disease, dementia with lewy bodies, hippocampal sclerosis, frontotemporal lobar dementia, and non-demented human controls

MicroRNAs (miRNAs) are small (20-22 nucleotides) regulatory non-coding RNAs that strongly influence gene expression. Most prior studies addressing the role of miRNAs in neurodegenerative diseases (NDs) have focused on individual diseases such as Alzheimer's disease (AD), making disease-to-disease comparisons impossible. Using RNA deep sequencing, we sought to analyze in detail the small RNAs (including miRNAs) in the temporal neocortex gray matter from non-demented controls (n = 2), AD (n = 5), dementia with Lewy bodies (n = 4), hippocampal sclerosis of aging (n = 4), and frontotemporal lobar dementia (FTLD) (n = 5) cases, together accounting for the most prevalent ND subtypes. All cases had short postmortem intervals, relatively high-quality RNA, and state-of-the-art neuropathological diagnoses. The resulting data (over 113 million reads in total, averaging 5.6 million reads per sample) and secondary expression analyses constitute an unprecedented look into the human cerebral cortical miRNome at a nucleotide resolution. While we find no apparent changes in isomiR or miRNA editing patterns in correlation with ND pathology, our results validate and extend previous miRNA profiling studies with regard to quantitative changes in NDs. In agreement with this idea, we provide independent cohort validation for changes in miR-132 expression levels in AD (n = 8) and FTLD (n = 14) cases when compared to controls (n = 8). The identification of common and ND-specific putative novel brain miRNAs and/or short-hairpin molecules is also presented. The challenge now is to better understand the impact of these and other alterations on neuronal gene expression networks and neuropathologies.

Divergent MicroRNA targetomes of closely related circulating strains of a polyomavirus

Hundreds of virus-encoded microRNAs (miRNAs) have been uncovered, but an in-depth functional understanding is lacking for most. A major challenge for the field is separating those miRNA targets that are biologically relevant from those that are not advantageous to the virus. Here, we show that miRNAs from related variants of the polyomavirus simian vacuolating virus 40 (SV40) have differing host target repertoires (targetomes) while their direct autoregulatory activity on virus-encoded early gene products is completely preserved. These results underscore the importance of miRNA-mediated viral gene autoregulation in some polyomavirus life cycles. More broadly, these findings imply that some host targets of virus-encoded miRNAs are likely to be of little selective advantage to the virus, and our approach provides a strategy for prioritizing relevant targets.

Interaction of small RNA-8105 and the intron of IbMYB1 RNA regulates IbMYB1 family genes through secondary siRNAs and DNA methylation after wounding

Small RNAs (sRNAs) play important roles in plants under stress conditions. However, limited research has been performed on the sRNAs involved in plant wound responses. In the present study, a novel wounding-induced sRNA, sRNA8105, was identified in sweet potato (Ipomoea batatas cv. Tainung 57) using microarray analysis. It was found that expression of sRNA8105 increased after mechanical wounding. Furthermore, Dicer-like 1 (DCL1) is required for the sRNA8105 precursor (pre-sRNA8105) to generate 22 and 24 nt mature sRNA8105. sRNA8105 targeted the first intron of IbMYB1 (MYB domain protein 1) before RNA splicing, and mediated RNA cleavage and DNA methylation of IbMYB1. The interaction between sRNA8105 and IbMYB1 was confirmed by cleavage site mapping, agro-infiltration analyses, and use of a transgenic sweet potato over-expressing pre-sRNA8105 gene. Induction of IbMYB1-siRNA was observed in the wild-type upon wounding and in transgenic sweet potato over-expressing pre-sRNA8105 gene without wounding, resulting in decreased expression of the whole IbMYB1 gene family, i.e. IbMYB1 and the IbMYB2 genes, and thus directing metabolic flux toward biosynthesis of lignin in the phenylpropanoid pathway. In conclusion, sRNA8105 induced by wounding binds to the first intron of IbMYB1 RNA to methylate IbMYB1, cleave IbMYB1 RNA, and trigger production of secondary siRNAs, further repressing the expression of the IbMYB1 family genes and regulating the phenylpropanoid pathway.

Functional specialization of the small interfering RNA pathway in response to virus infection

In Drosophila, post-transcriptional gene silencing occurs when exogenous or endogenous double stranded RNA (dsRNA) is processed into small interfering RNAs (siRNAs) by Dicer-2 (Dcr-2) in association with a dsRNA-binding protein (dsRBP) cofactor called Loquacious (Loqs-PD). siRNAs are then loaded onto Argonaute-2 (Ago2) by the action of Dcr-2 with another dsRBP cofactor called R2D2. Loaded Ago2 executes the destruction of target RNAs that have sequence complementarity to siRNAs. Although Dcr-2, R2D2, and Ago2 are essential for innate antiviral defense, the mechanism of virus-derived siRNA (vsiRNA) biogenesis and viral target inhibition remains unclear. Here, we characterize the response mechanism mediated by siRNAs against two different RNA viruses that infect Drosophila. In both cases, we show that vsiRNAs are generated by Dcr-2 processing of dsRNA formed during viral genome replication and, to a lesser extent, viral transcription. These vsiRNAs seem to preferentially target viral polyadenylated RNA to inhibit viral replication. Loqs-PD is completely dispensable for silencing of the viruses, in contrast to its role in silencing endogenous targets. Biogenesis of vsiRNAs is independent of both Loqs-PD and R2D2. R2D2, however, is required for sorting and loading of vsiRNAs onto Ago2 and inhibition of viral RNA expression. Direct injection of viral RNA into Drosophila results in replication that is also independent of Loqs-PD. This suggests that triggering of the antiviral pathway is not related to viral mode of entry but recognition of intrinsic features of virus RNA. Our results indicate the existence of a vsiRNA pathway that is separate from the endogenous siRNA pathway and is specifically triggered by virus RNA. We speculate that this unique framework might be necessary for a prompt and efficient antiviral response.

The RNA interference (RNAi) pathway utilizes small non-coding RNAs to silence gene expression. In insects, RNAi regulates endogenous genes and functions as an RNA-based immune system against viral infection. Here we have uncovered details of how RNAi is triggered by RNA viruses. Double-stranded RNA (dsRNA) generated as a replication intermediate or from transcription of the RNA virus can be used as substrate for the biogenesis of virus-derived small interfering RNAs (vsiRNAs). Unlike other dsRNAs, virus RNA processing involves Dicer but not its canonical partner protein Loqs-PD. Thus, vsiRNA biogenesis is mechanistically different from biogenesis of endogenous siRNAs or siRNAs derived from other exogenous RNA sources. Our results suggest a specialization of the pathway dedicated to silencing of RNA viruses versus other types of RNAi silencing. The understanding of RNAi mechanisms during viral infection could have implications for the control of insect-borne viruses and the use of siRNAs to treat viral infections in humans.

Small RNAs tackle large viruses: RNA interference-based antiviral defense against DNA viruses in insects

The antiviral RNA interference (RNAi) pathway processes viral double-stranded RNA (dsRNA) into viral small interfering RNAs (vsiRNA) that guide the recognition and cleavage of complementary viral target RNAs. In RNA virus infections, viral replication intermediates, dsRNA genomes or viral structured RNAs have been implicated as Dicer-2 substrates. In a recent publication, we demonstrated that a double-stranded DNA virus, Invertebrate iridescent virus 6, is a target of the Drosophila RNAi machinery, and we proposed that overlapping converging transcripts base pair to form the dsRNA substrates for vsiRNA biogenesis. Here, we discuss the role of RNAi in antiviral defense to DNA viruses in Drosophila and other invertebrate model systems.

Profiles of piRNA abundances at emerging or established piRNA loci are determined by local DNA sequences

Piwi-interacting RNAs (piRNAs) ensure transposable element silencing in Drosophila, thereby preserving genome integrity across generations. Primary piRNAs arise from the processing of long RNA transcripts produced in the germ line by a limited number of telomeric and pericentromeric loci. Primary piRNAs bound to the Argonaute protein Aubergine then drive the production of secondary piRNAs through the "ping-pong" amplification mechanism that involves an interplay with piRNAs bound to the Argonaute protein Argonaute-3. We recently discovered that clusters of P-element-derived transgenes produce piRNAs and mediate silencing of homologous target transgenes in the female germ line. We also demonstrated that some clusters are able to convert other homologous inactive transgene clusters into piRNA-producing loci, which then transmit their acquired silencing capacity over generations. This paramutation phenomenon is mediated by maternal inheritance of piRNAs homologous to the transgenes. Here we further mined our piRNA sequencing data sets generated from various strains carrying transgenes with partial sequence homology at distinct genomic sites. This analysis revealed that same sequences in different genomic contexts generate highly similar profiles of piRNA abundances. The strong tendency of piRNAs for bearing a U at their 5' end has long been recognized. Our observations support the notion that, in addition, the relative frequencies of Drosophila piRNAs are locally determined by the DNA sequence of piRNA loci.

Distinctive profile of IsomiR expression and novel microRNAs in rat heart left ventricle

MicroRNAs (miRNAs) are single-stranded non-coding RNAs that negatively regulate target gene expression through mRNA cleavage or translational repression. There is mounting evidence that they play critical roles in heart disease. The expression of known miRNAs in the heart has been studied at length by microarray and quantitative PCR but it is becoming evident that microRNA isoforms (isomiRs) are potentially physiologically important. It is well known that left ventricular (patho)physiology is influenced by transmural heterogeneity of cardiomyocyte phenotype, and this likely reflects underlying heterogeneity of gene expression. Given the significant role of miRNAs in regulating gene expression, knowledge of how the miRNA profile varies across the ventricular wall will be crucial to better understand the mechanisms governing transmural physiological heterogeneity. To determinine miRNA/isomiR expression profiles in the rat heart we investigated tissue from different locations across the left ventricular wall using deep sequencing. We detected significant quantities of 145 known rat miRNAs and 68 potential novel orthologs of known miRNAs, in mature, mature* and isomiR formation. Many isomiRs were detected at a higher frequency than their canonical sequence in miRBase and have different predicted targets. The most common miR-133a isomiR was more effective at targeting a construct containing a sequence from the gelsolin gene than was canonical miR-133a, as determined by dual-fluorescence assay. We identified a novel rat miR-1 homolog from a second miR-1 gene; and a novel rat miRNA similar to miR-676. We also cloned and sequenced the rat miR-486 gene which is not in miRBase (v18). Signalling pathways predicted to be targeted by the most highly detected miRNAs include Ubiquitin-mediated Proteolysis, Mitogen-Activated Protein Kinase, Regulation of Actin Cytoskeleton, Wnt signalling, Calcium Signalling, Gap junctions and Arrhythmogenic Right Ventricular Cardiomyopathy. Most miRNAs are not expressed in a gradient across the ventricular wall, with exceptions including miR-10b, miR-21, miR-99b and miR-486.

ShortStack: comprehensive annotation and quantification of small RNA genes

Small RNA sequencing allows genome-wide discovery, categorization, and quantification of genes producing regulatory small RNAs. Many tools have been described for annotation and quantification of microRNA loci (MIRNAs) from small RNA-seq data. However, in many organisms and tissue types, MIRNA genes comprise only a small fraction of all small RNA-producing genes. ShortStack is a stand-alone application that analyzes reference-aligned small RNA-seq data and performs comprehensive de novo annotation and quantification of the inferred small RNA genes. ShortStack's output reports multiple parameters of direct relevance to small RNA gene annotation, including RNA size distributions, repetitiveness, strandedness, hairpin-association, MIRNA annotation, and phasing. In this study, ShortStack is demonstrated to perform accurate annotations and useful descriptions of diverse small RNA genes from four plants (Arabidopsis, tomato, rice, and maize) and three animals (Drosophila, mice, and humans). ShortStack efficiently processes very large small RNA-seq data sets using modest computational resources, and its performance compares favorably to previously described tools. Annotation of MIRNA loci by ShortStack is highly specific in both plants and animals. ShortStack is freely available under a GNU General Public License.

RNA interference targets arbovirus replication in Culicoides cells

Arboviruses are transmitted to vertebrate hosts by biting arthropod vectors such as mosquitoes, ticks, and midges. These viruses replicate in both arthropods and vertebrates and are thus exposed to different antiviral responses in these organisms. RNA interference (RNAi) is a sequence-specific RNA degradation mechanism that has been shown to play a major role in the antiviral response against arboviruses in mosquitoes. Culicoides midges are important vectors of arboviruses, known to transmit pathogens of humans and livestock such as bluetongue virus (BTV) (Reoviridae), Oropouche virus (Bunyaviridae), and likely the recently discovered Schmallenberg virus (Bunyaviridae). In this study, we investigated whether Culicoides cells possess an antiviral RNAi response and whether this is effective against arboviruses, including those with double-stranded RNA (dsRNA) genomes, such as BTV. Using reporter gene-based assays, we established the presence of a functional RNAi response in Culicoides sonorensis-derived KC cells which is effective in inhibiting BTV infection. Sequencing of small RNAs from KC and Aedes aegypti-derived Aag2 cells infected with BTV or the unrelated Schmallenberg virus resulted in the production of virus-derived small interfering RNAs (viRNAs) of 21 nucleotides, similar to the viRNAs produced during arbovirus infections of mosquitoes. In addition, viRNA profiles strongly suggest that the BTV dsRNA genome is accessible to a Dicer-type nuclease. Thus, we show for the first time that midge cells target arbovirus replication by mounting an antiviral RNAi response mainly resembling that of other insect vectors of arboviruses.

Asymmetric purine-pyrimidine distribution in cellular small RNA population of papaya

Background

The small RNAs (sRNA) are a regulatory class of RNA mainly represented by the 21 and 24-nucleotide size classes. The cellular sRNAs are processed by RNase III family enzyme dicer (Dicer like in plant) from a self-complementary hairpin loop or other type of RNA duplexes. The papaya genome has been sequenced, but its microRNAs and other regulatory RNAs are yet to be analyzed.

Results

We analyzed the genomic features of the papaya sRNA population from three sRNA deep sequencing libraries made from leaves, flowers, and leaves infected with Papaya Ringspot Virus (PRSV). We also used the deep sequencing data to annotate the micro RNA (miRNA) in papaya. We identified 60 miRNAs, 24 of which were conserved in other species, and 36 of which were novel miRNAs specific to papaya. In contrast to the Chargaff’s purine-pyrimidine equilibrium, cellular sRNA was significantly biased towards a purine rich population. Of the two purine bases, higher frequency of adenine was present in 23nt or longer sRNAs, while 22nt or shorter sRNAs were over represented by guanine bases. However, this bias was not observed in the annotated miRNAs in plants. The 21nt species were expressed from fewer loci but expressed at higher levels relative to the 24nt species. The highly expressed 21nt species were clustered in a few isolated locations of the genome. The PRSV infected leaves showed higher accumulation of 21 and 22nt sRNA compared to uninfected leaves. We observed higher accumulation of miRNA* of seven annotated miRNAs in virus-infected tissue, indicating the potential function of miRNA* under stressed conditions.

Conclusions

We have identified 60 miRNAs in papaya. Our study revealed the asymmetric purine-pyrimidine distribution in cellular sRNA population. The 21nt species of sRNAs have higher expression levels than 24nt sRNA. The miRNA* of some miRNAs shows higher accumulation in PRSV infected tissues, suggesting that these strands are not totally functionally redundant. The findings open a new avenue for further investigation of the sRNA silencing pathway in plants.

Emerging Roles for Non-Coding RNAs in Male Reproductive Development in Flowering Plants

Knowledge of sexual reproduction systems in flowering plants is essential to humankind, with crop fertility vitally important for food security. Here, we review rapidly emerging new evidence for the key importance of non-coding RNAs in male reproductive development in flowering plants. From the commitment of somatic cells to initiating reproductive development through to meiosis and the development of pollen—containing the male gametes (sperm cells)—in the anther, there is now overwhelming data for a diversity of non-coding RNAs and emerging evidence for crucial roles for them in regulating cellular events at these developmental stages. A particularly exciting development has been the association of one example of cytoplasmic male sterility, which has become an unparalleled breeding tool for producing new crop hybrids, with a non-coding RNA locus.

Diverse correlation patterns between microRNAs and their targets during tomato fruit development indicates different modes of microRNA actions

MicroRNAs negatively regulate the accumulation of mRNAs therefore when they are expressed in the same cells their expression profiles show an inverse correlation. We previously described one positively correlated miRNA/target pair, but it is not known how widespread this phenomenon is. Here, we investigated the correlation between the expression profiles of differentially expressed miRNAs and their targets during tomato fruit development using deep sequencing, Northern blot and RT-qPCR. We found an equal number of positively and negatively correlated miRNA/target pairs indicating that positive correlation is more frequent than previously thought. We also found that the correlation between microRNA and target expression profiles can vary between mRNAs belonging to the same gene family and even for the same target mRNA at different developmental stages. Since microRNAs always negatively regulate their targets, the high number of positively correlated microRNA/target pairs suggests that mutual exclusion could be as widespread as temporal regulation. The change of correlation during development suggests that the type of regulatory circuit directed by a microRNA can change over time and can be different for individual gene family members. Our results also highlight potential problems for expression profiling-based microRNA target identification/validation.

Everybody wins! Poland hosts thrilling competitions of viruses, RNAi and football teams

The ESF-EMBO conference on 'Antiviral RNAi: From Molecular Biology towards Applications' took place in June 2012 in Pultusk, Poland. It brought together scientists working at the interface of RNAi and virus infections in different organisms, covering the complete range from basic mechanisms of RNA silencing to RNAi-based antiviral therapy.