Accumulation of Mn(II) in Deinococcus radiodurans facilitates gamma-radiation resistance.

Deinococcus radiodurans is extremely resistant to ionizing radiation. How this bacterium can grow under chronic gamma radiation [50 grays (Gy) per hour] or recover from acute doses greater than 10 kGy is unknown. We show that D. radiodurans accumulates very high intracellular manganese and low iron levels compared with radiation-sensitive bacteria and that resistance exhibits a concentration-dependent response to manganous chloride [Mn(II)]. Among the most radiation-resistant bacterial groups reported, Deinococcus, Enterococcus, Lactobacillus, and cyanobacteria accumulate Mn(II). In contrast, Shewanella oneidensis and Pseudomonas putida have high iron but low intracellular manganese concentrations and are very sensitive. We propose that Mn(II) accumulation facilitates recovery from radiation injury.

Genome-based bioinformatic selection of chromosomal Bacillus anthracis putative vaccine candidates coupled with proteomic identification of surface-associated antigens.

Bacillus anthracis (Ames strain) chromosome-derived open reading frames (ORFs), predicted to code for surface exposed or virulence related proteins, were selected as B. anthracis-specific vaccine candidates by a multistep computational screen of the entire draft chromosome sequence (February 2001 version, 460 contigs, The Institute for Genomic Research, Rockville, Md.). The selection procedure combined preliminary annotation (sequence similarity searches and domain assignments), prediction of cellular localization, taxonomical and functional screen and additional filtering criteria (size, number of paralogs). The reductive strategy, combined with manual curation, resulted in selection of 240 candidate ORFs encoding proteins with putative known function, as well as 280 proteins of unknown function. Proteomic analysis of two-dimensional gels of a B. anthracis membrane fraction, verified the expression of some gene products. Matrix-assisted laser desorption ionization-time-of-flight mass spectrometry analyses allowed identification of 38 spots cross-reacting with sera from B. anthracis immunized animals. These spots were found to represent eight in vivo immunogens, comprising of EA1, Sap, and 6 proteins whose expression and immunogenicity was not reported before. Five of these 8 immunogens were preselected by the bioinformatic analysis (EA1, Sap, 2 novel SLH proteins and peroxiredoxin/AhpC), as vaccine candidates. This study demonstrates that a combination of the bioinformatic and proteomic strategies may be useful in promoting the development of next generation anthrax vaccine.

Two C or not two C: recurrent disruption of Zn-ribbons, gene duplication, lineage-specific gene loss, and horizontal gene transfer in evolution of bacterial ribosomal proteins.

BACKGROUND: Ribosomal proteins are encoded in all genomes of cellular life forms and are, generally, well conserved during evolution. In prokaryotes, the genes for most ribosomal proteins are clustered in several highly conserved operons, which ensures efficient co-regulation of their expression. Duplications of ribosomal-protein genes are infrequent, and given their coordinated expression and functioning, it is generally assumed that ribosomal-protein genes are unlikely to undergo horizontal transfer. However, with the accumulation of numerous complete genome sequences of prokaryotes, several paralogous pairs of ribosomal protein genes have been identified. Here we analyze all such cases and attempt to reconstruct the evolutionary history of these ribosomal proteins. RESULTS: Complete bacterial genomes were searched for duplications of ribosomal proteins. Ribosomal proteins L36, L33, L31, S14 are each duplicated in several bacterial genomes and ribosomal proteins L11, L28, L7/L12, S1, S15, S18 are so far duplicated in only one genome each. Sequence analysis of the four ribosomal proteins, for which paralogs were detected in several genomes, two of the ribosomal proteins duplicated in one genome (L28 and S18), and the ribosomal protein L32 showed that each of them comes in two distinct versions. One form contains a predicted metal-binding Zn-ribbon that consists of four conserved cysteines (in some cases replaced by histidines), whereas, in the second form, these metal-chelating residues are completely or partially replaced. Typically, genomes containing paralogous genes for these ribosomal proteins encode both versions, designated C+ and C-, respectively. Analysis of phylogenetic trees for these seven ribosomal proteins, combined with comparison of genomic contexts for the respective genes, indicates that in most, if not all cases, their evolution involved a duplication of the ancestral C+ form early in bacterial evolution, with subsequent alternative loss of the C+ and C- forms in different lineages. Additionally, evidence was obtained for a role of horizontal gene transfer in the evolution of these ribosomal proteins, with multiple cases of gene displacement 'in situ', that is, without a change of the gene order in the recipient genome. CONCLUSIONS: A more complex picture of evolution of bacterial ribosomal proteins than previously suspected is emerging from these results, with major contributions of lineage-specific gene loss and horizontal gene transfer. The recurrent theme of emergence and disruption of Zn-ribbons in bacterial ribosomal proteins awaits a functional interpretation.

Horizontal gene transfer in prokaryotes: quantification and classification.

Comparative analysis of bacterial, archaeal, and eukaryotic genomes indicates that a significant fraction of the genes in the prokaryotic genomes have been subject to horizontal transfer. In some cases, the amount and source of horizontal gene transfer can be linked to an organism's lifestyle. For example, bacterial hyperthermophiles seem to have exchanged genes with archaea to a greater extent than other bacteria, whereas transfer of certain classes of eukaryotic genes is most common in parasitic and symbiotic bacteria. Horizontal transfer events can be classified into distinct categories of acquisition of new genes, acquisition of paralogs of existing genes, and xenologous gene displacement whereby a gene is displaced by a horizontally transferred ortholog from another lineage (xenolog). Each of these types of horizontal gene transfer is common among prokaryotes, but their relative contributions differ in different lineages. The fixation and long-term persistence of horizontally transferred genes suggests that they confer a selective advantage on the recipient organism. In most cases, the nature of this advantage remains unclear, but detailed examination of several cases of acquisition of eukaryotic genes by bacteria seems to reveal the evolutionary forces involved. Examples include isoleucyl-tRNA synthetases whose acquisition from eukaryotes by several bacteria is linked to antibiotic resistance, ATP/ADP translocases acquired by intracellular parasitic bacteria, Chlamydia and Rickettsia, apparently from plants, and proteases that may be implicated in chlamydial pathogenesis.

Genome sequence and comparative analysis of the solvent-producing bacterium Clostridium acetobutylicum.

The genome sequence of the solvent-producing bacterium Clostridium acetobutylicum ATCC 824 has been determined by the shotgun approach. The genome consists of a 3.94-Mb chromosome and a 192-kb megaplasmid that contains the majority of genes responsible for solvent production. Comparison of C. acetobutylicum to Bacillus subtilis reveals significant local conservation of gene order, which has not been seen in comparisons of other genomes with similar, or, in some cases closer, phylogenetic proximity. This conservation allows the prediction of many previously undetected operons in both bacteria. However, the C. acetobutylicum genome also contains a significant number of predicted operons that are shared with distantly related bacteria and archaea but not with B. subtilis. Phylogenetic analysis is compatible with the dissemination of such operons by horizontal transfer. The enzymes of the solventogenesis pathway and of the cellulosome of C. acetobutylicum comprise a new set of metabolic capacities not previously represented in the collection of complete genomes. These enzymes show a complex pattern of evolutionary affinities, emphasizing the role of lateral gene exchange in the evolution of the unique metabolic profile of the bacterium. Many of the sporulation genes identified in B. subtilis are missing in C. acetobutylicum, which suggests major differences in the sporulation process. Thus, comparative analysis reveals both significant conservation of the genome organization and pronounced differences in many systems that reflect unique adaptive strategies of the two gram-positive bacteria.

The structure and evolution of Penelope in the virilis species group of Drosophila: an ancient lineage of retroelements.

The Penelope element is the key element responsible for mobilization of other transposable elements in the course of hybrid dysgenesis in Drosophila virilis. Penelope has an unusually complex, highly variable organization in all studied species of the virlis group. Thc BRIDGE1 element from the fish Fugu rubripes is homologous to Penelope, and database searches detected additional homologous sequences among Expressed Sequence Tags from the flatworm Schistosoma mansonii and the nematode Ancylostoma caninum. Phylogenetic analysis shows that the reverse transcriptase of the Penelope group does not belong to any of the characterized major retroelement lineages, but apparently represents a novel branch of non-LTR retroelements. Sequence profile analysis results in the prediction that the C-terminal domain of the Penelope polyprotein is an active endonuclease related to intron-encoded endonucleases and the bacterial repair endonuclease UvrC, which could function as an integrase. No retroelements containing a predicted endonuclease of this family have been described previously. Phylogenetic analysis of Penelope copies isolated from several species of the virilis group reveals two subfamilies of Penelope elements, one of which includes full-length copies whose nucleotide sequences are almost identical, whereas the other one consists of highly diverged defective copies. Phylogenetic analysis of Penelope suggests both vertical transmission of the element and probable horizontal transfers. These findings support the notion that Penelope invasions occurred repeatedly in the evolution of the virilis group.

Lineage-specific gene expansions in bacterial and archaeal genomes.

Gene duplication is an important mechanistic antecedent to the evolution of new genes and novel biochemical functions. In an attempt to assess the contribution of gene duplication to genome evolution in archaea and bacteria, clusters of related genes that appear to have expanded subsequent to the diversification of the major prokaryotic lineages (lineage-specific expansions) were analyzed. Analysis of 21 completely sequenced prokaryotic genomes shows that lineage-specific expansions comprise a substantial fraction (approximately 5%-33%) of their coding capacities. A positive correlation exists between the fraction of the genes taken up by lineage-specific expansions and the total number of genes in a genome. Consistent with the notion that lineage-specific expansions are made up of relatively recently duplicated genes, >90% of the detected clusters consists of only two to four genes. The more common smaller clusters tend to include genes with higher pairwise similarity (as reflected by average score density) than larger clusters. Regardless of size, cluster members tend to be located more closely on bacterial chromosomes than expected by chance, which could reflect a history of tandem gene duplication. In addition to the small clusters, almost all genomes also contain rare large clusters of size > or =20. Several examples of the potential adaptive significance of these large clusters are explored. The presence or absence of clusters and their related genes was used as the basis for the construction of a similarity graph for completely sequenced prokaryotic genomes. The topology of the resulting graph seems to reflect a combined effect of common ancestry, horizontal transfer, and lineage-specific gene loss.

Conservation of the binding site for the arginine repressor in all bacterial lineages.

BACKGROUND: The arginine repressor ArgR/AhrC is a transcription factor universally conserved in bacterial genomes. Its recognition signal (the ARG box), a weak palindrome, is also conserved between genomes, despite a very low degree of similarity between individual sites within a genome. Thus, the arginine repressor is different from two other universal transcription factors - HrcA, whose recognition signal is very strongly conserved both within and between genomes, and LexA/DinR, whose signal is strongly conserved within, but not between, genomes. The arginine regulon is well studied in Escherichia coli and to some extent in Bacillus subtilis and some other genomes. Here, we apply the comparative genomic approach to the prediction of the ArgR-binding sites in all completely sequenced bacterial genomes. RESULTS: Orthologs of ArgR/AhrC were identified in the complete genomes of E. coli, Haemophilus influenzae, Vibrio choleras, B. subtilis, Mycobacterium tuberculosis, Thermotoga maritima, Chlamydia pneumoniae and Deinococcus radiodurans. Candidate arginine repressor binding sites were identified upstream of arginine transport and metabolism genes. CONCLUSIONS: We found that the ArgR/AhrC recognition signal is conserved in all genomes that contain genes encoding orthologous transcription factors of this family. All genomes studied except M. tuberculosis contain ABC transport cassettes (related to the Art system of E. coli) belonging to the candidate arginine regulons.

Genome of the extremely radiation-resistant bacterium Deinococcus radiodurans viewed from the perspective of comparative genomics.

The bacterium Deinococcus radiodurans shows remarkable resistance to a range of damage caused by ionizing radiation, desiccation, UV radiation, oxidizing agents, and electrophilic mutagens. D. radiodurans is best known for its extreme resistance to ionizing radiation; not only can it grow continuously in the presence of chronic radiation (6 kilorads/h), but also it can survive acute exposures to gamma radiation exceeding 1,500 kilorads without dying or undergoing induced mutation. These characteristics were the impetus for sequencing the genome of D. radiodurans and the ongoing development of its use for bioremediation of radioactive wastes. Although it is known that these multiple resistance phenotypes stem from efficient DNA repair processes, the mechanisms underlying these extraordinary repair capabilities remain poorly understood. In this work we present an extensive comparative sequence analysis of the Deinococcus genome. Deinococcus is the first representative with a completely sequenced genome from a distinct bacterial lineage of extremophiles, the Thermus-Deinococcus group. Phylogenetic tree analysis, combined with the identification of several synapomorphies between Thermus and Deinococcus, supports the hypothesis that it is an ancient group with no clear affinities to any of the other known bacterial lineages. Distinctive features of the Deinococcus genome as well as features shared with other free-living bacteria were revealed by comparison of its proteome to the collection of clusters of orthologous groups of proteins. Analysis of paralogs in Deinococcus has revealed several unique protein families. In addition, specific expansions of several other families including phosphatases, proteases, acyltransferases, and Nudix family pyrophosphohydrolases were detected. Genes that potentially affect DNA repair and recombination and stress responses were investigated in detail. Some proteins appear to have been horizontally transferred from eukaryotes and are not present in other bacteria. For example, three proteins homologous to plant desiccation resistance proteins were identified, and these are particularly interesting because of the correlation between desiccation and radiation resistance. Compared to other bacteria, the D. radiodurans genome is enriched in repetitive sequences, namely, IS-like transposons and small intergenic repeats. In combination, these observations suggest that several different biological mechanisms contribute to the multiple DNA repair-dependent phenotypes of this organism.

Gene conversions in genes encoding outer-membrane proteins in H. pylori and C. pneumoniae.

Helicobacter pylori and Chlamydia pneumoniae are both pathogenic to humans. Their genomes have recently been completed, allowing detailed study of their evolution and organization. Here we describe an evolutionary analysis of the H. pylori and C. pneumoniae genes that encode their outer-membrane proteins. By comparing complete genome sequences of two H. pylori strains and two C. pneumoniae strains, we identify multiple independent conversions among these genes. Such recombination events might provide a selective advantage for these bacterial pathogens.

SURVEY AND SUMMARY: holliday junction resolvases and related nucleases: identification of new families, phyletic distribution and evolutionary trajectories.

Holliday junction resolvases (HJRs) are key enzymes of DNA recombination. A detailed computer analysis of the structural and evolutionary relationships of HJRs and related nucleases suggests that the HJR function has evolved independently from at least four distinct structural folds, namely RNase H, endonuclease, endonuclease VII-colicin E and RusA. The endonuclease fold, whose structural prototypes are the phage lambda exonuclease, the very short patch repair nuclease (Vsr) and type II restriction enzymes, is shown to encompass by far a greater diversity of nucleases than previously suspected. This fold unifies archaeal HJRs, repair nucleases such as RecB and Vsr, restriction enzymes and a variety of predicted nucleases whose specific activities remain to be determined. Within the RNase H fold a new family of predicted HJRs, which is nearly ubiquitous in bacteria, was discovered, in addition to the previously characterized RuvC family. The proteins of this family, typified by Escherichia coli YqgF, are likely to function as an alternative to RuvC in most bacteria, but could be the principal HJRs in low-GC Gram-positive bacteria and AQUIFEX: Endonuclease VII of phage T4 is shown to serve as a structural template for many nucleases, including MCR:A and other type II restriction enzymes. Together with colicin E7, endonuclease VII defines a distinct metal-dependent nuclease fold. As a result of this analysis, the principal HJRs are now known or confidently predicted for all bacteria and archaea whose genomes have been completely sequenced, with many species encoding multiple potential HJRs. Horizontal gene transfer, lineage-specific gene loss and gene family expansion, and non-orthologous gene displacement seem to have been major forces in the evolution of HJRs and related nucleases. A remarkable case of displacement is seen in the Lyme disease spirochete Borrelia burgdorferi, which does not possess any of the typical HJRs, but instead encodes, in its chromosome and each of the linear plasmids, members of the lambda exonuclease family predicted to function as HJRs. The diversity of HJRs and related nucleases in bacteria and archaea contrasts with their near absence in eukaryotes. The few detected eukaryotic representatives of the endonuclease fold and the RNase H fold have probably been acquired from bacteria via horizontal gene transfer. The identity of the principal HJR(s) involved in recombination in eukaryotes remains uncertain; this function could be performed by topoisomerase IB or by a novel, so far undetected, class of enzymes. Likely HJRs and related nucleases were identified in the genomes of numerous bacterial and eukaryotic DNA viruses. Gene flow between viral and cellular genomes has probably played a major role in the evolution of this class of enzymes. This analysis resulted in the prediction of numerous previously unnoticed nucleases, some of which are likely to be new restriction enzymes.

Specific expansion of protein families in the radioresistant bacterium Deinococcus radiodurans.

Computer analysis of the complete genome of Deinococcus radiodurans R1 reveals a number of protein families, which are over-represented in this organism, compared to most other bacteria with known genome sequences. These families include both previously characterized and uncharacterized proteins. Most of the families whose functions are known or could be predicted seem to be related to stress-response and elimination of damage products (cell-cleaning). The two most prominent family expansions are the Nudix (MutT) family of pyrophosphohydrolases and a previously unnoticed family of proteins related to Bacillus subtilis DinB that could possess a metal-dependent enzymatic activity whose exact nature remains to be determined. Several proteins of the expanded families, particularly the Nudix family, are fused to other domains and form multidomain proteins that are so far unique for Deinococcus. The domain composition of some of these proteins indicates that they could be involved in novel DNA-repair pathways. Such unique proteins are good targets for knock-out and gene expression studies, which are aimed to shed light on the unusual features of this interesting bacterium.

Thermolysin and mitochondrial processing peptidase: how far structure-functional convergence goes.

The structure-functional convergence between two Zn-dependent proteases, namely thermolysin and mitochondrial processing peptidase (MPP), is described. These two families of nonhomologous enzymes show not only functional convergence of several active site residues as in chymotrypsin and subtilisin, but also structural convergence of overall molecular architectures including the beta-sheet arrangement and packing of the surrounding alpha-helices. The major functionally important structural elements are present in both enzymes with different topological connections and often in reverse main-chain orientation, but display similar packing. The structural comparison helps to rationalize sequence "inversion" of the HEXXH thermolysin consensus present as HXXEH in MPP. The described structural convergence may be due to a limited number of alternatives to build a Zn-protease that utilizes hydrogen bonding between a substrate main chain and the enzyme beta-sheet for substrate binding.

Genome sequence of the radioresistant bacterium Deinococcus radiodurans R1.

The complete genome sequence of the radiation-resistant bacterium Deinococcus radiodurans R1 is composed of two chromosomes (2,648,638 and 412,348 base pairs), a megaplasmid (177,466 base pairs), and a small plasmid (45,704 base pairs), yielding a total genome of 3,284, 156 base pairs. Multiple components distributed on the chromosomes and megaplasmid that contribute to the ability of D. radiodurans to survive under conditions of starvation, oxidative stress, and high amounts of DNA damage were identified. Deinococcus radiodurans represents an organism in which all systems for DNA repair, DNA damage export, desiccation and starvation recovery, and genetic redundancy are present in one cell.

The Zn-peptidase superfamily: functional convergence after evolutionary divergence.

Zn-dependent carboxypeptidases (ZnCP) cleave off the C-terminal amino acid residues from proteins and peptides. Here we describe a superfamily that unites classical ZnCP with other enzymes, most of which are known (or likely) to participate in metal-dependent peptide bond cleavage, but not necessarily in polypeptide substrates. It is demonstrated that aspartoacylase (ASP gene) and succinylglutamate desuccinylase (ASTE gene) are members of the ZnCP family. The Zn-binding site along with the structural core of the protein is shown to be conserved between ZnCP and another large family of hydrolases that includes mostly aminopeptidases (ZnAP). Both families (ZnCP and ZnAP) include not only proteases but also enzymes that perform N-deacylation, and enzymes that catalyze N-desuccinylation of amino acids. This is a result of functional convergence that apparently occurred after the divergence of the two families.

A superfamily of archaeal, bacterial, and eukaryotic proteins homologous to animal transglutaminases.

Computer analysis using profiles generated by the PSI-BLAST program identified a superfamily of proteins homologous to eukaryotic transglutaminases. The members of the new protein superfamily are found in all archaea, show a sporadic distribution among bacteria, and were detected also in eukaryotes, such as two yeast species and the nematode Caenorhabditis elegans. Sequence conservation in this superfamily primarily involves three motifs that center around conserved cysteine, histidine, and aspartate residues that form the catalytic triad in the structurally characterized transglutaminase, the human blood clotting factor XIIIa'. On the basis of the experimentally demonstrated activity of the Methanobacterium phage pseudomurein endoisopeptidase, it is proposed that many, if not all, microbial homologs of the transglutaminases are proteases and that the eukaryotic transglutaminases have evolved from an ancestral protease.

Comparative genomics of the Archaea (Euryarchaeota): evolution of conserved protein families, the stable core, and the variable shell.

Comparative analysis of the protein sequences encoded in the four euryarchaeal species whose genomes have been sequenced completely (Methanococcus jannaschii, Methanobacterium thermoautotrophicum, Archaeoglobus fulgidus, and Pyrococcus horikoshii) revealed 1326 orthologous sets, of which 543 are represented in all four species. The proteins that belong to these conserved euryarchaeal families comprise 31%-35% of the gene complement and may be considered the evolutionarily stable core of the archaeal genomes. The core gene set includes the great majority of genes coding for proteins involved in genome replication and expression, but only a relatively small subset of metabolic functions. For many gene families that are conserved in all euryarchaea, previously undetected orthologs in bacteria and eukaryotes were identified. A number of euryarchaeal synapomorphies (unique shared characters) were identified; these are protein families that possess sequence signatures or domain architectures that are conserved in all euryarchaea but are not found in bacteria or eukaryotes. In addition, euryarchaea-specific expansions of several protein and domain families were detected. In terms of their apparent phylogenetic affinities, the archaeal protein families split into bacterial and eukaryotic families. The majority of the proteins that have only eukaryotic orthologs or show the greatest similarity to their eukaryotic counterparts belong to the core set. The families of euryarchaeal genes that are conserved in only two or three species constitute a relatively mobile component of the genomes whose evolution should have involved multiple events of lineage-specific gene loss and horizontal gene transfer. Frequently these proteins have detectable orthologs only in bacteria or show the greatest similarity to the bacterial homologs, which might suggest a significant role of horizontal gene transfer from bacteria in the evolution of the euryarchaeota.

Short repeats and IS elements in the extremely radiation-resistant bacterium Deinococcus radiodurans and comparison to other bacterial species.

Computer analysis of the complete genome of Deinococcus radiodurans R1 has shown that the number of insertion sequences (ISs) and small noncoding repeats (SNRs) it contains is very high, and comparable with those of Escherichia coli. IS elements and several families of SNRs are described, together with their possible function in the D. radiodurans genome.