Intracellular Bacterial Symbionts of Aphids Possess Many Genomic Copies per Bacterium

Although Buchnera, the endosymbiotic bacteria of aphids, are close relatives of Escherichia coli, their genome size is only a seventh that of E. coli. In this study, we estimated the genomic copy number of Buchnera by dot-blot hybridization and fluorimetry using a video-intensified microscope photon-counting system and obtained convincing evidence that each cell of these bacteria contains an average of 120 genomic copies. Thus, the Buchnera symbiont, with many copies of a small-sized genome, is reminiscent of cell organelles such as mitochondria and chloroplasts. Received: 25 November 1998 / Accepted: 25 December 1998     

Accumulation of Species-Specific Amino Acid Replacements That Cause Loss of Particular Protein Functions in Buchnera, an Endocellular Bacterial Symbiont

Endosymbiotic bacteria live in animal cells and are transmitted vertically at the time of the host's reproduction. In view of their small and asexual populations with infrequent chances of recombination, these endocellular bacteria are expected to accumulate mildly deleterious mutations. Previous studies showed that the DNA sequences of these bacteria evolved faster than those of free-living bacteria. In this study, we compared all the ORFs of Buchnera, an endocellular bacterial symbiont of aphids, with those of 34 other prokaryotic organisms and estimated the effect of the accelerated evolution of Buchnera on the functions of its proteins. It was revealed that Buchnera proteins contain many mutations at the sites where sequences are conserved in their orthologues in many other organisms. In addition, amino acid replacements at the conserved sites are mostly changes to physicochemically different amino acids. These results suggest that functions and conformations of Buchnera proteins have been seriously impaired or strongly modified. Indeed, extensive loss of functional motifs was observed in some Buchnera proteins. In many Buchnera proteins mutations were not detected evenly throughout each molecule but tended to accumulate in some functional units, possibly leading to loss of specific functions. As Buchnera has an unusual and limited gene repertory, it is conceivable that the manner of interactions among its proteins has been changed, and thus, functional constraints over their amino acid residues have also been changed during evolution. This may account for the loss of some functional units only in the Buchnera proteins. We obtained evidence that amino acid replacements in Buchnera were not always deleterious, but neutral or, in some cases, even positively selected. Received: 14 December 2000 / Accepted: 12 March 2001     

Genetic Characterization of Plasmids Containing Genes Encoding Enzymes of Leucine Biosynthesis in Endosymbionts (Buchnera) of Aphids

The prokaryotic endosymbionts (Buchnera) of aphids are known to provision their hosts with amino acids that are limiting in the aphid diet. Buchnera from the aphids Schizaphis graminum and Diuraphis noxia have plasmids containing leuABCD, genes that encode enzymes of the leucine biosynthetic pathway, as well as genes encoding proteins probably involved in plasmid replication (repA1 and repA2) and an open reading frame (ORF1) of unknown function. The newly reported plasmids closely resemble a plasmid previously described in Buchnera of the aphid Rhopalosiphum padi [Bracho AM, Martínez-Torres D, Moya A, Latorre A (1995) J Mol Evol 41:67–73]. Nucleotide sequence comparisons indicate conserved regions which may correspond to an origin of replication and two promoters, as well as inverted repeats, one of which resembles a rho-independent terminator. Phylogenetic analyses based on amino acid sequences of leu gene products and ORF1 resulted in trees identical to those obtained from endosymbiont chromosomal genes and the plasmid-borne trpEG. These results are consistent with a single evolutionary origin of the leuABCD-containing plasmid in a common ancestor of Aphididae and the lack of plasmid exchange between endosymbionts of different aphid species. Trees for ORF1 and repA (based on both nucleotides and amino acids) are used to examine the basis for leu plasmid differences between Buchnera of Thelaxes suberi and Aphididae. The most plausible explanation is that a single transfer of the leu genes to an ancestral replicon was followed by rearrangements. The related replicon in Buchnera of Pemphigidae, which lacks leuABCD, appears to represent the ancestral condition, implying that the plasmid location of the leu genes arose after the Pemphigidae diverged from other aphid families. This conclusion parallels previously published observations for the unrelated trpEG plasmid, which is present in Aphididae and absent in Pemphigidae. Recruitment of amino acid biosynthetic genes to plasmids has been ongoing in Buchnera lineages after the infection of aphid hosts. Received: 9 March 1998 / Accepted: 18 May 1998.     

Construction of a Fibrobacter succinogenes Genomic Map and Demonstration of Diversity at the Genomic Level

The genomic cleavage map of the type strain Fibrobacter succinogenes S85 was constructed. The restriction enzymes AscI, AvrII, FseI, NotI, and SfiI generated DNA fragments of suitable size distribution that could be resolved by pulsed-field gel electrophoresis (PFGE). An average genome size of 3.6 Mb was obtained by summing the total fragment sizes. The linkages between the 15 AscI fragments of the genome were determined by combining two approaches: isolation of linking clones and cross-hybridization of restriction fragments. The genome of F. succinogenes was found to be represented by the single circular DNA molecule. Southern hybridization with specific probes allowed the eight genetic markers to be located on the restriction map. The genome of this bacterium contains at least three rRNA operons. PFGE of the other three strains of F. succinogenes gave estimated genome sizes close to that of the type strain. However, RFLP patterns of these strains generated by AscI digestion are completely different. Pairwise comparison of the genomic fragment distribution between the type strain and the three isolates showed a similarity level in the region of 14.3% to 31.3%. No fragment common to all of these F. succinogenes strains could be detected by PFGE. A marked degree of genomic heterogeneity among members of this species makes genomic RFLP a highly discriminatory and useful molecular typing tool for population studies. Received: 23 October 1996 / Accepted: 31 December 1996     

Limitations of compositional approach to identifying horizontally transferred genes

Genes with atypical G+C content and pattern of codon usage in a certain genome are possibly of exotic origin, and this idea has been applied to identify horizontal events. In this way, it was postulated that a total of 755 genes in the E. coli genome are relics of horizontal events after the divergence of E. coli from the Salmonella lineage 100 million years ago (Lawrence and Ochman, 1998). In this paper we propose a new way to study sequence composition more thoroughly. We found that although the 755 genes differ in composition from other genes in the E. coli genome, the difference is minor. If we accepted that these genes are horizontally transferred, then (1) it would be more likely that they were transferred from genomes evolutionarily closely related to E. coli; but (2) the dating method used by Lawrence and Ochman (1997, 1998) largely underestimated the average age of introduced sequences in the E. coli genome, in particular, most of the 755 genes should be introduced into E. coli before, instead of after, the divergence of E. coli from the Salmonella lineage. Our study reveals that atypical G+C content and pattern of codon usage are not reliable indicators of horizontal gene transfer events. Received: 27 September 2000 / Accepted: 9 April 2001     

Combined Physical and Genetic Map of the Pseudomonas putida KT2440 Chromosome

A combined physical and genetic map of the Pseudomonas putida KT2440 genome was constructed from data obtained by pulsed-field gel electrophoresis techniques (PFGE) and Southern hybridization. Circular genome size was estimated at 6.0 Mb by adding the sizes of 19 SwaI, 9 PmeI, 6 PacI, and 6 I-CeuI fragments. A complete physical map was achieved by combining the results of (i) analysis of PFGE of the DNA fragments resulting from digestion of the whole genome with PmeI, SwaI, I-CeuI, and PacI as well as double digestion with combinations of these enzymes and (ii) Southern hybridization analysis of the whole wild-type genome digested with different enzymes and hybridized against a series of probes obtained as cloned genes from different pseudomonads of rRNA group I and Escherichia coli, as P. putida DNA obtained by PCR amplification based on sequences deposited at the GenBank database, and by labeling of macrorestriction fragments of the P. putida genome eluted from agarose gels. As an alternative, 10 random mini-Tn5-Km mutants of P. putida KT2440 were used as a source of DNA, and the band carrying the mini-Tn5 in each mutant was identified after PFGE of a series of complete chromosomal digestions and hybridization with the kanamycin resistance gene of the mini-Tn5 as a probe. We established a circular genome map with an average resolution of 160 kb. Among the 63 genes located on the genetic map were key markers such as oriC, 6 rrn loci (rnnA to -F), recA, ftsZ, rpoS, rpoD, rpoN, and gyrB; auxotrophic markers; and catabolic genes for the metabolism of aromatic compounds. The genetic map of P. putida KT2440 was compared to those of Pseudomonas aeruginosa PAO1 and Pseudomonas fluorescens SBW25. The chromosomal backbone revealed some similarity in gene clustering among the three pseudomonads but differences in physical organization, probably as a result of intraspecific rearrangements.     

DNA Polymerase C of the Thermophilic Bacterium Thermus aquaticus: Classification and Phylogenetic Analysis of the Family C DNA Polymerases

Bacterial family C DNA polymerases (DNA pol IIIs), the major chromosomal replicative enzymes, have been provisionally classified based on primary sequences and domain structures into three classes: class I (Escherichia coli DNA pol C-type), class II (Bacillus subtilis DNA pol C-type), and class III (cyanobacterial DNA pol C-type), respectively. We have sequenced the structural gene encoding the DNA pol C catalytic subunit of the thermophilic bacterium Thermus aquaticus. This gene, designated the Taq DNA pol C gene, contains a 3660-bp open reading frame which specifies a polypeptide of molecular weight of 137,388 daltons. Comparative sequence analyses revealed that Taq DNA pol C is a class I family C DNA polymerase. The Taq DNA pol C is most closely related to the Deinococcus radiodurans DNA pol C. Although a phylogenetic tree based on the class I family C DNA pols is still in the provisional stage, some important conclusion can be drawn. First, the high-G+C and the low-G+C Gram-positive bacteria are not monophyletic. Second, the low-G+C Gram-positive bacteria contain multigenes of family C DNA pols (classes I and II). Third, the cyanobacterial family C DNA pol, classified as class III because it is encoded by a split gene, forms a group with the high-G+C Gram-positive bacteria. Received: 7 October 1998 / Accepted: 12 January 1999     

Comparative Molecular Evolution of Primary (<Emphasis Type="BoldItalic">Buchnera</Emphasis>) and Secondary Symbionts of Aphids Based on Two Protein-Coding Genes

A+T content, phylogenetic relationships, codon usage, evolutionary rates, and ratio of synonymous versus non-synonymous substitutions have been studied in partial sequences of the atpD and aroQ/pheA genes of primary (Buchnera) and secondary symbionts of aphids and a set of selected non-symbiotic bacteria, belonging to the five subdivisions of the Proteobacteria. Compared to the homologous genes of the last group, both genes belonging to Buchnera behave in a similar way, showing a higher A+T content, forming a monophyletic group, a loss in codon bias, especially in third base position, an evolutionary acceleration and an increase in the number of non-synonymous substitutions, confirming previous results reported elsewhere for other genes. When available, these properties have been partly observed with the secondary symbionts, but with values that are intermediate between Buchnera and free living Proteobacteria. They show high A+T content, but not as high as Buchnera, a non-solved phylogenetic position between Buchnera, and the other γ-Proteobacteria, a loss in codon bias, again not as high as in Buchnera and a significant evolutionary acceleration in the case of the three atpD genes, but not when considering aroQ/pheA genes. These results give support to the hypothesis that they are symbionts at different stages of the symbiotic accommodation to the host.     

Sequence Analysis of the Mitochondrial Genome of Sarcophyton glaucum: Conserved Gene Order Among Octocorals

The nucleotide sequence for an 11,715-bp segment of the mitochondrial genome of the octocoral Sarcophyton glaucum is presented, completing the analysis of the entire genome for this anthozoan member of the phylum Cnidaria. The genome contained the same 13 protein-coding and 2 ribosomal RNA genes as in other animals. However, it also included an unusual mismatch repair gene homologue reported previously and codes for only a single tRNA gene. Intermediate in length compared to two other cnidarians (17,443 and 18,911 bp), this organellar genome contained the smallest amount of noncoding DNA (428, compared to 1283 and 781 nt, respectively), making it the most compact one found for the phylum to date. The mitochondrial genes of S. glaucum exhibited an identical arrangement to that found in another octocoral, Renilla kolikeri, with five protein-coding genes in the same order as has been found in insect and vertebrate mitochondrial genomes. Although gene order appears to be highly conserved among octocorals, compared to the hexacoral, Metridium senile, few similarities were found. Like other metazoan mitochondrial genomes, the A + T composition was elevated and a general bias against codons ending in G or C was observed. However, an exception to this was the infrequent use of TGA compared to TGG to code for tryptophan. This divergent codon bias is unusual but appears to be a conserved feature among two rather distantly related anthozoans. Received: 27 January 1998 / Accepted: 25 May 1998     

The Mitochondrial Genome of Chlorogonium elongatum Inferred from the Complete Sequence

The 22,704-bp circular mitochondrial DNA (mtDNA) of the chlamydomonad alga Chlorogonium elongatum was completely cloned and sequenced. The genome encodes seven proteins of the respiratory electron transport chain, subunit 1 of the cytochrome oxidase complex (cox1), apocytochrome b (cob), five subunits of the NADH dehydrogenase complex (nad1, nad2, nad4, nad5, and nad6), a set of three tRNAs (Q, W, M), and the large (LSU)- and small (SSU)-subunit ribosomal RNAs. Six group-I introns were found, two each in the cox1, cob, and nad5 genes. In each intron an open reading frame (ORF) related to maturases or endonucleases was identified. Both the LSU and the SSU rRNA genes are split into fragments intermingled with each other and with other genes. Although the average A + T content is 62.2%, GC-rich clusters were detected in intergenic regions, in variable domains of the rRNA genes, and in introns and intron-encoded ORFs. A comparison of the genome maps reveals that C. elongatum and Chlamydomonas eugametos mtDNAs are more closely related to one another than either is to Chlamydomonas reinhardtii mtDNA. Received: 3 November 1997 / Accepted: 12 January 1998     

Complete Mitochondrial Genome of a Neotropical Fruit Bat, Artibeus jamaicensis, and a New Hypothesis of the Relationships of Bats to Other Eutherian Mammals

The complete mitochondrial genome was obtained from a microchiropteran bat, Artibeus jamaicensis. The presumptive amino acid sequence for the protein-coding genes was compared with predicted amino acid sequences from several representatives of other mammalian orders. Data were analyzed using maximum parsimony, maximum likelihood, and neighbor joining. All analyses placed bats as the sister group of carnivores, perissodactyls, artiodactyls, and cetaceans (e.g., 100% bootstrap value with both maximum parsimony and neighbor joining). The data strongly support a new hypothesis about the origin of bats, specifically a bat/ferungulate grouping. None of the analyses supported the superorder Archonta (bats, flying lemurs, primates, and tree shrews). Our hypothesis regarding the relationship of bats to other eutherian mammals is concordant with previous molecular studies and contrasts with hypotheses based solely on morphological criteria and an incomplete fossil record. The A. jamaicensis mitochondrial DNA control region has a complex pattern of tandem repeats that differs from previously reported chiropteran control regions. Received: 22 January 1998 / Accepted: 3 June 1998     

Conserved Gene Clusters in Bacterial Genomes Provide Further Support for the Primacy of RNA

Five complete bacterial genome sequences have been released to the scientific community. These include four (eu)Bacteria, Haemophilus influenzae, Mycoplasma genitalium, M. pneumoniae, and Synechocystis PCC 6803, as well as one Archaeon, Methanococcus jannaschii. Features of organization shared by these genomes are likely to have arisen very early in the history of the bacteria and thus can be expected to provide further insight into the nature of early ancestors. Results of a genome comparison of these five organisms confirm earlier observations that gene order is remarkably unpreserved. There are, nevertheless, at least 16 clusters of two or more genes whose order remains the same among the four (eu)Bacteria and these are presumed to reflect conserved elements of coordinated gene expression that require gene proximity. Eight of these gene orders are essentially conserved in the Archaea as well. Many of these clusters are known to be regulated by RNA-level mechanisms in Escherichia coli, which supports the earlier suggestion that this type of regulation of gene expression may have arisen very early. We conclude that although the last common ancestor may have had a DNA genome, it likely was preceded by progenotes with an RNA genome. Received: 10 March 1996 / Accepted: 20 May 1997     

Amelioration of Bacterial Genomes: Rates of Change and Exchange

Although bacterial species display wide variation in their overall GC contents, the genes within a particular species' genome are relatively similar in base composition. As a result, sequences that are novel to a bacterial genome—i.e., DNA introduced through recent horizontal transfer—often bear unusual sequence characteristics and can be distinguished from ancestral DNA. At the time of introgression, horizontally transferred genes reflect the base composition of the donor genome; but, over time, these sequences will ameliorate to reflect the DNA composition of the new genome because the introgressed genes are subject to the same mutational processes affecting all genes in the recipient genome. This process of amelioration is evident in a large group of genes involved in host-cell invasion by enteric bacteria and can be modeled to predict the amount of time required after transfer for foreign DNA to resemble native DNA. Furthermore, models of amelioration can be used to estimate the time of introgression of foreign genes in a chromosome. Applying this approach to a 1.43-megabase continuous sequence, we have calculated that the entire Escherichia coli chromosome contains more than 600 kb of horizontally transferred, protein-coding DNA. Estimates of amelioration times indicate that this DNA has accumulated at a rate of 31 kb per million years, which is on the order of the amount of variant DNA introduced by point mutations. This rate predicts that the E. coli and Salmonella enterica lineages have each gained and lost more than 3 megabases of novel DNA since their divergence. Received: 7 July 1996 / Accepted: 27 September 1996     

A physical and genetic map of theCorynebacterium glutamicum ATCC 13032 chromosome

A combined physical and genetic map of theCorynebacterium glutamicum ATCC 13032 chromosome was constructed using pulsed-field gel electrophoresis (PFGE) and hybridizations with cloned gene probes. Total genomic DNA was digested with the meganucleasesSwaI (5′-ATTTAAAT-3′),PacI (5′-TTAATTAA-3′), andPmeI (5′-GTTTAAAC-3′) yielding 26, 27, and 23 fragments, respectively. The chromosomal restriction fragments were then separated by PFGE. By summing up the lengths of the fragments generated with each of the three enzymes, a genome size of 3082 +/- 20 kb was determined. To identify adjacentSwaI fragments, a genomic cosmid library ofC. glutamicum was screened for chromosomal inserts containingSwaI sites. Southern blots of the PFGE gels were hybridized with these linking clones to connect theSwaI fragments in their natural order. By this method, about 90% of the genome could be ordered into three contigs. Two of the remaining gaps were closed by cross-hybridization of blottedSwaI digests using as probesPacI andPmeI fragments isolated from PFGE gels. The last gap in the chromosomal map was closed by hybridization experiments using partialSwaI digestions, thereby proving the circularity of the chromosome. By hybridization of gene probes toSwaI fragments separated by PFGE about 30 genes, including rRNA operons, IS element and transposon insertions were localized on the physical map.     

Horizontal Gene Transfer Involved in the Convergent Evolution of the Plasmid-Encoded Enantioselective 6-Hydroxynicotine Oxidases

The D- and L-specific nicotine oxidases are flavoproteins involved in the oxidative degradation of nicotine by the Gram-positive soil bacterium Arthrobacter nicotinovorans. Their structural genes are located on a 160-kbp plasmid together with those of other nicotine-degrading enzymes. They are structurally unrelated at the DNA as well as at the protein level. Each of these oxidases possesses a high degree of substrate specificity; their catalytic stereoselectivity is absolute, although they are able to bind both enantiomeric substrates with a similar affinity. It appears that the existence of these enzymes is the result of convergent evolution. The amino acid sequence of 6-hydroxy-l-nicotine oxidase (EC 1.5.3.6) as derived from the respective structural gene shows considerable structural similarity with eukaryotic monoamine oxidases (EC 1.4.3.4) but not with monoamine oxidases from prokaryotic bacteria including those of the genus Arthrobacter. These similarities are not confined to the nucleotide-binding sites. A 100-amino acid stretch at the N-terminal regions of 6-hydroxy-l-nicotine oxidase and human monoamine oxidases A possess a 35% homology. Overall, 27.0, 26.9, and 25.8% of the amino acid positions of the monoamine oxidases of Aspergillus niger (N), humans (A), and rainbow trout (Salmo gairdneri) are identical to those of 6-hydroxy-l-nicotine oxidase (Smith–Waterman algorithm). In addition, the G+C content of the latter enzyme is in the range of that of eukaryotic monoamine oxidases and definitely lower than that of the A. nicotinovorans DNA and even that of the pAO1 DNA. The primary structure of 6-hydroxy-d-nicotine oxidase (EC 1.5.3.5) does not reveal its evolutionary history as easily. Significant similarities are found with a mitomycin radical oxidase from Streptomyces lavendulae (23.3%) and a ``hypothetical protein' from Mycobacterium tuberculosis (26.0%). It is proposed that the plasmid-encoded gene of 6-hydroxy-l-nicotine oxidase evolved after horizontal transfer from an eukaryotic source. Received: 6 March 1998 / Accepted: 15 July 1998     

Base Composition Skews, Replication Orientation, and Gene Orientation in 12 Prokaryote Genomes

Variation in GC content, GC skew and AT skew along genomic regions was examined at third codon positions in completely sequenced prokaryotes. Eight out of nine eubacteria studied show GC and AT skews that change sign at the origin of replication. The leading strand in DNA replication is G-T rich at codon position 3 in six eubacteria, but C-T rich in two Mycoplasma species. In M. genitalium the AT and GC skews are symmetrical around the origin and terminus of replication, whereas its GC content variation has been shown to have a centre of symmetry elsewhere in the genome. Borrelia burgdorferi and Treponema pallidum show extraordinary extents of base composition skew correlated with direction of DNA replication. Base composition skews measured at third codon positions probably reflect mutational biases, whereas those measured over all bases in a sequence (or at codon positions 1 and 2) can be strongly affected by protein considerations due to the tendency in some bacteria for genes to be transcribed in the same direction that they are replicated. Consequently in some species the direction of skew for total genomic DNA is opposite to that for codon position 3. Received: 2 February 1998 / Accepted: 15 June 1998     

Synonymous Nucleotide Divergence and Saturation: Effects of Site-Specific Variations in Codon Bias and Mutation Rates

The synonymous divergence between Escherichia coli and Salmonella typhimurium is explained in a model where there is a large variation between mutation rates at different nucleotide sites in the genome. The model is based on the experimental observation that spontaneous mutation rates can vary over several orders of magnitude at different sites in a gene. Such site-specific variation must be taken into account when studying synonymous divergence and will result in an apparent saturation below the level expected from an assumption of uniform rates. Recently, it has been suggested that codon preference in enterobacteria has a very large site-specific variation and that the synonymous divergence between different species, e.g., E. coli and Salmonella, is saturated. In the present communication it is shown that when site-specific variation in mutation rates is introduced, there is no need to invoke assumptions of saturation and a large variability in codon preference. The same rate variation will also bring average mutation rates as estimated from synonymous sequence divergence into numerical agreement with experimental values. Received: 10 July 1998 / Accepted: 20 August 1998     

Identification of Sequences Encoding the Detoxification Metalloisomerase Glyoxalase I in Microbial Genomes from Several Pathogenic Organisms

The ubiquitous glyoxalase system, which is composed of two enzymes, removes cellular cytotoxic methylglyoxal (MG). In an effort to identify critical residues conserved in the evolution of the first enzyme in this system, glyoxalase I (GlxI), as well as the structural implications of sequence alterations in this enzyme, a search of the National Center for Biotechnology Information (NCBI) database of unfinished genomes was undertaken. Eleven putative GlxI sequences from pathogenic organisms were identified and analyses of these sequences in relation to the known and previously identified GlxI enzymes were performed. Several of these sequences show a very high similarity to the Escherichia coli GlxI sequence, most notably the 79% identity of the sequence identified from Yersinia pestis, the causative agent of bubonic plague. In addition to the conservation of residues critical to binding the catalytic metal in all of the proposed GlxI enzymes, four regions in the Homo sapiens GlxI enzyme are absent in all of the bacterial GlxI sequences, with the exception of Pseudomonas putida. Removal of these regions may alter the active-site conformation of the bacterial enzymes in relation to that of the H. sapiens. These differences may be targeted for the development of inhibitors selective to the bacterial enzymes. Received: 13 October 1999 / Accepted: 17 January 2000