Browse
Submit Data
Databases
API
Help

Dataset Information

34 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Model-Free RNA Sequence and Structure Alignment Informed by SHAPE Probing Reveals a Conserved Alternate Secondary Structure for 16S rRNA.

ABSTRACT: Discovery and characterization of functional RNA structures remains challenging due to deficiencies in de novo secondary structure modeling. Here we describe a dynamic programming approach for model-free sequence comparison that incorporates high-throughput chemical probing data. Based on SHAPE probing data alone, ribosomal RNAs (rRNAs) from three diverse organisms--the eubacteria E. coli and C. difficile and the archeon H. volcanii--could be aligned with accuracies comparable to alignments based on actual sequence identity. When both base sequence identity and chemical probing reactivities were considered together, accuracies improved further. Derived sequence alignments and chemical probing data from protein-free RNAs were then used as pseudo-free energy constraints to model consensus secondary structures for the 16S and 23S rRNAs. There are critical differences between these experimentally-informed models and currently accepted models, including in the functionally important neck and decoding regions of the 16S rRNA. We infer that the 16S rRNA has evolved to undergo large-scale changes in base pairing as part of ribosome function. As high-quality RNA probing data become widely available, structurally-informed sequence alignment will become broadly useful for de novo motif and function discovery.

SUBMITTER: Lavender CA

PROVIDER: S-EPMC4438973 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

How conserved are the conserved 16S-rRNA regions?

Project description:The 16S rRNA gene has been used as master key for studying prokaryotic diversity in almost every environment. Despite the claim of several researchers to have the best universal primers, the reality is that no primer has been demonstrated to be truly universal. This suggests that conserved regions of the gene may not be as conserved as expected. The aim of this study was to evaluate the conservation degree of the so-called conserved regions flanking the hypervariable regions of the 16S rRNA gene. Data contained in SILVA database (release 123) were used for the study. Primers reported as matches of each conserved region were assembled to form contigs; sequences sizing 12 nucleotides (12-mers) were extracted from these contigs and searched into the entire set of SILVA sequences. Frequency analysis shown that extreme regions, 1 and 10, registered the lowest frequencies. 12-mer frequencies revealed segments of contigs that were not as conserved as expected (≤90%). Fragments corresponding to the primer contigs 3, 4, 5b and 6a were recovered from all sequences in SILVA database. Nucleotide frequency analysis in each consensus demonstrated that only a small fraction of these so-called conserved regions is truly conserved in non-redundant sequences. It could be concluded that conserved regions of the 16S rRNA gene exhibit considerable variation that has to be considered when using this gene as biomarker.

| S-EPMC5333541 | biostudies-literature

TurboFold II: RNA structural alignment and secondary structure prediction informed by multiple homologs.

Project description:This paper presents TurboFold II, an extension of the TurboFold algorithm for predicting secondary structures for multiple RNA homologs. TurboFold II augments the structure prediction capabilities of TurboFold by additionally providing multiple sequence alignments. Probabilities for alignment of nucleotide positions between all pairs of input sequences are iteratively estimated in TurboFold II by incorporating information from both the sequence identity and secondary structures. A multiple sequence alignment is obtained from these probabilities by using a probabilistic consistency transformation and a hierarchically computed guide tree. To assess TurboFold II, its sequence alignment and structure predictions were compared with leading tools, including methods that focus on alignment alone and methods that provide both alignment and structure prediction. TurboFold II has comparable alignment accuracy with MAFFT and higher accuracy than other tools. TurboFold II also has comparable structure prediction accuracy as the original TurboFold algorithm, which is one of the most accurate methods. TurboFold II is part of the RNAstructure software package, which is freely available for download at http://rna.urmc.rochester.edu under a GPL license.

| S-EPMC5714223 | biostudies-literature

NAST: a multiple sequence alignment server for comparative analysis of 16S rRNA genes.

Project description:Microbiologists conducting surveys of bacterial and archaeal diversity often require comparative alignments of thousands of 16S rRNA genes collected from a sample. The computational resources and bioinformatics expertise required to construct such an alignment has inhibited high-throughput analysis. It was hypothesized that an online tool could be developed to efficiently align thousands of 16S rRNA genes via the NAST (Nearest Alignment Space Termination) algorithm for creating multiple sequence alignments (MSA). The tool was implemented with a web-interface at http://greengenes.lbl.gov/NAST. Each user-submitted sequence is compared with Greengenes' 'Core Set', comprising approximately 10,000 aligned non-chimeric sequences representative of the currently recognized diversity among bacteria and archaea. User sequences are oriented and paired with their closest match in the Core Set to serve as a template for inserting gap characters. Non-16S data (sequence from vector or surrounding genomic regions) are conveniently removed in the returned alignment. From the resulting MSA, distance matrices can be calculated for diversity estimates and organisms can be classified by taxonomy. The ability to align and categorize large sequence sets using a simple interface has enabled researchers with various experience levels to obtain bacterial and archaeal community profiles.

| S-EPMC1538769 | biostudies-literature

Identification and Quantification of Abundant Species from Pyrosequences of 16S rRNA by Consensus Alignment.

Project description:16S rRNA gene profiling has recently been boosted by the development of pyrosequencing methods. A common analysis is to group pyrosequences into Operational Taxonomic Units (OTUs), such that reads in an OTU are likely sampled from the same species. However, species diversity estimated from error-prone 16S rRNA pyrosequences may be inflated because the reads sampled from the same 16S rRNA gene may appear different, and current OTU inference approaches typically involve time-consuming pairwise/multiple distance calculation and clustering. I propose a novel approach AbundantOTU based on a Consensus Alignment (CA) algorithm, which infers consensus sequences, each representing an OTU, taking advantage of the sequence redundancy for abundant species. Pyrosequencing reads can then be recruited to the consensus sequences to give quantitative information for the corresponding species. As tested on 16S rRNA pyrosequence datasets from mock communities with known species, AbundantOTU rapidly reported identified sequences of the source 16S rRNAs and the abundances of the corresponding species. AbundantOTU was also applied to 16S rRNA pyrosequence datasets derived from real microbial communities and the results are in general agreement with previous studies.

| S-EPMC3217275 | biostudies-literature

Solution probing of metal ion binding by helix 27 from Escherichia coli 16S rRNA.

Project description:Helix (H)27 from Escherichia coli 16S ribosomal (r)RNA is centrally located within the small (30S) ribosomal subunit, immediately adjacent to the decoding center. Bacterial 30S subunit crystal structures depicting Mg(2+) binding sites resolve two magnesium ions within the vicinity of H27: one in the major groove of the G886-U911 wobble pair, and one within the GCAA tetraloop. Binding of such metal cations is generally thought to be crucial for RNA folding and function. To ask how metal ion-RNA interactions in crystals compare with those in solution, we have characterized, using solution NMR spectroscopy, Tb(3+) footprinting and time-resolved fluorescence resonance energy transfer (tr-FRET), location, and modes of metal ion binding in an isolated H27. NMR and Tb(3+) footprinting data indicate that solution secondary structure and Mg(2+) binding are generally consistent with the ribosomal crystal structures. However, our analyses also suggest that H27 is dynamic in solution and that metal ions localize within the narrow major groove formed by the juxtaposition of the loop E motif with the tandem G894-U905 and G895-U904 wobble pairs. In addition, tr-FRET studies provide evidence that Mg(2+) uptake by the H27 construct results in a global lengthening of the helix. We propose that only a subset of H27-metal ion interactions has been captured in the crystal structures of the 30S ribosomal subunit, and that small-scale structural dynamics afforded by solution conditions may contribute to these differences. Our studies thus highlight an example for differences between RNA-metal ion interactions observed in solution and in crystals.

| S-EPMC1370855 | biostudies-literature

SHAPE and DMS probing of SLNCR1 lncRNA conserved region 403-780

Project description:The goal of this study was to understand the underlying structure of nucleotides 403-780 of the lncRNA SLNCR1, in-cell and when extracted from nuclear and cytoplasmic fractions. SHAPE and DMS probing revealed that the region is largely unstructured inside and outside of the cell, and appears protein-bound in primary melanoma cells.

2019-11-07 | GSE140048 | GEO

ZEAL: protein structure alignment based on shape similarity.

Project description:MotivationMost protein-structure superimposition tools consider only Cartesian coordinates. Yet, much of biology happens on the surface of proteins, which is why proteins with shared ancestry and similar function often have comparable surface shapes. Superposition of proteins based on surface shape can enable comparison of highly divergent proteins, identify convergent evolution and enable detailed comparison of surface features and binding sites.ResultsWe present ZEAL, an interactive tool to superpose global and local protein structures based on their shape resemblance using 3D (Zernike-Canterakis) functions to represent the molecular surface. In a benchmark study of structures with the same fold, we show that ZEAL outperforms two other methods for shape-based superposition. In addition, alignments from ZEAL were of comparable quality to the coordinate-based superpositions provided by TM-align. For comparisons of proteins with limited sequence and backbone-fold similarity, where coordinate-based methods typically fail, ZEAL can often find alignments with substantial surface-shape correspondence. In combination with shape-based matching, ZEAL can be used as a general tool to study relationships between shape and protein function. We identify several categories of protein functions where global shape similarity is significantly more likely than expected by random chance, when comparing proteins with little similarity on the fold level. In particular, we find that global surface shape similarity is particular common among DNA binding proteins.Availability and implementationZEAL can be used online at https://andrelab.org/zeal or as a standalone program with command line or graphical user interface. Source files and installers are available at https://github.com/Andre-lab/ZEAL.Supplementary informationSupplementary data are available at Bioinformatics online.

| S-EPMC10262298 | biostudies-literature

Solution structure of an alternate conformation of helix27 from Escherichia coli16S rRNA.

Project description:Helix (H)27 of 16S ribosomal (r)RNA from Escherichia coli was dubbed the "switch helix" when mutagenesis suggested that two alternative base pair registers may have distinct functional roles in the bacterial ribosome. Although more recent genetic analyses suggest that H27 conformational switching is not required for translation, previous solution studies demonstrated that the isolated E. coli H27 can dynamically convert between the 885 and 888 conformations. Here, we have solved the nuclear magnetic resonance solution structure of a locked 888 conformation. NOE and residual dipolar coupling restraints reveal an architecture that markedly differs from that of the 885 conformation found in crystal structures of the bacterial ribosome. In place of the loop E motif that characterizes the 885 conformer and that the 888 conformer cannot adopt, we find evidence for an asymmetrical A-rich internal loop stabilized by stacking interactions among the unpaired A's. Comparison of the isolated H27 888 solution structure with the 885 crystal structure within the context of the ribosome suggests a difference in overall length of H27 that presents one plausible reason for the absence of H27 conformational switching within the sterically confining ribosome.

| S-EPMC3145048 | biostudies-literature

Multiple Streptomyces species with distinct secondary metabolomes have identical 16S rRNA gene sequences.

Project description:Microbial diversity studies using small subunit (SSU) rRNA gene sequences continue to advance our understanding of biological and ecological systems. Although a good predictor of overall diversity, using this gene to infer the presence of a species in a sample is more controversial. Here, we present a detailed polyphasic analysis of 10 bacterial strains isolated from three coastal lichens Lichina confinis, Lichina pygmaea and Roccella fuciformis with SSU rRNA gene sequences identical to the type strain of Streptomyces cyaneofuscatus. This analysis included phenotypic, microscopic, genetic and genomic comparisons and showed that despite their identical SSU rRNA sequences the strains had markedly different properties, and could be distinguished as 5 different species. Significantly, secondary metabolites profiles from these strains were also found to be different. It is thus clear that SSU rRNA based operational taxonomy units, even at the most stringent cut-off can represent multiple bacterial species, and that at least for the case of Streptomyces, strain de-replication based on SSU gene sequences prior to screening for bioactive molecules can miss potentially interesting novel molecules produced by this group that is notorious for the production of drug-leads.

| S-EPMC5593946 | biostudies-literature

Phylogenic inference using alignment-free methods for applications in microbial community surveys using 16s rRNA gene.

Project description:The diversity of microbiota is best explored by understanding the phylogenetic structure of the microbial communities. Traditionally, sequence alignment has been used for phylogenetic inference. However, alignment-based approaches come with significant challenges and limitations when massive amounts of data are analyzed. In the recent decade, alignment-free approaches have enabled genome-scale phylogenetic inference. Here we evaluate three alignment-free methods: ACS, CVTree, and Kr for phylogenetic inference with 16s rRNA gene data. We use a taxonomic gold standard to compare the accuracy of alignment-free phylogenetic inference with that of common microbiome-wide phylogenetic inference pipelines based on PyNAST and MUSCLE alignments with FastTree and RAxML. We re-simulate fecal communities from Human Microbiome Project data to evaluate the performance of the methods on datasets with properties of real data. Our comparisons show that alignment-free methods are not inferior to alignment-based methods in giving accurate and robust phylogenic trees. Moreover, consensus ensembles of alignment-free phylogenies are superior to those built from alignment-based methods in their ability to highlight community differences in low power settings. In addition, the overall running times of alignment-based and alignment-free phylogenetic inference are comparable. Taken together our empirical results suggest that alignment-free methods provide a viable approach for microbiome-wide phylogenetic inference.

| S-EPMC5685621 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data