The comparative RNA web (CRW) site: an online database of comparative sequence and structure information for ribosomal, intron, and other RNAs

Graphical representation of ribosomal RNA probe accessibility data using ARB software package

Background

Taxon specific hybridization probes in combination with a variety of commonly used hybridization formats nowadays are standard tools in microbial identification. A frequently applied technology, fluorescence in situ hybridization (FISH), besides single cell identification, allows the localization and functional studies of the microbial community composition. Careful in silico design and evaluation of potential oligonucleotide probe targets is therefore crucial for performing successful hybridization experiments.

Results

The PROBE Design tools of the ARB software package take into consideration several criteria such as number, position and quality of diagnostic sequence differences while designing oligonucleotide probes. Additionally, new visualization tools were developed to enable the user to easily examine further sequence associated criteria such as higher order structure, conservation, G+C content, transition-transversion profiles and in situ target accessibility patterns. The different types of sequence associated information (SAI) can be visualized by user defined background colors within the ARB primary and secondary structure editors as well as in the PROBE Match tool.

Conclusion

Using this tool, in silico probe design and evaluation can be performed with respect to in situ probe accessibility data. The evaluation of proposed probe targets with respect to higher-order rRNA structure is of importance for successful design and performance of in situ hybridization experiments. The entire ARB software package along with the probe accessibility data is available from the ARB home page .

Modular RNA architecture revealed by computational analysis of existing pseudoknots and ribosomal RNAs

Modular architecture is a hallmark of RNA structures, implying structural, and possibly functional, similarity among existing RNAs. To systematically delineate the existence of smaller topologies within larger structures, we develop and apply an efficient RNA secondary structure comparison algorithm using a newly developed two-dimensional RNA graphical representation. Our survey of similarity among 14 pseudoknots and subtopologies within ribosomal RNAs (rRNAs) uncovers eight pairs of structurally related pseudoknots with non-random sequence matches and reveals modular units in rRNAs. Significantly, three structurally related pseudoknot pairs have functional similarities not previously known: one pair involves the 3′ end of brome mosaic virus genomic RNA (PKB134) and the alternative hammerhead ribozyme pseudoknot (PKB173), both of which are replicase templates for viral RNA replication; the second pair involves structural elements for translation initiation and ribosome recruitment found in the viral internal ribosome entry site (PKB223) and the V4 domain of 18S rRNA (PKB205); the third pair involves 18S rRNA (PKB205) and viral tRNA-like pseudoknot (PKB134), which probably recruits ribosomes via structural mimicry and base complementarity. Additionally, we quantify the modularity of 16S and 23S rRNAs by showing that RNA motifs can be constructed from at least 210 building blocks. Interestingly, we find that the 5S rRNA and two tree modules within 16S and 23S rRNAs have similar topologies and tertiary shapes. These modules can be applied to design novel RNA motifs via build-up-like procedures for constructing sequences and folds.

Mechanistic approach to the problem of hybridization efficiency in fluorescent in situ hybridization

In fluorescent in situ hybridization (FISH), the efficiency of hybridization between the DNA probe and the rRNA has been related to the accessibility of the rRNA when ribosome content and cell permeability are not limiting. Published rRNA accessibility maps show that probe brightness is sensitive to the organism being hybridized and the exact location of the target site and, hence, it is highly unpredictable based on accessibility only. In this study, a model of FISH based on the thermodynamics of nucleic acid hybridization was developed. The model provides a mechanistic approach to calculate the affinity of the probe to the target site, which is defined as the overall Gibbs free energy change (DeltaG degrees overall) for a reaction scheme involving the DNA-rRNA and intramolecular DNA and rRNA interactions that take place during FISH. Probe data sets for the published accessibility maps and experiments targeting localized regions in the 16S rRNA of Escherichia coli were used to demonstrate that DeltaG degrees overall is a strong predictor of hybridization efficiency and superior to conventional estimates based on the dissociation temperature of the DNA/rRNA duplex. The use of the proposed model also allowed the development of mechanistic approaches to increase probe brightness, even in seemingly inaccessible regions of the 16S rRNA. Finally, a threshold DeltaG degrees overall of -13.0 kcal/mol was proposed as a goal in the design of FISH probes to maximize hybridization efficiency without compromising specificity.

Phylogeny and evolution of medical species of Candida and related taxa: a multigenic analysis

Hemiascomycetes are species of yeasts within the order Saccharomycetales. The order encompasses disparate genera with a variety of life styles, including opportunistic human pathogens (e.g., Candida albicans), plant pathogens (e.g., Eremothecium gossypii), and cosmopolitan yeasts associated with water and decaying vegetation. To analyze the phylogeny of medically important species of yeasts, we selected 38 human pathogenic and related strains in the order Saccharomycetales. The DNA sequences of six nuclear genes were analyzed by maximum likelihood and Bayesian phylogenetic methods. The maximum likelihood analysis of the combined data for all six genes resolved three major lineages with significant support according to Bayesian posterior probability. One clade was mostly comprised of pathogenic species of Candida. Another major group contained members of the family Metschnikowiaceae as a monophyletic group, three species of Debaryomyces, and strains of Candida guilliermondii. The third clade consisted exclusively of species of the family Saccharomycetaceae. Analysis of the evolution of key characters indicated that both codon reassignment and coenzyme Q(9) likely had single origins with multiple losses. Tests of correlated character evolution revealed that these two traits evolved independently.

Phylogeny of the New World true frogs (Rana)

Phylogenetic relationships among the species of true frogs (Rana) from North, South, and Central America were investigated based on the sequences of approximately 2 kb from the mitochondrial genome, sampled from most of the described species, as well as eight undescribed species. This analysis, combined with previous studies of the phylogeny of New World Rana, served as the basis for a revised classification of the group. The American species of Rana are not monophyletic; the western North American Amerana is more closely related to the R. temporaria group of Eurasia (together, these frogs form the group Laurasiarana). The remaining species from the Americas form the monophyletic group Novirana, which includes: R. sylvatica; Aquarana (the R. catesbeiana group); Ranula (the R. palmipes group, including the mostly upland Levirana species and the mostly lowland Lithobates species); Torrentirana (the R. tarahumarae group, or Zweifelia, plus R. sierramadrensis), Stertirana (the R. montezumae group, or Lacusirana, plus R. pipiens), Nenirana (the R. areolata group), and Scurrilirana (most of the southern and tropical leopard frogs). The mitochondrial sequences supported many of the previous hypotheses of relationships of New World Rana, although there were some differences involving the placement of the species R. pipiens, R. sierramadrensis, and R. sylvatica. Parametric bootstrap analyses indicated significant support for the relationships inferred from the mtDNA sequences, and rejected the previous hypotheses of relationships for these three species.

Dynamics inherent in helix 27 from Escherichia coli 16S ribosomal RNA

The original interpretation of a series of genetic studies suggested that the highly conserved Escherichia coli 16S ribosomal RNA helix 27 (H27) adopts two alternative secondary structure motifs, the 885 and 888 conformations, during each cycle of amino acid incorporation. Recent crystallographic and genetic evidence has called this hypothesis into question. To ask whether a slippery sequence such as that of H27 may harbor inherent conformational dynamics, we have designed a series of model RNAs based on E. coli H27 for in vitro physicochemical studies. One-dimensional (1)H NMR spectroscopy demonstrates that both the 885 and 888 conformations are occupied to approximately the same extent (f(888) = 0.427 +/- 0.04) in the native H27 sequence at low pH (6.4) and low ionic strength (50 mM NaCl). UV irradiation assays conducted under conditions analogous to those used for assays of ribosomal function (pH 7.5 and 20 mM MgCl(2)) suggest that nucleotides 892 and 905, which are too far apart in the known 885 crystal structures, can approach each other closely enough to form an efficient cross-link. The use of a fluorescence resonance energy transfer (FRET)-labeled RNA together with a partially complementary DNA oligonucleotide that induces a shift to the 888 conformation shows that H27 interchanges between the 885 and 888 conformations on the millisecond time scale, with an equilibrium constant of 0.33 +/-0.12. FRET assays also show that tetracycline interferes with the induced shift to the 888 conformation, a finding that is consistent with crystallographic localization of tetracycline bound to the 885 conformation of H27 in the 30S ribosomal subunit. Taken together, our data demonstrate the innate tendency of an isolated H27 to exist in a dynamic equilibrium between the 885 and 888 conformations. This begs the question of how these inherent structural dynamics are suppressed within the context of the ribosome.

Folding thermodynamics and kinetics of YNMG RNA hairpins: specific incorporation of 8-bromoguanosine leads to stabilization by enhancement of the folding rate

Modified nucleotides allow fundamental energetic and kinetic properties of nucleic acids to be probed. Here, we demonstrate that an RNA hairpin containing the nucleotide analogue 8-bromoguanosine (8BrG or G), gcUUCGgc, has enhanced stability relative to the unmodified hairpin, with DeltaDeltaG(37)(degrees)= -0.69 +/- 0.15 kcal mol(-1) and DeltaT(M) = +6.8 +/- 1.4 degrees C. NMR spectroscopic data suggest that the enhanced stability of gcUUCGgc does not arise from the native state; laser temperature-jump experiments support this notion, as gcUUCGgc and gcUUCGgc have similar unfolding rate constants, but the folding rate constant of gcUUCGgc is 4.1-fold faster at 37.5 degrees C and 2.8-fold faster under isoenergetic conditions. On the basis of these findings, we propose that 8BrG reduces the conformational entropy of the denatured state, resulting in an accelerated conformational search for the native state and enhanced stability.

The Ribosomal Database Project (RDP-II): sequences and tools for high-throughput rRNA analysis

The Ribosomal Database Project (RDP-II) provides the research community with aligned and annotated rRNA gene sequences, along with analysis services and a phylogenetically consistent taxonomic framework for these data. Updated monthly, these services are made available through the RDP-II website (http://rdp.cme.msu.edu/). RDP-II release 9.21 (August 2004) contains 101,632 bacterial small subunit rRNA gene sequences in aligned and annotated format. High-throughput tools for initial taxonomic placement, identification of related sequences, probe and primer testing, data navigation and subalignment download are provided. The RDP-II email address for questions or comments is rdpstaff@msu.edu.

Predicting RNA Structure Using Mutual Information

BACKGROUND

With the ever-increasing number of sequenced RNAs and the establishment of new RNA databases, such as the Comparative RNA Web Site and Rfam, there is a growing need for accurately and automatically predicting RNA structures from multiple alignments. Since RNA secondary structure is often conserved in evolution, the well known, but underused, mutual information measure for identifying covarying sites in an alignment can be useful for identifying structural elements. This article presents MIfold, a MATLAB toolbox that employs mutual information, or a related covariation measure, to display and predict conserved RNA secondary structure (including pseudoknots) from an alignment.

RESULTS

We show that MIfold can be used to predict simple pseudoknots, and that the performance can be adjusted to make it either more sensitive or more selective. We also demonstrate that the overall performance of MIfold improves with the number of aligned sequences for certain types of RNA sequences. In addition, we show that, for these sequences, MIfold is more sensitive but less selective than the related RNAalifold structure prediction program and is comparable with the COVE structure prediction package.

CONCLUSION

MIfold provides a useful supplementary tool to programs such as RNA Structure Logo, RNAalifold and COVE, and should be useful for automatically generating structural predictions for databases such as Rfam.

X-ray crystallography study on ribosome recycling: the mechanism of binding and action of RRF on the 50S ribosomal subunit

This study presents the crystal structure of domain I of the Escherichia coli ribosome recycling factor (RRF) bound to the Deinococcus radiodurans 50S subunit. The orientation of RRF is consistent with the position determined on a 70S-RRF complex by cryoelectron microscopy (cryo-EM). Alignment, however, requires a rotation of 7 degrees and a shift of the cryo-EM RRF by a complete turn of an alpha-helix, redefining the contacts established with ribosomal components. At 3.3 A resolution, RRF is seen to interact exclusively with ribosomal elements associated with tRNA binding and/or translocation. Furthermore, these results now provide a high-resolution structural description of the conformational changes that were suspected to occur on the 70S-RRF complex, which has implications for the synergistic action of RRF with elongation factor G (EF-G). Specifically, the tip of the universal bridge element H69 is shifted by 20 A toward h44 of the 30S subunit, suggesting that RRF primes the intersubunit bridge B2a for the action of EF-G. Collectively, our data enable a model to be proposed for the dual action of EF-G and RRF during ribosome recycling.

The identification of novel RNA structural motifs using COMPADRES: an automated approach to structural discovery

Recurring RNA structural motifs are important sites of tertiary interaction and as such, are integral to RNA macromolecular structure. Although numerous RNA motifs have been classified and characterized, the identification of new motifs is of great interest. In this study, we discovered four new conformationally recurring motifs: the pi-turn, the Omega-turn, the alpha-loop and the C2'-endo mediated flipped adenosine motif. Not only do they have complex and interesting structures, but they participate in contacts of high biological significance. In a first for the RNA field, new motifs were discovered by a fully automated algorithm. This algorithm, COMPADRES, utilized a reduced representation of the RNA backbone and was highly successful at discerning unique structural relationships. This study also shows that recurring RNA substructures are not necessarily accompanied by consistent primary or secondary structure.

A 16S rRNA-tRNA product containing a nucleotide phototrimer and specific for tRNA in the P/E hybrid state in the Escherichia coli ribosome

Ribosome complexes containing deacyl-tRNA1(Val) or biotinylvalyl-tRNA1(Val) and an mRNA analog have been irradiated with wavelengths specific for activation of the cmo5U nucleoside at position 34 in the tRNA1(Val) anticodon loop. The major product for both types of tRNA is the cross-link between 16S rRNA (C1400) and the tRNA (cmo5U34) characterized already by Ofengand and his collaborators [Prince et al. (1982) Proc. Natl Acad. Sci. USA, 79, 5450-5454]. However, in complexes containing deacyl-tRNA1(Val), an additional product is separated by denaturing PAGE and this is shown to involve C1400 and m5C967 of 16S rRNA and cmo5U34 of the tRNA. Puromycin treatment of the biotinylvalyl-tRNA1(Val) -70S complex followed by irradiation, results in the appearance of the unusual photoproduct, which indicates an immediate change in the tRNA interaction with the ribosome after peptide transfer. These results indicate an altered interaction between the tRNA anticodon and the 30S subunit for the tRNA in the P/E hybrid state compared with its interaction in the classic P/P state.

Structure of the Tetrahymena ribozyme: base triple sandwich and metal ion at the active site

The Tetrahymena intron is an RNA catalyst, or ribozyme. As part of its self-splicing reaction, this ribozyme catalyzes phosphoryl transfer between guanosine and a substrate RNA strand. Here we report the refined crystal structure of an active Tetrahymena ribozyme in the absence of its RNA substrate at 3.8 A resolution. The 3'-terminal guanosine (omegaG), which serves as the attacking group for RNA cleavage, forms a coplanar base triple with the G264-C311 base pair, and this base triple is sandwiched by three other base triples. In addition, a metal ion is present in the active site, contacting or positioned close to the ribose of the omegaG and five phosphates. All of these phosphates have been shown to be important for catalysis. Therefore, we provide a picture of how the ribozyme active site positions both a catalytic metal ion and the nucleophilic guanosine for catalysis prior to binding its RNA substrate.

RNA tertiary structure and cooperative assembly of a large ribonucleoprotein complex

The mechanisms that govern the ordered assembly of multiprotein ribonucleoprotein complexes are not well understood. The in vitro reconstitution of the small subunit of the bacterial ribosome provides a tractable system for the detailed study of ordered assembly. We present a quantitative thermodynamic description of the hierarchical binding of ribosomal proteins to 16S rRNA during assembly of the platform of the 30S ribosomal subunit. The binding of S8, S11, S15, and the S6:S18 heterodimer to the central domain of 16S rRNA has been measured both individually and in combination using isothermal titration calorimetry and gel mobility shift assays. Both enthalpy and free energy measurements demonstrate the cooperative binding of S15 and the S6:S18 heterodimer, but no cooperativity is observed for either S8 or S11. The results define a thermodynamic framework that describes cooperative platform assembly.

The consistent phylogenetic signal in genome trees revealed by reducing the impact of noise

Phylogenetic trees based on gene repertoires are remarkably similar to the current consensus of life history. Yet it has been argued that shared gene content is unreliable for phylogenetic reconstruction because of convergence in gene content due to horizontal gene transfer and parallel gene loss. Here we test this argument, by filtering out as noise those orthologous groups that have an inconsistent phylogenetic distribution, using two independent methods. The resulting phylogenies do indeed contain small but significant improvements. More importantly, we find that the majority of orthologous groups contain some phylogenetic signal and that the resulting phylogeny is the only detectable signal present in the gene distribution across genomes. Horizontal gene transfer or parallel gene loss does not cause systematic biases in the gene content tree.

Crystal structure of a phage Twort group I ribozyme–product complex

Group I introns are catalytic RNAs capable of orchestrating two sequential phosphotransesterification reactions that result in self-splicing. To understand how the group I intron active site facilitates catalysis, we have solved the structure of an active ribozyme derived from the orf142-I2 intron from phage Twort bound to a four-nucleotide product RNA at a resolution of 3.6 A. In addition to the three conserved domains characteristic of all group I introns, the Twort ribozyme has peripheral insertions characteristic of phage introns. These elements form a ring that completely envelops the active site, where a snug pocket for guanosine is formed by a series of stacked base triples. The structure of the active site reveals three potential binding sites for catalytic metals, and invokes a role for the 2' hydroxyl of the guanosine substrate in organization of the active site for catalysis.

Crystal structure of a group I intron splicing intermediate

A recently reported crystal structure of an intact bacterial group I self-splicing intron in complex with both its exons provided the first molecular view into the mechanism of RNA splicing. This intron structure, which was trapped in the state prior to the exon ligation reaction, also reveals the architecture of a complex RNA fold. The majority of the intron is contained within three internally stacked, but sequence discontinuous, helical domains. Here the tertiary hydrogen bonding and stacking interactions between the domains, and the single-stranded joiner segments that bridge between them, are fully described. Features of the structure include: (1) A pseudoknot belt that circumscribes the molecule at its longitudinal midpoint; (2) two tetraloop-tetraloop receptor motifs at the peripheral edges of the structure; (3) an extensive minor groove triplex between the paired and joiner segments, P6-J6/6a and P3-J3/4, which provides the major interaction interface between the intron's two primary domains (P4-P6 and P3-P9.0); (4) a six-nucleotide J8/7 single stranded element that adopts a mu-shaped structure and twists through the active site, making critical contacts to all three helical domains; and (5) an extensive base stacking architecture that realizes 90% of all possible stacking interactions. The intron structure was validated by hydroxyl radical footprinting, where strong correlation was observed between experimental and predicted solvent accessibility. Models of the pre-first and pre-second steps of intron splicing are proposed with full-sized tRNA exons. They suggest that the tRNA undergoes substantial angular motion relative to the intron between the two steps of splicing.

Assembly of the 30S ribosomal subunit: positioning ribosomal protein S13 in the S7 assembly branch

Studies of Escherichia coli 30S ribosomal subunit assembly have revealed a hierarchical and cooperative association of ribosomal proteins with 16S ribosomal RNA; these results have been used to compile an in vitro 30S subunit assembly map. In single protein addition and omission studies, ribosomal protein S13 was shown to be dependent on the prior association of ribosomal protein S20 for binding to the ribonucleoprotein particle. While the overwhelming majority of interactions revealed in the assembly map are consistent with additional data, the dependency of S13 on S20 is not. Structural studies position S13 in the head of the 30S subunit > 100 A away from S20, which resides near the bottom of the body of the 30S subunit. All of the proteins that reside in the head of the 30S subunit, except S13, have been shown to be part of the S7 assembly branch, that is, they all depend on S7 for association with the assembling 30S subunit. Given these observations, the assembly requirements for S13 were investigated using base-specific chemical footprinting and primer extension analysis. These studies reveal that S13 can bind to 16S rRNA in the presence of S7, but not S20. Additionally, interaction between S13 and other members of the S7 assembly branch have been observed. These results link S13 to the 3' major domain family of proteins, and the S7 assembly branch, placing S13 in a new location in the 30S subunit assembly map where its position is in accordance with much biochemical and structural data.

Ketolide antimicrobial activity persists after disruption of interactions with domain II of 23S rRNA

Ketolides are the latest derivatives developed from the macrolide erythromycin to improve antimicrobial activity. All macrolides and ketolides bind to the 50S ribosomal subunit, where they come into contact with adenosine 2058 (A2058) within domain V of the 23S rRNA and block protein synthesis. An additional interaction at nucleotide A752 in the rRNA domain II is made via the synthetic carbamate-alkyl-aryl substituent in the ketolides HMR3647 (telithromycin) and HMR3004, and this interaction contributes to their improved activities. Only a few macrolides, including tylosin, come into contact with domain II of the rRNA and do so via interactions with nucleotides G748 and A752. We have disrupted these macrolide-ketolide interaction sites in the rRNA to assess their relative importance for binding. Base substitutions at A752 were shown to confer low levels of resistance to telithromycin but not to HMR3004, while deletion of A752 confers low levels of resistance to both ketolides. Mutations at position 748 confer no resistance. Substitution of guanine at A2058 gives rise to the MLS(B) (macrolide, lincosamide, and streptogramin B) phenotype, which confers resistance to all the drugs. However, resistance to ketolides was abolished when the mutation at position 2058 was combined with a mutation in domain II of the same rRNA. In contrast, the same dual mutations in rRNAs conferred enhanced resistance to tylosin. Our results show that the domain II interactions of telithromycin and HMR3004 differ from each other and from those of tylosin. The data provide no indication that mutations within domain II, either alone or in combination with an A2058 mutation, can confer significant levels of telithromycin resistance.

RMF inactivates ribosomes by covering the peptidyl transferase centre and entrance of peptide exit tunnel

In gram-negative bacteria such as Escherichia coli, protein synthesis is suppressed by the formation of 100S ribosomes under stress conditions. The 100S ribosome, a dimer of 70S ribosomes, is formed by ribosome modulation factor (RMF) binding to the 70S ribosomes. During the stationary phase, most of the 70S ribosomes turn to 100S ribosomes, which have lost translational activity. This 100S formation is called the hibernation process in the ribosome cycle of the stationary phase. If stationary phase cells are transferred to fresh medium, the 100S ribosomes immediately go back to active 70S ribosomes, showing that inactive 100S <--> active 70S interconversion is a major system regulating translation activity in stationary phase cells. To elucidate the mechanisms of translational inactivation, the binding sites of RMF on 23S rRNA in 100S ribosome of E. coli were examined by a chemical probing method using dimethyl sulphate (DMS). As the results, the nine bases in 23S rRNA were protected from DMS modifications and the modification of one base was enhanced. Interestingly A2451 is included among the protected bases, which is thought to be directly involved in peptidyl transferase activity. We conclude that RMF inactivates ribosomes by covering the peptidyl transferase (PTase) centre and the entrance of peptide exit tunnel. It is surprising that the cell itself produces a protein that seems to inhibit protein synthesis in a similar manner to antibiotics and that it can reversibly bind to and release from the ribosome in response to environmental conditions.

The analysis of the complete mitochondrial genome of Lecanicillium muscarium (synonym Verticillium lecanii) suggests a minimum common gene organization in mtDNAs of Sordariomycetes: phylogenetic implications

The mitochondrial genome (mtDNA) of the entomopathogenic fungus Lecanicillium muscarium (synonym Verticillium lecanii) with a total size of 24,499-bp has been analyzed. So far, it is the smallest known mitochondrial genome among Pezizomycotina, with an extremely compact gene organization and only one group-I intron in its large ribosomal RNA (rnl) gene. It contains the 14 typical genes coding for proteins related to oxidative phosphorylation, the two rRNA genes, one intronic ORF coding for a possible ribosomal protein (rps), and a set of 25 tRNA genes which recognize codons for all amino acids, except alanine and cysteine. All genes are transcribed from the same DNA strand. Gene order comparison with all available complete fungal mtDNAs-representatives of all four Phyla are included-revealed some characteristic common features like uninterrupted gene pairs, overlapping genes, and extremely variable intergenic regions, that can all be exploited for the study of fungal mitochondrial genomes. Moreover, a minimum common mtDNA gene order could be detected, in two units, for all known Sordariomycetes namely nad1-nad4-atp8-atp6 and rns-cox3-rnl, which can be extended in Hypocreales, to nad4L-nad5-cob-cox1-nad1-nad4-atp8-atp6 and rns-cox3-rnl nad2-nad3, respectively. Phylogenetic analysis of all fungal mtDNA essential protein-coding genes as one unit, clearly demonstrated the superiority of small genome (mtDNA) over single gene comparisons.

Extensive intraindividual variation in plastid rDNA sequences from the holoparasite Cynomorium coccineum (Cynomoriaceae)

Ribosomal genes are considered to have a high degree of sequence conservation between species and also at higher taxonomic levels. In this paper we document a case where a single individual of Cynomorium coccineum (Cynomoriaceae), a nonphotosynthetic holoparasitic plant, contains highly divergent plastid ribosomal genes. PCR amplification a nearly complete ribosomal DNA cistron was performed using genomic DNA, the products cloned, and the 23S rDNA genes were sequenced from 19 colonies. Of these, five distinct types were identified. Fifteen of the sequences were nearly identical (11 or fewer differences) and these were designated Type I. The remaining types (II-V) were each represented by a single clone and differed from Type I by 93 to 255 changes. Compared with green vascular plants, we found that there are more substitutional differences in the 23S rDNA sequences within a single individual of Cynomorium than among all sequenced photosynthetic vascular plants. Several trends of molecular evolution observed in 16S rDNA from other holoparasitic angiosperms and heterotrophic green algae have been also observed in Cynomorium 23S rDNA. Higher-order structures were constructed for representatives of the five clone types, and in all cases these possessed complete complements of the major structural elements present in functional plastid 23S rRNAs. These data indicate that such molecules may be subject to purifying selection, thus providing indirect evidence that they have retained some degree of functionality. This intraindividual polymorphism is probably a case of plastid heteroplasmy but translocation of ribosomal cistrons to the nucleus or mitochondria has not been tested and therefore cannot be ruled out.

Testing the conservation of the translational machinery over evolution in diverse environments: assaying Thermus thermophilus ribosomes and initiation factors in a coupled transcription-translation system from Escherichia coli

Ribosomes from the extreme thermophile Thermus thermophilus are capable of translation in a coupled transcription-translation system derived from Escherichia coli. At 45 degrees C, T.thermophilus ribosomes translate at approximately 25-30% of the maximal rate of E.coli ribosomes, and synthesize full-length protein. T.thermophilus and E.coli subunits can be combined to effect translation, with the spectrum of proteins produced depending upon the source of the 30S subunit. In this system, T.thermophilus ribosomes function in concert with E.coli translational factors and tRNAs, with elongation and release factors being supplied from the E.coli extract, and purified initiation factors (IFs) being added exogenously. Cloned and purified T.thermophilus IF1, IF2 and IF3 supported the synthesis of the same products in vitro as the E.coli factors, although the relative levels of some polypeptides were factor dependent. We conclude that, at least between these two phylogenetically distant species, translational factor function and subunit-subunit interactions are conserved. This functional compatibility is remarkable given the extreme and highly divergent environments to which these species have adapted.

The molecular basis for A-site mutations conferring aminoglycoside resistance: relationship between ribosomal susceptibility and X-ray crystal structures

Aminoglycoside antibiotics target the 16S ribosomal RNA (rRNA) bacterial A site and induce misreading of the genetic code. Point mutations of the ribosomal A site may confer resistance to aminoglycoside antibiotics. The influence of bacterial mutations (introduced by site-directed mutagenesis) on ribosomal drug susceptibility was investigated in vivo by determination of minimal inhibitory concentrations. To determine the origin of the various resistance phenotypes at a molecular level, the in vivo results were compared with the previously published crystal structures of paromomycin, tobramycin, and geneticin bound to oligonucleotides containing the minimal A site. Two regions appear crucial for binding in the A site: the single adenine residue at position 1408 and the non-Watson-Crick U1406.U1495 pair. The effects of mutations at those positions are modulated by the nature of the substituent at position 6' (either hydroxy or ammonium group) on ring I, by the number of positive charges on the antibiotic, and by the linkage between rings I and III (either 4,5 or 4,6). In particular, the analysis demonstrates: 1) that the C1409-G1491 to A1409-U1491 polymorphism (observed in 15 % of bacteria) is not associated with resistance, which indicates that it does not affect the stacking of ring I on residue 1491, 2) that the high-level resistance to 6'-NH3+ aminoglycosides exhibited by the A1408G mutation most probably results from the inability of ring I forming a pseudo base pair with G1408, which prevents its insertion inside the A site helix, and 3) that mutations of the uracil residues forming the U1406.U1495 pair either to cytosine or to adenine residues mostly confer low to moderate levels of drug resistance, whereas the U1406C/U1495A double mutation confers high-level resistance (except for neomycin), which suggests that aminoglycoside binding to the wild-type A site and its functional consequences strongly depend on a particular geometry of the U1406.U1495 pair. The relationships between the resistance phenotypes observed in vivo and the interactions described at the molecular level define the biological importance of the different structural interactions observed by X-ray crystallography studies.

Spb1p-Directed Formation of Gm2922 in the Ribosome Catalytic Center Occurs at a Late Processing Stage

rRNA molecules undergo extensive posttranscriptional modification, predominantly 2'-O-ribose methylation and pseudouridine formation, both of which are guided by the numerous small nucleolar RNAs in eukaryotes. Here, we describe an exception to this rule. The essential yeast nucleolar protein Spb1p is a site-specific rRNA methyltransferase modifying the universally conserved G2922 that is located within the A loop of the catalytic center of the ribosome. The equivalent position in bacteria is the docking site for aminoacyl-tRNA, and it is critical for translation. In sharp contrast to other 2'-O-methylriboses that are formed on the primary transcript, Gm2922 appears at a late processing stage, during the maturation of the 27S pre-rRNA. Thus, eukaryotes have maintained a site-specific enzyme to catalyze the methylation of a nucleotide that plays a crucial role in ribosome biogenesis and translation.

The structural basis of macrolide-ribosome binding assessed using mutagenesis of 23S rRNA positions 2058 and 2059

Macrolides are a diverse group of antibiotics that inhibit bacterial growth by binding within the peptide tunnel of the 50S ribosomal subunit. There is good agreement about the architecture of the macrolide site from different crystallography studies of bacterial and archaeal 50S subunits. These structures show plainly that 23S rRNA nucleotides A2058 and A2059 are located accessibly on the surface of the tunnel wall where they act as key contact sites for macrolide binding. However, the molecular details of how macrolides fit into this site remain a matter of contention. Here, we have generated an isogenic set of single and dual substitutions at A2058 and A2059 in Mycobacterium smegmatis to investigate the effects of the rRNA mutations on macrolide binding. Resistances conferred to a comprehensive array of 11 macrolide compounds are used to assess models of macrolide binding predicted from the crystal structures. The data indicate that all macrolides and their derivatives bind at the same site in the tunnel with their C5 amino sugar in a similar orientation. Our data are compatible with the lactone rings of 14-membered and 16-membered macrolides adopting different conformations, enabling the latter compounds to avoid a steric clash with 2058G. This difference, together with interactions conveyed via substituents that are specific to certain ketolide and macrolide sub-classes, influences the binding to the large ribosomal subunit. Our genetic data show no support for a derivatized-macrolide binding site that has been proposed to be located further down the tunnel.

Compensatory evolution in the human malaria parasite Plasmodium ovale

The fixation of neutral compensatory mutations in a population depends on the effective population size of the species, which can fluctuate dramatically within a few generations, the mutation rate, and the selection intensity associated with the individual mutations. We observe compensatory mutations and intermediate states in populations of the malaria parasite Plasmodium ovale. The appearance of compensatory mutations and intermediate states in P. ovale raises interesting questions about population structure that could have considerable impact on the control of the associated disease.

The A2453-C2499 wobble base pair in Escherichia coli 23S ribosomal RNA is responsible for pH sensitivity of the peptidyltransferase active site conformation

Peptide bond formation, catalyzed by the ribosomal peptidyltransferase, has long been known to be sensitive to monovalent cation concentrations and pH. More recently, we and others have shown that residue A2451 in the peptidyltransferase center of the Escherichia coli 50S ribosomal subunit changes conformation in response to alterations in pH, depending on ionic conditions and temperature. Two wobble pairs, A2453-C2499 and A2450-C2063, have been proposed as potential candidates to convey pH-dependent flexibility to the peptidyltransferase center. Each is presumed to possess a near-neutral pK(a), and both lie in proximity to A2451. We show through mutagenesis and chemical probing that the identity of the A2453-C2499 base pair, but not the A2450-C2063 base pair, is critical for the pH-dependent structural rearrangement of A2451. We conclude that, while the A2453-C2499 base pair may be important for maintaining the structure of the active site in the E.coli peptidyltransferase center, its lack of conservation makes it, and consequently its near-neutral pK(a), unlikely to contribute to function during peptide bond formation.

A secondary structural model of the 28S rRNA expansion segments D2 and D3 from rootworms and related leaf beetles (Coleoptera: Chrysomelidae; Galerucinae)

We analysed the secondary structure of two expansion segments (D2, D3) of the 28S rRNA gene from 229 leaf beetles (Coleoptera: Chrysomelidae), the majority of which are in the subfamily Galerucinae. The sequences were compared in a multiple sequence alignment, with secondary structure inferred primarily from the compensatory base changes in the conserved helices of the rRNA molecules. This comparative approach yielded thirty helices comprised of base pairs with positional covariation. Based on these leaf beetle sequences, we report an annotated secondary structural model for the D2 and D3 expansion segments that will prove useful in assigning positional nucleotide homology for phylogeny reconstruction in these and closely related beetle taxa. This predicted structure, consisting of seven major compound helices, is mostly consistent with previously proposed models for the D2 and D3 expansion segments in insects. Despite a lack of conservation in the primary structure of these regions of insect 28S rRNA, the evolution of the secondary structure of these seven major motifs may be informative above the nucleotide level for higher-order phylogeny reconstruction of major insect lineages.

Intracellular folding of the Tetrahymena group I intron depends on exon sequence and promoter choice

The Tetrahymena group I intron splices 20 to 50 times faster in Tetrahymena than in vitro, implying that the intron rapidly adopts its active conformation in the cell. The importance of cotranscriptional folding and the contribution of the rRNA exons to the stability of the active pre-RNA structure were investigated by comparing the activity of minimal pre-RNAs expressed in Escherichia coli. Pre-RNAs containing exons derived from E. coli 23 S rRNA were three to four times more active than the wild-type Tetrahymena pre-RNA. E. coli transcripts of the chimeric E. coli pre-RNA were two to eight times more active than were T7 transcripts. However, the effect of cotranscriptional folding depends on exon sequences. Unexpectedly, the unspliced pre-RNA decays more slowly than predicted from the rate of splicing. This observation is best explained by partitioning of transcripts into active and inactive pools. We propose that the active pool splices within a few seconds, whereas the inactive pool is degraded without appreciable splicing.

A comprehensive comparison of comparative RNA structure prediction approaches

BACKGROUND

An increasing number of researchers have released novel RNA structure analysis and prediction algorithms for comparative approaches to structure prediction. Yet, independent benchmarking of these algorithms is rarely performed as is now common practice for protein-folding, gene-finding and multiple-sequence-alignment algorithms.

RESULTS

Here we evaluate a number of RNA folding algorithms using reliable RNA data-sets and compare their relative performance.

CONCLUSIONS

We conclude that comparative data can enhance structure prediction but structure-prediction-algorithms vary widely in terms of both sensitivity and selectivity across different lengths and homologies. Furthermore, we outline some directions for future research.

"Candidatus Endobugula glebosa," a specific bacterial symbiont of the marine bryozoan Bugula simplex

The bryozoans Bugula neritina and Bugula simplex harbor bacteria in the pallial sinuses of their larvae as seen by electron microscopy. In B. neritina, the bacterial symbiont has been characterized as a gamma-proteobacterium, "Candidatus Endobugula sertula." "Candidatus E. sertula" has been implicated as the source of the bryostatins, polyketides that provide chemical defense to the host and are also being tested for use in human cancer treatments. In this study, the bacterial symbiont in B. simplex larvae was identified by 16S rRNA-targeted PCR and sequencing as a gamma-proteobacterium closely related to and forming a monophyletic group with "Candidatus E. sertula." In a fluorescence in situ hybridization, a 16S ribosomal DNA probe specific to the B. simplex symbiont hybridized to long rod-shaped bacteria in the pallial sinus of a B. simplex larva. The taxonomic status "Candidatus Endobugula glebosa" is proposed for the B. simplex larval symbiont. Degenerate polyketide synthase (PKS) primers amplified a gene fragment from B. simplex that closely matched a PKS gene fragment from the bryostatin PKS cluster. PCR surveys show that the symbiont and this PKS gene fragment are consistently and uniquely associated with B. simplex. Bryostatin activity assays and chemical analyses of B. simplex extracts reveal the presence of compounds similar to bryostatins. Taken together, these findings demonstrate a symbiosis in B. simplex that is similar and evolutionarily related to that in B. neritina.

Factors Affecting Thermodynamic Stabilities of RNA 3 × 3 Internal Loops

Internal loops in RNA are important for folding and function. The 3 x 3 nucleotide internal loops are the smallest size symmetric loops with a potential noncanonical base pair (middle pair) flanked on both sides by a noncanonical base pair (loop-terminal pair). Thermodynamic and structural insights acquired for 3 x 3 loops should improve approximations for stabilities of 3 x 3 and larger internal loops. Most natural 3 x 3 internal loops are purine rich, which is also true of other internal loops. A series of oligoribonucleotides containing different 3 x 3 internal loops were studied by UV melting and imino proton NMR. Both loop-terminal and middle pairs contribute to the thermodynamic stabilities of 3 x 3 loops. Extra stabilization of -1.2 kcal/mol was found for a GA middle pair when flanked by at least one non-pyrimidine-pyrimidine loop-terminal pair. A penalty of approximately 1 kcal/mol was found for loops with a single loop-terminal GA pair that has a U 3' to the G of the GA pair. A revised model for predicting stabilities of 3 x 3 loops is derived by multiple linear regression.

Contribution of RPB2 to multilocus phylogenetic studies of the euascomycetes (Pezizomycotina, Fungi) with special emphasis on the lichen-forming Acarosporaceae and evolution of polyspory

Despite the recent progress in molecular phylogenetics, many of the deepest relationships among the main lineages of the largest fungal phylum, Ascomycota, remain unresolved. To increase both resolution and support on a large-scale phylogeny of lichenized and non-lichenized ascomycetes, we combined the protein coding-gene RPB2 with the traditionally used nuclear ribosomal genes SSU and LSU. Our analyses resulted in the naming of the new subclasses Acarosporomycetidae and Ostropomycetidae, and the new class Lichinomycetes, as well as the establishment of the phylogenetic placement and novel circumscription of the lichen-forming fungi family Acarosporaceae. The delimitation of this family has been problematic over the past century, because its main diagnostic feature, true polyspory (numerous spores issued from multiple post-meiosis mitoses) with over 100 spores per ascus, is probably not restricted to the Acarosporaceae. This observation was confirmed by our reconstruction of the origin and evolution of this form of true polyspory using maximum likelihood as the optimality criterion. The various phylogenetic analyses carried out on our data sets allowed us to conclude that: (1) the inclusion of phylogenetic signal from ambiguously aligned regions into the maximum parsimony analyses proved advantageous in reconstructing phylogeny; however, when more data become available, Bayesian analysis using different models of evolution is likely to be more efficient; (2) neighbor-joining bootstrap proportions seem to be more appropriate in detecting topological conflict between data partitions of large-scale phylogenies than posterior probabilities; and (3) Bayesian bootstrap proportion provides a compromise between posterior probability outcomes (i.e., higher accuracy, but with a higher number of significantly supported wrong internodes) vs. maximum likelihood bootstrap proportion outcomes (i.e., lower accuracy, with a lower number of significantly supported wrong internodes).

Diversity and evolution of Bdellovibrio-and-like organisms (BALOs), reclassification of Bacteriovorax starrii as Peredibacter starrii gen. nov., comb. nov., and description of the Bacteriovorax–Peredibacter clade as Bacteriovoracaceae fam. nov

A phylogenetic analysis of Bdellovibrio-and-like organisms (BALOs) was performed. It was based on the characterization of 71 strains and on all consequent 16S rRNA gene sequences available in databases, including clones identified by data-mining, totalling 120 strains from very varied biotopes. Amplified rDNA restriction analysis (ARDRA) accurately reflected the diversity and phylogenetic affiliation of BALOs, thereby providing an efficient screening tool. Extensive phylogenetic analysis of the 16S rRNA gene sequences revealed great diversity within the Bdellovibrio (>14 % divergence) and Bacteriovorax (>16 %) clades, which comprised nine and eight clusters, respectively, exhibiting more than 3 % intra-cluster divergence. The clades diverged by more than 20 %. The analysis of conserved 16S rRNA secondary structures showed that Bdellovibrio contained motifs atypical of the δ-Proteobacteria, suggesting that it is ancestral to Bacteriovorax. While none of the Bdellovibrio strains were of marine origin, Bacteriovorax included separate soil/freshwater and marine-specific groups. On the basis of their extensive diversity and the large distance separating the groups, it is proposed that Bacteriovorax starrii be placed into a new genus, Peredibacter gen. nov., with Peredibacter starrii A3.12T (=ATCC 15145T=NCCB 72004T) as its type strain. Also proposed is a redefinition of the Bdellovibrio and the Bacteriovorax–Peredibacter lineages as two different families, i.e. ‘Bdellovibrionaceae’ and a new family, Bacteriovoracaceae. Also, a re-evaluation of oligonucleotides targeting BALOs is presented, and the implications of the large diversity of these organisms and of their distribution in very different environments are discussed.

Co-transcriptional folding is encoded within RNA genes

Background

Most of the existing RNA structure prediction programs fold a completely synthesized RNA molecule. However, within the cell, RNA molecules emerge sequentially during the directed process of transcription. Dedicated experiments with individual RNA molecules have shown that RNA folds while it is being transcribed and that its correct folding can also depend on the proper speed of transcription.

Methods

The main aim of this work is to study if and how co-transcriptional folding is encoded within the primary and secondary structure of RNA genes. In order to achieve this, we study the known primary and secondary structures of a comprehensive data set of 361 RNA genes as well as a set of 48 RNA sequences that are known to differ from the originally transcribed sequence units. We detect co-transcriptional folding by defining two measures of directedness which quantify the extend of asymmetry between alternative helices that lie 5' and those that lie 3' of the known helices with which they compete.

Results

We show with statistical significance that co-transcriptional folding strongly influences RNA sequences in two ways: (1) alternative helices that would compete with the formation of the functional structure during co-transcriptional folding are suppressed and (2) the formation of transient structures which may serve as guidelines for the co-transcriptional folding pathway is encouraged.

Conclusions

These findings have a number of implications for RNA secondary structure prediction methods and the detection of RNA genes.

Evaluation of the suitability of free-energy minimization using nearest-neighbor energy parameters for RNA secondary structure prediction

BACKGROUND

A detailed understanding of an RNA's correct secondary and tertiary structure is crucial to understanding its function and mechanism in the cell. Free energy minimization with energy parameters based on the nearest-neighbor model and comparative analysis are the primary methods for predicting an RNA's secondary structure from its sequence. Version 3.1 of Mfold has been available since 1999. This version contains an expanded sequence dependence of energy parameters and the ability to incorporate coaxial stacking into free energy calculations. We test Mfold 3.1 by performing the largest and most phylogenetically diverse comparison of rRNA and tRNA structures predicted by comparative analysis and Mfold, and we use the results of our tests on 16S and 23S rRNA sequences to assess the improvement between Mfold 2.3 and Mfold 3.1.

RESULTS

The average prediction accuracy for a 16S or 23S rRNA sequence with Mfold 3.1 is 41%, while the prediction accuracies for the majority of 16S and 23S rRNA structures tested are between 20% and 60%, with some having less than 20% prediction accuracy. The average prediction accuracy was 71% for 5S rRNA and 69% for tRNA. The majority of the 5S rRNA and tRNA sequences have prediction accuracies greater than 60%. The prediction accuracy of 16S rRNA base-pairs decreases exponentially as the number of nucleotides intervening between the 5' and 3' halves of the base-pair increases.

CONCLUSION

Our analysis indicates that the current set of nearest-neighbor energy parameters in conjunction with the Mfold folding algorithm are unable to consistently and reliably predict an RNA's correct secondary structure. For 16S or 23S rRNA structure prediction, Mfold 3.1 offers little improvement over Mfold 2.3. However, the nearest-neighbor energy parameters do work well for shorter RNA sequences such as tRNA or 5S rRNA, or for larger rRNAs when the contact distance between the base-pairs is less than 100 nucleotides.

Using an RNA secondary structure partition function to determine confidence in base pairs predicted by free energy minimization

A partition function calculation for RNA secondary structure is presented that uses a current set of nearest neighbor parameters for conformational free energy at 37 degrees C, including coaxial stacking. For a diverse database of RNA sequences, base pairs in the predicted minimum free energy structure that are predicted by the partition function to have high base pairing probability have a significantly higher positive predictive value for known base pairs. For example, the average positive predictive value, 65.8%, is increased to 91.0% when only base pairs with probability of 0.99 or above are considered. The quality of base pair predictions can also be increased by the addition of experimentally determined constraints, including enzymatic cleavage, flavin mono-nucleotide cleavage, and chemical modification. Predicted secondary structures can be color annotated to demonstrate pairs with high probability that are therefore well determined as compared to base pairs with lower probability of pairing.

Phylogenetic relationships among 28 spirotrichous ciliates documented by rDNA

The contiguous sequence of the SSU rDNA, ITS 1, 5.8S, ITS 2, and approximately 1370 bp at the 5(') end of the LSU rDNA was determined in 25 stichotrichs, one oligotrich, and two hypotrichs. Maximum parsimony, neighbor-joining, and quartet-puzzling analyses were used to construct individual phylogenetic trees for SSU rDNA, for LSU rDNA, and ITS 1+5.8S+ITS 2, as well as for all these components combined. All trees were similar, with the greatest resolution obtained with the combined components. Phylogenetic relationships were largely consistent with classical taxonomy, with notable disagreements. DNA sequences indicate that Oxytricha granulifera and Oxytricha longa are rather distantly related. The oligotrich, Halteria grandinella, is placed well within the order Stichotrichida. Uroleptus pisces and Uroleptus gallina probably belong to different genera. Holosticha polystylata (family Holostichidae) and Urostyla grandis (family Urostylidae) are rather closely related. These rDNA sequence analyses imply the need for some modifications of classical taxonomic schemes.

Exploration of the conserved A+C wobble pair within the ribosomal peptidyl transferase center using affinity purified mutant ribosomes

Protein synthesis in the ribosome's large subunit occurs within an active site comprised exclusively of RNA. Mutational studies of rRNA active site residues could provide valuable insight into the mechanism of peptide bond formation, but many of these mutations cause a dominant lethal phenotype, which prevents production of the homogeneous mutant ribosomes needed for analysis. We report a general method to affinity purify in vivo assembled 50S ribosomal subunits containing lethal active site mutations via a U1A protein-binding tag inserted onto the 23S rRNA. The expected pH-dependent formation of the A2450+C2063 wobble pair has made it a potential candidate for the pH-dependent conformational change that occurs within the ribosomal active site. Using this approach, the active site A2450+C2063 pair was mutated to the isosteric, but pH-independent, G2450*U2063 wobble pair, and 50S subunits containing the mutations were affinity purified. The G*U mutation caused the adjacent A2451 to become hyper-reactive to dimethylsulfate (DMS) modification in a pH-independent manner. Furthermore, the G*U mutation decreased both the rate of peptide bond formation and the affinity of the post-translocation complex for puromycin. The reaction rate (k(pep)) was reduced approximately 200-fold for both puromycin and the natural aminoacyl-tRNA A-site substrate. The mutations also substantially altered the pH dependence of the reaction. Mutation of this base pair has significant deleterious effects upon peptidyl transferase activity, but because G*U mutation disrupts several tertiary contacts with the wobble pair, the assignment of A2450 as the active site residue with the neutral pK(a) important for the peptidyl transferase reaction cannot be fully supported or excluded based upon these data.

Crystal structure of KsgA, a universally conserved rRNA adenine dimethyltransferase in Escherichia coli

The bacterial enzyme KsgA catalyzes the transfer of a total of four methyl groups from S-adenosyl-l-methionine (S-AdoMet) to two adjacent adenosine bases in 16S rRNA. This enzyme and the resulting modified adenosine bases appear to be conserved in all species of eubacteria, eukaryotes, and archaebacteria, and in eukaryotic organelles. Bacterial resistance to the aminoglycoside antibiotic kasugamycin involves inactivation of KsgA and resulting loss of the dimethylations, with modest consequences to the overall fitness of the organism. In contrast, the yeast ortholog, Dim1, is essential. In yeast, and presumably in other eukaryotes, the enzyme performs a vital role in pre-rRNA processing in addition to its methylating activity. Another ortholog has been discovered recently, h-mtTFB in human mitochondria, which has a second function; this enzyme is a nuclear-encoded mitochondrial transcription factor. The KsgA enzymes are homologous to another family of RNA methyltransferases, the Erm enzymes, which methylate a single adenosine base in 23S rRNA and confer resistance to the MLS-B group of antibiotics. Despite their sequence similarity, the two enzyme families have strikingly different levels of regulation that remain to be elucidated. We have crystallized KsgA from Escherichia coli and solved its structure to a resolution of 2.1A. The structure bears a strong similarity to the crystal structure of ErmC' from Bacillus stearothermophilus and a lesser similarity to sc-mtTFB, the Saccharomyces cerevisiae version of h-mtTFB. Comparison of the three crystal structures and further study of the KsgA protein will provide insight into this interesting group of enzymes.

RAG: RNA-As-Graphs web resource

BACKGROUND

The proliferation of structural and functional studies of RNA has revealed an increasing range of RNA's structural repertoire. Toward the objective of systematic cataloguing of RNA's structural repertoire, we have recently described the basis of a graphical approach for organizing RNA secondary structures, including existing and hypothetical motifs.

DESCRIPTION

We now present an RNA motif database based on graph theory, termed RAG for RNA-As-Graphs, to catalogue and rank all theoretically possible, including existing, candidate and hypothetical, RNA secondary motifs. The candidate motifs are predicted using a clustering algorithm that classifies RNA graphs into RNA-like and non-RNA groups. All RNA motifs are filed according to their graph vertex number (RNA length) and ranked by topological complexity.

CONCLUSIONS

RAG's quantitative cataloguing allows facile retrieval of all classes of RNA secondary motifs, assists identification of structural and functional properties of user-supplied RNA sequences, and helps stimulate the search for novel RNAs based on predicted candidate motifs.

High overall diversity and dominance of microdiverse relationships in salt marsh sulphate-reducing bacteria

The biogeochemistry of North Atlantic salt marshes is characterized by the interplay between the marsh grass Spartina and sulphate-reducing bacteria (SRB), which mineralize the diverse carbon substrates provided by the plants. It was hypothesized that SRB populations display high diversity within the sediment as a result of the rich spatial and chemical structuring provided by Spartina roots. A 2000-member 16S rRNA gene library, prepared with delta-proteobacterial SRB-selective primers, was analysed for diversity patterns and phylogenetic relationships. Sequence clustering detected 348 16S rRNA sequence types (ribotypes) related to delta-proteobacterial SRB, and it was estimated that a total of 623 ribotypes were present in the library. Similarity clustering showed that approximately 46% of these sequences fell into groups with < 1% divergence; thus, microheterogeneity accounts for a large portion of the observable genetic diversity. Phylogenetic comparison revealed that sequences most frequently recovered were associated with the Desulfobacteriaceae and Desulfobulbaceae families. Sequences from the Desulfovibrionaceae family were also observed, but were infrequent. Over 80% of the delta-proteobacterial ribotypes clustered with cultured representatives of Desulfosarcina, Desulfococcus and Desulfobacterium genera, suggesting that complete oxidizers with high substrate versatility dominate. The large-scale approach demonstrates the co-existence of numerous SRB-like sequences and reveals an unexpected amount of microdiversity.

Long-term evolution of the S788 fungal nuclear small subunit rRNA group I introns

More than 1000 group I introns have been identified in fungal rDNA. Little is known, however, of the splicing and secondary structure evolution of these ribozymes. Here, we use a combination of comparative and biochemical methods to address the evolution and splicing of a vertically inherited group I intron found at position 788 in the fungal small subunit (S) rRNA. The ancestral state of the S788 intron contains a highly conserved core and an extended P5 domain typical of IC1 introns. In contrast, the more derived introns have lost most of P5, and have an accelerated divergence rate within the core region with three functionally important substitutions that unambiguously separate them from the ancestral pool. Of 14 S788 group I introns that were tested for splicing, five, all of the ancestral type, were able to self-splice and produced intron RNA circles in vitro. The more derived S788 introns did not self-splice, and potentially rely on fungal-specific factors to facilitate splicing. In summary, we demonstrate one possible fate of vertically inherited group I introns, the loss of secondary structure elements, lessened selective constraints in the intron core, and ultimately, dependence on host-mediated splicing.