The Value of Next-Generation Sequencing in Diagnosis and Therapy of Critically Ill Patients with Suspected Bloodstream Infections: A Retrospective Cohort Study

Bloodstream infection (BSI), a frequent cause of severe sepsis, is a life-threatening complication in critically ill patients and still associated with a high mortality rate. Rapid pathogen identification from blood is crucial for an early diagnosis and the treatment of patients with suspected BSI. For this purpose, novel diagnostic tools on the base of genetic analysis have emerged for clinical application. The aim of this study was to assess the diagnostic value of additional next-generation sequencing (NGS) pathogen test for patients with suspected BSI in a surgical ICU and its potential impact on antimicrobial therapy. In this retrospective single-centre study, clinical data and results from blood culture (BC) and NGS pathogen diagnostics were analysed for ICU patients with suspected BSI. Consecutive changes in antimicrobial therapy and diagnostic procedures were evaluated. Results: 41 cases with simultaneous NGS and BC sampling were assessed. NGS showed a statistically non-significant higher positivity rate than BC (NGS: 58.5% (24/41 samples) vs. BC: 21.9% (9/41); p = 0.056). NGS detected eight different potentially relevant bacterial species, one fungus and six different viruses, whereas BC detected four different bacterial species and one fungus. NGS results affected antimicrobial treatment in 7.3% of cases. Conclusions: NGS-based diagnostics have the potential to offer a higher positivity rate than conventional culture-based methods in patients with suspected BSI. Regarding the high cost, their impact on anti-infective therapy is currently limited. Larger randomized prospective clinical multicentre studies are required to assess the clinical benefit of this novel diagnostic technology.

Effectiveness of CPR in Hypogravity Conditions-A Systematic Review

(1) Background: Cardiopulmonary resuscitation (CPR), as a form of basic life support, is critical for maintaining cardiac and cerebral perfusion during cardiac arrest, a medical condition with high expected mortality. Current guidelines emphasize the importance of rapid recognition and prompt initiation of high-quality CPR, including appropriate cardiac compression depth and rate. As space agencies plan missions to the Moon or even to explore Mars, the duration of missions will increase and with it the chance of life-threatening conditions requiring CPR. The objective of this review was to examine the effectiveness and feasibility of chest compressions as part of CPR following current terrestrial guidelines under hypogravity conditions such as those encountered on planetary or lunar surfaces; (2) Methods: A systematic literature search was conducted by two independent reviewers (PubMed, Cochrane Register of Controlled Trials, ResearchGate, National Aeronautics and Space Administration (NASA)). Only controlled trials conducting CPR following guidelines from 2010 and after with advised compression depths of 50 mm and above were included; (3) Results: Four different publications were identified. All studies examined CPR feasibility in 0.38 G simulating the gravitational force on Mars. Two studies also simulated hypogravity on the Moon with a force of 0.17 G/0,16 G. All CPR protocols consisted of chest compressions only without ventilation. A compression rate above 100/s could be maintained in all studies and hypogravity conditions. Two studies showed a significant reduction of compression depth in 0.38 G (-7.2 mm/-8.71 mm) and 0.17 G (-12.6 mm/-9.85 mm), respectively, with nearly similar heart rates, compared to 1 G conditions. In the other two studies, participants with higher body weight could maintain a nearly adequate mean depth while effort measured by heart rate (+23/+13.85 bpm) and VO_2max (+5.4 mL·kg^-1·min^-1) increased significantly; (4) Conclusions: Adequate CPR quality in hypogravity can only be achieved under increased physical stress to compensate for functional weight loss. Without this extra effort, the depth of compression quickly falls below the guideline level, especially for light-weight rescuers. This means faster fatigue during resuscitation and the need for more frequent changes of the resuscitator than advised in terrestrial guidelines. Alternative techniques in the straddling position should be further investigated in hypogravity.

Experimental determination and system level analysis of essential genes in Escherichia coli MG1655

Defining the gene products that play an essential role in an organism's functional repertoire is vital to understanding the system level organization of living cells. We used a genetic footprinting technique for a genome-wide assessment of genes required for robust aerobic growth of Escherichia coli in rich media. We identified 620 genes as essential and 3,126 genes as dispensable for growth under these conditions. Functional context analysis of these data allows individual functional assignments to be refined. Evolutionary context analysis demonstrates a significant tendency of essential E. coli genes to be preserved throughout the bacterial kingdom. Projection of these data over metabolic subsystems reveals topologic modules with essential and evolutionarily preserved enzymes with reduced capacity for error tolerance.

Use of contiguity on the chromosome to predict functional coupling

The availability of a growing number of completely sequenced genomes opens new opportunities for understanding of complex biological systems. Success of genome-based biology will, to a large extent, depend on the development of new approaches and tools for efficient comparative analysis of the genomes and their organization. We have developed a technique for detecting possible functional coupling between genes based on detection of potential operons. The approach involves computation of "pairs of close bidirectional best hits", which are pairs of genes that apparently occur within operons in multiple genomes. Using these pairs, one can compose evidence (based on the number of distinct genomes and the phylogenetic distance between the orthologous pairs) that a pair of genes is potentially functionally coupled. The technique has revealed a surprisingly rich and apparently accurate set of functionally coupled genes. The approach depends on the use of a relatively large number of genomes, and the amount of detected coupling grows dramatically as the number of genomes increases.

Archaeal shikimate kinase, a new member of the GHMP-kinase family

Shikimate kinase (EC 2.7.1.71) is a committed enzyme in the seven-step biosynthesis of chorismate, a major precursor of aromatic amino acids and many other aromatic compounds. Genes for all enzymes of the chorismate pathway except shikimate kinase are found in archaeal genomes by sequence homology to their bacterial counterparts. In this study, a conserved archaeal gene (gi1500322 in Methanococcus jannaschii) was identified as the best candidate for the missing shikimate kinase gene by the analysis of chromosomal clustering of chorismate biosynthetic genes. The encoded hypothetical protein, with no sequence similarity to bacterial and eukaryotic shikimate kinases, is distantly related to homoserine kinases (EC 2.7.1.39) of the GHMP-kinase superfamily. The latter functionality in M. jannaschii is assigned to another gene (gi591748), in agreement with sequence similarity and chromosomal clustering analysis. Both archaeal proteins, overexpressed in Escherichia coli and purified to homogeneity, displayed activity of the predicted type, with steady-state kinetic parameters similar to those of the corresponding bacterial kinases: K(m,shikimate) = 414 +/- 33 microM, K(m,ATP) = 48 +/- 4 microM, and k(cat) = 57 +/- 2 s(-1) for the predicted shikimate kinase and K(m,homoserine) = 188 +/- 37 microM, K(m,ATP) = 101 +/- 7 microM, and k(cat) = 28 +/- 1 s(-1) for the homoserine kinase. No overlapping activity could be detected between shikimate kinase and homoserine kinase, both revealing a >1,000-fold preference for their own specific substrates. The case of archaeal shikimate kinase illustrates the efficacy of techniques based on reconstruction of metabolism from genomic data and analysis of gene clustering on chromosomes in finding missing genes.

Analysis of the Thermotoga maritima genome combining a variety of sequence similarity and genome context tools

The proliferation of genome sequence data has led to the development of a number of tools and strategies that facilitate computational analysis. These methods include the identification of motif patterns, membership of the query sequences in family databases, metabolic pathway involvement and gene proximity. We re-examined the completely sequenced genome of Thermotoga maritima by employing the combined use of the above methods. By analyzing all 1877 proteins encoded in this genome, we identified 193 cases of conflicting annotations (10%), of which 164 are new function predictions and 29 are amendments of previously proposed assignments. These results suggest that the combined use of existing computational tools can resolve inconclusive sequence similarities and significantly improve the prediction of protein function from genome sequence.

An archaeal genomic signature

Comparisons of complete genome sequences allow the most objective and comprehensive descriptions possible of a lineage's evolution. This communication uses the completed genomes from four major euryarchaeal taxa to define a genomic signature for the Euryarchaeota and, by extension, the Archaea as a whole. The signature is defined in terms of the set of protein-encoding genes found in at least two diverse members of the euryarchaeal taxa that function uniquely within the Archaea; most signature proteins have no recognizable bacterial or eukaryal homologs. By this definition, 351 clusters of signature proteins have been identified. Functions of most proteins in this signature set are currently unknown. At least 70% of the clusters that contain proteins from all the euryarchaeal genomes also have crenarchaeal homologs. This conservative set, which appears refractory to horizontal gene transfer to the Bacteria or the Eukarya, would seem to reflect the significant innovations that were unique and fundamental to the archaeal "design fabric." Genomic protein signature analysis methods may be extended to characterize the evolution of any phylogenetically defined lineage. The complete set of protein clusters for the archaeal genomic signature is presented as supplementary material (see the PNAS web site, www.pnas.org).

Functional analysis of gapped microbial genomes: amino acid metabolism of Thiobacillus ferrooxidans

A gapped genome sequence of the biomining bacterium Thiobacillus ferrooxidans strain ATCC23270 was assembled from sheared DNA fragments (3.2-times coverage) into 1,912 contigs. A total of 2,712 potential genes (ORFs) were identified in 2.6 Mbp (megabase pairs) of Thiobacillus genomic sequence. Of these genes, 2,159 could be assigned functions by using the WIT-Pro/EMP genome analysis system, most with a high degree of certainty. Nine hundred of the genes have been assigned roles in metabolic pathways, producing an overview of cellular biosynthesis, bioenergetics, and catabolism. Sequence similarities, relative gene positions on the chromosome, and metabolic reconstruction (placement of gene products in metabolic pathways) were all used to aid gene assignments and for development of a functional overview. Amino acid biosynthesis was chosen to demonstrate the analytical capabilities of this approach. Only 10 expected enzymatic activities, of the nearly 150 involved in the biosynthesis of all 20 amino acids, are currently unassigned in the Thiobacillus genome. This result compares favorably with 10 missing genes for amino acid biosynthesis in the complete Escherichia coli genome. Gapped genome analysis can therefore give a decent picture of the central metabolism of a microorganism, equivalent to that of a complete sequence, at significantly lower cost.

WIT: integrated system for high-throughput genome sequence analysis and metabolic reconstruction

The WIT (What Is There) (http://wit.mcs.anl.gov/WIT2/) system has been designed to support comparative analysis of sequenced genomes and to generate metabolic reconstructions based on chromosomal sequences and metabolic modules from the EMP/MPW family of databases. This system contains data derived from about 40 completed or nearly completed genomes. Sequence homologies, various ORF-clustering algorithms, relative gene positions on the chromosome and placement of gene products in metabolic pathways (metabolic reconstruction) can be used for the assignment of gene functions and for development of overviews of genomes within WIT. The integration of a large number of phylogenetically diverse genomes in WIT facilitates the understanding of the physiology of different organisms.

Universal protein families and the functional content of the last universal common ancestor

The phylogenetic distribution of Methanococcus jannaschii proteins can provide, for the first time, an estimate of the genome content of the last common ancestor of the three domains of life. Relying on annotation and comparison with reference to the species distribution of sequence similarities results in 324 proteins forming the universal family set. This set is very well characterized and relatively small and nonredundant, containing 301 biochemical functions, of which 246 are unique. This universal function set contains mostly genes coding for energy metabolism or information processing. It appears that the Last Universal Common Ancestor was an organism with metabolic networks and genetic machinery similar to those of extant unicellular organisms.

The use of gene clusters to infer functional coupling

Previously, we presented evidence that it is possible to predict functional coupling between genes based on conservation of gene clusters between genomes. With the rapid increase in the availability of prokaryotic sequence data, it has become possible to verify and apply the technique. In this paper, we extend our characterization of the parameters that determine the utility of the approach, and we generalize the approach in a way that supports detection of common classes of functionally coupled genes (e.g., transport and signal transduction clusters). Now that the analysis includes over 30 complete or nearly complete genomes, it has become clear that this approach will play a significant role in supporting efforts to assign functionality to the remaining uncharacterized genes in sequenced genomes.

A new version of the RDP (Ribosomal Database Project)

The Ribosomal Database Project (RDP-II), previously described by Maidak et al. [ Nucleic Acids Res. (1997), 25, 109-111], is now hosted by the Center for Microbial Ecology at Michigan State University. RDP-II is a curated database that offers ribosomal RNA (rRNA) nucleotide sequence data in aligned and unaligned forms, analysis services, and associated computer programs. During the past two years, data alignments have been updated and now include >9700 small subunit rRNA sequences. The recent development of an ObjectStore database will provide more rapid updating of data, better data accuracy and increased user access. RDP-II includes phylogenetically ordered alignments of rRNA sequences, derived phylogenetic trees, rRNA secondary structure diagrams, and various software programs for handling, analyzing and displaying alignments and trees. The data are available via anonymous ftp (ftp.cme.msu. edu) and WWW (http://www.cme.msu.edu/RDP). The WWW server provides ribosomal probe checking, approximate phylogenetic placement of user-submitted sequences, screening for possible chimeric rRNA sequences, automated alignment, and a suggested placement of an unknown sequence on an existing phylogenetic tree. Additional utilities also exist at RDP-II, including distance matrix, T-RFLP, and a Java-based viewer of the phylogenetic trees that can be used to create subtrees.

The complete genome of the hyperthermophilic bacterium Aquifex aeolicus

Aquifex aeolicus was one of the earliest diverging, and is one of the most thermophilic, bacteria known. It can grow on hydrogen, oxygen, carbon dioxide, and mineral salts. The complex metabolic machinery needed for A. aeolicus to function as a chemolithoautotroph (an organism which uses an inorganic carbon source for biosynthesis and an inorganic chemical energy source) is encoded within a genome that is only one-third the size of the E. coli genome. Metabolic flexibility seems to be reduced as a result of the limited genome size. The use of oxygen (albeit at very low concentrations) as an electron acceptor is allowed by the presence of a complex respiratory apparatus. Although this organism grows at 95 degrees C, the extreme thermal limit of the Bacteria, only a few specific indications of thermophily are apparent from the genome. Here we describe the complete genome sequence of 1,551,335 base pairs of this evolutionarily and physiologically interesting organism.

The complete genome sequence of the hyperthermophilic, sulphate-reducing archaeon Archaeoglobus fulgidus

Archaeoglobus fulgidus is the first sulphur-metabolizing organism to have its genome sequence determined. Its genome of 2,178,400 base pairs contains 2,436 open reading frames (ORFs). The information processing systems and the biosynthetic pathways for essential components (nucleotides, amino acids and cofactors) have extensive correlation with their counterparts in the archaeon Methanococcus jannaschii. The genomes of these two Archaea indicate dramatic differences in the way these organisms sense their environment, perform regulatory and transport functions, and gain energy. In contrast to M. jannaschii, A. fulgidus has fewer restriction-modification systems, and none of its genes appears to contain inteins. A quarter (651 ORFs) of the A. fulgidus genome encodes functionally uncharacterized yet conserved proteins, two-thirds of which are shared with M. jannaschii (428 ORFs). Another quarter of the genome encodes new proteins indicating substantial archaeal gene diversity.

A reconstruction of the metabolism of Methanococcus jannaschii from sequence data

The interpretation of the Methanococcus jannaschii genome will inevitably require many years of effort. This initial attempt to connect the sequence data to aspects of known biochemistry and to provide an overview of what is already apparent from the sequence data will be refined. Numerous issues remain that can be resolved only by direct biochemical analysis. Let us draw the reader's attention to just a few that might be considered central: (1) We are still missing key enzymes from the glycolytic pathway, and the conjecture is that this is due to ADP-dependency. The existence of glycolytic activity in the cell-free extract should be tested. (2) The issue of whether the Calvin cycle is present needs to be examined. (3) We need to determine whether the 2-oxoglutarate synthase (ferredoxin-dependent) (EC 1.2.7.3) activity is present. (4) The issue of whether cyclic 2,3-bisphosphate is detectable in the cell-free extracts needs to be checked. If it is, this result would confirm our assertion of the two pathways controlling synthesis and degradation of cyclic 2,3-bisphosphate.

The metabolic pathway collection: an update

The Metabolic Pathway Collection from EMP is an extraction of data from the larger Enzymes and Metabolic Pathways database (EMP). This extraction has been made publicly available in the hope that others will find it useful for a variety of purposes. The original release in October 1995 contained 1814 distinct pathways. The current collection contains 2180. Metabolic reconstructions for the first completely sequenced organisms-Haemophilus influenzae,Mycoplasma genitalium,Saccharomyces cerevisiaeandMethanococcus janaschii-are all included in the current release. All of the pathways in the collections are available as ASCII files in the form generated by the main curator, Evgeni Selkov. In addition, we are offering a more structured encoding of a subset of the collection; our initial release of this subcollection includes all of the pathways inMycoplasma genitalium, and we ultimately intend to offer the entire collection in this form as well.

The RDP (Ribosomal Database Project)

The Ribosomal Database Project (RDP) is a curated database that offers ribosome-related data, analysis services and associated computer programs. The offerings include phylogenetically ordered alignments of ribosomal RNA (rRNA) sequences, derived phylogenetic trees, rRNA secondary structure diagrams, and various software for handling, analyzing and displaying alignments and trees. The data are available via anonymous FTP (rdp.life.uiuc.edu), electronic mail (server@rdp.life.uiuc.edu), gopher (rdpgopher.life.uiuc.edu) and WWW (http://rdpwww.life.uiuc.edu/ ). The electronic mail and WWW servers provide ribosomal probe checking, approximate phylogenetic placement of user-submitted sequences, screening for possible chimeric rRNA sequences, automated alignment, and a suggested placement of an unknown sequence on an existing phylogenetic tree.

Complete genome sequence of the methanogenic archaeon, Methanococcus jannaschii

The complete 1.66-megabase pair genome sequence of an autotrophic archaeon, Methanococcus jannaschii, and its 58- and 16-kilobase pair extrachromosomal elements have been determined by whole-genome random sequencing. A total of 1738 predicted protein-coding genes were identified; however, only a minority of these (38 percent) could be assigned a putative cellular role with high confidence. Although the majority of genes related to energy production, cell division, and metabolism in M. jannaschii are most similar to those found in Bacteria, most of the genes involved in transcription, translation, and replication in M. jannaschii are more similar to those found in Eukaryotes.

The Ribosomal Database Project (RDP)

The Ribosomal Database Project (RDP) is a curated database that offers ribosome-related data, analysis services and associated computer programs. The offerings include phylogenetically ordered alignments of ribosomal RNA (rRNA) sequences, derived phylogenetic trees, rRNA secondary structure diagrams and various software for handling, analyzing and displaying alignments and trees. The data are available via anonymous ftp (rdp.life.uiuc.edu), electronic mail (server@rdp.life.uiuc.edu), gopher (rdpgopher.life.uiuc.edu) and World Wide Web (WWW)(http://rdpwww.life.uiuc.edu/). The electronic mail and WWW servers provide ribosomal probe checking, screening for possible chimeric rRNA sequences, automated alignment and approximate phylogenetic placement of user-submitted sequences on an existing phylogenetic tree.

The metabolic pathway collection from EMP: the enzymes and metabolic pathways database

The Enzymes and Metabolic Pathways database (EMP) is an encoding of the contents of over 10 000 original publications on the topics of enzymology and metabolism. This large body of information has been transformed into a queryable database. An extraction of over 1800 pictorial representations of metabolic pathways from this collection is freely available on the World Wide Web. We believe that this collection will play an important role in the interpretation of genetic sequence data, as well as offering a meaningful framework for the integration of many other forms of biological data.

The 60 kDa heat shock proteins in the hyperthermophilic archaeon Sulfolobus shibatae

One of the most abundant proteins in the hyperthermophilic archaeon Sulfolobus shibatae is the 59 kDa heat shock protein (TF55) that is believed to form a homo-oligomeric double ring complex structurally similar to the bacterial chaperonins. We discovered a second protein subunit in the S. shibatae ring complex (referred to as alpha) that is stoichiometric with TF55 (renamed beta). The gene and flanking regions of alpha were cloned and sequenced and its inferred amino acid sequence has 54.4% identity and 74.4% similarity to beta. Transcription start sites for both alpha and beta were mapped and three potential transcription regulatory regions were identified. Northern analyses of cultures shifted from normal growth temperatures (70 to 75 degrees C) to heat shock temperatures (85 to 90 degrees C) indicated that the levels of alpha and beta mRNAs increased during heat shock, but at all temperatures their relative proportions remained constant. Monitoring protein synthesis by autoradiography of total proteins from cultures pulse labeled with L(-)[35S]methionine at normal and heat shock temperatures indicated significant increases in alpha and beta synthesis during heat shock. Under extreme heat shock conditions (> or = 90 degrees C) alpha and beta appeared to be the only two proteins synthesized. The purified alpha and beta subunits combined to form high molecular mass complexes with similar mobilities on native polyacrylamide gels to the complexes isolated directly from cells. Equal proportions of the two subunits gave the greatest yield of the complex, which we refer to as a "rosettasome". It is argued that the rosettasome consists of two homo-oligomeric rings; one of alpha and the other of beta. Polyclonal antibodies against alpha and beta from S. shibatae cross-reacted with proteins of similar molecular mass in 10 out of the 17 archaeal species tested, suggesting that the two rosettasome proteins are highly conserved among the archaea. The archaeal sequences were aligned with bacterial and eukaryotic chaperonins to generate a phylogenetic tree. The tree reveals the close relationship between the archaeal rosettasomes and the eukaryotic TCP1 protein family and the distant relationship to the bacterial GroEL/HSP60 proteins.

The genetic data environment an expandable GUI for multiple sequence analysis

An X-Windows-based graphic user interface is presented which allows the seamless integration of numerous existing biomolecular programs into a single analysis environment. This environment is based on a core multiple sequence editor that is linked to external programs by a user-expandable menu system and is supported on Sun and DEC workstations. There is no limitation to the number of external functions that can be linked to the interface. The length and number of sequences that can be handled are limited only by the size of virtual memory present on the workstation. The sequence data itself is used as the reference point from which analysis is done, and scalable graphic views are supported. It is suggested that future software development utilizing this expandable, user-defined menu system and the I/O linkage of external programs will allow biologists to easily integrate expertise from disparate fields into a single environment.

The Ribosomal Database Project

The Ribosomal Database Project (RDP) is a curated database that offers ribosome-related data, analysis services, and associated computer programs. The offerings include phylogenetically ordered alignments of ribosomal RNA (rRNA) sequences, derived phylogenetic trees, rRNA secondary structure diagrams, and various software for handling, analyzing and displaying alignments and trees. The data are available via anonymous ftp (rdp.life.uiuc.edu), electronic mail (server/rdp.life.uiuc.edu) and gopher (rdpgopher.life.uiuc.edu). The electronic mail server also provides ribosomal probe checking, approximate phylogenetic placement of user-submitted sequences, screening for chimeric nature of newly sequenced rRNAs, and automated alignment.

fastDNAmL: a tool for construction of phylogenetic trees of DNA sequences using maximum likelihood

We have developed a new tool, called fastDNAml, for constructing phylogenetic trees from DNA sequences. The program can be run on a wide variety of computers ranging from Unix workstations to massively parallel systems, and is available from the Ribosomal Database Project (RDP) by anonymous FTP. Our program uses a maximum likelihood approach and is based on version 3.3 of Felsenstein's dnaml program. Several enhancements, including algorithmic changes, significantly improve performance and reduce memory usage, making it feasible to construct even very large trees. Trees containing 40-100 taxa have been easily generated, and phylogenetic estimates are possible even when hundreds of sequences exist. We are currently using the tool to construct a phylogenetic tree based on 473 small subunit rRNA sequences from prokaryotes.

The ribosomal database project

The Ribosomal Database Project (RDP) is a curated database that offers ribosome data along with related programs and services. The offerings include phylogenetically ordered alignments of ribosomal RNA (rRNA) sequences, derived phylogenetic trees, rRNA secondary structure diagrams and various software packages for handling, analyzing and displaying alignments and trees. The data are available via ftp and electronic mail. Certain analytic services are also provided by the electronic mail server.

The Ribosomal Database Project

The Ribosomal Database Project (RDP) complies ribosomal sequences and related data, and redistributes them in aligned and phylogenetically ordered form to its user community. It also offers various software packages for handling, analyzing and displaying sequences. In addition, the RDP offers (or will offer) certain analytic services. At present the project is in an intermediate stage of development.

Structure detection through automated covariance search

This paper summarizes our investigations into the computational detection of secondary and tertiary structure of ribosomal RNA. We have developed a new automated procedure that not only identifies potential secondary and tertiary structural interactions, but also provides the covariation evidence that supports the proposed bondings, and any counterevidence that can be detected in the known sequences. A small number of previously unknown higher-order structural features have been detected in individual RNA molecules (16S rRNA and 7S RNA) through the use of our automated procedure. We are systematically studying mitochondrial rRNA, seeking tertiary structure within 16S rRNA and quaternary structure between 16S and 23S rRNA. To test hypotheses suggested by an examination of our program's output, our colleagues in biology are sequencing key portions of the 23S ribosomal RNA for species in which the known 16S ribosomal RNA exhibits variation (from the dominant pattern) at the site of a proposed bonding. Our ultimate hope is that automated covariation analysis will contribute significantly to a refined picture of ribosomal structure.