Characterization of 29 newly isolated bacteriophages as a potential therapeutic agent against IMP-6-producing Klebsiella pneumoniae from clinical specimens.

Carbapenemase-producing Enterobacteriaceae (CPE) are one of the most detrimental species of antibiotic-resistant bacteria globally. Phage therapy has emerged as an effective strategy for the treatment of CPE infections. In western Japan, the rise of Klebsiella pneumoniae strains harboring the pKPI-6 plasmid encoding bla IMP-6 is of increasing concern. To address this challenge, we isolated 29 phages from Japanese sewage, specifically targeting 31 K. pneumoniae strains and one Escherichia coli strain harboring the pKPI-6 plasmid. Electron microscopy analysis revealed that among the 29 isolated phages, 21 (72.4%), 5 (17.2%), and 3 (10.3%) phages belonged to myovirus, siphovirus, and podovirus morphotypes, respectively. Host range analysis showed that 18 Slopekvirus strains within the isolated phages infected 25-26 K. pneumoniae strains, indicating that most of the isolated phages have a broad host range. Notably, K. pneumoniae strain Kp21 was exclusively susceptible to phage øKp_21, whereas Kp22 exhibited susceptibility to over 20 phages. Upon administering a phage cocktail composed of 10 phages, we observed delayed emergence of phage-resistant bacteria in Kp21 but not in Kp22. Intriguingly, phage-resistant Kp21 exhibited heightened sensitivity to other bacteriophages, indicating a "trade-off" for resistance to phage øKp_21. Our proposed phage set has an adequate number of phages to combat the K. pneumoniae strain prevalent in Japan, underscoring the potential of a well-designed phage cocktail in mitigating the occurrence of phage-resistant bacteria. IMPORTANCE The emergence of Klebsiella pneumoniae harboring the bla IMP-6 plasmid poses an escalating threat in Japan. In this study, we found 29 newly isolated bacteriophages that infect K. pneumoniae strains carrying the pKPI-6 plasmid from clinical settings in western Japan. Our phages exhibited a broad host range. We applied a phage cocktail treatment composed of 10 phages against two host strains, Kp21 and Kp22, which displayed varying phage susceptibility patterns. Although the phage cocktail delayed the emergence of phage-resistant Kp21, it was unable to hinder the emergence of phage-resistant Kp22. Moreover, the phage-resistant Kp21 became sensitive to other phages that were originally non-infective to the wild-type Kp21 strains. Our study highlights the potential of a well-tailored phage cocktail in reducing the occurrence of phage-resistant bacteria.

Comparison of mucosa-associated microbiota in Crohn's disease patients with and without anti-tumor necrosis factor-α therapy.

Most studies on the gut microbiome of Crohn's disease have been conducted using feces, instead of intestinal mucus to analyze the mucosa-associated microbiota. To investigate the characteristics of mucosa-associated microbiota in Crohn's disease patients and the effect of anti-tumor necrosis factor (TNF)-α therapy on mucosa-associated microbiota, we analyzed microbiota in Crohn's disease patients using brushing samples taken from terminal ileum. The recruited subjects were 18 Crohn's disease patients and 13 controls. There were 10 patients with anti-TNF-α therapy in Crohn's disease group. Crohn's disease patients had significantly reduced α-diversity in Shannon index compared to the controls. The comparative analysis of the taxonomic composition at the genus level between the Crohn's disease group and the controls indicated that butyrate-producing bacteria were less abundant in the Crohn's disease group compared to the controls. There were no differences in the diversity between the patients taking anti-TNF-α therapy and the patients without. The comparative analysis of the taxonomic composition at the genus level between the two groups indicated that some of anti-inflammatory bacteria were less abundant in the anti-TNF-α therapy group than the other. Reduction of specific bacteria producing anti-inflammatory molecules, especially butyrate-producing bacteria may play important roles in the pathophysiology of Crohn's disease.

Distribution of Antimicrobial Resistance and Virulence Genes within the Prophage-Associated Regions in Nosocomial Pathogens.

Prophages are often involved in host survival strategies and contribute toward increasing the genetic diversity of the host genome. Prophages also drive horizontal propagation of various genes as vehicles. However, there are few retrospective studies contributing to the propagation of antimicrobial resistance (AMR) and virulence factor (VF) genes by prophage. We extracted the complete genome sequences of seven pathogens, including ESKAPE bacteria and Escherichia coli from a public database, and examined the distribution of both the AMR and VF genes in prophage-like regions. We found that the ratios of AMR and VF genes greatly varied among the seven species. More than 70% of Enterobacter cloacae strains had VF genes, but only 1.2% of Klebsiella pneumoniae strains had VF genes from prophages. AMR and VF genes are unlikely to exist together in the same prophage region except in E. coli and Staphylococcus aureus, and the distribution patterns of prophage types containing AMR genes are distinct from those of VF gene-carrying prophage types. AMR genes in the prophage were located near transposase and/or integrase. The prophage containing class 1 integrase possessed a significantly greater number of AMR genes than did prophages with no class 1 integrase. The results of this study present a comprehensive picture of AMR and VF genes present within, or close to, prophage-like elements and different prophage patterns between AMR- or VF-encoding prophage-like elements. IMPORTANCE Although we believe phages play an important role in horizontal gene transfer in exchanging genetic material, we do not know the distribution of the antimicrobial resistance (AMR) and/or virulence factor (VF) genes in prophages. We collected different prophage elements from the complete genome sequences of seven species-Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, Enterobacter cloacae, and Escherichia coli-and characterized the distribution of antimicrobial resistance and virulence genes located in the prophage region. While virulence genes in prophage were species specific, antimicrobial resistance genes in prophages were highly conserved in various species. An integron structure was detected within specific prophage regions such as P1-like prophage element. Maximum of 10 antimicrobial resistance genes were found in a single prophage region, suggesting that prophages act as a reservoir for antimicrobial resistance genes. The results of this study show the different characteristic structures between AMR- or VF-encoding prophages.

Development of multiplexing gene silencing system using conditionally induced polycistronic synthetic antisense RNAs in Escherichia coli.

Although efficient methods of gene silencing have been established in eukaryotes, many different techniques are still used in bacteria due to the lack of a standardized tool. Here, we developed a convenient and efficient method to downregulate the expression of a specific gene using ∼140 nucleotide RNA with a 24-nucleotide antisense region from an arabinose-inducible expression plasmid by taking Escherichia coli lacZ and phoA genes encoding ß-galactosidase and alkaline phosphatase, respectively, as target genes to evaluate the model. We examined the antisense RNA (asRNA) design, including targeting position, uORF stability elements at the 5'-end, and Hfq-binding module at the 3'-end, and inducer amount required to obtain effective experimental conditions for gene silencing. Furthermore, we constructed multiplexed dual-acting asRNA genes in the plasmid, which were transcribed as polycistronic RNA and were able to knockdown multiple target genes simultaneously. We observed the highest inhibition level of 98.6% when lacZ was targeted using the pMKN104 asRNA expression plasmid, containing a five times stronger PBAD -10 promoter sequence with no requirement of the Hfq protein for repression. These features allow the system to be utilized as an asRNA expression platform in many bacteria, besides E. coli, for gene regulation.

Comprehensive analysis of liver and blood miRNA in precancerous conditions.

Streptozotocin administration to mice (STZ-mice) induces type I diabetes and hepatocellular carcinoma (HCC). We attempted to elucidate the carcinogenic mechanism and the miRNA expression status in the liver and blood during the precancerous state. Serum and liver tissues were collected from STZ-mice and non-treated mice (CTL-mice) at 6, 10, and 12 W. The exosome enriched fraction extracted from serum was used. Hepatic histological examination and hepatic and exosomal miRNA expression analysis were serially performed using next-generation sequencing (NGS). Human miRNA expression analysis of chronic hepatitis liver tissue and exosomes, which were collected before starting the antiviral treatment, were also performed. No inflammation or fibrosis was found in the liver of CTL-mice during the observation period. In STZ-mice, regeneration and inflammation of hepatocytes was found at 6 W and nodules of atypical hepatocytes were found at 10 and 12 W. In the liver tissue, during 6-12 W, the expression levels of let-7f-5p, miR-143-3p, 148a-3p, 191-5p, 192-5p, 21a-5p, 22-3p, 26a-5p, and 92a-3p was significantly increased in STZ-mice, and anti-oncogenes of their target gene candidates were down-regulated. miR-122-5p was also significantly down-regulated in STZ-mice. Fifteen exosomal miRNAs were upregulated in STZ-mice. Six miRNAs (let-7f-5p, miR-10b-5p, 143-3p, 191-5p, 21a-5p, and 26a-5p) were upregulated, similarly to human HCC cases. From the precancerous state, aberrant expression of hepatic miRNAs has already occurred, and then, it can promote carcinogenesis. In exosomes, the expression pattern of common miRNAs between mice and humans before carcinogenesis was observed and can be expected to be developed as a cancer predictive marker.

Antimicrobial antisense RNA delivery to F-pili producing multidrug-resistant bacteria via a genetically engineered bacteriophage.

Multidrug-resistant bacteria are a growing issue worldwide. This study developed a convenient and effective method to downregulate the expression of a specific gene to produce a novel antimicrobial tool using a small (140 nucleotide) RNA with a 24-nucleotide antisense (as) region from an arabinose-inducible expression phagemid vector in Escherichia coli. Knockdown effects of rpoS encoding RNA polymerase sigma factor were observed using this inducible artificial asRNA approach. asRNAs targeting several essential E. coli genes produced significant growth defects, especially when targeted to acpP and ribosomal protein coding genes rplN, rplL, and rpsM. Growth inhibited phenotypes were facilitated in hfq- conditions. Phage lysates were prepared from cells harboring phagemids as a lethal-agent delivery tool. Targeting the rpsM gene by phagemid-derived M13 phage infection of E. coli containing a carbapenem-producing F-plasmid and multidrug-resistant Klebsiella pneumoniae containing an F-plasmid resulted in the death of over 99.99% of infected bacteria. This study provides a possible strategy for treating bacterial infection and can be applied to any F-pilus producing bacterial species.

Artificial small RNA-mediated growth inhibition in Escherichia coli and Salmonella enterica serovar Typhimurium.

We developed a synthetic RNA approach to identify growth inhibition sequences by cloning random 24-nucleotide (nt) sequences into an arabinose-inducible expression vector. This vector expressed a small RNA (sRNA) of ∼140 nt containing a 24 nt random sequence insert. After transforming Escherichia coli with the vector, 10 out of 954 transformants showed strong growth defect phenotypes and two clones caused cell lysis. We then examined growth inhibition phenotypes in the Salmonella Typhimurium LT2 strain using the twelve sRNAs that exerted an inhibitory effect on E. coli growth. Three of these clones showed strong growth inhibition phenotypes in S. Typhimurium LT2. The most effective sRNA contained the same insert (N1) in both bacteria. The 24 nt random sequence insert of N1 was abundant in guanine residues (ten out of 24 nt), and other random sequences causing growth defects were also highly enriched for guanine (G) nucleotides. We, therefore, generated clones that express sRNAs containing a stretch of 16 to 24 continuous guanine sequences (poly-G16, -G18, -G20, -G22, and -G24). All of these clones induced growth inhibition in both liquid and agar plate media and the poly-G20 clone showed the strongest effect in E. coli. These results demonstrate that our sRNA expression system can be used to identify nucleotide sequences that are potential candidates for oligonucleotide antimicrobial drugs.

Functional characterization of the type I toxin Lpt from Lactobacillus rhamnosus by fluorescence and atomic force microscopy.

Lpt is a 29 amino acid long type I toxin identified in the plasmid DNA of wild Lactobacillus rhamnosus strains isolated from food. We previously reported that transcription of the encoding gene was upregulated under nutritional starvation conditions mimicking cheese ripening environment. The heterologous expression of the Lpt peptide in E. coli resulted in cell growth inhibition, nucleoid condensation and compromised integrity of the cell membrane. Fusion of the Lpt peptide with the fluorescent protein mCherry allowed to visualize the accumulation of the peptide into the membrane, while mutagenesis experiments showed that either the insertion of a negatively charged amino acid into the hydrophobic α-helix or deletion of the hydrophilic C-terminal region, leads to a non-toxic peptide. AFM imaging of Lpt expressing E. coli cells has revealed the presence of surface defects that are compatible with the loss of portions of the outer membrane bilayer. This observation provides support for the so-called "carpet" model, by which the Lpt peptide is supposed to destabilize the phospholipid packing through a detergent-like mechanism leading to the removal of small patches of bilayer through micellization.

Circulating osteocyte-derived exosomes contain miRNAs which are enriched in exosomes from MLO-Y4 cells.

Signaling molecules produced by osteocytes have been proposed to serve as soluble factors that contribute to bone remodeling, as well as to homeostasis of other organs. However, to the best of our knowledge, there are currently no studies investigating the role of osteocyte-secreted exosomes. In the present study, ablation of osteocytes in mice [osteocyte-less (OL)] was used to examine the microRNA (miRNA) levels of plasma-circulating exosomes. In order to investigate the function of osteocyte-secreted exosomes, exosomes derived from MLO-Y4 cells were extracted and their miRNA expression levels were examined using miRNA array analysis and deep sequencing. Comparison of miRNA expression levels between plasma exosomes from OL mouse plasma and MLO-Y4-derived exosomes revealed that decreases in the number of miRNAs from exosomes circulating in the OL mouse plasma may be caused by a decrease in secretion of exosomes from osteocytes. These results suggest that osteocytes secrete exosomes containing characterized miRNAs and then circulate in the blood, and may thus transfer their components, including miRNAs, to recipient cells where they function as signaling molecules in other organs and/or tissues to regulate biological responses.

Involvement of an intracellular vesicular transport process in naked-sgRNA-mediated TRUE gene silencing.

tRNase ZL-utilizing efficacious gene silencing (TRUE gene silencing) is an RNA-mediated gene expression control technology with therapeutic potential. Recently, our group demonstrated that a heptamer, mh1 (Bcl­2), targeting human Bcl-2 mRNA, can be taken up by cells without the use of any transfection reagents and can induce the apoptosis of leukemia cells. However, little is known regarding the mechanism of naked small guide (sg)RNA uptake by cultured cells. Therefore, in the present study the effects of various inhibitors on the induction of apoptosis by naked sgRNA treatment were investigated in order to identify the uptake pathway required for sgRNA function in cultured cells. Addition of the endocytosis inhibitors chlorpromazine, nystatin or methyl­ß­cyclodextrin together with naked effective sgRNA was unable to diminish the apoptosis­inducing effects of naked sgRNA or the reduction in target mRNA, suggesting that functional uptake of sgRNA by cells is clathrin­, caveolae­ and raft­independent. Next, chloroquine, an inhibitor of lysosome acidification, and brefeldin A, an inhibitor that blocks protein transport from the Golgi apparatus to the endoplasmic reticulum were administered. In the presence of these compounds, the apoptosis­inducing effects of naked sgRNA were reduced. These results suggest that a vesicular transport process is involved in sgRNA­mediated TRUE gene silencing. A greater understanding of how naked sgRNAs enter cells and how they reach their target RNAs may aid in the design of more specifically­targeted and potent sgRNA drugs.

Potential small guide RNAs for tRNase ZL from human plasma, peripheral blood mononuclear cells, and cultured cell lines.

Several pieces of evidence suggest that small RNA degradation products together with tRNase ZL appear to form another layer of the whole gene regulatory network. The degraded RNA such as a 5'-half-tRNA and an rRNA fragment function as small guide RNA (sgRNA) to guide the enzyme to target RNA. We were curious whether there exist RNAs in plasma that can function as sgRNAs for tRNase ZL, whether these RNAs are working as signaling molecules between cells to fulfill physiological roles, and whether there are any differences in plasma sgRNA species and levels between normal and pathological conditions. Here, we analyzed small plasma RNAs from three healthy persons and three multiple myeloma patients for potential sgRNAs by deep sequencing. We also examined small RNAs from peripheral blood mononuclear cells (PBMC) of three healthy persons and three myeloma patients and from various cultured human cell lines for sgRNAs. We found that read-number distribution patterns of plasma and PBMC RNAs differ between persons in the range of 5-40 nt and that there are many RNA species that exist significantly more or less abundantly in the plasma or PBMC of the myeloma patients than those of the healthy persons. Furthermore, we found that there are many potential sgRNAs in the 5-40-nt RNAs and that, among them, a 31-nt RNA fragment derived from 94-nt Y4-RNA, which can function as a 5'-half-tRNA-type sgRNA, is overwhelmingly abundant in the plasma of 2/3 of the examinees. These observations suggest that the gene regulatory network via tRNase ZL and sgRNA may be extended intercellularly.

Screening of a heptamer-type sgRNA library for potential therapeutic agents against hematological malignancies.

tRNase-Z(L)-utilizing efficacious (TRUE) gene silencing is an RNA-mediated gene expression control technology that has therapeutic potential. This technology is based on the property of tRNase Z(L) that it can cleave any target RNA at any desired site under the direction of an appropriate artificial small guide RNA (sgRNA). To search for novel potential therapeutic sgRNAs for hematological malignancies, we screened a library composed of 156 sgRNAs, and found that 20 sgRNAs can efficiently induce apoptosis in leukemia and/or myeloma cells. Furthermore, we demonstrated that 4 of the 20 sgRNAs can reduce growth rates of HL60 cells in mouse xenograft models.

Comparison of CAGE and RNA-seq transcriptome profiling using clonally amplified and single-molecule next-generation sequencing.

CAGE (cap analysis gene expression) and RNA-seq are two major technologies used to identify transcript abundances as well as structures. They measure expression by sequencing from either the 5' end of capped molecules (CAGE) or tags randomly distributed along the length of a transcript (RNA-seq). Library protocols for clonally amplified (Illumina, SOLiD, 454 Life Sciences [Roche], Ion Torrent), second-generation sequencing platforms typically employ PCR preamplification prior to clonal amplification, while third-generation, single-molecule sequencers can sequence unamplified libraries. Although these transcriptome profiling platforms have been demonstrated to be individually reproducible, no systematic comparison has been carried out between them. Here we compare CAGE, using both second- and third-generation sequencers, and RNA-seq, using a second-generation sequencer based on a panel of RNA mixtures from two human cell lines to examine power in the discrimination of biological states, detection of differentially expressed genes, linearity of measurements, and quantification reproducibility. We found that the quantified levels of gene expression are largely comparable across platforms and conclude that CAGE and RNA-seq are complementary technologies that can be used to improve incomplete gene models. We also found systematic bias in the second- and third-generation platforms, which is likely due to steps such as linker ligation, cleavage by restriction enzymes, and PCR amplification. This study provides a perspective on the performance of these platforms, which will be a baseline in the design of further experiments to tackle complex transcriptomes uncovered in a wide range of cell types.

Temporal dynamics and transcriptional control using single-cell gene expression analysis.

BACKGROUND: Changes in environmental conditions lead to expression variation that manifest at the level of gene regulatory networks. Despite a strong understanding of the role noise plays in synthetic biological systems, it remains unclear how propagation of expression heterogeneity in an endogenous regulatory network is distributed and utilized by cells transitioning through a key developmental event. RESULTS: Here we investigate the temporal dynamics of a single-cell transcriptional network of 45 transcription factors in THP-1 human myeloid monocytic leukemia cells undergoing differentiation to macrophages. We systematically measure temporal regulation of expression and variation by profiling 120 single cells at eight distinct time points, and infer highly controlled regulatory modules through which signaling operates with stochastic effects. This reveals dynamic and specific rewiring as a cellular strategy for differentiation. The integration of both positive and negative co-expression networks further identifies the proto-oncogene MYB as a network hinge to modulate both the pro- and anti-differentiation pathways. CONCLUSIONS: Compared to averaged cell populations, temporal single-cell expression profiling provides a much more powerful technique to probe for mechanistic insights underlying cellular differentiation. We believe that our approach will form the basis of novel strategies to study the regulation of transcription at a single-cell level.

Divergently overlapping cis-encoded antisense RNA regulating toxin-antitoxin systems from E. coli: hok/sok, ldr/rdl, symE/symR.

Toxin-antitoxin (TA) systems are categorized into three classes based on the type of antitoxin. In type I TA systems, the antitoxin is a small antisense RNA that inhibits translation of small toxic proteins by binding to the corresponding mRNAs. Those type I TA systems were originally identified as plasmid stabilization modules rendering a post-segregational killing (PSK) effect on the host cells. The type I TA loci also exist on the Escherichia coli chromosome but their biological functions are less clear. Genetic organization and regulatory elements of hok/sok and ldr/rdl families are very similar and the toxins are predicted to contain a transmembrane domain, but otherwise share no detectable sequence similarity. This review will give an overview of the type I TA modules of E. coli K-12, especially hok/sok, ldr/rdl and SOS-inducible symE/symR systems, which are regulated by divergently overlapping cis-encoded antisense RNAs.

Novel small noncoding RNAs in mouse spermatozoa, zygotes and early embryos.

The recent discovery of a significant amount of RNA in spermatozoa contradicted the previously held belief that paternal contribution was limited to one copy of the genome. Furthermore, detection of RNA in sperm raised the intriguing question of its possible role in embryonic development. The possibility that RNAs may serve as epigenetic determinants was supported by experiments showing inheritance of epigenetic traits in mice mediated by RNA. We used high-throughput, large-scale sequencing technology to analyze sperm RNA. The RNA sequences generated were diverse in terms of length and included mRNAs, rRNAs, piRNAs, and miRNAs. We studied two small noncoding RNAs enriched in mature sperm, designated sperm RNAs (spR) -12 and -13. They are both encoded in a piRNA locus on chromosome 17, but neither their length (20-21 nt), nor their sequences correspond to known piRNAs or miRNAs. They are resistant to periodate-oxidation-mediated reaction, implying that they undergo terminal post-transcriptional modification. Both were detected in sperm and ovulated unfertilized oocytes, present in one-cell embryos and maintained in preimplantation stages, but not at later differentiation stages. These findings offer a new perspective regarding a possibly important role for gamete-specific small RNAs in early embryogenesis.

pre-miRNA profiles obtained through application of locked nucleic acids and deep sequencing reveals complex 5'/3' arm variation including concomitant cleavage and polyuridylation patterns.

Recent research hints at an underappreciated complexity in pre-miRNA processing and regulation. Global profiling of pre-miRNA and its potential to increase understanding of the pre-miRNA landscape is impeded by overlap with highly expressed classes of other non coding (nc) RNA. Here, we present a data set excluding these RNA before sequencing through locked nucleic acids (LNA), greatly increasing pre-miRNA sequence counts with no discernable effect on pre-miRNA or mature miRNA sequencing. Analysis of profiles generated in total, nuclear and cytoplasmic cell fractions reveals that pre-miRNAs are subject to a wide range of regulatory processes involving loci-specific 3'- and 5'-end variation entailing complex cleavage patterns with co-occurring polyuridylation. Additionally, examination of nuclear-enriched flanking sequences of pre-miRNA, particularly those derived from polycistronic miRNA transcripts, provides insight into miRNA and miRNA-offset (moRNA) production, specifically identifying novel classes of RNA potentially functioning as moRNA precursors. Our findings point to particularly intricate regulation of the let-7 family in many ways reminiscent of DICER1-independent, pre-mir-451-like processing, introduce novel and unify known forms of pre-miRNA regulation and processing, and shed new light on overlooked products of miRNA processing pathways.

miRNA profiling of bilateral rat hippocampal CA3 by deep sequencing.

MicroRNAs (miRNAs) have been demonstrated to be potent post-trascriptional modulators of protein expression. miRNA expression was profiled in the left and right dorsal hippocampal CA3 of mature rats by high-throughput deep sequencing. Among the sequenced and cross-mapped small RNAs, 88% belonged to the miRNAs annotated in the miRBase 15 database. Nearly half of the small RNAs belonged to the let-7 family miRNA. Seven percent of the sequenced small RNAs were not annotated in miRBase 15. Bioinformatic analysis of the unannotated small RNA sequences suggested seventeen novel miRNA candidates with relatively high expression levels (>100 tags per million). The left:right expression ratios were similar for all highly expressed miRNAs with less than 10% differences. These results provide a basic idea of the relative expression strengths of known and unknown miRNAs in the dorsal hippocampal CA3.

Reduction of non-insert sequence reads by dimer eliminator LNA oligonucleotide for small RNA deep sequencing.

Here we describe a method for constructing small RNA libraries for high-throughput sequencing in which we have made a significant improvement to commonly available standard protocols. We added a locked nucleic acid (LNA) oligonucleotide--named dimer eliminator--that is complementary to the adapter-dimer ligation products during the reverse transcription reaction. It reduces adapter-dimers, which often contaminate standard libraries and increase the number of non-insert sequence reads. This simple technology can be used for simultaneous multiplex sequencing of various barcoded samples as well as nonbarcoded small RNA library sequencing. In this study we also evaluated the reproducibility and quantitative design of the eight barcoded tags by comparing the Pearson's correlation values in the expression analysis between each barcoded sample. This method improves the sequencing yield and efficiency, while simplifying library construction, and makes it easier to perform large-scale small RNA analysis under multiple conditions with next-generation sequencers.

Cross-mapping and the identification of editing sites in mature microRNAs in high-throughput sequencing libraries.

MicroRNAs (miRNAs) are short (20-23 nt) RNAs that are sequence-specific mediators of transcriptional and post-transcriptional regulation of gene expression. Modern high-throughput technologies enable deep sequencing of such RNA species on an unprecedented scale. We find that the analysis of small RNA deep-sequencing libraries can be affected by cross-mapping, in which RNA sequences originating from one locus are inadvertently mapped to another. Similar to cross-hybridization on microarrays, cross-mapping is prevalent among miRNAs, as they tend to occur in families, are similar or derived from repeat or structural RNAs, or are post-transcriptionally modified. Here, we develop a strategy to correct for cross-mapping, and apply it to the analysis of RNA editing in mature miRNAs. In contrast to previous reports, our analysis suggests that RNA editing in mature miRNAs is rare in animals.