Conservative Fragments in Bacterial 16S rRNA Genes and Primer Design for 16S Ribosomal DNA Amplicons in Metagenomic Studies

Bacterial 16S ribosomal DNA (rDNA) amplicons have been widely used in the classification of uncultured bacteria inhabiting environmental niches. Primers targeting conservative regions of the rDNAs are used to generate amplicons of variant regions that are informative in taxonomic assignment. One problem is that the percentage coverage and application scope of the primers used in previous studies are largely unknown. In this study, conservative fragments of available rDNA sequences were first mined and then used to search for candidate primers within the fragments by measuring the coverage rate defined as the percentage of bacterial sequences containing the target. Thirty predicted primers with a high coverage rate (>90%) were identified, which were basically located in the same conservative regions as known primers in previous reports, whereas 30% of the known primers were associated with a coverage rate of <90%. The application scope of the primers was also examined by calculating the percentages of failed detections in bacterial phyla. Primers A519–539, E969–983, E1063–1081, U515 and E517, are highly recommended because of their high coverage in almost all phyla. As expected, the three predominant phyla, Firmicutes, Gemmatimonadetes and Proteobacteria, are best covered by the predicted primers. The primers recommended in this report shall facilitate a comprehensive and reliable survey of bacterial diversity in metagenomic studies

A Modern Mode of Activation for Nucleic Acid Enzymes

Through evolution, enzymes have developed subtle modes of activation in order to ensure the sufficiently high substrate specificity required by modern cellular metabolism. One of these modes is the use of a target-dependent module (i.e. a docking domain) such as those found in signalling kinases. Upon the binding of the target to a docking domain, the substrate is positioned within the catalytic site. The prodomain acts as a target-dependent module switching the kinase from an off state to an on state. As compared to the allosteric mode of activation, there is no need for the presence of a third partner. None of the ribozymes discovered to date have such a mode of activation, nor does any other known RNA. Starting from a specific on/off adaptor for the hepatitis delta virus ribozyme, that differs but has a mechanism reminiscent of this signalling kinase, we have adapted this mode of activation, using the techniques of molecular engineering, to both catalytic RNAs and DNAs exhibiting various activities. Specifically, we adapted three cleaving ribozymes (hepatitis delta virus, hammerhead and hairpin ribozymes), a cleaving 10-23 deoxyribozyme, a ligating hairpin ribozyme and an artificially selected capping ribozyme. In each case, there was a significant gain in terms of substrate specificity. Even if this mode of control is unreported for natural catalytic nucleic acids, its use needs not be limited to proteinous enzymes. We suggest that the complexity of the modern cellular metabolism might have been an important selective pressure in this evolutionary process

Structural Constraints Identified with Covariation Analysis in Ribosomal RNA

Covariation analysis is used to identify those positions with similar patterns of sequence variation in an alignment of RNA sequences. These constraints on the evolution of two positions are usually associated with a base pair in a helix. While mutual information (MI) has been used to accurately predict an RNA secondary structure and a few of its tertiary interactions, early studies revealed that phylogenetic event counting methods are more sensitive and provide extra confidence in the prediction of base pairs. We developed a novel and powerful phylogenetic events counting method (PEC) for quantifying positional covariation with the Gutell lab’s new RNA Comparative Analysis Database (rCAD). The PEC and MI-based methods each identify unique base pairs, and jointly identify many other base pairs. In total, both methods in combination with an N-best and helix-extension strategy identify the maximal number of base pairs. While covariation methods have effectively and accurately predicted RNAs secondary structure, only a few tertiary structure base pairs have been identified. Analysis presented herein and at the Gutell lab’s Comparative RNA Web (CRW) Site reveal that the majority of these latter base pairs do not covary with one another. However, covariation analysis does reveal a weaker although significant covariation between sets of nucleotides that are in proximity in the three-dimensional RNA structure. This reveals that covariation analysis identifies other types of structural constraints beyond the two nucleotides that form a base pair

Novel high-rank phylogenetic lineages within a sulfur spring (Zodletone Spring, Oklahoma), revealed using a combined pyrosequencing-Sanger approach

The utilization of high-throughput sequencing technologies in 16S rRNA gene-based diversity surveys has indicated that within most ecosystems, a significant fraction of the community could not be assigned to known microbial phyla. Accurate determination of the phylogenetic affiliation of such sequences is difficult due to the short-read-length output of currently available high-throughput technologies. This fraction could harbor multiple novel phylogenetic lineages that have so far escaped detection. Here we describe our efforts in accurate assessment of the novelty and phylogenetic affiliation of selected unclassified lineages within a pyrosequencing data set generated from source sediments of Zodletone Spring, a sulfide- and sulfur-rich spring in southwestern Oklahoma. Lineage-specific forward primers were designed for 78 putatively novel lineages identified within the pyrosequencing data set, and representative nearly full-length small-subunit (SSU) rRNA gene sequences were obtained by pairing those primers with reverse universal bacterial primers. Of the 78 lineages tested, amplifiable products were obtained for 52, 32 of which had at least one nearly full-length sequence that was representative of the lineage targeted. Analysis of phylogenetic affiliation of the obtained Sanger sequences identified 5 novel candidate phyla and 10 novel candidate classes (within Fibrobacteres, Planctomycetes, and candidate phyla BRC1, GN12, TM6, TM7, LD1, WS2, and GN06) in the data set, in addition to multiple novel orders and families. The discovery of multiple novel phyla within a pilot study of a single ecosystem clearly shows the potential of the approach in identifying novel diversities within the rare biosphere.Peer reviewedMicrobiology and Molecular Genetic

Mother and Adolescent Reports of Associations Between Child Behavior Problems and Mother-Child Relationship Qualities: Separating Shared Variance from Individual Variance

This study contrasts results from different correlational methods for examining links between mother and child (N = 72 dyads) reports of early adolescent (M = 11.5 years) behavior problems and relationship negativity and support. Simple (Pearson) correlations revealed a consistent pattern of statistically significant associations, regardless of whether scores came from the same reporter or from different reporters. When correlations between behavior problems and relationship quality differed, within-reporter correlations were always greater in magnitude than between-reporter correlations. Dyadic (common fate) analyses designed for interdependent data decomposed within-reporter correlations into variance shared across reporters (dyadic correlations) and variance unique to specific reporters (individual correlations). Dyadic correlations were responsible for most associations between adolescent behavior problems and relationship negativity; after partitioning variance shared across reporters, no individual correlations emerged as statistically significant. In contrast, adolescent behavior problems were linked to relationship support via both shared variance and variance unique to maternal perceptions. Dyadic analyses provide a parsimonious alternative to multiple contrasts in instances when identical measures have been collected from multiple reporters. Findings from these analyses indicate that same-reporter variance bias should not be assumed in the absence of dyadic statistical analyses

Insights into substrate stabilization from snapshots of the peptidyl transferase center of the intact 70S ribosome

Protein synthesis is catalyzed in the peptidyl transferase center (PTC), located in the large (50S) subunit of the ribosome. No high-resolution structure of the intact ribosome has contained a complete active site including both A- and P-site tRNAs. In addition, although past structures of the 50S subunit have found no ordered proteins at the PTC, biochemical evidence suggests that specific proteins are capable of interacting with the 3′ ends of tRNA ligands. Here we present structures, at 3.6-Å and 3.5-Å resolution respectively, of the 70S ribosome in complex with A- and P-site tRNAs that mimic pre- and post-peptidyl-transfer states. These structures demonstrate that the PTC is very similar between the 50S subunit and the intact ribosome. They also reveal interactions between the ribosomal proteins L16 and L27 and the tRNA substrates, helping to elucidate the role of these proteins in peptidyl transfer

Loss of Ribosomal Protein L11 Affects Zebrafish Embryonic Development through a p53-Dependent Apoptotic Response

Ribosome is responsible for protein synthesis in all organisms and ribosomal proteins (RPs) play important roles in the formation of a functional ribosome. L11 was recently shown to regulate p53 activity through a direct binding with MDM2 and abrogating the MDM2-induced p53 degradation in response to ribosomal stress. However, the studies were performed in cell lines and the significance of this tumor suppressor function of L11 has yet to be explored in animal models. To investigate the effects of the deletion of L11 and its physiological relevance to p53 activity, we knocked down the rpl11 gene in zebrafish and analyzed the p53 response. Contrary to the cell line-based results, our data indicate that an L11 deficiency in a model organism activates the p53 pathway. The L11-deficient embryos (morphants) displayed developmental abnormalities primarily in the brain, leading to embryonic lethality within 6–7 days post fertilization. Extensive apoptosis was observed in the head region of the morphants, thus correlating the morphological defects with apparent cell death. A decrease in total abundance of genes involved in neural patterning of the brain was observed in the morphants, suggesting a reduction in neural progenitor cells. Upregulation of the genes involved in the p53 pathway were observed in the morphants. Simultaneous knockdown of the p53 gene rescued the developmental defects and apoptosis in the morphants. These results suggest that ribosomal dysfunction due to the loss of L11 activates a p53-dependent checkpoint response to prevent improper embryonic development

Arginine Cofactors on the Polymerase Ribozyme

The RNA world hypothesis states that the early evolution of life went through a stage in which RNA served both as genome and as catalyst. The central catalyst in an RNA world organism would have been a ribozyme that catalyzed RNA polymerization to facilitate self-replication. An RNA polymerase ribozyme was developed previously in the lab but it is not efficient enough for self-replication. The factor that limits its polymerization efficiency is its weak sequence-independent binding of the primer/template substrate. Here we tested whether RNA polymerization could be improved by a cationic arginine cofactor, to improve the interaction with the substrate. In an RNA world, amino acid-nucleic acid conjugates could have facilitated the emergence of the translation apparatus and the transition to an RNP world. We chose the amino acid arginine for our study because this is the amino acid most adept to interact with RNA. An arginine cofactor was positioned at ten different sites on the ribozyme, using conjugates of arginine with short DNA or RNA oligonucleotides. However, polymerization efficiency was not increased in any of the ten positions. In five of the ten positions the arginine reduced or modulated polymerization efficiency, which gives insight into the substrate-binding site on the ribozyme. These results suggest that the existing polymerase ribozyme is not well suited to using an arginine cofactor

Disease-Associated Mutations That Alter the RNA Structural Ensemble

Genome-wide association studies (GWAS) often identify disease-associated mutations in intergenic and non-coding regions of the genome. Given the high percentage of the human genome that is transcribed, we postulate that for some observed associations the disease phenotype is caused by a structural rearrangement in a regulatory region of the RNA transcript. To identify such mutations, we have performed a genome-wide analysis of all known disease-associated Single Nucleotide Polymorphisms (SNPs) from the Human Gene Mutation Database (HGMD) that map to the untranslated regions (UTRs) of a gene. Rather than using minimum free energy approaches (e.g. mFold), we use a partition function calculation that takes into consideration the ensemble of possible RNA conformations for a given sequence. We identified in the human genome disease-associated SNPs that significantly alter the global conformation of the UTR to which they map. For six disease-states (Hyperferritinemia Cataract Syndrome, β-Thalassemia, Cartilage-Hair Hypoplasia, Retinoblastoma, Chronic Obstructive Pulmonary Disease (COPD), and Hypertension), we identified multiple SNPs in UTRs that alter the mRNA structural ensemble of the associated genes. Using a Boltzmann sampling procedure for sub-optimal RNA structures, we are able to characterize and visualize the nature of the conformational changes induced by the disease-associated mutations in the structural ensemble. We observe in several cases (specifically the 5′ UTRs of FTL and RB1) SNP–induced conformational changes analogous to those observed in bacterial regulatory Riboswitches when specific ligands bind. We propose that the UTR and SNP combinations we identify constitute a “RiboSNitch,” that is a regulatory RNA in which a specific SNP has a structural consequence that results in a disease phenotype. Our SNPfold algorithm can help identify RiboSNitches by leveraging GWAS data and an analysis of the mRNA structural ensemble

Global Regulatory Functions of the Staphylococcus aureus Endoribonuclease III in Gene Expression

RNA turnover plays an important role in both virulence and adaptation to stress in the Gram-positive human pathogen Staphylococcus aureus. However, the molecular players and mechanisms involved in these processes are poorly understood. Here, we explored the functions of S. aureus endoribonuclease III (RNase III), a member of the ubiquitous family of double-strand-specific endoribonucleases. To define genomic transcripts that are bound and processed by RNase III, we performed deep sequencing on cDNA libraries generated from RNAs that were co-immunoprecipitated with wild-type RNase III or two different cleavage-defective mutant variants in vivo. Several newly identified RNase III targets were validated by independent experimental methods. We identified various classes of structured RNAs as RNase III substrates and demonstrated that this enzyme is involved in the maturation of rRNAs and tRNAs, regulates the turnover of mRNAs and non-coding RNAs, and autoregulates its synthesis by cleaving within the coding region of its own mRNA. Moreover, we identified a positive effect of RNase III on protein synthesis based on novel mechanisms. RNase III–mediated cleavage in the 5′ untranslated region (5′UTR) enhanced the stability and translation of cspA mRNA, which encodes the major cold-shock protein. Furthermore, RNase III cleaved overlapping 5′UTRs of divergently transcribed genes to generate leaderless mRNAs, which constitutes a novel way to co-regulate neighboring genes. In agreement with recent findings, low abundance antisense RNAs covering 44% of the annotated genes were captured by co-immunoprecipitation with RNase III mutant proteins. Thus, in addition to gene regulation, RNase III is associated with RNA quality control of pervasive transcription. Overall, this study illustrates the complexity of post-transcriptional regulation mediated by RNase III