Copy number of ribosomal operons in prokaryotes and its effect on phylogenic analyses

Different aspects of the presence of multiple copies of ribosomal operons in prokaryotic genomes are reviewed. Structure of prokaryotic ribosomal operons is briefly described. The available data are summarized regarding the copy number of ribosomal genes in various prokaryotic genomes, the degree of polymorphism of their individual copies, physiological and evolutionary aspects of the presence of the multiple copies of ribosomal genes. The review also considers the influence of the presence of multiple copies of ribosomal genes on the results of identification of prokaryotic isolates and of the studies of prokaryotic diversity in environmental samples based on phylogenetic analysis of 16S rRNA gene sequences.     

Conservation of ribosomal protein gene ordering in 16 complete genomes

The organization of ribosomal proteins in 16 prokaryotic genomes was studied as an example of comparative genome analyses of gene systems. Hypothetical ribosomal protein-containing operons were constructed. These operons also contained putative genes and other non-ribosomal genes. The correspondences among these genes across different organisms were clarified by sequence homology computations. In this way a cross tabulation of 70 ribosomal proteins genes was constructed. On average, these were organized into 9-14 operons in each genome. There were also 25 non-ribosomal or putative genes in these mainly ribosomal protein operons. Hence the table contains 95 genes in total. It was found that: (i) the conservation of the block of about 20 r-proteins in the L3 and L4 operons across almost the entire eubacteria and ar-chaebacteria is remarkable; (ii) some operons only belong to eubacteria or archaebacte-ria; (iii) although the ribosomal protein operons are highly conserved within domain, there are fine variat     

Conservation of ribosomal protein gene ordering in 16 complete genomes

The organization of ribosomal proteins in 16 prokaryotic genomes was studied as an example of comparative genome analyses of gene systems. Hypothetical ribosomal protein-containing operons were constructed. These operons also contained putative genes and other non-ribosomal genes. The correspondences among these genes across different organisms were clarified by sequence homology computations. In this way a cross tabulation of 70 ribosomal proteins genes was constructed. On average, these were organized into 9-14 operons in each genome. There were also 25 non-ribosomal or putative genes in these mainly ribosomal protein operons. Hence the table contains 95 genes in total. It was found that: (i) the conservation of the block of about 20 r-proteins in the L3 and L4 operons across almost the entire eubacteria and archaebacteria is remarkable; (ii) some operons only belong to eubacteria or archaebacteria; (iii) although the ribosomal protein operons are highly conserved within domain, there are fine variations in some operons across different organisms within each domain, and these variations are informative on the evolutionary relations among the organisms. This method provides a new potential for studying the origin and evolution of old species.     

Organisation of the S10, spc and alpha ribosomal protein gene clusters in prokaryotic genomes

Although it is well known that there is no long range colinearity in gene order in bacterial genomes, it is thought that there are several regions that are under strong structural constraints during evolution, in which gene order is extremely conserved. One such region is the str locus, containing the S10-spc-alpha operons. These operons contain genes coding for ribosomal proteins and for a number of housekeeping genes. We compared the organisation of these gene clusters in 111 sequenced prokaryotic genomes (99 bacterial and 12 archaeal genomes). We also compared the organisation to the phylogeny based on 16S ribosomal RNA gene sequences and the sequences of the ribosomal proteins L22, L16 and S14. Our data indicate that there is much variation in gene order and content in these gene clusters, both in bacterial as well as in archaeal genomes. Our data indicate that differential gene loss has occurred on multiple occasions during evolution. We also noted several discrepancies between phylogenetic trees based on 16S rRNA gene sequences and sequences of ribosomal proteins L16, L22 and S14, suggesting that horizontal gene transfer did play a significant role in the evolution of the S10-spc-alpha gene clusters.     

Computational approaches for the analysis of gene neighbourhoods in prokaryotic genomes

Gene order in prokaryotes is conserved to a much lesser extent than protein sequences. Only some operons, primarily those that encode physically interacting proteins, are conserved in all or most of the bacterial and archaeal genomes. Nevertheless, even the limited conservation of operon organisation that is observed provides valuable evolutionary and functional clues through multiple genome comparisons. With the rapid growth in the number and diversity of sequenced prokaryotic genomes, functional inferences for uncharacterized genes located in the same conserved gene neighborhood with well-studied genes are becoming increasingly important. In this review, we discuss various computational approaches for identification of conserved gene strings and construction of local alignments of gene orders in prokaryotic genomes.     

System to study horizontal gene exchange among microorganisms without cultivation of recipients

Ribosomal RNA genes are characterized by highly conserved sequences and are present in multiple copies in most prokaryotic chromosomes. In principle, therefore, they might serve as sites for homologous recombination between unrelated microorganisms. Plasmids containing 23S ribosomal gene sequences, from different bacteria, which had been interrupted by insertion of a kanamycin-resistance gene, were used to transform Acinetobacter sp. DSM587 (former name: Acinetobacter calcoaceticus BD413-ivl10). In all cases, homologies between the 23S rRNA genes of phylogenetically distant bacteria and Acinetobac-ter sp. DSM587 were sufficient for replacement recombination events. The integration events, resulting in inactivation of any one of the seven rrn operons of Acinetobacter sp. DSM587, had no observable influence on cell growth. These results suggest the possibility of rRNA genes serving as natural vehicles for horizontal gene transfer. They also provide the basis of a novel strategy to analyse gene transfer without selection or cultivation of recipient cells. Because of the highly conserved structure of bacterial rrn operons, recombination events subsequent to gene transfer can be readily identified by polymerase chain reaction amplification of the recombinant sequence using a universal forward primer for the 16S rRNA gene and a reverse primer specific for the integrated marker gene.     

Detection of operons

Operons are clusters of genes that are transcribed as a single message, and regulated by the same gene expression machinery. They are found primarily in prokaryotic genomes. Because genes in the same operon are likely to have related functions, identification of the operon structure is potentially useful for assigning gene function. We report the development and benchmarking of two different methods for detecting operons, based on an analysis of 42 fully sequenced prokaryotic organisms. The Gene Neighbor method (GNM) utilizes the relatively high conservation of gene order in operons, compared with genes in general. The Gene Gap Method (GGM) makes use of the relatively short gap between genes in operons compared with that otherwise found between adjacent genes. The methods have been benchmarked using KEGG pathway data and RegulonDB Escherichia coli operon data. With optimum parameters, the specificity of the GNM is 93% and the sensitivity is 70%. For the GGM, the specificity is 95% and the sensitivity is 68%. Together, the two methods have a sensitivity of 87.2%, while joint predictions have a sensitivity of 50% and a specificity of 98%. The methods are used to infer possible functions for some hypothetical genes in prokaryotic genomes. The methods have proven a useful addition to structure information in deriving protein function in a structural genomics project.     

Loss of DNA recombinational repair enzymes in the initial stages of genome degeneration

Many obligate intracellular pathogens and symbionts undergo genome degeneration during long-term association with eukaryotic hosts; however, very little is known about genome changes that occur in the initial stages of such intracellular associations. By focusing on a clade of bacteria that have recently established symbiotic associations with insect hosts, we have identified events that may contribute to the reduction and degeneration of symbiont genomes. Unlike virtually all other bacteria, the obligate symbionts of maize and rice weevils each display substantial sequence divergence between multiple copies of their rDNA genes, resulting from a reduction in the efficacy of recombinational gene conversion, coincident with the inactivation of the recombinational repair gene recF in the common ancestor of both symbionts. The maize weevil endosymbiont also lacks a functional recA, resulting in further reduction in the efficacy of gene conversion between paralogous rDNAs and in a novel IS-mediated deletion in a 23S rDNA gene. Similar events may be pervasive during the evolution of symbiosis because symbiont genomes typically lack recombinational repair genes and have reduced numbers of ribosomal operons.     

Intra- and Intergenomic Variation of Ribosomal RNA Operons in Concurrent Alteromonas macleodii Strains

Biodiversity estimates based on ribosomal operon sequence diversity rely on the premise that a sequence is characteristic of a single specific taxon or operational taxonomic unit (OTU). Here, we have studied the sequence diversity of 14 ribosomal RNA operons (rrn) contained in the genomes of two isolates (five operons in each genome) and four metagenomic fosmids, all from the same seawater sample. Complete sequencing of the isolate genomes and the fosmids establish that they represent strains of the same species, Alteromonas macleodii, with average nucleotide identity (ANI) values >97 %. Nonetheless, we observed high levels of intragenomic heterogeneity (i.e., variability between operons of a single genome) affecting multiple regions of the 16S and 23S rRNA genes as well as the internally transcribed spacer 1 (ITS-1) region. Furthermore, the ribosomal operons exhibited intergenomic heterogeneity (i.e., variability between operons located in separate genomes) in each of these regions, compounding the variability. Our data reveal the extensive heterogeneity observed in natural populations of A. macleodii at a single point in time and support the idea that distinct lineages of A. macleodii exist in the deep Mediterranean. These findings highlight the potential of rRNA fingerprinting methods to misrepresent species diversity while simultaneously failing to recognize the ecological significance of individual strains.     

Predicting regulons and their cis-regulatory motifs by comparative genomics

We have combined and compared three techniques for predicting functional interactions based on comparative genomics (methods based on conserved operons, protein fusions and correlated evolution) and optimized these methods to predict coregulated sets of genes in 24 complete genomes, including Saccharomyces cerevisiae, Caernorhabditis elegans and 22 prokaryotes. The method based on conserved operons was the most useful for this purpose. Upstream regions of the genes comprising these predicted regulons were then used to search for regulatory motifs in 22 prokaryotic genomes using the motif-discovery program AlignACE. Many significant upstream motifs, including five known Escherichia coli regulatory motifs, were identified in this manner. The presence of a significant regulatory motif was used to refine the members of the predicted regulons to generate a final set of predicted regulons that share significant regulatory elements.     

The life-cycle of operons

Operons are a major feature of all prokaryotic genomes, but how and why operon structures vary is not well understood. To elucidate the life-cycle of operons, we compared gene order between Escherichia coli K12 and its relatives and identified the recently formed and destroyed operons in E. coli. This allowed us to determine how operons form, how they become closely spaced, and how they die. Our findings suggest that operon evolution may be driven by selection on gene expression patterns. First, both operon creation and operon destruction lead to large changes in gene expression patterns. For example, the removal of lysA and ruvA from ancestral operons that contained essential genes allowed their expression to respond to lysine levels and DNA damage, respectively. Second, some operons have undergone accelerated evolution, with multiple new genes being added during a brief period. Third, although genes within operons are usually closely spaced because of a neutral bias toward deletion and because of selection against large overlaps, genes in highly expressed operons tend to be widely spaced because of regulatory fine-tuning by intervening sequences. Although operon evolution may be adaptive, it need not be optimal: new operons often comprise functionally unrelated genes that were already in proximity before the operon formed.     

Gene Conversion Drives Within Genic Sequences: Concerted Evolution of Ribosomal RNA Genes in Bacteria and Archaea

Multiple copies of a given ribosomal RNA gene family undergo concerted evolution such that sequences of all gene copies are virtually identical within a species although they diverge normally between species. In eukaryotes, gene conversion and unequal crossing over are the proposed mechanisms for concerted evolution of tandemly repeated sequences, whereas dispersed genes are homogenized by gene conversion. However, the homogenization mechanisms for multiple-copy, normally dispersed, prokaryotic rRNA genes are not well understood. Here we compared the sequences of multiple paralogous rRNA genes within a genome in 12 prokaryotic organisms that have multiple copies of the rRNA genes. Within a genome, putative sequence conversion tracts were found throughout the entire length of each individual rRNA genes and their immediate flanks. Individual conversion events convert only a short sequence tract, and the conversion partners can be any paralogous genes within the genome. Interestingly, the genic sequences undergo much slower divergence than their flanking sequences. Moreover, genomic context and operon organization do not affect rRNA gene homogenization. Thus, gene conversion underlies concerted evolution of bacterial rRNA genes, which normally occurs within genic sequences, and homogenization of flanking regions may result from co-conversion with the genic sequence. Received: 31 March 2000 / Accepted: 15 June 2000     

Phylogenomic analysis identifies a sodium-translocating decarboxylating oxidoreductase in thermotogae

Bacterial sodium-dependent decarboxylases were the first enzymes exemplifying sodium-dependent bioenergetics. These enzyme complexes couple decarboxylation of organic acids with the export of sodium ions via a special membrane subunit. In 711 representative prokaryotic genomes, we have analyzed genomic neighborhoods of the genes that code the membrane subunit of sodium decarboxylases. In representatives of Thermotogae, the operons with the gene of this subunit lack the genes of subunits that perform non-oxidative decarboxylation. Instead, these operons contain the genes of alphaand delta-subunits of decarboxylating oxidoreductases of alpha-ketoacids. The genes of betaand gamma-subunits of the decarboxylating oxidoreductases were found within the genomes of respective Thermotogae species as separate, twogene operons. We suggest that the described two operons code together for sodium-translocating decarboxylating oxidoreductases capable of coupling oxidative decarboxylation of alpha-ketoacids with the export of sodium ions, which is a novel type of bioenergetic coupling.