A proposal for the secondary structure of a variable area of eukaryotic small ribosomal subunit RNA involving the existence of a pseudoknot

About the Analysis of 18S rDNA Sequence Data from Trypanosomes in Barcoding and Phylogenetics: Tracing a Continuation Error Occurring in the Literature

The variable regions (V1-V9) of the 18S rDNA are routinely used in barcoding and phylogenetics. In handling these data for trypanosomes, we have noticed a misunderstanding that has apparently taken a life of its own in the literature over the years. In particular, in recent years, when studying the phylogenetic relationship of trypanosomes, the use of V7/V8 was systematically established. However, considering the current numbering system for all other organisms (including other Euglenozoa), V7/V8 was never used. In Maia da Silva et al. [Parasitology 2004, 129, 549-561], V7/V8 was promoted for the first time for trypanosome phylogenetics, and since then, more than 70 publications have replicated this nomenclature and even discussed the benefits of the use of this region in comparison to V4. However, the primers used to amplify the variable region of trypanosomes have actually amplified V4 (concerning the current 18S rDNA numbering system).

Metabolism of Milk Oligosaccharides in Preterm Pigs Sensitive to Necrotizing Enterocolitis

Human milk oligosaccharides (HMO) are major components of breast milk that may have local effects in the gastrointestinal tract and systemic functions after being absorbed, both depending on their metabolism. Using preterm pigs, we investigated the metabolic fate of HMO in three experiments with two different HMO blends. In addition, we examined effects on the colonic microbiota in the presence or absence of necrotizing enterocolitis (NEC). Thus, preterm pigs (n = 112) were fed formula without or with HMO supplementation (5-10) g/L of a mixture of 4 (4-HMO) or >25 HMO (25-HMO) for 5 (Experiment 1 and 2) or 11 days (Experiment 3). Individual HMO were quantified in colon contents and urine using MALDI-TOF-MS (matrix-assisted laser desorption ionization mass spectrometry) and HPAEC-PAD (high-performance anion-exchange chromatography with pulsed amperometric detection). Microbial colonization was analyzed by sequencing of 16S rRNA gene tags. Intestinal permeability was measured by lactulose to mannitol ratio in urine. HMO supplemented to formula were detected in urine and colon contents in preterm piglets after 5 and 11 days in all three experiments. The amount of HMO excreted via the gut or the kidneys showed large individual variations. Microbial diversity in the colon changed from high levels of Firmicutes (dominated by Clostridium) at day 5 (Exp 2) to high levels of Proteobacteria dominated by Helicobacter and Campylobacter at day 11 (Exp 3). Colonic microbiota composition as well as HMO excretion pattern varied greatly among piglets. Interestingly, the 5-day supplementation of the complex 25-HMO blend led to low concentrations of 3-fucosyllactose (FL) and lacto-N-fucopentaose (LNFP) I in colonic contents, indicating a preferred utilization of these two HMO. Although the interpretation of the data from our piglet study is difficult due to the large individual variation, the presence of Bifidobacteria, although low in total numbers, was correlated with total HMO contents, and specifically with 2'FL levels in colonic content. However, early supplementation of formula with HMO did not affect NEC incidence.

Stability and resilience of the intestinal microbiota in children in daycare - a 12 month cohort study

BACKGROUND

We performed a 12-month cohort study of the stability and resilience of the intestinal microbiota of healthy children in daycare in Denmark in relation to diarrheal events and exposure to known risk factors for gastrointestinal health such as travelling and antibiotic use. In addition, we analyzed how gut microbiota recover from such exposures.

RESULTS

We monitored 32 children in daycare aged 1-6 years. Fecal samples were submitted every second month during a one-year observational period. Information regarding exposures and diarrheal episodes was obtained through questionnaires. Bacterial communities were identified using 16S rRNA gene sequencing. The core microbiota (mean abundance > 95%) dominated the intestinal microbiota, and none of the tested exposures (diarrheal events, travel, antibiotic use) were associated with decreases in the relative abundance of the core microbiota. Samples exhibited lower intra-individual variation than inter-individual variation. Half of all the variation between samples was explained by which child a sample originated from. Age explained 7.6-9.6% of the variation, while traveling, diarrheal events, and antibiotic use explained minor parts of the beta diversity. We found an age-dependent increase of alpha diversity in children aged 1-3 years, and while diarrheal events caused a decrease in alpha diversity, a recovery time of 40-45 days was observed. Among children having had a diarrheal event, we observed a 10x higher relative abundance of Prevotella. After travelling, a higher abundance of two Bacteroides species and 40% less Lachnospiraceae were seen. Antibiotic use did not correlate with changes in the abundance of any bacteria.

CONCLUSION

We present data showing that Danish children in daycare have stable intestinal microbiota, resilient to the exposures investigated. An early age-dependent increase in the diversity was demonstrated. Diarrheal episodes decreased alpha diversity with an estimated recovery time of 40-45 days.

Maturation of the gut microbiome and risk of asthma in childhood

The composition of the human gut microbiome matures within the first years of life. It has been hypothesized that microbial compositions in this period can cause immune dysregulations and potentially cause asthma. Here we show, by associating gut microbial composition from 16S rRNA gene amplicon sequencing during the first year of life with subsequent risk of asthma in 690 participants, that 1-year-old children with an immature microbial composition have an increased risk of asthma at age 5 years. This association is only apparent among children born to asthmatic mothers, suggesting that lacking microbial stimulation during the first year of life can trigger their inherited asthma risk. Conversely, adequate maturation of the gut microbiome in this period may protect these pre-disposed children.

Colonization of commensal bacteria is thought to impact immune development, especially in the earliest years of life. Here, the authors show, by analyzing the development of the gut microbiome of 690 children, that microbial composition at the age of 1 year is associated with asthma diagnosed in the first 5 years of life.

The developing hypopharyngeal microbiota in early life

BACKGROUND

The airways of healthy humans harbor a distinct microbial community. Perturbations in the microbial community have been associated with disease, yet little is known about the formation and development of a healthy airway microbiota in early life. Our goal was to understand the establishment of the airway microbiota within the first 3 months of life. We investigated the hypopharyngeal microbiota in the unselected COPSAC₂₀₁₀ cohort of 700 infants, using 16S rRNA gene sequencing of hypopharyngeal aspirates from 1 week, 1 month, and 3 months of age.

RESULTS

Our analysis shows that majority of the hypopharyngeal microbiota of healthy infants belong to each individual's core microbiota and we demonstrate five distinct community pneumotypes. Four of these pneumotypes are dominated by the genera Staphylococcus, Streptococcus, Moraxella, and Corynebacterium, respectively. Furthermore, we show temporal pneumotype changes suggesting a rapid development towards maturation of the hypopharyngeal microbiota and a significant effect from older siblings. Despite an overall common trajectory towards maturation, individual infants' microbiota are more similar to their own, than to others, over time.

CONCLUSIONS

Our findings demonstrate a consolidation of the population of indigenous bacteria in healthy airways and indicate distinct trajectories in the early development of the hypopharyngeal microbiota.

Evolution of helix formation in the ribosomal Internal Transcribed Spacer 2 (ITS2) and its significance for RNA secondary structures

Helices are the most common elements of RNA secondary structure. Despite intensive investigations of various types of RNAs, the evolutionary history of the formation of new helices (novel helical structures) remains largely elusive. Here, by studying the nuclear ribosomal Internal Transcribed Spacer 2 (ITS2), a fast-evolving part of the eukaryotic nuclear ribosomal operon, we identify two possible types of helix formation: one type is "dichotomous helix formation"--transition from one large helix to two smaller helices by invagination of the apical part of a helix, which significantly changes the shape of the original secondary structure but does not increase its complexity (i.e., the total length of the RNA). An alternative type is "lateral helix formation"--origin of an extra helical region by the extension of a bulge loop or a spacer in a multi-helix loop of the original helix, which does not disrupt the pre-existing structure but increases RNA size. Moreover, we present examples from the RNA sequence literature indicating that both types of helix formation may have implications for RNA evolution beyond ITS2.

Molecular diversity of phototrophic biofilms on building stone

Composition and diversity of aeroterrestrial phototrophic microbial communities are up to now poorly understood. Here, we present a comparative study addressing the composition of algal communities on sandstone substrata based upon the analysis of rRNA gene clone libraries from environmental samples and crude cultures. From a west-facing, shaded wall area of the mediaeval castle ruin Gleichen (Thuringia, Germany), sequences mainly related to the green algae Prasiococcus and Trebouxia (Trebouxiophyceae) were retrieved. A south-west-facing, sun-exposed wall area was mainly colonized by Apatococcus and a Phyllosiphon-related alga. Just a few species, in particular Stichococcus-related strains, were ubiquitous in both areas. Samples from a basement vault exposed to low irradiance exhibited Chlorophyceae like Chromochloris and Bracteacoccus. Thus, most green algae on the daylight-exposed walls were affiliated to Trebouxiophyceae, whereas Chlorophyceae were dominant in samples taken from the site kept under low irradiance. Accordingly, cyanobacterial communities were different: the sun-exposed area was dominated by Synechococcus-related organisms, while on the shaded wall area, cyanobacteria were almost absent. The filamentous Leptolyngbya dominated samples from the basement vault. Scanning electron microscopy revealed endolithic algal morphotypes (coccoid algae and diatoms) dominant in open pores between mineral particles. Here, the organisms may be also involved in biogenic weathering of stone.

Molecular Evidence for the Wide Distribution of Two Lineages of Terrestrial Green Algae (Chlorophyta) over Tropics to Temperate Zone

Phylogenetic analyses of 18S rDNA sequences from environmental clones and culture strains revealed a widespread distribution of two subaerial green algal lineages, Jenufa and Xylochloris, recently described from rainforests in southeast Asia. A new lineage of Jenufa (Chlorophyceae), most closely related to or even conspecific with J. minuta, was formed by sequences of European origin. Two more lineages of Jenufa were formed by three additional sequences from Ecuador and Panama. The other lineage was a close relative of Xylochloris irregularis (Trebouxiophyceae), probably representing a new species of the genus and distinct from the only so far described species, X. irregularis. It comprised two distinct clades each containing almost identical sequences from Germany and Ecuador. Analyses of the new sequences for both genera allowed to presume a preference of J. minuta to subaerial growth on rock or artificial hard substrates combined with a remarkable adaptation to extended periods of darkness, whereas Xylochloris may preferably occur on tree bark or in the soil.

A comparison of the crystal structures of eukaryotic and bacterial SSU ribosomal RNAs reveals common structural features in the hypervariable regions

While the majority of the ribosomal RNA structure is conserved in the three major domains of life--archaea, bacteria, and eukaryotes, specific regions of the rRNA structure are unique to at least one of these three primary forms of life. In particular, the comparative secondary structure for the eukaryotic SSU rRNA contains several regions that are different from the analogous regions in the bacteria. Our detailed analysis of two recently determined eukaryotic 40S ribosomal crystal structures, Tetrahymena thermophila and Saccharomyces cerevisiae, and the comparison of these results with the bacterial Thermus thermophilus 30S ribosomal crystal structure: (1) revealed that the vast majority of the comparative structure model for the eukaryotic SSU rRNA is substantiated, including the secondary structure that is similar to both bacteria and archaea as well as specific for the eukaryotes, (2) resolved the secondary structure for regions of the eukaryotic SSU rRNA that were not determined with comparative methods, (3) identified eukaryotic helices that are equivalent to the bacterial helices in several of the hypervariable regions, (4) revealed that, while the coaxially stacked compound helix in the 540 region in the central domain maintains the constant length of 10 base pairs, its two constituent helices contain 5+5 bp rather than the 6+4 bp predicted with comparative analysis of archaeal and eukaryotic SSU rRNAs.

Bacterial community dynamics during the application of a Myxococcus xanthus-inoculated culture medium used for consolidation of ornamental limestone

In this study, we investigated under laboratory conditions the bacterial communities inhabiting quarry and decayed ornamental carbonate stones before and after the application of a Myxococcus xanthus-inoculated culture medium used for consolidation of the stones. The dynamics of the community structure and the prevalence of the inoculated bacterium, M. xanthus, were monitored during the time course of the consolidation treatment (30 days). For this purpose, we selected a molecular strategy combining fingerprinting by denaturing gradient gel electrophoresis (DGGE) with the screening of eubacterial 16S rDNA clone libraries by DGGE and sequencing. Quantification of the inoculated strain was performed by quantitative real-time PCR (qPCR) using M. xanthus-specific primers designed in this work. Results derived from DGGE and sequencing analysis showed that, irrespective of the origin of the stone,the same carbonatogenic microorganisms were activated by the application of a M. xanthus culture. Those microorganisms were Pseudomonas sp., Bacillus sp., and Brevibacillus sp. The monitoring of M. xanthus in the culture media of treated stones during the time course experiment showed disparate results depending on the applied technique. By culture-dependent methods, the detection of this bacterium was only possible in the first day of the treatment, showing the limitation of these conventional techniques. By PCR-DGGE analysis, M. xanthus was detected during the first 3-6 days of the experiment. At this time, the population of this bacterium in the culture media varied between 108-106 cells ml-1, as showed by qPCR analyses. Thereafter, DGGE analyses showed to be not suitable for the detection of M. xanthus in a mixed culture. Nevertheless, qPCR analysis using specific primers for M. xanthus showed to bea more sensitive technique for the detection of thisbacterium, revealing a population of 104 cells ml-1 in the culture media of both treated stones at the end of the consolidation treatment. The molecular strategy used in this study is proposed as an effective monitoring system to evaluate the impact of the application of a bacterially induced carbonate mineralization as restoration/conservation treatment for ornamental stones.

Molecular and microscopical investigation of the microflora inhabiting a deteriorated Italian manuscript dated from the thirteenth century

This case study shows the application of nontraditional diagnostic methods to investigate the microbial consortia inhabiting an ancient manuscript. The manuscript was suspected to be biologically deteriorated and SEM observations showed the presence of fungal spores attached to fibers, but classic culturing methods did not succeed in isolating microbial contaminants. Therefore, molecular methods, including PCR, denaturing gradient gel electrophoresis (DGGE), and clone libraries, were used as a sensitive alternative to conventional cultivation techniques. DGGE fingerprints revealed a high biodiversity of both bacteria and fungi inhabiting the manuscript. DNA sequence analysis confirmed the existence of fungi and bacteria in manuscript samples. A number of fungal clones identified on the manuscript showed similarity to fungal species inhabiting dry or saline environments, suggesting that the manuscript environment selects for osmophilic or xerophilic fungal species. Most of the bacterial sequences retrieved from the manuscript belong to phylotypes with cellulolytic activities.

In situ study of tetrachloroethylene bioremediation with different microbial community shifting

In this study, we characterized the microbial community in groundwater contaminated with tetrachloroethylene (PCE) in order to evaluate the intrinsic and enhanced bioremediation of PCE. Variable behaviour of microbes was observed between natural attenuation and biostimulation, where the latter was mediated by the addition of nutrients. Results of denaturing gradient gel electrophoresis (DGGE) of amplified bacterial 16S rDNA in the case of biostimulation showed that the microbial community was dominated by species phylogenetically related to the beta-proteobacteria. With regards to natural attenuation, sequences were found belonging to multiple species of different phyla. Interestingly, we found sequences that matched the species belonging to the Firmicutes, which contains bacteria capable of reductive dehalogenation. These results suggest the possibility of the presence of some Clostridium-like PCE degraders within the microbial community when using bioremediation or biostimulation.

Comprehensive molecular structure of the eukaryotic ribosome

Despite the emergence of a large number of X-ray crystallographic models of the bacterial 70S ribosome over the past decade, an accurate atomic model of the eukaryotic 80S ribosome is still not available. Eukaryotic ribosomes possess more ribosomal proteins and ribosomal RNA than do bacterial ribosomes, which are implicated in extraribosomal functions in the eukaryotic cells. By combining cryo-EM with RNA and protein homology modeling, we obtained an atomic model of the yeast 80S ribosome complete with all ribosomal RNA expansion segments and all ribosomal proteins for which a structural homolog can be identified. Mutation or deletion of 80S ribosomal proteins can abrogate maturation of the ribosome, leading to several human diseases. We have localized one such protein unique to eukaryotes, rpS19e, whose mutations are associated with Diamond-Blackfan anemia in humans. Additionally, we characterize crucial interactions between the dynamic stalk base of the ribosome with eukaryotic elongation factor 2.

X-cell parasites in the European dab Limanda limanda are related to other X-cell organisms: a discussion on the potential identity of this new group of parasites

Unusual tumour-like pathologies caused by mysterious cells termed ‘X-cells’ have been reported from numerous fish groups worldwide. After nearly 100 years of research, the tumour-like growths have recently been shown to be caused by a protozoan parasite. In the present study, histopathology and small subunit ribosomal DNA (SSU rDNA) sequences are used to assess whether the X-cell parasite infecting Atlantic dab Limanda limanda L. is distinct from the X-cell parasite infecting Japanese flounder and goby, and to determine their systematic position within the protists. SSU rDNA from Scottish dab was 89·3% and 86·7% similar to Japanese X-cell sequences from flounder and goby respectively, indicating that the parasite infecting dab in the Atlantic is distinct from the Pacific species. Histological studies revealed significant gill pathology and demonstrated the precise location of the parasites within the gill tissues using specific in situ hybridization probes. Phylogenetic analyses showed that the X-cell parasites from Scotland and Japan form a monophyletic group within the Myzozoa, and are basal alveolates. However, ultrastructure of X-cells from dab fails to confirm this systematic placement.

The secondary structure of the unusually long 18S ribosomal RNA of the myxozoan Sphaerospora truttae and structural evolutionary trends in the Myxozoa

The nearly complete 18S rRNA sequence of the myxozoan parasite Sphaerospora truttae shows an extraordinary length (2,552bp) in comparison with other myxozoans and with metazoans in general (average 1,800-1,900bp). The sequence shows nucleotide insertions in most variable regions of the 18S rRNA (V2, V4, V5 and V7), with especially large expansion segments in V4 and V7. In the myxozoans, nucleotide insertions and specific secondary structures in these regions of the gene were found to be strongly related to large scale phylogenetic clustering and thus with the invertebrate host type. Whereas expansion segments were generally found to be absent in the malacasporeans and the clade of primary marine myxozoan species, they occur in all taxa of the clade containing freshwater species, where they showed a consistent secondary structure throughout. The longest expansion segments occur in S. truttae, Sphaerospora elegans and Leptotheca ranae, which represent a clade that has emerged after the malacosporeans and before the radiation of all other myxozoan genera. These three species demonstrate structural links to the malacosporeans as well as other unique features. A smaller number of nucleotide insertions in different subhelices and specific secondary structures appear to have evolved independently in two marine genera, i.e. Ceratomyxa and Parvicapsula. The secondary structural elements of V4 and V7 of the myxozoan 18S rRNAs were found to be highly informative and revealed evolutionary trends of various regions of the gene hitherto unknown, since previous analyses have been based on primary sequence data excluding these regions. Furthermore, the unique features of the V4 region in S. truttae allowed for the design of a highly specific PCR assay for this species.

High local and global diversity of Flavobacteria in marine plankton

Members of the phylum Bacteroidetes are among the most abundant microbes in coastal marine waters, but it is unclear to which extent the diversity within this phylum is covered by currently available 16S rRNA gene sequence information. We, thus, obtained a comprehensive collection of sequence types affiliated with Bacteroidetes in coastal North Sea surface waters and we compared this local diversity with the available sequences of marine planktonic and other aquatic Bacteroidetes. Approximately 15% of > 600 clones from two libraries (August 2000, June 2001) were related to Bacteroidetes, specifically to the Flavobacteria. Local diversity appeared to be almost exhaustively sampled. However, the diversity of the two libraries virtually did not overlap, indicating a pronounced temporal variability of the planktonic Flavobacteria assemblage. The majority of sequence types represented novel phylogenetic lineages, adding 6-7% to the currently known genera and species of Bacteroidetes in marine waters. Different diversity estimators suggested that so far only approximately half of the global diversity of planktonic marine Bacteroidetes has been described. The data set moreover indicated that cultivation-independent techniques and isolation approaches have recovered almost equally sized and virtually non-overlapping fractions of the currently known diversity within this phylum. Interestingly, only 15% of genera of Bacteroidetes from various aquatic environments appear to occur in more than one habitat type.

Characteristics of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) rRNA genes of Apis mellifera (Insecta: Hymenoptera): structure, organization, and retrotransposable elements

As an accompanying manuscript to the release of the honey bee genome, we report the entire sequence of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) ribosomal RNA (rRNA)-encoding gene sequences (rDNA) and related internally and externally transcribed spacer regions of Apis mellifera (Insecta: Hymenoptera: Apocrita). Additionally, we predict secondary structures for the mature rRNA molecules based on comparative sequence analyses with other arthropod taxa and reference to recently published crystal structures of the ribosome. In general, the structures of honey bee rRNAs are in agreement with previously predicted rRNA models from other arthropods in core regions of the rRNA, with little additional expansion in non-conserved regions. Our multiple sequence alignments are made available on several public databases and provide a preliminary establishment of a global structural model of all rRNAs from the insects. Additionally, we provide conserved stretches of sequences flanking the rDNA cistrons that comprise the externally transcribed spacer regions (ETS) and part of the intergenic spacer region (IGS), including several repetitive motifs. Finally, we report the occurrence of retrotransposition in the nuclear large subunit rDNA, as R2 elements are present in the usual insertion points found in other arthropods. Interestingly, functional R1 elements usually present in the genomes of insects were not detected in the honey bee rRNA genes. The reverse transcriptase products of the R2 elements are deduced from their putative open reading frames and structurally aligned with those from another hymenopteran insect, the jewel wasp Nasonia (Pteromalidae). Stretches of conserved amino acids shared between Apis and Nasonia are illustrated and serve as potential sites for primer design, as target amplicons within these R2 elements may serve as novel phylogenetic markers for Hymenoptera. Given the impending completion of the sequencing of the Nasonia genome, we expect our report eventually to shed light on the evolution of the hymenopteran genome within higher insects, particularly regarding the relative maintenance of conserved rDNA genes, related variable spacer regions and retrotransposable elements.

Chlamydomonas pitschmannii Ettl, a little known species from thermoacidic environments

Three Chlamydomonas strains were isolated from the soils of a hot spring located in the Campi Flegrei Caldera (Naples, Italy). Ecophysiological, morpho-cytological and molecular features were used to characterize these isolates and to compare them with chlamydomonax acidophila strains from algal culture collections. The strains were collected from three points of the volcanic site, differing in their physico-chemical conditions. Among the examined Chlamydomonas strains, only the isolates from Campi Flegrei could grow optimally at pH values < or =3.0. These isolates also showed a high tolerance to desiccation and high temperatures, not evidenced by the other Chlamydomonas strains included in the study. 18S rDNA phylogeny indicates that the isolates from Campi Flegrei are closely related to Chlamydomonas pitschmannii and two strains isolated in Canada and Europe, that have been designated as Chlamydomonas acidophila. A Chlamydomonas acidophila strain isolated from the type locality in Japan is less closely related according to its molecular phylogeny, and can also be discerned by light and electron microscopy. Moreover, vegetative cells and sporangia of Chlamydomonas acidophila from Japan showed a median trilaminar structure not observed in the other strains. Our results show that Chlamydomonas pitschmannii could represent a hitherto unknown extremophilic Chlamydomonas species.

New Rapid Polymerase Chain Reaction-Immunochromatographic Assay forPorphyromonas gingivalis

BACKGROUND

A simple and rapid method for Porphyromonas gingivalis detection in clinical samples has been developed using polymerase chain reaction (PCR) and an immunochromatographic assay (ICA) with a lateral-flow device (strip) to detect species-specific 16S rRNA genes.

METHODS

The PCR used a pair of primer sets labeled with fluorescein isothiocyanate (FITC) or biotin at each 5' terminus. The strip used a nitrocellulose membrane containing streptavidin conjugated to gold particles and anti-FITC line.

RESULTS

PCR and ICA detected as few as 1 and 10 cells of P. gingivalis, respectively. ICA required 5 to 10 minutes more than the initial PCR. The amplifications were not observed in other oral black-pigmented bacteria at concentrations of 10(6) colony forming unit (CFU). The ICA strips showed bands at more than 10(4) CFU/ml equivalents in clinical samples from periodontitis.

CONCLUSIONS

A diagnostic assay based on PCR-ICA was developed for the detection of P. gingivalis, and results were obtained visually in 3 hours. PCR-ICA will be a valuable tool for the rapid detection of target bacteria by chair side.

Autoclave method for rapid preparation of bacterial PCR-template DNA

An autoclave method for preparing bacterial DNA for PCR template is presented, it eliminates the use of detergents, organic solvents, and mechanical cellular disruption approaches, thereby significantly reducing processing time and costs while increasing reproducibility. Bacteria are lysed by rapid heating and depressurization in an autoclave. The lysate, cleared by microcentrifugation, was either used directly in the PCR reaction, or concentrated by ultrafiltration. This approach was compared with seven established methods of DNA template preparation from four bacterial sources which included boiling Triton X-100 and SDS, bead beating, lysozyme/proteinase K, and CTAB lysis method components. Bacteria examined were Enterococcus and Escherichia coli, a natural marine bacterial community and an Antarctic cyanobacterial-mat. DNAs were tested for their suitability as PCR templates by repetitive element random amplified polymorphic DNA (RAPD) and denaturing gradient gel electrophoresis (DGGE) analysis. The autoclave method produced PCR amplifiable template comparable or superior to the other methods, with greater reproducibility, much shorter processing time, and at a significantly lower cost.

Secondary structure of two regions in expansion segments ES3 and ES6 with the potential of forming a tertiary interaction in eukaryotic 40S ribosomal subunits

The 18S rRNA of the small eukaryotic ribosomal subunit contains several expansion segments. Electron microscopy data indicate that two of the largest expansion segments are juxtaposed in intact 40S subunits, and data from phylogenetic sequence comparisons indicate that these two expansion segments contain complementary sequences that could form a direct tertiary interaction on the ribosome. We have investigated the secondary structure of the two expansion segments in the region around the putative tertiary interaction. Ribosomes from yeast, wheat, and mouse-three organisms representing separate eukaryotic kingdoms-were isolated, and the structure of ES3 and part of the ES6 region were analyzed using the single-strand-specific chemical reagents CMCT and DMS and the double-strand-specific ribonuclease V1. The modification patterns were analyzed by primer extension and gel electrophoresis on an ABI 377 automated DNA sequencer. The investigated sequences were relatively exposed to chemical and enzymatic modification. This is in line with their indicated location on the surface at the solvent side of the subunit. The complementary ES3 and ES6 sequences were clearly inaccessible to single-strand modification, but available for cleavage by double-strand-specific RNase V1. The results are compatible with a direct helical interaction between bases in ES3 and ES6. Almost identical results were obtained with ribosomes from the three organisms investigated.

A structural and primary sequence comparison of the viral RNA-dependent RNA polymerases

A systematic bioinformatic approach to identifying the evolutionarily conserved regions of proteins has verified the universality of a newly described conserved motif in RNA-dependent RNA polymerases (motif F). In combination with structural comparisons, this approach has defined two regions that may be involved in unwinding double-stranded RNA (dsRNA) for transcription. One of these is the N-terminal portion of motif F and the second is a large insertion in motif F present in the RNA-dependent RNA polymerases of some dsRNA viruses.

A possible tertiary rRNA interaction between expansion segments ES3 and ES6 in eukaryotic 40S ribosomal subunits

Eukaryotic 16S-like ribosomal RNAs contain 12 so-called expansion segments, i.e., sequences not included in the RNA secondary structure core common to eubacteria, archaea, and eukarya. Two of these expansion segments, ES3 and ES6, are juxtaposed in the recent three-dimensional model of the eukaryotic 40S ribosomal subunit. We have analyzed ES3 and ES6 sequences from more than 2900 discrete eukaryotic species, for possible sequence complementarity between the two expansion segments. The data show that ES3 and ES6 could interact by forming a helix consisting of seven to nine contiguous base pairs in almost all analyzed species. We, therefore, suggest that ES3 and ES6 form a direct RNA-RNA contact in the ribosome.

Comparative analysis of more than 3000 sequences reveals the existence of two pseudoknots in area V4 of eukaryotic small subunit ribosomal RNA

The secondary structure of V4, the largest variable area of eukaryotic small subunit ribosomal RNA, was re-examined by comparative analysis of 3253 nucleotide sequences distributed over the animal, plant and fungal kingdoms and a diverse set of protist taxa. An extensive search for compensating base pair substitutions and for base covariation revealed that in most eukaryotes the secondary structure of the area consists of 11 helices and includes two pseudoknots. In one of the pseudoknots, exchange of base pairs between the two stems seems to occur, and covariation analysis points to the presence of a base triple. The area also contains three potential insertion points where additional hairpins or branched structures are present in a number of taxa scattered throughout the eukaryotic domain.

Coupled nucleotide covariations reveal dynamic RNA interaction patterns

Evolutionarily conserved structures in related RNA molecules contain coordinated variations (covariations) of paired nucleotides. Analysis of covariations is a very powerful approach to deduce phylogenetically conserved (i.e., functional) conformations, including tertiary interactions. Here we discuss conserved RNA folding pathways that are revealed by covariation patterns. In such pathways, structural requirements for alternative pairings cause some nucleotides to covary with two different partners. Such "coupled" covariations between three or more nucleotides were found in various types of RNAs. The analysis of coupled covariations can unravel important features of RNA folding dynamics and improve phylogeny reconstruction in some cases. Importantly, it is necessary to distinguish between multiple covariations determined by mutually exclusive structures and those determined by tertiary contacts.

Hyaloraphidium curvatum is not a green alga, but a lower fungus; Amoebidium parasiticum is not a fungus, but a member of the DRIPs

The unicellular heterotrophic protist Hyaloraphidium is classified with a family of green algae, the Ankistrodesmaceae. The only species that exists in pure culture and that is available for taxonomic studies is H. curvatum. Comparative 18S ribosomal RNA sequence analyses showed that H. curvatum belongs to the fungi rather than to the algae. Within the fungi, H. curvatum preferentially clustered with Chytridiomycetes. Unlike Chytridiomycetes, H. curvatum propagates by autosporulation, and the presence of flagella has never been reported. Transmission electron microscopy indicated that H. curvatum in some respects resembles Chytridiomycetes, but no elements of a flagellar apparatus were detected. The habitus of H. curvatum is unlike that of other fungi except the trichomycete Amoebidium parasiticum. The cell wall sugar composition of H. curvatum was unique, but to some extent resembled that of A. parasiticum. However, H. curvatum and A. parasiticum are not closely related to each other according to 18S rRNA sequence data. Moreover, A. parasiticum clustered with protistan animals, the Mesomycetozoa (DRIPs). Combined molecular, ultrastructural and chemical data do not allow assignment of H. curvatum to any recognized clade of fungi. This suggests that H. curvatum may represent an independent evolutionary lineage within the fungi.

Cirripede phylogeny using a novel approach: molecular morphometrics

We present a new method using nucleic acid secondary structure to assess phylogenetic relationships among species. In this method, which we term "molecular morphometrics," the measurable structural parameters of the molecules (geometrical features, bond energies, base composition, etc.) are used as specific characters to construct a phylogenetic tree. This method relies both on traditional morphological comparison and on molecular sequence comparison. Applied to the phylogenetic analysis of Cirripedia, molecular morphometrics supports the most recent morphological analyses arguing for the monophyly of Cirripedia sensu stricto (Thoracica + Rhizocephala + Acrothoracica). As a proof, a classical multiple alignment was also performed, either using or not using the structural information to realign the sequence segments considered in the molecular morphometrics analysis. These methods yielded the same tree topology as the direct use of structural characters as a phylogenetic signal. By taking into account the secondary structure of nucleic acids, the new method allows investigators to use the regions in which multiple alignments are barely reliable because of a large number of insertions and deletions. It thus appears to be complementary to classical primary sequence analysis in phylogenetic studies.

Phylogenetic relationships among mytilidae (Bivalvia): 18S rRNA data suggest convergence in mytilid body plans

Nearly complete sequences were determined for small-subunit (18S) rRNA genes from seven species representative of four subfamilies of Mytilidae: Modiolus modiolus and M. auriculatus (Modiolinae); Lithophaga lithophaga and L. nigra (Lithophaginae); Musculus senhousie and M. discors (Crenellinae); and Hormomya domingensis (Mytilinae). Small-subunit rRNA gene sequences were also determined for Solemya reidi (Subclass Protobranchia), Mya arenaria (Subclass Heterodonta), and Elliptio complanata (Subclass Paleoheterodonta) as outgroup taxa. Phylogenetic analyses including these and other nearly complete bivalve small-subunit rRNA sequences demonstrate support for the monophyly of the family Mytilidae and the subfamilies Crenellinae and Lithophaginae. However, the subfamilies Mytilinae and Modiolinae appear polyphyletic. Likelihood, parsimony, and distance analyses support the placement of H. domingensis (Mytilinae) in a clade with G. demissa (Modiolinae). This clade is distinct from those containing other species traditionally assigned to these two subfamilies. Kishino-Hasegawa tests support these nontraditional relationships, suggesting that the mytiliform and/or modioliform body plans have evolved independently in at least two mytilid lineages.

Long branch attraction effects and the status of "basal eukaryotes": phylogeny and structural analysis of the ribosomal RNA gene cluster of the free-living diplomonad Trepomonas agilis

The three taxa emerging at the base of the eukaryotic ribosomal RNA phylogenetic tree (Diplomonadida, Microspora, and Parabasalia) include a wide array of parasitic species. and some free-living organisms that appear to be derived from a parasitic ancestry. The basal position of these taxa, which lack mitochondria, has recently been questioned. I sequenced most of the ribosomal RNA gene cluster of the free-living diplomonad Trepomonas agilis and a secondary structure model was reconstructed for the SSU rRNA. I conducted a RASA matrix analysis to identify, independently from tree reconstruction, putative long branch attraction effects in the data matrix. The results show that each of the basal clades and the euglenozoan clade act, indeed, as long branches and may have been engaged in a process of accelerated rate of evolution. A nucleotide signature analysis was conducted in the conserved regions for positions defining the three great domains of life (Eubacteria, Archea, and Eukaryota). For the three basal taxa, this analysis showed the presence of a significant number of different non-eukaryotic nucleotides. A precise study of the nature and location of these nucleotides led to conclusions supporting the results of the RASA analysis. Altogether, these findings suggest that the basal placement of these taxa in the SSU ribosomal RNA phylogenetic tree is artifactual, and flawed by long branch attraction effects.

Modelling the secondary structures of slippage-prone hypervariable RNA regions: the example of the tiger beetle 18S rRNA variable region V4

Variable regions within ribosomal RNAs frequently vary in length as a result of incorporating products of slippage. This makes constructing secondary structure models problematic because base homology is difficult or impossible to establish between species. Here, we model such a region by comparing the results of the MFOLD suboptimal folding algorithm for different species to identify conserved structures. Based on the reconstruction of base change on a phylogenetic tree of the species and comparison against null models of character change, we devise a statistical analysis to assess support of these structures from compensatory and semi-compensatory (i.e. G.C to G.U or A.U to G.U) mutations. As a model system we have used variable region V4 from cicindelid (tiger beetle) small subunit ribosomal RNAs (SSU rRNAs). This consists of a mixture of conserved and highly variable subregions and has been subject to extensive comparative analysis in the past. The model that results is similar to a previously described model of this variable region derived from a different set of species and contains a novel structure in the central, highly variable part. The method we describe may be useful in modelling other RNA regions that are subject to slippage.

RAGA: RNA sequence alignment by genetic algorithm

We describe a new approach for accurately aligning two homologous RNA sequences when the secondary structure of one of them is known. To do so we developed two software packages, called RAGA and PRAGA, which use a genetic algorithm approach to optimize the alignments. RAGA is mainly an extension of SAGA, an earlier package for multiple protein sequence alignment. In PRAGA several genetic algorithms run in parallel and exchange individual solutions. This method allows us to optimize an objective function that describes the quality of a RNA pairwise alignment, taking into account both primary and secondary structure, including pseudoknots. We report results obtained using PRAGA on nine test cases of pairs of eukaryotic small subunit rRNA sequence (nuclear and mitochondrial).

Database on the structure of small ribosomal subunit RNA

The Antwerp database on small ribosomal subunit RNA now offers more than 6000 nucleotide sequences (August 1996). All these sequences are stored in the form of an alignment based on the adopted secondary structure model, which is corroborated by the observation of compensating substitutions in the alignment. Besides the primary and secondary structure information, literature references, accession numbers and detailed taxonomic information are also compiled. For ease of use, the complete database is made available to the scientific community via World Wide Web at URL http://rrna.uia.ac.be/ssu/ .

Chemical and computer probing of RNA structure

Ribonucleic acids (RNAs) are one of the most important types of biopolymers. RNAs play key roles in the storage and multiplication of genetic information. They are important in catalysis and RNA splicing and are the most important steps of translation. This chapter describes experimental methods for probing RNA structure and theoretical methods allowing the prediction of thermodynamically favorable RNA folding. These methods are complementary and together they provide a powerful approach to determine the structure of RNAs. The three-dimensional (tertiary) structure of RNA is formed by hydrogen-bonding among functional groups of nucleosides in different regions of the molecule, by coordination of polyvalent cations, and by stacking between the double-stranded regions present in the RNA. The tertiary structures of only some small RNAs have been determined by high-resolution X-ray crystallographic analysis and nuclear magnetic resonance analysis. The most widely used approach for the investigation of RNA structure is chemical and enzymatic probing, in combination with theoretical methods and phylogenetic studies allowing the prediction of variants of RNA folding. Investigations of RNA structures with different enzymatic and chemical probes can provide detailed data allowing the identification of double-stranded regions of the molecules and nucleotides involved in tertiary interactions.

Database on the structure of small ribosomal subunit RNA

The Antwerp database on small ribosomal subunit RNA offers over 4300 nucleotide sequences (August 1995). All these sequences are stored in the form of an alignment based on the adopted secondary structure model, which in turn is corroborated by the observation of compensating substitutions in the alignment. Besides the primary and secondary structure information, literature references, accession numbers and detailed taxonomic information are also compiled. The complete database is made available to the scientific community through anonymous ftp and World Wide Web(WWW).

Antisense snoRNAs: a family of nucleolar RNAs with long complementarities to rRNA

A growing subset of small nucleolar RNAs (snoRNAs) contains long stretches of sequence complementarity to conserved sequences in mature ribosomal RNA (rRNA). This article reviews current knowledge about these complementarities and proposes that these antisense snoRNAs might function in pre-rRNA folding, base modification and ribosomal ribonucleoprotein assembly, in some cases acting as RNA chaperones.

The origin of land plants: phylogenetic relationships among charophytes, bryophytes, and vascular plants inferred from complete small-subunit ribosomal RNA gene sequences

Complete nuclear-encoded small-subunit 18S rRNA (= SSU rRNA) gene sequences were determined for the prasinophyte green alga Mantoniella squamata; the charophycean green algae Chara foetida, Coleochaete scutata, Klebsormidium flaccidum, and Mougeotia scalaris; the bryophytes Marchantia polymorpha, Fossombronia pusilla, and Funaria hygrometrica; and the lycopod Selaginella galleottii to get a better insight into the sequential evolution from green algae to land plants. The sequences were aligned with several previously published SSU rRNA sequences from chlorophytic and charophytic algae as well as from land plants to infer the evolutionary relationships for major evolutionary lineages within the Chlorobionta by distance matrix, maximum parsimony, and maximum likelihood analyses. Phylogenetic trees created by the different methods consistently placed the Charophyceae on the branch leading to the land plants. The Charophyceae were shown to be polyphyletic with the Charales ("charalean" algae) diverging earlier than the Coleochaetales, Klebsormidiales, Chlorokybales, and Zygnematales ("charophycean" algae) which branch from a point closer to the land plants in most analyses. Maximum parsimony and maximum likelihood analyses imply a successive evolution from "charophycean" algae, particularly Coleochaetales, to bryophytes, lycopods, and seed plants. In contrast, distance matrix methods group the bryophytes together with the "charophycean" algae, suggesting a separate evolution of these organisms compared with the club moss and the seed plants.

Probing the structure of mouse Ehrlich ascites cell 5.8S, 18S and 28S ribosomal RNA in situ

The secondary structure of mouse Ehrlich ascites 18S, 5.8S and 28S ribosomal RNA in situ was investigated by chemical modification using dimethyl sulphate and 1-cyclohexyl-3-(morpholinoethyl) carbodiimide metho-p-toluene sulphonate. These reagents specifically modify unpaired bases in the RNA. The reactive bases were localized by primer extension followed by gel electrophoresis. The three rRNA species were equally accessible for modification i.e. approximately 10% of the nucleotides were reactive. The experimental data support the theoretical secondary structure models proposed for 18S and 5.8/28S rRNA as almost all modified bases were located in putative single-strand regions of the rRNAs or in helical regions that could be expected to undergo dynamic breathing. However, deviations from the suggested models were found in both 18S and 28S rRNA. In 18S rRNA some putative helices in the 5'-domain were extensively modified by the single-strand specific reagents as was one of the suggested helices in domain III of 28S rRNA. Of the four eukaryote specific expansion segments present in mouse Ehrlich ascites cell 28S rRNA, segments I and III were only partly available for modification while segments II and IV showed average to high modification.

Mosquito large subunit ribosomal RNA: simultaneous alignment of primary and secondary structure

We report the sequence and propose a secondary structure for the cytoplasmic large subunit (5.8S and 28S) ribosomal RNA of the mosquito, Aedes albopictus, in an aligned format that incorporates secondary structure comparisons with Homo sapiens, Drosophila melanogaster, and Escherichia coli ribosomal RNAs. This format facilitates comparison of subtle differences between models, allowing nucleotide by nucleotide analysis at each position of discrepancy. Comparison of the A. albopictus large subunit ribosomal RNA gene with those from other species revealed new compensatory base changes. The aligned format focuses attention to the specific contribution of the A. albopictus sequence by facilitating comparison with the sequence of another dipteran, D. melanogaster. This is the second report of a complete large subunit rRNA sequence from an arthropod, and the first 28S rRNA sequence for a member of the lower Diptera (Nematocera).

Reconstructing evolution from eukaryotic small-ribosomal-subunit RNA sequences: calibration of the molecular clock

The detailed descriptions now available for the secondary structure of small-ribosomal-subunit RNA, including areas of highly variable primary structure, facilitate the alignment of nucleotide sequences. However, for optimal exploitation of the information contained in the alignment, a method must be available that takes into account the local sequence variability in the computation of evolutionary distance. A quantitative definition for the variability of an alignment position is proposed in this study. It is a parameter in an equation which expresses the probability that the alignment position contains a different nucleotide in two sequences, as a function of the distance separating these sequences, i.e., the number of substitutions per nucleotide that occurred during their divergence. This parameter can be estimated from the distance matrix resulting from the conversion of pairwise sequence dissimilarities into pairwise distances. Alignment positions can then be subdivided into a number of sets of matching variability, and the average variability of each set can be derived. Next, the conversion of dissimilarity into distance can be recalculated for each set of alignment positions separately, using a modified version of the equation that corrects for multiple substitutions and changing for each set the parameter that reflects its average variability. The distances computed for each set are finally averaged, giving a more precise distance estimation. Trees constructed by the algorithm based on variability calibration have a topology markedly different from that of trees constructed from the same alignments in the absence of calibration. This is illustrated by means of trees constructed from small-ribosomal-subunit RNA sequences of Metazoa. A reconstruction of vertebrate evolution based on calibrated alignments matches the consensus view of paleontologists, contrary to trees based on uncalibrated alignments. In trees derived from sequences covering several metazoan phyla, artefacts in topology that are probably due to a high clock rate in certain lineages are avoided.

Structural arrangement of the codon-anticodon interaction area in human placenta ribosomes. Affinity labelling of the 40S subunits by derivatives of oligoribonucleotides containing the AUG codon

Using the derivatives of the oligoribonucleotides pAUGUn and AUGUnC (n = 0; 3) bearing an alkylating group at either the 5' or 3' end, respectively (mRNA analogues), the affinity labelling of the human placenta 40S ribosomal subunits has been investigated in model initiation complexes obtained in the presence of the ternary complex eIF-2.GTP.Met-tRNA(fMet). The regions of 18S rRNA and ribosomal proteins labelled with these mRNA analogues were identified. The sites of covalent attachment of the pAUGUn derivatives with a reactive group at the 5' end were located between 18S rRNA positions 976 and 1164. The derivative of AUGU3C with an alkylating group at the 3' end modified 18S rRNA mainly at the 593-673 region. The main targets of the 3' end derivative of AUGC were located between positions 1610 and 1869. The proteins S3/S3a, S6, S7 and S14/S15 were modified by both types of the oligoribonucleotide derivatives regardless of the point of the reactive group attachment to the oligonucleotide moiety. The proteins S2 and S4 were modified by both the 3' end derivative of AUGC and 5' end derivative of pAUGU3; and the protein S8 was modified by the 3' end derivative of AUGC. The proteins S5 and S9 were labelled by the 5' end derivative of pAUGU3, and the protein S17 was modified by the 5' end derivative of pAUG.