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

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

Broadly applicable oligonucleotide mass spectrometry for the analysis of RNA writers and erasers in vitro.

RNAs are post-transcriptionally modified by dedicated writer or eraser enzymes that add or remove specific modifications, respectively. Mass spectrometry (MS) of RNA is a useful tool to study the modification state of an oligonucleotide (ON) in a sensitive manner. Here, we developed an ion-pairing reagent free chromatography for positive ion detection of ONs by low- and high-resolution MS, which does not interfere with other types of small compound analyses done on the same instrument. We apply ON-MS to determine the ONs from an RNase T1 digest of in vitro transcribed tRNA, which are purified after ribozyme-fusion transcription by automated size exclusion chromatography. The thus produced tRNAValAAC is substrate of the human tRNA ADAT2/3 enzyme and we confirm the deamination of adenosine to inosine and the formation of tRNAValIACin vitro by ON-MS. Furthermore, low resolution ON-MS is used to monitor the demethylation of ONs containing 1-methyladenosine by bacterial AlkB in vitro. The power of high-resolution ON-MS is demonstrated by the detection and mapping of modified ONs from native total tRNA digested with RNase T1. Overall, we present an oligonucleotide MS method which is broadly applicable to monitor in vitro RNA (de-)modification processes and native RNA.

Purification, characterization and cytotoxic activities of individual tRNAs from Escherichia coli.

Transfer RNAs (tRNAs) are the most abundant class in small non-coding RNAs which have been proved to be pharmacologically active. In the present study, we evaluated the potential anticancer activities of tRNAs from Escherichia coli MRE 600 to investigate the relationship between non-pathogenic Escherichia coli strain and colorectal cancer. To purify individual tRNAs, we firstly developed a two-dimensional liquid chromatography (2D-LC) and successfully obtained two pure tRNAs. Nuclease mediated base-specific digestions coupled with UHPLC-MS/MS techniques led to an identification of these two tRNAs as tRNA-Val(UAC) and tRNA-Leu(CAG) with typical cloverleaf-like secondary structure. MTT assay demonstrated that both tRNA-1 and tRNA-2 exhibit strong cytotoxicity with IC50 of 113.0 nM and 124.8 nM on HCT-8 cells in a dose-dependent manner. Further clonogenic assay revealed that the purified tRNAs exhibit significant inhibition in colony formation with survival percentage of 79.0 ±â€¯1.6 and 71.2 ±â€¯2.2 at the concentration of 100 nM. These findings provided evidences of anticancer activities of tRNAs from non-pathogenic Escherichia coli strain, indicating that the pharmacological effects of these neglected biomacromolecules from microorganisms should be emphasized. This study put new insights into the therapeutic effects of intestinal microorganism on human diseases, therefore broadened our knowledge of the biological functions of gut microbiota.

The complete mitochondrial genome sequence of an endometrial cancer inbred Donryu rat model.

The Donryu rat strain is a commonly used model for endometrial cancer disease study. We sequenced this rat strain mitochondrial genome for the first time (GenBank Accession No. KM114605). Its mitogenome was 16,307 bp and coding 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes. A total of 96 SNPs were examined when compared to reference BN sequence.

Complete mitochondrial genome sequence and mutations of the prostate cancer model inbred Sprague-Dawley strain.

In the present work we undertook the complete mitochondrial genome sequencing of a important prostate cancer model inbred Sprague-Dawley strain for the first time. The total length of the mitogenome was 16,308 bp. It harbored 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes and one non-coding control region (D-loop region). The mutation events were also reported.

Transfer RNA detection by small RNA deep sequencing and disease association with myelodysplastic syndromes.

BACKGROUND: Although advances in sequencing technologies have popularized the use of microRNA (miRNA) sequencing (miRNA-seq) for the quantification of miRNA expression, questions remain concerning the optimal methodologies for analysis and utilization of the data. The construction of a miRNA sequencing library selects RNA by length rather than type. However, as we have previously described, miRNAs represent only a subset of the species obtained by size selection. Consequently, the libraries obtained for miRNA sequencing also contain a variety of additional species of small RNAs. This study looks at the prevalence of these other species obtained from bone marrow aspirate specimens and explores the predictive value of these small RNAs in the determination of response to therapy in myelodysplastic syndromes (MDS). METHODS: Paired pre and post treatment bone marrow aspirate specimens were obtained from patients with MDS who were treated with either azacytidine or decitabine (24 pre-treatment specimens, 23 post-treatment specimens) with 22 additional non-MDS control specimens. Total RNA was extracted from these specimens and submitted for next generation sequencing after an additional size exclusion step to enrich for small RNAs. The species of small RNAs were enumerated, single nucleotide variants (SNVs) identified, and finally the differential expression of tRNA-derived species (tDRs) in the specimens correlated with diseasestatus and response to therapy. RESULTS: Using miRNA sequencing data generated from bone marrow aspirate samples of patients with known MDS (N = 47) and controls (N = 23), we demonstrated that transfer RNA (tRNA) fragments (specifically tRNA halves, tRHs) are one of the most common species of small RNA isolated from size selection. Using tRNA expression values extracted from miRNA sequencing data, we identified six tRNA fragments that are differentially expressed between MDS and normal samples. Using the elastic net method, we identified four tRNAs-derived small RNAs (tDRs) that together can explain 67 % of the variation in treatment response for MDS patients. Similar analysis of specifically mitochondrial tDRs (mt-tDRs) identified 13 mt-tDRs which distinguished disease status in the samples and a single mt-tDR which predited response. Finally, 14 SNVs within the tDRs were found in at least 20 % of the MDS samples and were not observed in any of the control specimens. DISCUSSION: This study highlights the prevalence of tDRs in RNA-seq studies focused on small RNAs. The potential etiologies of these species, both technical and biologic, are discussed as well as important challenges in the interpretation of tDR data. CONCLUSIONS: Our analysis results suggest that tRNA fragments can be accurately detected through miRNA sequencing data and that the expression of these species may be useful in the diagnosis of MDS and the prediction of response to therapy.

Identification and functional characterization of tRNA-derived RNA fragments (tRFs) in respiratory syncytial virus infection.

The discovery of small noncoding RNAs (sncRNAs) with regulatory functions is a recent breakthrough in biology. Among sncRNAs, microRNA (miRNA), derived from host or virus, has emerged as elements with high importance in control of viral replication and host responses. However, the expression pattern and functional aspects of other types of sncRNAs, following viral infection, are unexplored. In order to define expression patterns of sncRNAs, as well as to discover novel regulatory sncRNAs in response to viral infection, we applied deep sequencing to cells infected with human respiratory syncytial virus (RSV), the most common cause of bronchiolitis and pneumonia in babies. RSV infection leads to abundant production of transfer RNA (tRNA)-derived RNA Fragments (tRFs) that are ~30 nucleotides (nts) long and correspond to the 5'-half of mature tRNAs. At least one tRF, which is derived from tRNA-Glu-CTC, represses target mRNA in the cytoplasm and promotes RSV replication. This demonstrates that this tRF is not a random by-product of tRNA degradation but a functional molecule. The biogenesis of this tRF is also specific, as it is mediated by the endonuclease angiogenin (ANG), not by other nucleases. In summary, our study presents novel information on the induction of a functional tRF by viral infection.

Biosynthesis of 4-thiouridine in tRNA in the methanogenic archaeon Methanococcus maripaludis.

4-Thiouridine (s(4)U) is a conserved modified nucleotide at position 8 of bacterial and archaeal tRNAs and plays a role in protecting cells from near-UV killing. Escherichia coli employs the following two enzymes for its synthesis: the cysteine desulfurase IscS, which forms a Cys persulfide enzyme adduct from free Cys; and ThiI, which adenylates U8 and transfers sulfur from IscS to form s(4)U. The C-terminal rhodanese-like domain (RLD) of ThiI is responsible for the sulfurtransferase activity. The mechanism of s(4)U biosynthesis in archaea is not known as many archaea lack cysteine desulfurase and an RLD of the putative ThiI. Using the methanogenic archaeon Methanococcus maripaludis, we show that deletion of ThiI (MMP1354) abolished the biosynthesis of s(4)U but not of thiamine. MMP1354 complements an Escherichia coli ΔthiI mutant for s(4)U formation, indicating that MMP1354 is sufficient for sulfur incorporation into s(4)U. In the absence of an RLD, MMP1354 uses Cys(265) and Cys(268) located in the PP-loop pyrophosphatase domain to generate persulfide and disulfide intermediates for sulfur transfer. In vitro assays suggest that S(2-) is a physiologically relevant sulfur donor for s(4)U formation catalyzed by MMP1354 (K(m) for Na(2)S is ∼1 mm). Thus, methanogenic archaea developed a strategy for sulfur incorporation into s(4)U that differs from bacteria; this may be an adaptation to life in sulfide-rich environments.

Identification of N6,N6-dimethyladenosine in transfer RNA from Mycobacterium bovis Bacille Calmette-Guérin.

There are more than 100 different ribonucleoside structures incorporated as post-transcriptional modifications, mainly in tRNA and rRNA of both prokaryotes and eukaryotes, and emerging evidence suggests that these modifications function as a system in the translational control of cellular responses. However, our understanding of this system is hampered by the paucity of information about the complete set of RNA modifications present in individual organisms. To this end, we have employed a chromatography-coupled mass spectrometric approach to define the spectrum of modified ribonucleosides in microbial species, starting with Mycobacterium bovis BCG. This approach revealed a variety of ribonucleoside candidates in tRNA from BCG, of which 12 were definitively identified based on comparisons to synthetic standards and 5 were tentatively identified by exact mass comparisons to RNA modification databases. Among the ribonucleosides observed in BCG tRNA was one not previously described in tRNA, which we have now characterized as N6,N6-dimethyladenosine.

Deep sequencing of human nuclear and cytoplasmic small RNAs reveals an unexpectedly complex subcellular distribution of miRNAs and tRNA 3' trailers.

BACKGROUND: MicroRNAs (miRNAs) are approximately 22-nt small non-coding regulatory RNAs that have generally been considered to regulate gene expression at the post-transcriptional level in the cytoplasm. However, recent studies have reported that some miRNAs localize to and function in the nucleus. METHODOLOGY/PRINCIPAL FINDINGS: To determine the number of miRNAs localized to the nucleus, we systematically investigated the subcellular distribution of small RNAs (sRNAs) by independent deep sequencing sequenced of the nuclear and cytoplasmic pools of 18- to 30-nucleotide sRNAs from human cells. We identified 339 nuclear and 324 cytoplasmic known miRNAs, 300 of which overlap, suggesting that the majority of miRNAs are imported into the nucleus. With the exception of a few miRNAs evidently enriched in the nuclear pool, such as the mir-29b, the ratio of miRNA abundances in the nuclear fraction versus in the cytoplasmic fraction vary to some extent. Moreover, our results revealed that a large number of tRNA 3' trailers are exported from the nucleus and accumulate in the cytoplasm. These tRNA 3' trailers accumulate in a variety of cell types, implying that the biogenesis of tRNA 3' trailers is conserved and that they have a potential functional role in vertebrate cells. CONCLUSION/SIGNIFICANCE: Our results provide the first comprehensive view of the subcellular distribution of diverse sRNAs and new insights into the roles of miRNAs and tRNA 3' trailers in the cell.

Aquifex aeolicus tRNA (N2,N2-guanine)-dimethyltransferase (Trm1) catalyzes transfer of methyl groups not only to guanine 26 but also to guanine 27 in tRNA.

Transfer RNA (N2,N2-guanine)-dimethyltransferase (Trm1) catalyzes N2,N2-dimethylguanine formation at position 26 (m(2)(2)G26) in tRNA. In the reaction, N2-guanine at position 26 (m(2)G26) is generated as an intermediate. The trm1 genes are found only in archaea and eukaryotes, although it has been reported that Aquifex aeolicus, a hyper-thermophilic eubacterium, has a putative trm1 gene. To confirm whether A. aeolicus Trm1 has tRNA methyltransferase activity, we purified recombinant Trm1 protein. In vitro methyl transfer assay revealed that the protein has a strong tRNA methyltransferase activity. We confirmed that this gene product is expressed in living A. aeolicus cells and that the enzymatic activity exists in cell extract. By preparing 22 tRNA transcripts and testing their methyl group acceptance activities, it was demonstrated that this Trm1 protein has a novel tRNA specificity. Mass spectrometry analysis revealed that it catalyzes methyl transfers not only to G26 but also to G27 in substrate tRNA. Furthermore, it was confirmed that native tRNA(Cys) has an m(2)(2)G26m(2)G27 or m(2)(2)G26m(2)(2)G27 sequence, demonstrating that these modifications occur in living cells. Kinetic studies reveal that the m2G26 formation is faster than the m(2)G27 formation and that disruption of the G27-C43 base pair accelerates velocity of the G27 modification. Moreover, we prepared an additional 22 mutant tRNA transcripts and clarified that the recognition sites exist in the T-arm structure. This long distance recognition results in multisite recognition by the enzyme.

A modified DNA isolation protocol for obtaining pure RT-PCR grade RNA.

We provide a simple but very efficient method for RNA preparation from Saccharomyces cerevisiae based on a standard chromosomal DNA isolation protocol. The method yields DNA-free total RNA, including mRNA, rRNA, and tRNA but can easily be adjusted to considerably enrich low molecular weight RNAs, such as tRNAs and the small rRNA species (5S and 5.8S). The procedure was proven and validated by verification of cDNAs belonging to four different genes, two of which encoding polypeptides and two tRNA genes. Besides its simplicity, the method is further advantageous in terms of safety (omitting hazardous phenol) and cost efficiency.

The many applications of acid urea polyacrylamide gel electrophoresis to studies of tRNAs and aminoacyl-tRNA synthetases.

Here we describe the many applications of acid urea polyacrylamide gel electrophoresis (acid urea PAGE) followed by Northern blot analysis to studies of tRNAs and aminoacyl-tRNA synthetases. Acid urea PAGE allows the electrophoretic separation of different forms of a tRNA, discriminated by changes in bulk, charge, and/or conformation that are brought about by aminoacylation, formylation, or modification of a tRNA. Among the examples described are (i) analysis of the effect of mutations in the Escherichia coli initiator tRNA on its aminoacylation and formylation; (ii) evidence of orthogonality of suppressor tRNAs in mammalian cells and yeast; (iii) analysis of aminoacylation specificity of an archaeal prolyl-tRNA synthetase that can aminoacylate archaeal tRNA(Pro) with cysteine, but does not aminoacylate archaeal tRNA(Cys) with cysteine; (iv) identification and characterization of the AUA-decoding minor tRNA(Ile) in archaea; and (v) evidence that the archaeal minor tRNA(Ile) contains a modified base in the wobble position different from lysidine found in the corresponding eubacterial tRNA.

Diet-dependent depletion of queuosine in tRNAs in Caenorhabditis elegans does not lead to a developmental block.

Queuosine (Q), a hypermodified nucleoside,occurs at the wobble position of transfer RNAs (tRNAs)with GUN anticodons. In eubacteria, absence of Q affects messenger RNA (mRNA) translation and reduces the virulence of certain pathogenic strains. In animal cells,changes in the abundance of Q have been shown to correlate with diverse phenomena including stress tolerance, cell proliferation and tumour growth but the function of Q in animals is poorly understood. Animals are thought to obtain Q (or its analogues) as a micronutrient from dietary sources such as gut micro flora. However,the difficulty of maintaining animals under bacteria-free conditions on Q-deficient diets has severely hampered the study of Q metabolism and function in animals. In this study,we show that as in higher animals, tRNAs in the nematode Caenorhabditis elegans are modified by Q and its sugar derivatives. When the worms were fed on Q-deficient Escherichia coli, Q modification was absent from the worm tRNAs suggesting that C.elegans lacks a de novo pathway of Q biosynthesis. The inherent advantages of C.elegans as a model organism, and the simplicity of conferring a Q-deficient phenotype on it make it an ideal system to investigate the function of Q modification in tRNA.

Early gene expression profiles during intraoperative myocardial ischemia-reperfusion in cardiac surgery.

OBJECTIVE: The effects of cold cardioplegic arrest and reperfusion on human ventricular gene expression are unknown. We tested the hypothesis that intraoperative ischemia-reperfusion under conditions of blood cardioplegic arrest would induce a unique myocardial genomic profile indicative of a cardioprotective response. METHODS: Right ventricular samples were serially acquired during surgical repair of ventricular septal defect. RESULTS: Expression profiling revealed 3 patterns of gene expression: (1) increased expression above control levels within 1 hour of cardioplegic arrest, with further amplification during early reperfusion; (2) increased expression limited to the reperfusion phase; and (3) reduced expression during reperfusion. Functional annotation and network mapping of differentially expressed genes indicated activation of multiple signaling pathways regulated by phosphatidylinositide 3'-OH kinase convergent on cellular growth and reparative programs. Also observed was increased expression of genes regulating hemoglobin synthesis, suggesting a novel cardioprotective pathway evoked during ischemia-reperfusion. CONCLUSION: Reversible myocardial ischemia-reperfusion during cardiac surgery is associated with an immediate genomic response that predicts a net cardioprotective phenotype.

Determination of queuosine derivatives by reverse-phase liquid chromatography for the hypomodification study of Q-bearing tRNAs from various mammal liver cells.

Three queuosine derivatives (Q-derivatives) have been found at position 34 of four mammalian so-called Q-tRNAs: queuosine (Q) in tRNA(Asn) and tRNA(His), mannosyl-queuosine (manQ) in tRNA(Asp), and galactosyl-queuosine (galQ) in tRNA(Tyr). An analytical procedure based on the combined means of purified tRNA isolation from liver cells and ribonucleoside analysis by reverse-phase high performance liquid chromatography coupled with real-time UV-spectrometry (RPLC-UV) was developed for the quantitative analysis of the three Q-derivatives present in total tRNA from liver tissues and liver cell cultures. Using this analytical procedure, the rates of Q-tRNA modification were studied in total tRNAs from various mammalian hepatic cells. Our results show that the four Q-tRNAs are fully modified in liver tissues from adult mammals, regardless of the mammal species. However, a lack in the Q-modification level was observed in Q-tRNAs from newborn rat liver, as well in Q-tRNAs from normal rat liver cell cultures growing in a low queuine content medium, and from a rat hepatoma cell line. It is noteworthy that in all cases of Q-tRNA hypomodification, our analytical procedure showed that tRNA(Asp) is always the least affected by the hypomodification. The biological significance of this phenomenon is discussed.

Identification and sequencing of two isopentenyladenosine-modified transfer RNAs from Chinese hamster ovary cells.

To determine the presence and identity of isopentenyladenosine-containing transfer RNAs (tRNAs) in a mammalian cell line, we adopted a novel method to isolate, clone and sequence these RNAs. This method was based on 3' polyadenylation of the tRNA prior to cDNA synthesis, PCR amplification, cloning and DNA sequencing. Using this unique procedure, we report the cloning and sequencing of the selenocysteine-tRNA and mitochondrial tryptophan-tRNA from Chinese hamster ovary cells which contain this specific tRNA modification. This new method will be useful in the identification of other tRNAs and other small RNAs where the primary sequence is unknown.

Identification of the tRNAs which up-regulate agrostin, barley RIP and PAP-S, three ribosome-inactivating proteins of plant origin.

Ribosome-inactivating proteins (RIP) are RNA-N-glycosidases widely diffused in plants which depurinate ribosomal RNA at a specific universally conserved position, A4324 in rat ribosomes. A small group of RIPs (cofactor-dependent RIPs) require ATP and tRNA to reach maximal activity on isolated ribosomes. The tRNA which stimulates gelonin was identified as tRNA(Trp). The present paper reports the identification of three other tRNAs which stimulate agrostin (tRNA(Ala)), barley RIP (tRNA(Ala), tRNA(Val)) and PAP-S (tRNA(Gly)), while for tritin-S no particular stimulating tRNA emerged. The sequences of tRNA(Val) and tRNA(Gly) correspond to the already known ones (rabbit and man, respectively). The tRNA(Ala) (anticodon IGC) identifies a new isoacceptor. Only the stimulating activity of the tRNA(Ala) for agrostin approaches the specificity previously observed for the couple gelonin-tRNA(Trp).

Inhibition and substrate competition kinetics in analysis of porcine thyroid alkaline ribonuclease's specificity toward synthetic RNA's and tRNA.

Inhibition and substrate competition kinetics demonstrated that tRNA is a highly preferred substrate of thyroid alkaline RNase. The pyrimidine-specific RNase cleaved poly(C) 2.8 x 10(5) faster than poly(U). kcat:K(M) ratios for tRNA and poly(C) based on molecular weights failed to predict preference when both were present. Competition experiments between poly(C) and tRNA revealed tRNA was a tight-binding competing substrate and the cytidylate residues in the 3'-CCA terminus to tRNA were preferred about 280:1 over those in poly(C). Poly(U) was competitive with tRNA. When poly(C) was the substrate, inhibition type by poly(G) depended on poly(G) concentration. Neither tRNA lacking its 3' terminal cytidylyl(3'-5')adenosine and terminating in a 2':3' cCMP residue, tRNA lacking its 3' terminal 5'AMP residue, guanosine, nor guanylyl(3'-5')guanylyl(3'-5')guanosine were inhibitors. Product inhibition by adenosine and 2':3' cCMP showed the kinetic mechanism for cleavage of tRNA was ordered uni bi.

Nucleotides of tRNA (Glu) involved in recognition by barley chloroplast glutamyl-tRNA synthetase and glutamyl-tRNA reductase.

The biosynthesis of delta-aminolevulinate (ALA), via the C-5 pathway, requires tRNA(Glu) as a cofactor for the glutamyl tRNA(Glu) synthetase and the glutamyl tRNA(Glu) reductase which are the first two enzymes in this three step pathway. These two enzymes form a ternary complex with the tRNA(Glu) in Chlamydomonas reinhardtii suggesting that the recognition elements on the tRNA cofactor are different for each enzyme. Chemical modification and comparative studies with tRNA(Glu)s from a number of species were used to determine the nucleotides involved in the recognition of the barley chloroplast tRNA(Glu) by the barley enzymes. The barley chloroplast tRNA(Glu) is chemically modified both before and after ligation to glutamate with monobromobimane or CNBr. The chemically modified tRNA(Glu) is a poor substrate for the glutamyl-tRNA synthetase and the chemically modified glutamyl-tRNA(Glu) is used as a substrate for glutamyl-tRNA(Glu) reductase. The tRNA(Glu) from the chloroplasts if barley, Chlamydomonas reinhardtii, tobacco, cucumber, wheat and spinach and tRNA(Glu) from Synechocystis PCC6803, Escherichia coli, barley germ and bakers yeast and the barley chloroplast tRNA(Gln) are all effective substrates for the barley chloroplast glutamyl-tRNA synthetase. A comparison of the sequences of these tRNAs shows 19 conserved bases and five of these bases, G10, A26, U34, U35 and A37 are suggested as recognition elements of barley glutamyl tRNA(Glu) synthetase by assuming a similar binding orientation as in the crystal structure of the E. coli tRNA(Gln) GlnRS complex. The glutamyl-tRNA(Glu) from E. coli, bakers yeast and barley germ and the barley chloroplast glutamyl-tRNA(Gln) are not effective substrates for the barley chloroplast glutamyl-tRNA(Glu) reductase. A comparison of the sequences of these four tRNA species with the sequences of the tRNA(Glu) species that can be used as substrate by the glutamyl-tRNA(Glu) reductase yields seven common differences in the primary sequence. These 7 nucleotides, A7-U66, U29-A41, A53-U61, and U72 are expected to be required for recognition by the barley chloroplast glutamyl-tRNA(Glu) reductase.