Simultaneous characterization of cellular RNA structure and function with in-cell SHAPE-Seq

Many non-coding RNAs form structures that interact with cellular machinery to control gene expression. A central goal of molecular and synthetic biology is to uncover design principles linking RNA structure to function to understand and engineer this relationship. Here we report a simple, high-throughput method called in-cell SHAPE-Seq that combines in-cell probing of RNA structure with a measurement of gene expression to simultaneously characterize RNA structure and function in bacterial cells. We use in-cell SHAPE-Seq to study the structure–function relationship of two RNA mechanisms that regulate translation in Escherichia coli. We find that nucleotides that participate in RNA–RNA interactions are highly accessible when their binding partner is absent and that changes in RNA structure due to RNA–RNA interactions can be quantitatively correlated to changes in gene expression. We also characterize the cellular structures of three endogenously expressed non-coding RNAs: 5S rRNA, RNase P and the btuB riboswitch. Finally, a comparison between in-cell and in vitro folded RNA structures revealed remarkable similarities for synthetic RNAs, but significant differences for RNAs that participate in complex cellular interactions. Thus, in-cell SHAPE-Seq represents an easily approachable tool for biologists and engineers to uncover relationships between sequence, structure and function of RNAs in the cell.     

Bacterial community structure analysis of a hot spring soil by next generation sequencing of ribosomal RNA

In the present study the whole bacterial community structure of Tapovan hot spring soil located in the state of Uttarakhand, India was analysed through next generation sequencing. The hot spring soil is slightly alkaline in nature with abundance of sulphur. The spring soil was rich in various metallic and non-metallic elements required for bacterial survival. The community was found to comprise of 14 bacterial phyla with representation of members belonging to Firmicutes, Proteobacteria, Thermi, Bacteroidetes, Aquificae, Actinobacteria, chloroflexi, TM7, Fusobacteria etc. At the genus level Bacillus, Pseudomonas, Symbiobacterium, Thermus, Geobacillus, Anoxybacillus were found in abundance as compared to other genera like Flavobacterium, Ureibacillus, Clostridium, Meiothermus, Acinetobacter, Desulfotomaculum and Rheinheimera.     

High-throughput RNA interference screening using pooled shRNA libraries and next generation sequencing

RNA interference (RNAi) screening is a state-of-the-art technology that enables the dissection of biological processes and disease-related phenotypes. The commercial availability of genome-wide, short hairpin RNA (shRNA) libraries has fueled interest in this area but the generation and analysis of these complex data remain a challenge. Here, we describe complete experimental protocols and novel open source computational methodologies, shALIGN and shRNAseq, that allow RNAi screens to be rapidly deconvoluted using next generation sequencing. Our computational pipeline offers efficient screen analysis and the flexibility and scalability to quickly incorporate future developments in shRNA library technology.     

RNA folding with soft constraints: reconciliation of probing data and thermodynamic secondary structure prediction

Thermodynamic folding algorithms and structure probing experiments are commonly used to determine the secondary structure of RNAs. Here we propose a formal framework to reconcile information from both prediction algorithms and probing experiments. The thermodynamic energy parameters are adjusted using 'pseudo-energies' to minimize the discrepancy between prediction and experiment. Our framework differs from related approaches that used pseudo-energies in several key aspects. (i) The energy model is only changed when necessary and no adjustments are made if prediction and experiment are consistent. (ii) Pseudo-energies remain biophysically interpretable and hold positional information where experiment and model disagree. (iii) The whole thermodynamic ensemble of structures is considered thus allowing to reconstruct mixtures of suboptimal structures from seemingly contradicting data. (iv) The noise of the energy model and the experimental data is explicitly modeled leading to an intuitive weighting factor through which the problem can be seen as folding with 'soft' constraints of different strength. We present an efficient algorithm to iteratively calculate pseudo-energies within this framework and demonstrate how this approach can be used in combination with SHAPE chemical probing data to improve secondary structure prediction. We further demonstrate that the pseudo-energies correlate with biophysical effects that are known to affect RNA folding such as chemical nucleotide modifications and protein binding.     

The VSGB 2.0 model: a next generation energy model for high resolution protein structure modeling

A novel energy model (VSGB 2.0) for high resolution protein structure modeling is described, which features an optimized implicit solvent model as well as physics‐based corrections for hydrogen bonding, π–π interactions, self‐contact interactions, and hydrophobic interactions. Parameters of the VSGB 2.0 model were fit to a crystallographic database of 2239 single side chain and 100 11–13 residue loop predictions. Combined with an advanced method of sampling and a robust algorithm for protonation state assignment, the VSGB 2.0 model was validated by predicting 115 super long loops up to 20 residues. Despite the dramatically increasing difficulty in reconstructing longer loops, a high accuracy was achieved: all of the lowest energy conformations have global backbone RMSDs better than 2.0 Å from the native conformations. Average global backbone RMSDs of the predictions are 0.51, 0.63, 0.70, 0.62, 0.80, 1.41, and 1.59 Å for 14, 15, 16, 17, 18, 19, and 20 residue loop predictions, respectively. When these results are corrected for possible statistical bias as explained in the text, the average global backbone RMSDs are 0.61, 0.71, 0.86, 0.62, 1.06, 1.67, and 1.59 Å. Given the precision and robustness of the calculations, we believe that the VSGB 2.0 model is suitable to tackle “real” problems, such as biological function modeling and structure‐based drug discovery. Proteins 2011; © 2011 Wiley‐Liss, Inc.     

Temperature dependent chemical and enzymatic probing of the tRNA-like structure of TYMV RNA

In this paper we report on the thermal unfolding of the tRNA-like structure present at the 3' end of turnip yellow mosaic virus (TYMV) RNA. Diethyl pyrocarbonate (DEP), sodium bisulphite, nuclease S1 and ribonuclease T1 were used as structure probes at a broad range of temperatures. In this way most of the nucleotides present in the tRNA-like moiety were analysed. The melting behaviour of both secondary and tertiary interactions could be followed on the basis of the temperature dependent accessibility of the individual nucleotides or bases towards the various probes. The three-dimensional model of the tRNA-like domain (Dumas et al., J. Biomol. Struct. and Dyn. 4, 707 (1987] was supported by the results to a large extent. The interactions occurring between the T- and D-loop appear to be more complex than proposed in the latter model. Additional evidence for the presence of the RNA pseudoknot (Rietveld et al., Nucleic Acids Res. 10, 1929 (1982] was derived from the fact that the three coaxially stacked helical segments in the aminoacylacceptor arm displayed different melting transitions under certain experimental conditions. Aspects of melting behaviour and thermal stability of double helical regions within the tRNA-like structure are discussed, as well as the applicability of nucleases and modifying reagents at various temperatures in the analysis of RNA structure.     

Characteristic chemical probing patterns of loop motifs improve prediction accuracy of RNA secondary structures

RNA structures play a fundamental role in nearly every aspect of cellular physiology and pathology. Gaining insights into the functions of RNA molecules requires accurate predictions of RNA secondary structures. However, the existing thermodynamic folding models remain less accurate than desired, even when chemical probing data, such as selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE) reactivities, are used as restraints. Unlike most SHAPE-directed algorithms that only consider SHAPE restraints for base pairing, we extract two-dimensional structural features encoded in SHAPE data and establish robust relationships between characteristic SHAPE patterns and loop motifs of various types (hairpin, internal, and bulge) and lengths (2–11 nucleotides). Such characteristic SHAPE patterns are closely related to the sugar pucker conformations of loop residues. Based on these patterns, we propose a computational method, SHAPELoop, which refines the predicted results of the existing methods, thereby further improving their prediction accuracy. In addition, SHAPELoop can provide information about local or global structural rearrangements (including pseudoknots) and help researchers to easily test their hypothesized secondary structures.     

Dynamic interactions between RNA viruses and human hosts unravelled by a decade of next generation sequencing

Background

Next generation sequencing (NGS) methods have significantly contributed to a paradigm shift in genomic research for nearly a decade now. These methods have been useful in studying the dynamic interactions between RNA viruses and human hosts.

NG6: Integrated next generation sequencing storage and processing environment

ABSTRACT: BACKGROUND: Next generation sequencing platforms are now well implanted in sequencing centres and some laboratories. Upcoming smaller scale machines such as the 454 junior from Roche or the MiSeq from Illumina will increase the number of laboratories hosting a sequencer. In such a context, it is important to provide these teams with an easily manageable environment to store and process the produced reads. RESULTS: We describe a user-friendly information system able to manage large sets of sequencing data. It includes, on one hand, a workflow environment already containing pipelines adapted to different input formats (sff, fasta, fastq and qseq), different sequencers (Roche 454, Illumina HiSeq) and various analyses (quality control, assembly, alignment, diversity studies,...) and, on the other hand, a secured web site giving access to the results. The connected user will be able to download raw and processed data and browse through the analysis result statistics. The provided workflows can easily be modified or extended and new ones can be added. Ergatis is used as a workflow building, running and monitoring system. The analyses can be run locally or in a cluster environment using Sun Grid Engine. CONCLUSIONS: NG6 is a complete information system designed to answer the needs of a sequencing platform. It provides a user-friendly interface to process, store and download high-throughput sequencing data.     

Rapid evaluation and quality control of next generation sequencing data with FaQCs

Background

Next generation sequencing (NGS) technologies that parallelize the sequencing process and produce thousands to millions, or even hundreds of millions of sequences in a single sequencing run, have revolutionized genomic and genetic research. Because of the vagaries of any platform’s sequencing chemistry, the experimental processing, machine failure, and so on, the quality of sequencing reads is never perfect, and often declines as the read is extended. These errors invariably affect downstream analysis/application and should therefore be identified early on to mitigate any unforeseen effects.

Results

Here we present a novel FastQ Quality Control Software (FaQCs) that can rapidly process large volumes of data, and which improves upon previous solutions to monitor the quality and remove poor quality data from sequencing runs. Both the speed of processing and the memory footprint of storing all required information have been optimized via algorithmic and parallel processing solutions. The trimmed output compared side-by-side with the original data is part of the automated PDF output. We show how this tool can help data analysis by providing a few examples, including an increased percentage of reads recruited to references, improved single nucleotide polymorphism identification as well as de novo sequence assembly metrics.

Conclusion

FaQCs combines several features of currently available applications into a single, user-friendly process, and includes additional unique capabilities such as filtering the PhiX control sequences, conversion of FASTQ formats, and multi-threading. The original data and trimmed summaries are reported within a variety of graphics and reports, providing a simple way to do data quality control and assurance.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0366-2) contains supplementary material, which is available to authorized users.     

Ad 2.0: a novel recombineering platform for high-throughput generation of tailored adenoviruses

Recombinant adenoviruses containing a double-stranded DNA genome of 26–45 kb were broadly explored in basic virology, for vaccination purposes, for treatment of tumors based on oncolytic virotherapy, or simply as a tool for efficient gene transfer. However, the majority of recombinant adenoviral vectors (AdVs) is based on a small fraction of adenovirus types and their genetic modification. Recombineering techniques provide powerful tools for arbitrary engineering of recombinant DNA. Here, we adopted a seamless recombineering technology for high-throughput and arbitrary genetic engineering of recombinant adenoviral DNA molecules. Our cloning platform which also includes a novel recombination pipeline is based on bacterial artificial chromosomes (BACs). It enables generation of novel recombinant adenoviruses from different sources and switching between commonly used early generation AdVs and the last generation high-capacity AdVs lacking all viral coding sequences making them attractive candidates for clinical use. In combination with a novel recombination pipeline allowing cloning of AdVs containing large and complex transgenes and the possibility to generate arbitrary chimeric capsid-modified adenoviruses, these techniques allow generation of tailored AdVs with distinct features. Our technologies will pave the way toward broader applications of AdVs in molecular medicine including gene therapy and vaccination studies.