Reproducibility of high-throughput mRNA and small RNA sequencing across laboratories.

RNA sequencing is an increasingly popular technology for genome-wide analysis of transcript sequence and abundance. However, understanding of the sources of technical and interlaboratory variation is still limited. To address this, the GEUVADIS consortium sequenced mRNAs and small RNAs of lymphoblastoid cell lines of 465 individuals in seven sequencing centers, with a large number of replicates. The variation between laboratories appeared to be considerably smaller than the already limited biological variation. Laboratory effects were mainly seen in differences in insert size and GC content and could be adequately corrected for. In small-RNA sequencing, the microRNA (miRNA) content differed widely between samples owing to competitive sequencing of rRNA fragments. This did not affect relative quantification of miRNAs. We conclude that distributing RNA sequencing among different laboratories is feasible, given proper standardization and randomization procedures. We provide a set of quality measures and guidelines for assessing technical biases in RNA-seq data.

Genetic profiling of osteoblast-like cells cultured on a novel bone reconstructive material, consisting of poly-L-lactide, carbon nanotubes and microhydroxyapatite, in the presence of bone morphogenetic protein-2.

In bone tissue engineering composite materials have been introduced, combining a degradable polymer matrix with, for instance, carbon nanotubes (CNTs) to improve mechanical properties or with microhydroxyapatite (µHA) to improve osteoconduction. The addition of bone morphogenetic protein-2 (BMP-2) can further improve the biological response to the material. However, the influence of such an elaborate composite formation on osteoprogenitor cells is unknown. To examine this, rat bone marrow (RBM) cells were cultured on porous poly-L-lactic acid and composite scaffolds, with or without added BMP-2. Cell proliferation and differentiation were studied using DNA, alkaline phosphatase and scanning electron microscopic analysis. Further, genetic profiles were examined by microarray investigation. Results showed that the composite scaffold had no significant effect on the proliferation of RBM cells, but indicated a negative effect on cell differentiation. The addition of BMP-2 also had no significant effect on the proliferation of RBM cells, but differentiation towards the osteogenic lineage was confirmed. In the arrays results, the addition of BMP-2 alone led to the expression of genes involved in (minor) inflammation. The composite scaffold, and even more distinctly the combination of the composite scaffold with BMP-2, led to the expression of genes, based on gene ontology, connected to tumorigenesis. Therefore, CNT- and µHA-containing composite materials are not recommended as a bone restorative material.

Evidence for induced fit in bacterial RNase P RNA-mediated cleavage.

RNase P with its catalytic RNA subunit is involved in the processing of a number of RNA precursors with different structures. However, precursor tRNAs are the most abundant substrates for RNase P. Available data suggest that a tRNA is folded into its characteristic structure already at the precursor state and that RNase P recognizes this structure. The tRNA D-/T-loop domain (TSL-region) is suggested to interact with the specificity domain of RNase P RNA while residues in the catalytic domain interact with the cleavage site. Here, we have studied the consequences of a productive interaction between the TSL-region and its binding site (TBS) in the specificity domain using tRNA precursors and various hairpin-loop model substrates. The different substrates were analyzed with respect to cleavage site recognition, ground-state binding, cleavage as a function of the concentration of Mg(2+) and the rate of cleavage under conditions where chemistry is suggested to be rate limiting using wild-type Escherichia coli RNase P RNA, M1 RNA, and M1 RNA variants with structural changes in the TBS-region. On the basis of our data, we conclude that a productive TSL/TBS interaction results in a conformational change in the M1 RNA substrate complex that has an effect on catalysis. Moreover, it is likely that this conformational change comprises positioning of chemical groups (and Mg(2+)) at and in the vicinity of the cleavage site. Hence, our findings are consistent with an induced-fit mechanism in RNase P RNA-mediated cleavage.

Synthesis of RNA using 2'-O-DTM protection.

tert-Butyldithiomethyl (DTM), a novel hydroxyl protecting group, cleavable under reductive conditions, was developed and applied for the protection of 2'-OH during solid-phase RNA synthesis. This function is compatible with all standard protecting groups used in oligonucleotide synthesis, and allows for fast and high-yield synthesis of RNA. Oligonucleotides containing the 2'-O-DTM groups can be easily deprotected under the mildest possible aqueous and homogeneous conditions. The preserved 5'-O-DMTr function can be used for high-throughput cartridge RNA purification.

The exocyclic amine at the RNase P cleavage site contributes to substrate binding and catalysis.

Most tRNAs carry a G at their 5' termini, i.e. at position +1. This position corresponds to the position immediately downstream of the site of cleavage in tRNA precursors. Here we studied RNase P RNA-mediated cleavage of substrates carrying substitutions/modifications at position +1 in the absence of the RNase P protein, C5, to investigate the role of G at the RNase P cleavage site. We present data suggesting that the exocyclic amine (2NH2) of G+1 contributes to cleavage site recognition, ground state binding and catalysis by affecting the rate of cleavage. This is in contrast to O6, N7 and 2'OH that are suggested to affect ground state binding and rate of cleavage to significantly lesser extent. We also provide evidence that the effects caused by the absence of 2NH2 at position +1 influenced the charge distribution and conceivably Mg2+ binding at the RNase P cleavage site. These findings are consistent with models where the 2NH2 at the cleavage site (when present) interacts with RNase P RNA and/or influences the positioning of Mg2+ in the vicinity of the cleavage site. Moreover, our data suggest that the presence of the base at +1 is not essential for cleavage but its presence suppresses miscleavage and dramatically increases the rate of cleavage. Together our findings provide reasons why most tRNAs carry a guanosine at their 5' end.

Lead(II) cleavage analysis of RNase P RNA in vivo.

The overall conformation of M1 RNA, the catalytic RNA subunit of RNase P in Escherichia coli, was analyzed in vivo and, in the presence of the C5 protein subunit, in vitro by lead(II) acetate probing. The partial cleavage patterns obtained are congruent with previous structure mapping performed in vitro. Most of the known major and minor cleavages in M1 RNA were supported and could be mapped onto a secondary structure model. The data obtained indicate that C5 has only minor effects on the overall structure of the RNA subunit. The similar cleavage patterns obtained in vitro and in vivo furthermore suggest that the intracellular environment does not greatly alter the overall conformation of M1 RNA within the holoenzyme complex. Moreover, our data indicate that M1 RNA in vivo is present in at least two states-the major fraction is bound to tRNA substrates and a minor fraction is substrate free. Finally, both in this and previous work we found that lead(II) probing data from in vivo experiments conducted on longer RNAs (tmRNA and M1 RNA) generally gives superior resolution compared to parallel in vitro experiments. This may reflect the absence of alternative conformers present in vitro and the more natural state of these RNAs in the cell due to proper, co-transcriptional folding pathways and possibly the presence of RNA chaperones.

Complexity in orchestration of chemical groups near different cleavage sites in RNase P RNA mediated cleavage.

RNase P mediated cleavage of the tRNA(His) precursor does not rely on the formation of the "+73/294 interaction" to give the correct cleavage product, i.e. cleavage at -1, while other tRNA precursors that are cleaved at the canonical site +1 do. A previous model, here referred to as the "2'OH-model", predicts that the 2'OH at the canonical cleavage site would affect cleavage at -1. Here we used model RNA hairpin substrates mimicking the structural architecture of the tRNA(His) precursor cleavage site to investigate the role of 2'OH with respect to ground state binding and rate of cleavage in the presence and absence of the +73/294 interaction. Our data emphasize the importance of the 2'OH in the immediate vicinity of the scissile bond. Moreover, introduction of 2'H at the cleavage site did not affect cleavage at an alternative cleavage site to any significant extent. Our findings are therefore inconsistent with the 2'OH model. We favor a model where the 2'OH at the cleavage site influence Mg2+ binding in its vicinity, however we do not exclude the possibility that the 2'OH at the cleavage site interacts with RNase P RNA. Studying the importance of the 2'OH at different cleavage sites also indicated a higher dependence on the 2'OH at the cleavage site in the absence of the +73/294 interaction than in its presence. Finally, we provide data suggesting that N3 of U at position -1 in the substrate is most likely not involved in an interaction with RNase P RNA.

Substrate discrimination in RNase P RNA-mediated cleavage: importance of the structural environment of the RNase P cleavage site.

Like the translational elongation factor EF-Tu, RNase P interacts with a large number of substrates where RNase P with its RNA subunit generates tRNAs with matured 5' termini by cleaving tRNA precursors immediately 5' of the residue at +1, i.e. at the position that corresponds to the first residue in tRNA. Most tRNAs carry a G+1C+72 base pair at the end of the aminoacyl acceptor-stem whereas in tRNA(Gln) G+1C+72 is replaced with U+1A+72. Here, we investigated RNase P RNA-mediated cleavage as a function of having G+1C+72 versus U+1A+72 in various substrate backgrounds, two full-size tRNA precursors (pre-tRNA(Gln) and pre-tRNA(Tyr)Su3) and a model RNA hairpin substrate (pATSer). Our data showed that replacement of G+1C+72 with U+1A+72 influenced ground state binding, cleavage efficiency under multiple and single turnover conditions in a substrate-dependent manner. Interestingly, we observed differences both in ground state binding and rate of cleavage comparing two full-size tRNA precursors, pre-tRNA(Gln) and pre-tRNA(Tyr)Su3. These findings provide evidence for substrate discrimination in RNase P RNA-mediated cleavage both at the level of binding, as previously observed for EF-Tu, as well as at the catalytic step. In our experiments where we used model substrate derivatives further indicated the importance of the +1/+72 base pair in substrate discrimination by RNase P RNA. Finally, we provide evidence that the structural architecture influences Mg2+ binding, most likely in its vicinity.

Cross talk between the +73/294 interaction and the cleavage site in RNase P RNA mediated cleavage.

To monitor functionally important metal ions and possible cross talk in RNase P RNA mediated cleavage we studied cleavage of substrates, where the 2'OH at the RNase P cleavage site (at -1) and/or at position +73 had been replaced with a 2' amino group (or 2'H). Our data showed that the presence of 2' modifications at these positions affected cleavage site recognition, ground state binding of substrate and/or rate of cleavage. Cleavage of 2' amino substituted substrates at different pH showed that substitution of Mg2+ by Mn2+ (or Ca2+), identity of residues at and near the cleavage site, and addition of C5 protein influenced the frequency of miscleavage at -1 (cleavage at the correct site is referred to as +1). From this we infer that these findings point at effects mediated by protonation/deprotonation of the 2' amino group, i.e. an altered charge distribution, at the site of cleavage. Moreover, our data suggested that the structural architecture of the interaction between the 3' end of the substrate and RNase P RNA influence the charge distribution at the cleavage site as well as the rate of cleavage under conditions where the chemistry is suggested to be rate limiting. Thus, these data provide evidence for cross talk between the +73/294 interaction and the cleavage site in RNase P RNA mediated cleavage. We discuss the role metal ions might play in this cross talk and the likelihood that at least one functionally important metal ion is positioned in the vicinity of, and use the 2'OH at the cleavage site as an inner or outer sphere ligand.

Importance of the +73/294 interaction in Escherichia coli RNase P RNA substrate complexes for cleavage and metal ion coordination.

We have studied an interaction, the "73/294-interaction", between residues 294 in M1 RNA (the catalytic subunit of Escherichia coli RNase P) and +73 in the tRNA precursor substrate. The 73/294-interaction is part of the "RCCA-RNase P RNA interaction", which anchors the 3' R(+73)CCA-motif of the substrate to M1 RNA (interacting residues underlined). Considering that in a large fraction of tRNA precursors residue +73 is base-paired to nucleotide -1 immediately 5' of the cleavage site, formation of the 73/294-interaction results in exposure of the cleavage site. We show that the nature/orientation of the 73/294-interaction is important for cleavage site recognition and cleavage efficiency. Our data further suggest that this interaction is part of a metal ion-binding site and that specific chemical groups are likely to act as ligands in binding of Mg(2+) or other divalent cations important for function. We argue that this Mg(2+) is involved in metal ion cooperativity in M1 RNA-mediated cleavage. Moreover, we suggest that the 73/294-interaction operates in concert with displacement of residue -1 in the substrate to ensure efficient and correct cleavage. The possibility that the residue at -1 binds to a specific binding surface/pocket in M1 RNA is discussed. Our data finally rationalize why the preferred residue at position 294 in M1 RNA is U.

The residue immediately upstream of the RNase P cleavage site is a positive determinant.

We have studied the importance of the residue at the position immediately upstream of the RNase P RNA cleavage site using model substrates that mimic the structure at and near the cleavage site of the tRNA(His) precursor. The various model substrates were studied with respect to cleavage site recognition as well as the kinetics of cleavage using M1 RNA, the catalytic subunit of Escherichia coli RNase P. Our studies showed that the identity of the residue immediately upstream of the cleavage site critically influences both these aspects. Among the ones tested, U is the preferred nucleotide at this position. Hence, these findings rationalize why most bacterial tRNA(His) genes/transcripts harbor a U immediately upstream of the RNase P cleavage site and extend our understanding of the cleavage site recognition process in general and the unusual cleavage of the tRNA(His) precursor in particular. Based on our as well as the data of others, we suggest that the nucleotide immediately upstream of the cleavage site is a positive determinant for cleavage by RNase P in general and the expression of tRNA genes is influenced by structural elements localized outside the promoter region i.e. in the leader and spacer regions of tRNA transcripts.