5S Ribosomal RNA Data Bank

Uniparental silencing of 5S rRNA genes in plant allopolyploids - insights from Cardamine (Brassicaceae)

While the phenomenon of uniparental silencing of 35S rDNA in interspecific hybrids and allopolyploids is well documented, there is a notable absence of information regarding whether such silencing extends to the 5S RNA component of ribosomes. To address this gap in knowledge, we analyzed the 5S and 35S rDNA expression in Cardamine (Brassicaceae) allopolyploids, namely C. × insueta (2n = 3x = 24, genome composition RRA), C. flexuosa (2n = 4x = 32, AAHH), and C. scutata (2n = 4x = 32, PPAA) which share a common diploid ancestor (AA). We employed high-throughput sequencing of transcriptomes and genomes and phylogenetic analyses of 5S rRNA variants. The genomic organization of rDNA was further scrutinized through clustering and fluorescence in situ hybridization. In the C. × insueta allotriploid, we observed uniparental dominant expression of 5S and 35S rDNA loci. In the C. flexuosa and C. scutata allotetraploids, the expression pattern differed, with the 35S rDNA being expressed from the A subgenome, whereas the 5S rDNA was expressed from the partner subgenome. Both C. flexuosa and C. scutata but not C. × insueta showed copy and locus number changes. We conclude that in stabilized allopolyploids, transcription of ribosomal RNA components occurs from different subgenomes. This phenomenon appears to result in the formation of chimeric ribosomes comprising rRNA molecules derived from distinct parental origins. We speculate that the interplay of epigenetic silencing and rDNA rearrangements introduces an additional layer of variation in multimolecule ribosomal complexes, potentially contributing to the evolutionary success of allopolyploids.

Metabolic Survival Adaptations of Plasmodium falciparum Exposed to Sublethal Doses of Fosmidomycin

The malaria parasite Plasmodium falciparum contains the apicoplast organelle that synthesizes isoprenoids, which are metabolites necessary for posttranslational modification of Plasmodium proteins. We used fosmidomycin, an antibiotic that inhibits isoprenoid biosynthesis, to identify mechanisms that underlie the development of the parasite's adaptation to the drug at sublethal concentrations. We first determined a concentration of fosmidomycin that reduced parasite growth by ∼50% over one intraerythrocytic developmental cycle (IDC). At this dose, we maintained synchronous parasite cultures for one full IDC and collected metabolomic and transcriptomic data at multiple time points to capture global and stage-specific alterations. We integrated the data with a genome-scale metabolic model of P. falciparum to characterize the metabolic adaptations of the parasite in response to fosmidomycin treatment. Our simulations showed that, in treated parasites, the synthesis of purine-based nucleotides increased, whereas the synthesis of phosphatidylcholine during the trophozoite and schizont stages decreased. Specifically, the increased polyamine synthesis led to increased nucleotide synthesis, while the reduced methyl-group cycling led to reduced phospholipid synthesis and methyltransferase activities. These results indicate that fosmidomycin-treated parasites compensate for the loss of prenylation modifications by directly altering processes that affect nucleotide synthesis and ribosomal biogenesis to control the rate of RNA translation during the IDC. This also suggests that combination therapies with antibiotics that target the compensatory response of the parasite, such as nucleotide synthesis or ribosomal biogenesis, may be more effective than treating the parasite with fosmidomycin alone.

Expansion segments in bacterial and archaeal 5S ribosomal RNAs

The large ribosomal RNAs of eukaryotes frequently contain expansion sequences that add to the size of the rRNAs but do not affect their overall structural layout and are compatible with major ribosomal function as an mRNA translation machine. The expansion of prokaryotic ribosomal RNAs is much less explored. In order to obtain more insight into the structural variability of these conserved molecules, we herein report the results of a comprehensive search for the expansion sequences in prokaryotic 5S rRNAs. Overall, 89 expanded 5S rRNAs of 15 structural types were identified in 15 archaeal and 36 bacterial genomes. Expansion segments ranging in length from 13 to 109 residues were found to be distributed among 17 insertion sites. The strains harboring the expanded 5S rRNAs belong to the bacterial orders Clostridiales, Halanaerobiales, Thermoanaerobacterales, and Alteromonadales as well as the archael order Halobacterales When several copies of a 5S rRNA gene are present in a genome, the expanded versions may coexist with normal 5S rRNA genes. The insertion sequences are typically capable of forming extended helices, which do not seemingly interfere with folding of the conserved core. The expanded 5S rRNAs have largely been overlooked in 5S rRNA databases.

Planktonic Euryarchaeota are a significant source of archaeal tetraether lipids in the ocean

Archaea are ubiquitous in marine plankton, and fossil forms of archaeal tetraether membrane lipids in sedimentary rocks document their participation in marine biogeochemical cycles for >100 million years. Ribosomal RNA surveys have identified four major clades of planktonic archaea but, to date, tetraether lipids have been characterized in only one, the Marine Group I Thaumarchaeota. The membrane lipid composition of the other planktonic archaeal groups--all uncultured Euryarchaeota--is currently unknown. Using integrated nucleic acid and lipid analyses, we found that Marine Group II Euryarchaeota (MG-II) contributed significantly to the tetraether lipid pool in the North Pacific Subtropical Gyre at shallow to intermediate depths. Our data strongly suggested that MG-II also synthesize crenarchaeol, a tetraether lipid previously considered to be a unique biomarker for Thaumarchaeota. Metagenomic datasets spanning 5 y indicated that depth stratification of planktonic archaeal groups was a stable feature in the North Pacific Subtropical Gyre. The consistent prevalence of MG-II at depths where the bulk of exported organic matter originates, together with their ubiquitous distribution over diverse oceanic provinces, suggests that this clade is a significant source of tetraether lipids to marine sediments. Our results are relevant to archaeal lipid biomarker applications in the modern oceans and the interpretation of these compounds in the geologic record.

Multilign: an algorithm to predict secondary structures conserved in multiple RNA sequences

MOTIVATION

With recent advances in sequencing, structural and functional studies of RNA lag behind the discovery of sequences. Computational analysis of RNA is increasingly important to reveal structure-function relationships with low cost and speed. The purpose of this study is to use multiple homologous sequences to infer a conserved RNA structure.

RESULTS

A new algorithm, called Multilign, is presented to find the lowest free energy RNA secondary structure common to multiple sequences. Multilign is based on Dynalign, which is a program that simultaneously aligns and folds two sequences to find the lowest free energy conserved structure. For Multilign, Dynalign is used to progressively construct a conserved structure from multiple pairwise calculations, with one sequence used in all pairwise calculations. A base pair is predicted only if it is contained in the set of low free energy structures predicted by all Dynalign calculations. In this way, Multilign improves prediction accuracy by keeping the genuine base pairs and excluding competing false base pairs. Multilign has computational complexity that scales linearly in the number of sequences. Multilign was tested on extensive datasets of sequences with known structure and its prediction accuracy is among the best of available algorithms. Multilign can run on long sequences (> 1500 nt) and an arbitrarily large number of sequences.

AVAILABILITY

The algorithm is implemented in ANSI C++ and can be downloaded as part of the RNAstructure package at: http://rna.urmc.rochester.edu.

An algorithm for searching RNA motifs in genomic sequences

RNA molecules, which are found in all living cells, fold into characteristic structures that account for their diverse functional activities. Many of these RNA structures consist of a collection of fundamental RNA motifs. The various combinations of RNA basic components form different RNA classes and define their unique structural and functional properties. The availability of many genome sequences makes it possible to search computationally for functional RNAs. Biological experiments indicate that functional RNAs have characteristic RNA structural motifs represented by specific combinations of base pairings and conserved nucleotides in the loop regions. The searching for those well-ordered RNA structures and their homologues in genomic sequences is very helpful for the understanding of RNA-based gene regulation. In this paper, we consider the following problem: given an RNA sequence with a known secondary structure, efficiently determine candidate segments in genomic sequences that can potentially form RNA secondary structures similar to the given RNA secondary structure. Our new bottom-up approach searches all potential stem-loops similar to ones of the given RNA secondary structure first, and then based on located stem-loops, detects potential homologous structural RNAs in genomic sequences.

Approximate matching of structured motifs in DNA sequences

Several methods have been developed for identifying more or less complex RNA structures in a genome. All these methods are based on the search for conserved primary and secondary sub-structures. In this paper, we present a simple formal representation of a helix, which is a combination of sequence and folding constraints, as a constrained regular expression. This representation allows us to develop a well-founded algorithm that searches for all approximate matches of a helix in a genome. The algorithm is based on an alignment graph constructed from several copies of a pushdown automaton, arranged one on top of another. This is a first attempt to take advantage of the possibilities of pushdown automata in the context of approximate matching. The worst time complexity is O(krpn), where k is the error threshold, n the size of the genome, p the size of the secondary expression, and r its number of union symbols. We then extend the algorithm to search for pseudo-knots and secondary structures containing an arbitrary number of helices.

Fast protocols for the 5S rDNA and ITS-2 based identification ofOenococcus oeni

To identify specific marker sequences for the rapid identification of Oenococcus oeni, we sequenced the 23S-5S internal transcribed spacer (ITS-2) region and the 5S rDNA of five different O. oeni strains and three phylogenetically related lactic acid bacteria (LAB). Comparative analysis revealed 100% identity among the ITS-2 region of the O. oeni strains and remarkable differences in length and sequence compared to related LAB. These results enabled us to develop a primer set for a rapid PCR-identification of O. oeni within three hours. Moreover, the comparison of the 5S rDNA sequences and the highly conserved secondary structure provided the template for the design of three fluorescence-labeled specific oligonucleotides for fluorescence in situ hybridization (FISH). These probes are partial complementary to each other. This feature promotes the accessibility to the target sequence within the ribosome and enhances the fluorescence signal. For the rapid identification of Oenococci both the 5S rRNA gene and the ITS-2 region are useful targets.

Novel perspectives in wood certification and forensics: dry wood as a source of DNA

The importance of wood for human societies can hardly be understated. If dry wood were amenable to molecular genetic investigations, this could lead to major applications in wood forensics, certification, archaeology and palaeobotany. To evaluate the potential of wood for molecular genetic investigations, we have attempted to isolate and amplify, by PCR, DNA fragments of increasing size corresponding to all three plant genomes from different regions of 10 oak logs. Stringent procedures to avoid contamination with external DNA were used in order to demonstrate the authenticity of the fragments amplified. This authenticity was further confirmed by demonstrating genetic uniformity within each log using both nuclear and chloroplast microsatellites. For most wood samples DNA was degraded, and the sequences that gave the best results were those of small size and present in high copy number (chloroplast, mitochondrial, or repeated nuclear sequences). Both storage conditions and storage duration play a role in DNA conservation. Overall, this work demonstrates that molecular markers from all three plant genomes can be used for genetic analysis on dry oak wood, but outlines some limitations and the need for further evaluation of the potential of wood for DNA analysis.

The specific hydrolysis of HIV-1 TAR RNA element with the anti-TAR hammerhead ribozyme: structural and functional implications

The main transcriptional regulator of the human immunodeficiency virus is the Tat protein, which recognises and binds to a fragment RNA at the 5' end of viral mRNA, named transactivation response element (TAR) RNA. Extensive mutagenesis studies have shown that a region of TAR RNA important for Tat binding involves a set of nucleotides surrounding a characteristic UCU nucleotide bulge. The specific Tat-TAR complex formation enhances the rate of transcription elongation but inhibition of that interaction prevents the human immunodeficiency virus type 1 (HIV-1) replication. If so, a possibility of virus inactivation would be a site specific degradation of the TAR RNA element. To break down and inactivate TAR RNA, we designated the anti-hammerhead (HH) ribozyme to cleave nucleosides within the bulge. We showed for the first time the new type of the AUC hammerhead ribozyme, which hydrolyses specifically the TAR RNA element at C8 nucleotide in the bulge (C24 in the standard TAR RNA numbering). The cleavage reaction has broad magnesium requirements. Mn and particularly Ca are less efficient. Argininamide interferes with the cleavage of TAR RNA induced by the ribozyme. These results have two implications; (i) structural, where the HIV-1 TAR RNA element in solution occurs in equilibrium of only two forms, one of which, a double stranded RNA, meets structural requirements for ribozyme pairing and cleavage, and (ii) functional, the HH ribozyme can be explored for an inactivation of HIV-1 through the TAR RNA element deintegration.

Two crystal forms of helix II of Xenopus laevis 5S rRNA with a cytosine bulge

The crystal structure of r(GCCACCCUG).r(CAGGGUCGGC), helix II of the Xenopus laevis 5S rRNA with a cytosine bulge (underlined), has been determined in two forms at 2.2 A (Form I, space group P4(2)2(1)2, a = b = 57.15 A and c = 43.54 A) and 1.7 A (Form II, space group P4(3)2(1)2, a = b = 32.78 A and c = 102.5 A). The helical regions of the nonamers are found in the standard A-RNA conformations and the two forms have an RMS deviation of 0.75 A. However, the cytosine bulge adopts two significantly different conformations with an RMS deviation of 3.9 A. In Form I, the cytosine bulge forms an intermolecular C+*G.C triple in the major groove of a symmetry-related duplex with intermolecular hydrogen bonds between N4C and O6G, and between protonated N3+C and N7G. In contrast, a minor groove C*G.C triple is formed in Form II with intermolecular hydrogen bonds between O2C and N2G, and between N3C and N3G with a water bridge. A partial major groove opening was observed in Form I structure at the bulge site. Two Ca2+ ions were found in Form I helix whereas there were none in Form II. The structural comparison of these two forms indicates that bulged residues can adopt a variety of conformations with little perturbation to the global helix structure. This suggests that bulged residues could function as flexible latches in bridging double helical motifs and facilitate the folding of large RNA molecules.

Re-annotating the Mycoplasma pneumoniae genome sequence: adding value, function and reading frames

Four years after the original sequence submission, we have re-annotated the genome of Mycoplasma pneumoniae to incorporate novel data. The total number of ORFss has been increased from 677 to 688 (10 new proteins were predicted in intergenic regions, two further were newly identified by mass spectrometry and one protein ORF was dismissed) and the number of RNAs from 39 to 42 genes. For 19 of the now 35 tRNAs and for six other functional RNAs the exact genome positions were re-annotated and two new tRNA(Leu) and a small 200 nt RNA were identified. Sixteen protein reading frames were extended and eight shortened. For each ORF a consistent annotation vocabulary has been introduced. Annotation reasoning, annotation categories and comparisons to other published data on M.pneumoniae functional assignments are given. Experimental evidence includes 2-dimensional gel electrophoresis in combination with mass spectrometry as well as gene expression data from this study. Compared to the original annotation, we increased the number of proteins with predicted functional features from 349 to 458. The increase includes 36 new predictions and 73 protein assignments confirmed by the published literature. Furthermore, there are 23 reductions and 30 additions with respect to the previous annotation. mRNA expression data support transcription of 184 of the functionally unassigned reading frames.

Conformational dynamics of a 5S rRNA hairpin domain containing loop D and a single nucleotide bulge

Molecular modeling and molecular dynamics have been employed to study the conformation and flexibility of a 15-nucleotide fragment of the plant 5S rRNA containing loop D and a single uridine bulge. Two different model built initial structures were used: one with the bulge localized inside the helical stem and another with the bulge pointing out from the helix. Several independent 700-ps-long trajectories in aqueous solution with Na(+) conterions were produced for each starting structure. The bulge nucleotide inside the helix stayed in two main conformations, both of which affected the geometry of the stem part opposite the bulge. When the bulge nucleotide was located outside the helix, we found high base mobility and local backbone flexibility. The dynamics of the hydrogen bond network and conformational changes from a direct to a water mediated hydrogen bond in the sheared G-A basepair in the tetraloop was described. Our results correlate with lead ion induced cleavage patterns in 5S rRNA. Sites resistant to nonspecific lead cleavage appeared in all our simulations as the most rigid fragments independent of the localization of the bulge nucleotide.

Organization, structure, and variability of the rRNA operon of the Whipple's disease bacterium (Tropheryma whippelii)

Whipple's disease is a systemic disorder associated with a cultivation-resistant, poorly characterized actinomycete, Tropheryma whippelii. We determined a nearly complete rRNA operon sequence of T. whippelii from specimens from 3 patients with Whipple's disease, as well as partial operon sequences from 43 patients. Variability was observed in the 16S-23S rRNA spacer sequences, leading to the description of five distinct sequence types. One specimen contained two spacer sequence types, raising the possibility of a double infection. Secondary structure models for the primary rRNA transcript and mature rRNAs revealed rare or unique features.

5S ribosomal RNA database Y2K

This paper presents the updated version (Y2K) of the database of ribosomal 5S ribonucleic acids (5S rRNA) and their genes (5S rDNA), http://rose.man/poznan.pl/5SData/index.html. This edition of the database contains 1985primary structures of 5S rRNA and 5S rDNA. They include 60 archaebacterial, 470 eubacterial, 63 plastid, nine mitochondrial and 1383 eukaryotic sequences. The nucleotide sequences of the 5S rRNAs or 5S rDNAs are divided according to the taxonomic position of the source organisms.