Identification of microRNA-Size, Small RNAs in Escherichia coli

Noncoding small regulatory RNA molecules control gene expression and microRNAs provide one of the best examples in eukaryotes. However, bacterial RNAs of comparable size to eukaryotic microRNAs have received little attention. Here, we demonstrate the existence of microRNA-size, small RNAs (msRNAs) in the model bacterium Escherichia coli. We examined the small RNAs in E. coli using a deep sequencing approach, and analyzed 33.2 million small RNA clone reads after size fractionation. Bioinformatic analysis of the whole set revealed more than 400 individual msRNA species. The cellular contents of selected highly expressed msRNAs were verified by quantitative RT-PCR and northern blotting. Although, the functional significance of these RNAs is unclear, their high abundance suggests that they may play specialized roles in bacteria, analogous to miRNAs in eukaryotes.     

Detection of 5'- and 3'-UTR-derived small RNAs and cis-encoded antisense RNAs in Escherichia coli

Evidence is accumulating that small, noncoding RNAs are important regulatory molecules. Computational and experimental searches have led to the identification of ~60 small RNA genes in Escherichia coli. However, most of these studies focused on the intergenic regions and assumed that small RNAs were >50 nt. Thus, the previous screens missed small RNAs encoded on the antisense strand of protein-coding genes and small RNAs of <50 nt. To identify additional small RNAs, we carried out a cloning-based screen focused on RNAs of 30–65 nt. In this screen, we identified RNA species corresponding to fragments of rRNAs, tRNAs and known small RNAs. Several of the small RNAs also corresponded to 5′- and 3′-untranslated regions (UTRs) and internal fragments of mRNAs. Four of the 3′-UTR-derived RNAs were highly abundant and two showed expression patterns that differed from the corresponding mRNAs, suggesting independent functions for the 3′-UTR-derived small RNAs. We also detected three previously unidentified RNAs encoded in intergenic regions and RNAs from the long direct repeat and hok/sok elements. In addition, we identified a few small RNAs that are expressed opposite protein-coding genes and could base pair with 5′ or 3′ ends of the mRNAs with perfect complementarity.     

Discovery of Replicating Circular RNAs by RNA-Seq and Computational Algorithms

Replicating circular RNAs are independent plant pathogens known as viroids, or act to modulate the pathogenesis of plant and animal viruses as their satellite RNAs. The rate of discovery of these subviral pathogens was low over the past 40 years because the classical approaches are technical demanding and time-consuming. We previously described an approach for homology-independent discovery of replicating circular RNAs by analysing the total small RNA populations from samples of diseased tissues with a computational program known as progressive filtering of overlapping small RNAs (PFOR). However, PFOR written in PERL language is extremely slow and is unable to discover those subviral pathogens that do not trigger in vivo accumulation of extensively overlapping small RNAs. Moreover, PFOR is yet to identify a new viroid capable of initiating independent infection. Here we report the development of PFOR2 that adopted parallel programming in the C++ language and was 3 to 8 times faster than PFOR. A new computational program was further developed and incorporated into PFOR2 to allow the identification of circular RNAs by deep sequencing of long RNAs instead of small RNAs. PFOR2 analysis of the small RNA libraries from grapevine and apple plants led to the discovery of Grapevine latent viroid (GLVd) and Apple hammerhead viroid-like RNA (AHVd-like RNA), respectively. GLVd was proposed as a new species in the genus Apscaviroid, because it contained the typical structural elements found in this group of viroids and initiated independent infection in grapevine seedlings. AHVd-like RNA encoded a biologically active hammerhead ribozyme in both polarities, and was not specifically associated with any of the viruses found in apple plants. We propose that these computational algorithms have the potential to discover novel circular RNAs in plants, invertebrates and vertebrates regardless of whether they replicate and/or induce the in vivo accumulation of small RNAs.     

Small RNAs of Rous Sarcoma Virus: Characterization by Two-Dimensional Polyacrylamide Gel Electrophoresis and Fingerprint Analysis

Approximately 15 to 20 different species of small (4 to 7S) RNAs have been purified by two-dimensional polyacrylamide gel electrophoresis of RNA isolated from virions of Schmidt-Ruppin D strain of Rous sarcoma virus. Each species of small RNA has been isolated free of 70S RNA; nine of them, including 5S and 7S RNAs, were also found associated with the 70S genomic RNA. Most of the 4S RNAs are present at an average of less than one copy per virion. The 4S RNAs have T1 RNase (EC 2.7.7.26) fingerprints, which are very similar to those of tRNAs. One of the smallest 4S RNAs, which can act as a primer for initiation of RNA-directed DNA synthesis, is associated with the 70S RNA in 1 to 2 copies per complex, whereas an additional 6 to 8 copies of this molecule are free.     

RNA-Seq analyses reveal the order of tRNA processing events and the maturation of C/D box and CRISPR RNAs in the hyperthermophile Methanopyrus kandleri

The methanogenic archaeon Methanopyrus kandleri grows near the upper temperature limit for life. Genome analyses revealed strategies to adapt to these harsh conditions and elucidated a unique transfer RNA (tRNA) C-to-U editing mechanism at base 8 for 30 different tRNA species. Here, RNA-Seq deep sequencing methodology was combined with computational analyses to characterize the small RNome of this hyperthermophilic organism and to obtain insights into the RNA metabolism at extreme temperatures. A large number of 132 small RNAs were identified that guide RNA modifications, which are expected to stabilize structured RNA molecules. The C/D box guide RNAs were shown to exist as circular RNA molecules. In addition, clustered regularly interspaced short palindromic repeats RNA processing and potential regulatory RNAs were identified. Finally, the identification of tRNA precursors before and after the unique C8-to-U8 editing activity enabled the determination of the order of tRNA processing events with termini truncation preceding intron removal. This order of tRNA maturation follows the compartmentalized tRNA processing order found in Eukaryotes and suggests its conservation during evolution.     

The Complex Exogenous RNA Spectra in Human Plasma: An Interface with Human Gut Biota?

Human plasma has long been a rich source for biomarker discovery. It has recently become clear that plasma RNA molecules, such as microRNA, in addition to proteins are common and can serve as biomarkers. Surveying human plasma for microRNA biomarkers using next generation sequencing technology, we observed that a significant fraction of the circulating RNA appear to originate from exogenous species. With careful analysis of sequence error statistics and other controls, we demonstrated that there is a wide range of RNA from many different organisms, including bacteria and fungi as well as from other species. These RNAs may be associated with protein, lipid or other molecules protecting them from RNase activity in plasma. Some of these RNAs are detected in intracellular complexes and may be able to influence cellular activities under in vitro conditions. These findings raise the possibility that plasma RNAs of exogenous origin may serve as signaling molecules mediating for example the human-microbiome interaction and may affect and/or indicate the state of human health.     

Biosynthesis of transfer RNA: isolation and characterization of precursors to transfer RNA in the posterior silkgland of Bombyx mori.

Distinct low molecular weight RNA species that have properties expected for the precursor to tRNA have been isolated from the posterior silkglands of the silkworm Bombyx mori. These RNAs migrate between 4 S and 5 S markers on acrylamide gels and are labeled preferentially in vivo in relation to tRNA. The precursor RNAs can be converted specifically into molecules indistinguishable in size from tRNA upon incubation with “cleavage” enzymes isolated from the silkgland ribosomes. Two of the three low molecular weight RNAs contain the modified residues, pseudouridine, dihydrouridine and ribothymidine, and are methylated in vivo, suggesting that these base modifications occur while the tRNA is still in its precursor stage.     

Circular RNAs are abundant,conserved, and associated with ALU repeats

Circular RNAs composed of exonic sequence have been described in a small number of genes. Thought to result from splicing errors, circular RNA species possess no known function. To delineate the universe of endogenous circular RNAs, we performed high-throughput sequencing (RNA-seq) of libraries prepared from ribosome-depleted RNA with or without digestion with the RNA exonuclease, RNase R. We identified >25,000 distinct RNA species in human fibroblasts that contained non-colinear exons (a “backsplice”) and were reproducibly enriched by exonuclease degradation of linear RNA. These RNAs were validated as circular RNA (ecircRNA), rather than linear RNA, and were more stable than associated linear mRNAs in vivo. In some cases, the abundance of circular molecules exceeded that of associated linear mRNA by >10-fold. By conservative estimate, we identified ecircRNAs from 14.4% of actively transcribed genes in human fibroblasts. Application of this method to murine testis RNA identified 69 ecircRNAs in precisely orthologous locations to human circular RNAs. Of note, paralogous kinases HIPK2 and HIPK3 produce abundant ecircRNA from their second exon in both humans and mice. Though HIPK3 circular RNAs contain an AUG translation start, it and other ecircRNAs were not bound to ribosomes. Circular RNAs could be degraded by siRNAs and, therefore, may act as competing endogenous RNAs. Bioinformatic analysis revealed shared features of circularized exons, including long bordering introns that contained complementary ALU repeats. These data show that ecircRNAs are abundant, stable, conserved and nonrandom products of RNA splicing that could be involved in control of gene expression.     

Cloning and Expression Studies of Novel Small RNAs in Tetraploid Cotton

Small RNAs are a group of non-coding RNAs that downregulate gene expression in a sequence-specific manner to control plant growth and development. The objective of the present study was to clone and characterize several small RNAs in cotton. To identify small RNAs that are involved in the development of cotton bolls and fibers, we generated cDNA libraries from cotton bolls at 13?days post-anthesis from two cotton cultivars, Pima Phy 76 (Gossypium bardadense) and Acala 1517?C99 (Gossypium hirsutum). Screening of these libraries identified eight small RNAs, seven of which have not been reported in other plant species and appear to be absent in the known sequences of other plant species. Their predicted target genes are known to be involved in cotton fiber development. The cloned small RNAs displayed lower and differential expression in the examined boll developmental stages using RT-PCR and quantitative RT-PCR. The genetic polymorphism of the small RNAs at the DNA level was evaluated by miRNA-amplified fragment length polymorphism (AFLP) analysis using primers designed from the small miRNA genes in combination with AFLP primers. Homologous small RNA gene sequences were further isolated using this homology-based genotyping approach, and potential hairpin structures were identified. The results represent a novel method to isolate small including miRNA genes at the RNA and DNA levels in many plant species where genome sequences are not available or expressed sequence tags are limited.     

The small RNA content of human sperm reveals pseudogene-derived piRNAs complementary to protein-coding genes

At the end of mammalian sperm development, sperm cells expel most of their cytoplasm and dispose of the majority of their RNA. Yet, hundreds of RNA molecules remain in mature sperm. The biological significance of the vast majority of these molecules is unclear. To better understand the processes that generate sperm small RNAs and what roles they may have, we sequenced and characterized the small RNA content of sperm samples from two human fertile individuals. We detected 182 microRNAs, some of which are highly abundant. The most abundant microRNA in sperm is miR-1246 with predicted targets among sperm-specific genes. The most abundant class of small noncoding RNAs in sperm are PIWI-interacting RNAs (piRNAs). Surprisingly, we found that human sperm cells contain piRNAs processed from pseudogenes. Clusters of piRNAs from human testes contain pseudogenes transcribed in the antisense strand and processed into small RNAs. Several human protein-coding genes contain antisense predicted targets of pseudogene-derived piRNAs in the male germline and these piRNAs are still found in mature sperm. Our study provides the most extensive data set and annotation of human sperm small RNAs to date and is a resource for further functional studies on the roles of sperm small RNAs. In addition, we propose that some of the pseudogene-derived human piRNAs may regulate expression of their parent gene in the male germline.     

Evidence for a role of RNA in eukaryotic chromosome structure. Metabolically stable, small nuclear RNA species are covalently linked to chromosomal DNA in HeLa cells

An investigation of metabolically stable, chromatin-associated RNA in HeLa cells has revealed that three small RNA species, 193, 171 and 127 nucleotides in length, are covalently linked to double-stranded chromosomal DNA through phosphodiester bonds. These DNA-linked RNAs appear to be members of the small nuclear RNA species that have been identified in a wide variety of eukaryotic cells, and they are tentatively identified as species C, D and G′, in the nomenclature system currently employed for HeLa cell small nuclear RNAs. These DNA-linked RNAs do not appear to be involved in priming DNA replication, since they are of relatively high metabolic stability ($\text{t}{}_{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{= 19}$ hours in HeLa cells with a 21·5-hour cell generation time) and since their covalently contiguous DNA stretches are not enriched in newly replicated material. They lack saturated pyrimidine bases (level of detection = 0·15 mol %) and are therefore not “chromosomal RNA”, as defined by its proponents. The covalent linkage of these small RNA species with chromosomal DNA was discovered by virtue of the fact that when highly purified HeLa cell chromatin is dissociated by chaotropic solutes, these RNAs are released in association with small pieces of double-stranded DNA (approx. 475 nucleotide pairs). These DNA-RNA complexes can then be purified by removing the bulk, high molecular weight DNA by ultra-centrifugation. The resulting DNA-RNA complexes are shown to be covalently joined by several criteria, including equilibrium density-gradient centrifugation in either Cs₂SO₄/dimethylsulfoxide or aqueous Cs₂SO₄/formaldehyde after thermal denaturation (90 °C in 50% formamide, which is 55 deg. C above the melting temperature of this DNA), by the chromat ographicbehavior of the complexes on hydroxylapatite before and after thermal denaturation, and by the demonstration of alkali-resistant ribonucleotides flanking the 3′ hydroxyl termini of the DNA, the latter criterion providing evidence for 3′ to 5′ DNA-RNA phosphodiester bonds. Reconstruction experiments involving addition of the purified RNAs to nuclei or chromatin demonstrate that the covalent DNA-RNA linkages do not arise by ligation events during cell fractionation. Further experiments indicate the existence of a dynamic equilibrium of these small nuclear RNA species between chromosomal and nucleoplasmic loci in vivo, and other considerations suggest that this equilibrium may be cell cycle-dependent. The DNA adjacent to these covalently linked RNAs has the same melting temperature as total HeLa chromosomal DNA and its reassociation kinetics reveal the presence of both repeated and non-repeated sequences, implying that the DNA-linked RNAs are widely distributed throughout the HeLa cell genome. It is proposed that these DNA-linked RNAs are involved in the tertiary structure of chromatin, particularly in relation to the cell cycle.     

siRNAs from miRNA sites mediate DNA methylation of target genes

Arabidopsis microRNA (miRNA) genes (MIR) give rise to 20- to 22-nt miRNAs that are generated predominantly by the type III endoribonuclease Dicer-like 1 (DCL1) but do not require any RNA-dependent RNA Polymerases (RDRs) or RNA Polymerase IV (Pol IV). Here, we identify a novel class of non-conserved MIR genes that give rise to two small RNA species, a 20- to 22-nt species and a 23- to 27-nt species, at the same site. Genetic analysis using small RNA pathway mutants reveals that the 20- to 22-nt small RNAs are typical miRNAs generated by DCL1 and are associated with Argonaute 1 (AGO1). In contrast, the accumulation of the 23- to 27-nt small RNAs from the miRNA-generating sites is dependent on DCL3, RDR2 and Pol IV, components of the typical heterochromatic small interfering RNA (hc-siRNA) pathway. We further demonstrate that these MIR-derived siRNAs associate with AGO4 and direct DNA methylation at some of their target loci in trans. In addition, we find that at the miRNA-generating sites, some conserved canonical MIR genes also produce siRNAs, which also induce DNA methylation at some of their target sites. Our systematic examination of published small RNA deep sequencing datasets of rice and moss suggests that this type of dual functional MIRs exist broadly in plants.     

Characterisation of the Small RNAs in the Biomedically Important Green-Bottle Blowfly Lucilia sericata

Background

The green bottle fly maggot, Lucilia sericata, is a species with importance in medicine, agriculture and forensics. Improved understanding of this species’ biology is of great potential benefit to many research communities. MicroRNAs (miRNA) are a short non-protein coding regulatory RNA, which directly regulate a host of protein coding genes at the translational level. They have been shown to have developmental and tissue specific distributions where they impact directly on gene regulation. In order to improve understanding of the biology of L. sericata maggots we have performed small RNA-sequencing of their secretions and tissue at different developmental stages.

Results

We have successfully isolated RNA from the secretions of L. sericata maggots. Illumina small RNA-sequencing of these secretions and the three tissues (crop, salivary gland, gut) revealed that the most common small RNA fragments were derived from ribosomal RNA and transfer RNAs of both insect and bacterial origins. These RNA fragments were highly specific, with the most common tRNAs, such as GlyGCC, predominantly represented by reads derived from the 5’ end of the mature maggot tRNA. Each library also had a unique profile of miRNAs with a high abundance of miR-10-5p in the maggot secretions and gut and miR-8 in the food storage organ the crop and salivary glands. The pattern of small RNAs in the bioactive maggot secretions suggests they originate from a combination of saliva, foregut and hindgut tissues. Droplet digital RT-PCR validation of the RNA-sequencing data shows that not only are there differences in the tissue profiles for miRNAs and small RNA fragments but that these are also modulated through developmental stages of the insect.

Conclusions

We have identified the small-RNAome of the medicinal maggots L. sericata and shown that there are distinct subsets of miRNAs expressed in specific tissues that also alter during the development of the insect. Furthermore there are very specific RNA fragments derived from other non-coding RNAs present in tissues and in the secretions. This new knowledge has applicability in diverse research fields including wound healing, agriculture and forensics.