Comparative genome analysis of the pathogenic spirochetes Borrelia burgdorferi and Treponema pallidum

Broadly conserved FlgV controls flagellar assembly and Borrelia burgdorferi dissemination in mice

Flagella propel pathogens through their environments, yet are expensive to synthesize and are immunogenic. Thus, complex hierarchical regulatory networks control flagellar gene expression. Spirochetes are highly motile bacteria, but peculiarly, the archetypal flagellar regulator σ28 is absent in the Lyme spirochete Borrelia burgdorferi. Here, we show that gene bb0268 (flgV) in B. burgdorferi, previously and incorrectly annotated to encode the RNA-binding protein Hfq, is instead a structural flagellar component that modulates flagellar assembly. The flgV gene is broadly conserved in the flagellar superoperon alongside σ28 in many Spirochaetae, Firmicutes and other phyla, with distant homologs in Epsilonproteobacteria. We find that B. burgdorferi FlgV is localized within flagellar basal bodies, and strains lacking flgV produce fewer and shorter flagellar filaments and are defective in cell division and motility. During the enzootic cycle, flgV-deficient B. burgdorferi survive and replicate in Ixodes ticks but are attenuated for infection and dissemination in mice. Our work defines infection timepoints when spirochete motility is most crucial and implicates FlgV as a broadly distributed structural flagellar component that modulates flagellar assembly.

Broadly conserved FlgV controls flagellar assembly and Borrelia burgdorferi dissemination in mice

Flagella propel pathogens through their environments yet are expensive to synthesize and are immunogenic. Thus, complex hierarchical regulatory networks control flagellar gene expression. Spirochetes are highly motile bacteria, but peculiarly in the Lyme spirochete Borrelia burgdorferi, the archetypal flagellar regulator σ28 is absent. We rediscovered gene bb0268 in B. burgdorferi as flgV, a broadly-conserved gene in the flagellar superoperon alongside σ28 in many Spirochaetes, Firmicutes and other phyla, with distant homologs in Epsilonproteobacteria. We found that B. burgdorferi FlgV is localized within flagellar motors. B. burgdorferi lacking flgV construct fewer and shorter flagellar filaments and are defective in cell division and motility. During the enzootic cycle, B. burgdorferi lacking flgV survive and replicate in Ixodes ticks but are attenuated for dissemination and infection in mice. Our work defines infection timepoints when spirochete motility is most crucial and implicates FlgV as a broadly distributed structural flagellar component that modulates flagellar assembly.

Genomic insights into the c-di-GMP signaling and biofilm development in the saprophytic spirochete Leptospira biflexa

C-di-GMP is a bacterial second messenger with central role in biofilm formation. Spirochete bacteria from Leptospira genus present a wide diversity, with species of medical importance and environmental species, named as saprophytic. Leptospira form biofilms in the rat's reservoir kidneys and in the environment. Here, we performed genomic analyses to identify enzymatic and effector c-di-GMP proteins in the saprophytic biofilm-forming species Leptospira biflexa serovar Patoc. We identified 40 proteins through local alignments. Amongst them, 16 proteins are potentially functional diguanylate cyclases, phosphodiesterases, or hybrid proteins. We also identified nine effectors, including PilZ proteins. Enrichment analyses suggested that c-di-GMP interacts with cAMP signaling system, CsrA system, and flagella assembly regulation during biofilm development of L. biflexa. Finally, we identified eight proteins in the pathogen Leptospira interrogans serovar Copenhageni that share high similarity with L. biflexa c-di-GMP-related proteins. This work revealed proteins related to c-di-GMP turnover and cellular response in Leptospira and their potential roles during biofilm development.

Investigation of the immune escape mechanism of Treponema pallidum

BACKGROUND

Syphilis is a chronic sexually transmitted disease caused by Treponema pallidum subspecies pallidum (T. pallidum), which is a public health problem that seriously affects human health worldwide. T. pallidum is characterized by early transmission and immune escape and is therefore termed an "invisible pathogen".

METHODS

This review systematically summarizes the host's innate and adaptive immune responses to T. pallidum infection as well as the escape mechanisms of T. pallidum.

PURPOSE

To lay the foundation for assessing the pathogenic mechanism and the systematic prevention and treatment of syphilis.

CONCLUSION

The immune escape mechanism of T. pallidum plays an important role in its survival. Exploring the occurrence and development of these mechanisms has laid the foundation for the development of syphilis vaccine.

Treponema pallidum genome sequencing from six continents reveals variability in vaccine candidate genes and dominance of Nichols clade strains in Madagascar

In spite of its immutable susceptibility to penicillin, Treponema pallidum (T. pallidum) subsp. pallidum continues to cause millions of cases of syphilis each year worldwide, resulting in significant morbidity and mortality and underscoring the urgency of developing an effective vaccine to curtail the spread of the infection. Several technical challenges, including absence of an in vitro culture system until very recently, have hampered efforts to catalog the diversity of strains collected worldwide. Here, we provide near-complete genomes from 196 T. pallidum strains-including 191 T. pallidum subsp. pallidum-sequenced directly from patient samples collected from 8 countries and 6 continents. Maximum likelihood phylogeny revealed that samples from most sites were predominantly SS14 clade. However, 99% (84/85) of the samples from Madagascar formed two of the five distinct Nichols subclades. Although recombination was uncommon in the evolution of modern circulating strains, we found multiple putative recombination events between T. pallidum subsp. pallidum and subsp. endemicum, shaping the genomes of several subclades. Temporal analysis dated the most recent common ancestor of Nichols and SS14 clades to 1717 (95% HPD: 1543-1869), in agreement with other recent studies. Rates of SNP accumulation varied significantly among subclades, particularly among different Nichols subclades, and was associated in the Nichols A subclade with a C394F substitution in TP0380, a ERCC3-like DNA repair helicase. Our data highlight the role played by variation in genes encoding putative surface-exposed outer membrane proteins in defining separate lineages, and provide a critical resource for the design of broadly protective syphilis vaccines targeting surface antigens.

Identification and Functional Assessment of the First Placental Adhesin of Treponema pallidum That May Play Critical Role in Congenital Syphilis

Syphilis is a global, re-emerging sexually transmitted infection and congenital syphilis remains a major cause of adverse pregnancy outcomes due to bacterial infection in developing nations with a high rate of fetus loss. The molecular mechanisms involved in pathogenesis of the causative agent, Treponema pallidum subsp. pallidum remain poorly understood due to the difficulties of working with this pathogen, including the inability to grow it in pure culture. To reduce the spread of syphilis, we must first increase our knowledge of the virulence factors of T. pallidum and their contribution to syphilis manifestations. Tp0954 was predicted to be a surface lipoprotein of T. pallidum. Therefore, we experimentally demonstrated that Tp0954 is indeed a surface protein and further investigated its role in mediating bacterial attachment to various mammalian host cells. We found that expression of Tp0954 in a poorly adherent, but physiologically related derivative strain of the Lyme disease causing spirochete Borrelia burgdorferi B314 strain promotes its binding to epithelial as well as non-epithelial cells including glioma and placental cell lines. We also found that Tp0954 expression facilitates binding of this strain to purified dermatan sulfate and heparin, and also that bacterial binding to mammalian cell lines is mediated by the presence of heparan sulfate and dermatan sulfate in the extracellular matrix of the specific cell lines. These results suggest that Tp0954 may be involved not only in initiating T. pallidum infection by colonizing skin epithelium, but it may also contribute to disseminated infection and colonization of distal tissues. Significantly, we found that Tp0954 promotes binding to the human placental choriocarcinoma BeWo cell line, which is of trophoblastic endocrine cell type, as well as human placental tissue sections, suggesting its role in placental colonization and possible contribution to transplacental transmission of T. pallidum. Altogether, these novel findings offer an important step toward unraveling syphilis pathogenesis, including placental colonization and T. pallidum vertical transmission from mother to fetus during pregnancy.

Leptospira interrogans Bat proteins impair host hemostasis by fibrinogen cleavage and platelet aggregation inhibition

Leptospirosis is a worldwide spread zoonosis, caused by pathogenic Leptospira. Evidences suggest that compromised hemostasis might be involved in the leptospirosis pathophysiology. In the genome of L. interrogans serovar Copenhageni, we found two genes coding for proteins which comprise von Willebrand factor (VWF) A domains (BatA and BatB). As VWF A domains exhibit multiple binding sites which contributes to human VWF hemostatic functions, we hypothesized that the L. interrogans BatA and BatB proteins could be involved in the hemostatic impairment during leptospirosis. We have cloned, expressed in Escherichia coli, and purified recombinant BatA and BatB. The influence of recombinant BatA and BatB on different in vitro hemostatic assays evaluating the enzymatic activity, platelet aggregation and fibrinogen integrity was investigated. We describe BatB as a new serine protease which is able to cleave thrombin chromogenic substrate, fibrin, fibrinogen, gelatin and casein; while BatA is active only towards fibrinogen. BatA and BatB interfere with the platelet aggregation induced by VWF/ristocetin and thrombin. Our results suggest an important role of the L. interrogans serovar Copenhageni Bat proteins in the hemostasis dysfunction observed during leptospirosis and contribute to the understanding of the leptospirosis pathophysiological mechanisms.

Evolutionary convergence and divergence in archaeal chromosomal proteins and Chromo-like domains from bacteria and eukaryotes

SH3-fold-β-barrel domains of the chromo-like superfamily recognize epigenetic marks in eukaryotic proteins. Their provenance has been placed either in archaea, based on apparent structural similarity to chromatin-compacting Sul7d and Cren7 proteins, or in bacteria based on the presence of sequence homologs. Using sequence and structural evidence we establish that the archaeal Cren7/Sul7 proteins emerged from a zinc ribbon (ZnR) ancestor. Further, we show that the ancestral eukaryotic chromo-like domains evolved from bacterial versions, likely acquired from early endosymbioses, which already possessed an aromatic cage for recognition of modified amino-groups. These bacterial versions are part of a radiation of secreted SH3-fold domains, which spawned both chromo-like domains and classical SH3 domains in the context of peptide-recognition in the peptidoglycan or the extracellular matrix. This establishes that Cren7/Sul7 converged to a “SH3”-like state from a ZnR precursor via the loss of metal-chelation and acquisition of stronger hydrophobic interactions; it is unlikely to have participated in the evolution of the chromo-like domains. We show that archaea possess several Cren7/Sul7-related proteins with intact Zn-chelating ligands, which we predict to play previously unstudied roles in chromosome segregation during cell-division comparable to the PRC barrel and CdvA domain proteins.

Characterization of the inner membrane protein BB0173 from Borrelia burgdorferi

BACKGROUND

The bacterial spirochete Borrelia burgdorferi is the causative agent of the most commonly reported arthropod-borne illness in the United States, Lyme disease. A family of proteins containing von Willebrand Factor A (VWFA) domains adjacent to a MoxR AAA+ ATPase have been found to be highly conserved in the genus Borrelia. Previously, a VWFA domain containing protein of B. burgdorferi, BB0172, was determined to be an outer membrane protein capable of binding integrin α3β1. In this study, the characterization of a new VWFA domain containing membrane protein, BB0173, is evaluated in order to define the location and topology of this multi-spanning membrane protein. In addition, functional predictions are made.

RESULTS

Our results show that BB0173, in contrast to BB0172, is an inner membrane protein, in which the VWFA domain is exposed to the periplasmic space. Further, BB0173 was predicted to have an aerotolerance regulator domain, and expression of BB0173 and the surrounding genes was evaluated under aerobic and microaerophilic conditions, revealing that these genes are downregulated under aerobic conditions. Since the VWFA domain containing proteins of B. burgdorferi are highly conserved, they are likely required for survival of the pathogen through sensing diverse environmental oxygen conditions.

CONCLUSIONS

Presently, the complex mechanisms that B. burgdorferi uses to detect and respond to environmental changes are not completely understood. However, studying the mechanisms that allow B. burgdorferi to survive in the highly disparate environments of the tick vector and mammalian host could allow for the development of novel methods of preventing acquisition, survival, or transmission of the spirochete. In this regard, a putative membrane protein, BB0173, was characterized. BB0173 was found to be highly conserved across pathogenic Borrelia, and additionally contains several truly transmembrane domains, and a Bacteroides aerotolerance-like domain. The presence of these functional domains and the highly conserved nature of this protein, strongly suggests a required function of BB0173 in the survival of B. burgdorferi.

1,2-Diacylglycerol choline phosphotransferase catalyzes the final step in the unique Treponema denticola phosphatidylcholine biosynthesis pathway

Treponema denticola synthesizes phosphatidylcholine through a licCA-dependent CDP-choline pathway identified only in the genus Treponema. However, the mechanism of conversion of CDP-choline to phosphatidylcholine remained unclear. We report here characterization of TDE0021 (herein designated cpt) encoding a 1,2-diacylglycerol choline phosphotransferase homologous to choline phosphotransferases that catalyze the final step of the highly conserved Kennedy pathway for phosphatidylcholine synthesis in eukaryotes. T. denticola Cpt catalyzed in vitro phosphatidylcholine formation from CDP-choline and diacylglycerol, and full activity required divalent manganese. Allelic replacement mutagenesis of cpt in T. denticola resulted in abrogation of phosphatidylcholine synthesis. T. denticola Cpt complemented a Saccharomyces cerevisiae CPT1 mutant, and expression of the entire T. denticola LicCA-Cpt pathway in E. coli resulted in phosphatidylcholine biosynthesis. Our findings show that T. denticola possesses a unique phosphatidylcholine synthesis pathway combining conserved prokaryotic choline kinase and CTP:phosphocholine cytidylyltransferase activities with a 1,2-diacylglycerol choline phosphotransferase that is common in eukaryotes. Other than in a subset of mammalian host-associated Treponema that includes T. pallidum, this pathway is found in neither bacteria nor Archaea. Molecular dating analysis of the Cpt gene family suggests that a horizontal gene transfer event introduced this gene into an ancestral Treponema well after its divergence from other spirochetes.

Impact of genomics on the understanding of microbial evolution and classification: the importance of Darwin's views on classification

Analyses of genome sequences, by some approaches, suggest that the widespread occurrence of horizontal gene transfers (HGTs) in prokaryotes disguises their evolutionary relationships and have led to questioning of the Darwinian model of evolution for prokaryotes. These inferences are critically examined in the light of comparative genome analysis, characteristic synapomorphies, phylogenetic trees and Darwin's views on examining evolutionary relationships. Genome sequences are enabling discovery of numerous molecular markers (synapomorphies) such as conserved signature indels (CSIs) and conserved signature proteins (CSPs), which are distinctive characteristics of different prokaryotic taxa. Based on these molecular markers, exhibiting high degree of specificity and predictive ability, numerous prokaryotic taxa of different ranks, currently identified based on the 16S rRNA gene trees, can now be reliably demarcated in molecular terms. Within all studied groups, multiple CSIs and CSPs have been identified for successive nested clades providing reliable information regarding their hierarchical relationships and these inferences are not affected by HGTs. These results strongly support Darwin's views on evolution and classification and supplement the current phylogenetic framework based on 16S rRNA in important respects. The identified molecular markers provide important means for developing novel diagnostics, therapeutics and for functional studies providing important insights regarding prokaryotic taxa.

Transcription factors, chromatin proteins and the diversification of Hemiptera

Availability of complete genomes provides a means to explore the evolution of enormous developmental, morphological, and behavioral diversity among insects. Hemipterans in particular show great diversity of both morphology and life history within a single order. To better understand the role of transcription regulators in the diversification of hemipterans, using sequence profile searches and hidden Markov models we computationally analyzed transcription factors (TFs) and chromatin proteins (CPs) in the recently available Rhodnius prolixus genome along with 13 other insect and 4 non-insect arthropod genomes. We generated a comprehensive collection of TFs and CPs across arthropods including 303 distinct types of domains in TFs and 139 in CPs. This, along with the availability of two hemipteran genomes, R. prolixus and Acyrthosiphon pisum, helped us identify possible determinants for their dramatic morphological and behavioral divergence. We identified five domain families (i.e. Pipsqueak, SAZ/MADF, THAP, FLYWCH and BED finger) as having undergone differential patterns of lineage-specific expansion in hemipterans or within hemipterans relative to other insects. These expansions appear to be at least in part driven by transposons, with the DNA-binding domains of transposases having provided the raw material for emergence of new TFs. Our analysis suggests that while R. prolixus probably retains a state closer to the ancestral hemipteran, A. pisum represents a highly derived state, with the emergence of asexual reproduction potentially favoring genome duplication and transposon expansion. Both hemipterans are predicted to possess active DNA methylation systems. However, in the course of their divergence, aphids seem to have expanded the ancestral hemipteran DNA methylation along with a distinctive linkage to the histone methylation system, as suggested by expansion of SET domain methylases, including those fused to methylated CpG recognition domains. Thus, differential use of DNA methylation and histone methylation might have played a role in emergence of polyphenism and cyclic parthenogenesis from the ancestral hemipteran.

Phosphoenolpyruvate Phosphotransferase System Components Modulate Gene Transcription and Virulence of Borrelia burgdorferi

The phosphoenolpyruvate phosphotransferase system (PEP-PTS) and adenylate cyclase (AC) IV (encoded by BB0723 [cyaB]) are well conserved in different species of Borrelia. However, the functional roles of PEP-PTS and AC in the infectious cycle of Borrelia have not been characterized previously. We examined 12 PEP-PTS transporter component mutants by needle inoculation of mice to assess their ability to cause mouse infection. Transposon mutants with mutations in the EIIBC components (ptsG) (BB0645, thought to be involved in glucose-specific transport) were unable to cause infection in mice, while all other tested PEP-PTS mutants retained infectivity. Infectivity was partially restored in an in trans-complemented strain of the ptsG mutant. While the ptsG mutant survived normally in unfed as well as fed ticks, it was unable to cause infection in mice by tick transmission, suggesting that the function of ptsG is essential to establish infection by either needle inoculation or tick transmission. In Gram-negative organisms, the regulatory effects of the PEP-PTS are mediated by adenylate cyclase and cyclic AMP (cAMP) levels. A recombinant protein encoded by B. burgdorferi BB0723 (a putative cyaB homolog) was shown to have adenylate cyclase activity in vitro; however, mutants with mutations in this gene were fully infectious in the tick-mouse infection cycle, indicating that its function is not required in this process. By transcriptome analysis, we demonstrated that the ptsG gene may directly or indirectly modulate gene expression of Borrelia burgdorferi. Overall, the PEP-PTS glucose transporter PtsG appears to play important roles in the pathogenesis of B. burgdorferi that extend beyond its transport functions.

Genomic survey and expression analysis of DNA repair genes in the genus Leptospira

Leptospirosis is an emerging zoonosis with important economic and public health consequences and is caused by pathogenic leptospires. The genus Leptospira belongs to the order Spirochaetales and comprises saprophytic (L. biflexa), pathogenic (L. interrogans) and host-dependent (L. borgpetersenii) members. Here, we present an in silico search for DNA repair pathways in Leptospira spp. The relevance of such DNA repair pathways was assessed through the identification of mRNA levels of some genes during infection in animal model and after exposition to spleen cells. The search was performed by comparison of available Leptospira spp. genomes in public databases with known DNA repair-related genes. Leptospires exhibit some distinct and unexpected characteristics, for instance the existence of a redundant mechanism for repairing a chemically diverse spectrum of alkylated nucleobases, a new mutS-like gene and a new shorter version of uvrD. Leptospira spp. shares some characteristics from Gram-positive, as the presence of PcrA, two RecQ paralogs and two SSB proteins; the latter is considered a feature shared by naturally competent bacteria. We did not find a significant reduction in the number of DNA repair-related genes in both pathogenic and host-dependent species. Pathogenic leptospires were enriched for genes dedicated to base excision repair and non-homologous end joining. Their evolutionary history reveals a remarkable importance of lateral gene transfer events for the evolution of the genus. Up-regulation of specific DNA repair genes, including components of SOS regulon, during infection in animal model validates the critical role of DNA repair mechanisms for the complex interplay between host/pathogen.

Elucidating the function of the RpoS regulon

ABSTRACT: 

Bacterial adaptation to suboptimal nutrient environments, including host and/or extreme environments, is subject to complex, coordinated control involving many proteins and RNAs. Among the γ-proteobacteria, which includes many pathogens, the RpoS regulon has been a key focus for many years. Although the RpoS regulator was first identified as a growth phase-dependent regulator, our current understanding of RpoS is now more nuanced as this central regulator also has roles in exponential phase, biofilm development, bacterial virulence and bacterial persistence, as well as in stress adaptation. Induction of RpoS can also exert substantial metabolic effects by negatively regulating key systems including flagella biosynthesis, cryptic phage gene expression and the tricarboxylic acid cycle. Although core RpoS-controlled metabolic functions are conserved, there are substantial differences in RpoS regulation even among closely related bacteria, indicating that regulatory plasticity may be an important aspect of RpoS regulation, which is important in evolutionary adaptation to specialized environments.

Complete Genome of Bacillus megaterium Siphophage Slash

The complete annotated genome sequence of Bacillus megaterium bacteriophage Slash is described here. Several key features related to morphogenesis, replication/recombination, DNA metabolism, and lysis are described. Slash also encodes a homolog of SleB, a germination-specific cell wall amidase.

Borrelia burgdorferi and tick proteins supporting pathogen persistence in the vector

Borrelia burgdorferi, a pathogen transmitted by Ixodes ticks, is responsible for a prevalent illness known as Lyme disease, and a vaccine for human use is unavailable. Recently, genome sequences of several B. burgdorferi strains and Ixodes scapularis ticks have been determined. In addition, remarkable progress has been made in developing molecular genetic tools to study the pathogen and vector, including their intricate relationship. These developments are helping unravel the mechanisms by which Lyme disease pathogens survive in a complex enzootic infection cycle. Notable discoveries have already contributed to understanding the spirochete gene regulation accounting for the temporal and spatial expression of B. burgdorferi genes during distinct phases of the lifecycle. A number of pathogen and vector gene products have also been identified that contribute to microbial virulence and/or persistence. These research directions will enrich our knowledge of vector-borne infections and contribute towards the development of preventative strategies against Lyme disease.

BB0172, a Borrelia burgdorferi outer membrane protein that binds integrin α3β1

Lyme disease is a multisystemic disorder caused by Borrelia burgdorferi infection. Upon infection, some B. burgdorferi genes are upregulated, including members of the microbial surface components recognizing adhesive matrix molecule (MSCRAMM) protein family, which facilitate B. burgdorferi adherence to extracellular matrix components of the host. Comparative genome analysis has revealed a new family of B. burgdorferi proteins containing the von Willebrand factor A (vWFA) domain. In the present study, we characterized the expression and membrane association of the vWFA domain-containing protein BB0172 by using in vitro transcription/translation systems in the presence of microsomal membranes and with detergent phase separation assays. Our results showed evidence of BB0172 localization in the outer membrane, the orientation of the vWFA domain to the extracellular environment, and its function as a metal ion-dependent integrin-binding protein. This is the first report of a borrelial adhesin with a metal ion-dependent adhesion site (MIDAS) motif that is similar to those observed in eukaryotic integrins and has a similar function.

Insights from the architecture of the bacterial transcription apparatus

We provide a portrait of the bacterial transcription apparatus in light of the data emerging from structural studies, sequence analysis and comparative genomics to bring out important but underappreciated features. We first describe the key structural highlights and evolutionary implications emerging from comparison of the cellular RNA polymerase subunits with the RNA-dependent RNA polymerase involved in RNAi in eukaryotes and their homologs from newly identified bacterial selfish elements. We describe some previously unnoticed domains and the possible evolutionary stages leading to the RNA polymerases of extant life forms. We then present the case for the ancient orthology of the basal transcription factors, the sigma factor and TFIIB, in the bacterial and the archaeo-eukaryotic lineages. We also present a synopsis of the structural and architectural taxonomy of specific transcription factors and their genome-scale demography. In this context, we present certain notable deviations from the otherwise invariant proteome-wide trends in transcription factor distribution and use it to predict the presence of an unusual lineage-specifically expanded signaling system in certain firmicutes like Paenibacillus. We then discuss the intersection between functional properties of transcription factors and the organization of transcriptional networks. Finally, we present some of the interesting evolutionary conundrums posed by our newly gained understanding of the bacterial transcription apparatus and potential areas for future explorations.

Scoring protein relationships in functional interaction networks predicted from sequence data

The abundance of diverse biological data from various sources constitutes a rich source of knowledge, which has the power to advance our understanding of organisms. This requires computational methods in order to integrate and exploit these data effectively and elucidate local and genome wide functional connections between protein pairs, thus enabling functional inferences for uncharacterized proteins. These biological data are primarily in the form of sequences, which determine functions, although functional properties of a protein can often be predicted from just the domains it contains. Thus, protein sequences and domains can be used to predict protein pair-wise functional relationships, and thus contribute to the function prediction process of uncharacterized proteins in order to ensure that knowledge is gained from sequencing efforts. In this work, we introduce information-theoretic based approaches to score protein-protein functional interaction pairs predicted from protein sequence similarity and conserved protein signature matches. The proposed schemes are effective for data-driven scoring of connections between protein pairs. We applied these schemes to the Mycobacterium tuberculosis proteome to produce a homology-based functional network of the organism with a high confidence and coverage. We use the network for predicting functions of uncharacterised proteins.

AVAILABILITY

Protein pair-wise functional relationship scores for Mycobacterium tuberculosis strain CDC1551 sequence data and python scripts to compute these scores are available at http://web.cbio.uct.ac.za/~gmazandu/scoringschemes.

Treponema pallidum, the stealth pathogen, changes, but how?

Treponema pallidum rapidly disseminates from a genital site of inoculation to diverse organs where it establishes persistent infection. T. pallidum has long been regarded as a stealth pathogen because of its poorly antigenic and non-inflammatory surface. There is now increasing evidence that antigenic variation also contributes to the ability of the spirochaete to evade host defences. Among the small number of proteins encoded by the T. pallidum genome with sequence similarity to well-characterized transcription factors is TP0262, an orthologue for cAMP regulatory protein (CRP) of Escherichia coli. Giacani and co-workers identified sequences matching the CRP consensus-binding motif upstream of the promoters of tprE, tprG and tprJ, three members of the T. pallidum repeat (tpr) gene family (subfamily II). Using electrophoretic mobility shift assay, DNaseI footprinting and an E. coli-based reporter system, they demonstrated that TP0262 specifically recognizes the putative binding sequences and that DNA binding is cAMP-dependent. Their report, a major methodological advance for syphilis research, suggests that T. pallidum has appropriated a paradigmatic global regulator of metabolic processes in heterotrophic bacteria to further its capacity for immune evasion in its obligate human host.

Borrelia burgdorferi vlsE antigenic variation is not mediated by RecA

RecA is a key protein linking genetic recombination to DNA replication and repair in bacteria. Previous functional characterization of Borrelia burgdorferi RecA indicated that the protein is mainly involved in genetic recombination rather than DNA repair. Genetic recombination may play a role in B. burgdorferi persistence by generation of antigenic variation. We report here the isolation of a recA null mutant in an infectious B. burgdorferi strain. Comparison of the in vitro growth characteristics of the mutant with those of the wild-type strain under various conditions showed no significant differences. While the RecA mutant was moderately more sensitive to UV irradiation and mitomycin C than the wild-type strain, the lack of RecA abolished allelic exchange in the mutant. Absence of RecA did not affect the ability of the mutant to infect mice. However, the RecA mutant was attenuated for joint infection in competitive-infection assays with the wild-type strain. vlsE sequence variation in mice was observed in both wild-type and RecA mutant spirochetes, indicating that the mechanism of antigenic variation is not homologous genetic recombination.

New metrics for comparative genomics

The availability of genome sequences from a variety of organisms presents an opportunity to apply this sequence information to solving the key problems of molecular biology. One of the principal roadblocks on this path is the lack of appropriate descriptors and metrics that could succinctly represent the new knowledge stemming from the genomic data. Several new metrics have recently been used in comparative genome analysis, yet challenges remain in finding an appropriate language for the emerging discipline of systems biology.

TprK sequence diversity accumulates during infection of rabbits with Treponema pallidum subsp. pallidum Nichols strain

The tprK gene in Treponema pallidum undergoes antigenic variation. In all T. pallidum isolates examined to date, except the Nichols type strain, heterogeneous tprK sequences have been identified. This heterogeneity is localized to seven variable (V) regions, and tprK sequence diversity accumulates with serial passage in naïve rabbits. The T. pallidum Nichols genome described a single tprK sequence, and after decades of independent passage, only minor tprK sequence diversity is seen among the Nichols strains from different laboratories. We hypothesized that T. pallidum Nichols is capable of only limited tprK diversification. To address this hypothesis, we passaged the T. pallidum Nichols strain in naïve rabbits at the peak of infection (rapid passage) or after the adaptive immune response had cleared most organisms in vivo (slow passage). After 22 rapid passages (9- to 10-day intervals), no tprK V region sequence changes were observed. In contrast, after two slow passages (30- to 35-day intervals), three V regions had sequences that were completely different from that of the original inoculum. New sequences were observed in all seven V regions by the fifth slow passage. In contrast to the rapid-passaged Nichols strain, rapid-passaged Chicago C, a clonal strain isolated from the highly diverse parent Chicago strain, developed significant tprK diversification. These findings suggest that tprK variation can occur, but at a lower rate, in Nichols and that immune pressure may be required for accumulation of bacteria with diverse tprK sequences. Adaptation to growth in rabbits may explain the limited repertoire of V region sequences seen in the Nichols strain.

Paths of lateral gene transfer of lysyl-aminoacyl-tRNA synthetases with a unique evolutionary transition stage of prokaryotes coding for class I and II varieties by the same organisms

Background

While the premise that lateral gene transfer (LGT) is a dominant evolutionary force is still in considerable dispute, the case for widespread LGT in the family of aminoacyl-tRNA synthetases (aaRS) is no longer contentious. aaRSs are ancient enzymes, guarding the fidelity of the genetic code. They are clustered in two structurally unrelated classes. Only lysine aminoacyl-tRNA synthetase (LysRS) is found both as a class 1 and a class 2 enzyme (LysRS1-2). Remarkably, in several extant prokaryotes both classes of the enzyme coexist, a unique phenomenon that has yet to receive its due attention.

Results

We applied a phylogenetic approach for determining the extent and origin of LGT in prokaryotic LysRS. Reconstructing species trees for Archaea and Bacteria, and inferring that their last common ancestors encoded LysRS1 and LysRS2, respectively, we studied the gains and losses of both classes. A complex pattern of LGT events emerged. In specific groups of organisms LysRS1 was replaced by LysRS2 (and vice versa). In one occasion, within the alpha proteobacteria, a LysRS2 to LysRS1 LGT was followed by reversal to LysRS2. After establishing the most likely LGT paths, we studied the possible origins of the laterally transferred genes. To this end, we reconstructed LysRS gene trees and evaluated the likely origins of the laterally transferred genes. While the sources of LysRS1 LGTs were readily identified, those for LysRS2 remain, for now, uncertain. The replacement of one LysRS by another apparently transits through a stage simultaneously coding for both synthetases, probably conferring a selective advantage to the affected organisms.

Conclusion

The family of LysRSs features complex LGT events. The currently available data were sufficient for identifying unambiguously the origins of LysRS1 but not of LysRS2 gene transfers. A selective advantage is suggested to organisms encoding simultaneously LysRS1-2.

Weighted genome trees: refinements and applications

There are many ways to group completed genome sequences in hierarchical patterns (trees) reflecting relationships between their genes. Such groupings help us organize biological information and bear crucially on underlying processes of genome and organismal evolution. Genome trees make use of all comparable genes but can variously weight the contributions of these genes according to similarity, congruent patterns of similarity, or prevalence among genomes. Here we explore such possible weighting strategies, in an analysis of 142 prokaryotic and 5 eukaryotic genomes. We demonstrate that alternate weighting strategies have different advantages, and we propose that each may have its specific uses in systematic or evolutionary biology. Comparisons of results obtained with different methods can provide further clues to major events and processes in genome evolution.

Comparative analysis of apicomplexa and genomic diversity in eukaryotes

The apicomplexans Plasmodium and Cryptosporidium have developed distinctive adaptations via lineage-specific gene loss and gene innovation in the process of diverging from a common parasitic ancestor. The two lineages have acquired distinct but overlapping sets of surface protein adhesion domains typical of animal proteins, but in no case do they share multidomain architectures identical to animals. Cryptosporidium, but not Plasmodium, possesses an animal-type O-linked glycosylation pathway, along with >30 predicted surface proteins having mucin-like segments. The two parasites have notable qualitative differences in conserved protein architectures associated with chromatin dynamics and transcription. Cryptosporidium shows considerable reduction in the number of introns and a concomitant loss of spliceosomal machinery components. We also describe additional molecular characteristics distinguishing Apicomplexa from other eukaryotes for which complete genome sequences are available.

Cross-species hybridization of a Borrelia burgdorferi DNA array reveals infection- and culture-associated genes of the unsequenced genome of the relapsing fever agent Borrelia hermsii

The known genome sequence of Borrelia burgdorferi, an agent of Lyme borreliosis, was used to study the genetic content and gene expression in B. hermsii, another spirochete pathogen and a cause of relapsing fever. Cross-species hybridization of a DNA array representing 1628 open reading frames (ORF) of B. burgdorferi with genomic DNA of B. hermsii indicated that the latter organism has at least 81% of the chromosomal genes and 43% of the plasmid genes of B. burgdorferi. We then carried out quantitative hybridization of the arrays with multiple replicates of cDNA produced from B. hermsii cells growing in the blood of infected mice or in culture medium that was adjusted to the same pH, temperature and a spirochete density as infected blood. Of 642 B. burgdorferi ORFs hybridized by all replicates under both conditions, 12 (1.9%) demonstrated differential expression by a regularized t-test and stringent criteria. BBP07 and BBG30, two plasmid-borne ORFs with the greatest measurable difference in expression between in vivo and in vitro conditions, putatively encode proteins of unknown function. Orthologues of BBP07 in B. hermsii were identified, and increased expression in infected mice was demonstrated by quantitative reverse-transcriptase polymerase chain reaction.

Functional properties of Borrelia burgdorferi recA

Functions of the Borrelia burgdorferi RecA protein were investigated in Escherichia coli recA null mutants. Complementation with B. burgdorferi recA increased survival of E. coli recA mutants by 3 orders of magnitude at a UV dose of 2,000 microJ/cm(2). The viability at this UV dose was about 10% that provided by the homologous recA gene. Expression of B. burgdorferi recA resulted in survival of E. coli at levels of mitomycin C that were lethal to noncomplemented hosts. B. burgdorferi RecA was as effective as E. coli RecA in mediating homologous recombination in E. coli. Furthermore, E. coli lambda phage lysogens complemented with B. burgdorferi recA produced phage even in the absence of UV irradiation. The level of phage induction was 55-fold higher than the level in cells complemented with the homologous recA gene, suggesting that B. burgdorferi RecA may possess an enhanced coprotease activity. This study indicates that B. burgdorferi RecA mediates the same functions in E. coli as the homologous E. coli protein mediates. However, the rapid loss of viability and the absence of induction in recA expression after UV irradiation in B. burgdorferi suggest that recA is not involved in the repair of UV-induced damage in B. burgdorferi. The primary role of RecA in B. burgdorferi is likely to be a role in some aspect of recombination.

Genome features of Leptospira interrogans serovar Copenhageni

We report novel features of the genome sequence of Leptospira interrogans serovar Copenhageni, a highly invasive spirochete. Leptospira species colonize a significant proportion of rodent populations worldwide and produce life-threatening infections in mammals. Genomic sequence analysis reveals the presence of a competent transport system with 13 families of genes encoding for major transporters including a three-member component efflux system compatible with the long-term survival of this organism. The leptospiral genome contains a broad array of genes encoding regulatory system, signal transduction and methyl-accepting chemotaxis proteins, reflecting the organism's ability to respond to diverse environmental stimuli. The identification of a complete set of genes encoding the enzymes for the cobalamin biosynthetic pathway and the novel coding genes related to lipopolysaccharide biosynthesis should bring new light to the study of Leptospira physiology. Genes related to toxins, lipoproteins and several surface-exposed proteins may facilitate a better understanding of the Leptospira pathogenesis and may serve as potential candidates for vaccine.

Identification of novel virulence-associated genes via genome analysis of hypothetical genes

The sequencing of bacterial genomes has opened new perspectives for identification of targets for treatment of infectious diseases. We have identified a set of novel virulence-associated genes (vag genes) by comparing the genome sequences of six human pathogens that are known to cause persistent or chronic infections in humans: Yersinia pestis, Neisseria gonorrhoeae, Helicobacter pylori, Borrelia burgdorferi, Streptococcus pneumoniae, and Treponema pallidum. This comparison was limited to genes annotated as hypothetical in the T. pallidum genome project. Seventeen genes with unknown functions were found to be conserved among these pathogens. Insertional inactivation of 14 of these genes generated nine mutants that were attenuated for virulence in a mouse infection model. Out of these nine genes, five were found to be specifically associated with virulence in mice as demonstrated by infection with Yersinia pseudotuberculosis in-frame deletion mutants. In addition, these five vag genes were essential only in vivo, since all the mutants were able to grow in vitro. These genes are broadly conserved among bacteria. Therefore, we propose that the corresponding vag gene products may constitute novel targets for antimicrobial therapy and that some vag mutants could serve as carrier strains for live vaccines.

The response regulator Rrp2 is essential for the expression of major membrane lipoproteins in Borrelia burgdorferi

Borrelia burgdorferi (Bb), the agent of Lyme disease, exists in nature through a complex enzootic life cycle that involves both ticks and mammals. As Bb transitions between its two diverse niches, profound adaptive changes occur that are reflected in differential patterns of gene expression, particularly involving lipoprotein genes. Using a mutagenesis approach, we show that Rrp2 (gene BB0763), one of the proteins predicted by the Bb genome (www.tigr.org) to be a response regulator of a two-component sensory transduction system, is a pivotal regulator governing the expression of major membrane lipoproteins such as OspC, DbpA, and Mlp8, as well as many other mammalian infection-associated immunogens of Bb. Sequence analysis additionally suggested that Rrp2 is a bacterial enhancer-binding protein, essential for sigma54-dependent gene activation. Mutagenesis of a key amino acid residue within a putative activation domain revealed that Rrp2 controlled lipoprotein expression by governing the expression of the alternative sigma-factor sigmas in a sigma54-dependent manner. We therefore propose a signal transduction pathway involving Rrp2, sigma54, and sigmas, which in concert control the expression of key lipoproteins and other infection-associated immunogens in Bb.

New knowledge from old: in silico discovery of novel protein domains in Streptomyces coelicolor

BACKGROUND

Streptomyces coelicolor has long been considered a remarkable bacterium with a complex life-cycle, ubiquitous environmental distribution, linear chromosomes and plasmids, and a huge range of pharmaceutically useful secondary metabolites. Completion of the genome sequence demonstrated that this diversity carried through to the genetic level, with over 7000 genes identified. We sought to expand our understanding of this organism at the molecular level through identification and annotation of novel protein domains. Protein domains are the evolutionary conserved units from which proteins are formed.

RESULTS

Two automated methods were employed to rapidly generate an optimised set of targets, which were subsequently analysed manually. A final set of 37 domains or structural repeats, represented 204 times in the genome, was developed. Using these families enabled us to correlate items of information from many different resources. Several immediately enhance our understanding both of S. coelicolor and also general bacterial molecular mechanisms, including cell wall biosynthesis regulation and streptomycete telomere maintenance.

DISCUSSION

Delineation of protein domain families enables detailed analysis of protein function, as well as identification of likely regions or residues of particular interest. Hence this kind of prior approach can increase the rate of discovery in the laboratory. Furthermore we demonstrate that using this type of in silico method it is possible to fairly rapidly generate new biological information from previously uncorrelated data.

Expression of spoT in Borrelia burgdorferi during serum starvation

Borrelia burgdorferi, the causative agent of Lyme disease, is transmitted by the tick Ixodes scapularis. A 2.9-kb fragment containing a putative spoT gene was isolated from B. burgdorferi genomic DNA by PCR amplification and cloned into a pBAD24 vector. The cloned gene complemented Escherichia coli mutant strain CF1693, which contains deletions of both the relA and spoT genes. The spoT gene in E. coli encodes a bifunctional enzyme capable of synthesizing and degrading (p)ppGpp, which mediates the stringent response during carbon source starvation. B. burgdorferi has been reported to have a stress response to serum starvation. Thin-layer chromatography was used to detect (p)ppGpp extracted from H(3)(32)PO(4)-labeled B. burgdorferi cells starved for serum in RPMI. B. burgdorferi spoT gene expression was characterized during fatty acid starvation. Northern analysis of spoT revealed detectable message at 2.5 min of starvation in RPMI. Expression of spoT during serum starvation increased approximately 6-fold during the 30 min that starvation conditions were maintained. Further, expression of spoT decreased when serum was added to serum-starved cells. Reverse transcriptase PCR (RT-PCR) was used to detect spoT mRNA from approximately 10(6) cells starved for serum in RPMI for 2.5 to 30 min or incubated in tick saliva for 15 min. Northern blot analysis suggests that spoT transcript was approximately 900 nucleotides in length. RT-PCR amplification of the transcript using several sets of primers confirmed this finding. Additionally, a truncated clone containing only the first 950 bp of the 2,001-bp spoT open reading frame was able to complement E. coli CF1693. The data suggest that B. burgdorferi exhibits a stringent response to serum starvation and during incubation in tick saliva.

Rational identification of new antibacterial drug targets that are essential for viability using a genomics-based approach

In the last two decades, the search for completely novel antibacterial agents has acquired a new sense of urgency due to the remarkable rise of antibiotic resistance among key bacterial pathogens. More recently, the advent of bacterial genomics has provided investigators with the data and bioinformatic tools to rationally identify novel antibacterial targets and the genome-scaled methodologies to validate them. Only 6 years have elapsed since the publication of the first complete bacterial genome sequence, but more than 50 complete microbial genome sequences are now available. This review will discuss the advantages and limitations of the existing bacterial genome dataset for the rational identification of novel antibacterial targets. Since the ability to rapidly identify essential genes where loss of function is coincident with loss of viability is the most important task of genomics-based target validation, essentiality testing methodologies (in which molecular genetic techniques are used to determine whether or not a gene product is required for viability of the parent cell) will be surveyed and their amenability to genome-scaled analysis assessed. Finally, we will discuss the impact of bacterial genomics to date on the development of novel and effective antibiotics.

Examination of the Borrelia burgdorferi transcriptome in Ixodes scapularis during feeding

Borrelia burgdorferi gene expression within the guts of engorging Ixodes scapularis ticks was examined by use of differential immunoscreening and differential expression with a customized amplified library. Fourteen chromosomal genes involved in energy metabolism, substrate transport, and signal transduction and 10 (4 chromosomal and 6 plasmid) genes encoding putative lipoproteins and periplasmic proteins were preferentially expressed in engorging ticks. These data demonstrate a new approach to the global analysis of B. burgdorferi genes that are preferentially expressed within the vector during feeding.

Sodium ion cycle in bacterial pathogens: evidence from cross-genome comparisons

Analysis of the bacterial genome sequences shows that many human and animal pathogens encode primary membrane Na+ pumps, Na+-transporting dicarboxylate decarboxylases or Na+ translocating NADH:ubiquinone oxidoreductase, and a number of Na+ -dependent permeases. This indicates that these bacteria can utilize Na+ as a coupling ion instead of or in addition to the H+ cycle. This capability to use a Na+ cycle might be an important virulence factor for such pathogens as Vibrio cholerae, Neisseria meningitidis, Salmonella enterica serovar Typhi, and Yersinia pestis. In Treponema pallidum, Chlamydia trachomatis, and Chlamydia pneumoniae, the Na+ gradient may well be the only energy source for secondary transport. A survey of preliminary genome sequences of Porphyromonas gingivalis, Actinobacillus actinomycetemcomitans, and Treponema denticola indicates that these oral pathogens also rely on the Na+ cycle for at least part of their energy metabolism. The possible roles of the Na+ cycling in the energy metabolism and pathogenicity of these organisms are reviewed. The recent discovery of an effective natural antibiotic, korormicin, targeted against the Na+ -translocating NADH:ubiquinone oxidoreductase, suggests a potential use of Na+ pumps as drug targets and/or vaccine candidates. The antimicrobial potential of other inhibitors of the Na+ cycle, such as monensin, Li+ and Ag+ ions, and amiloride derivatives, is discussed.

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.

Gene re-arrangement and fusion in Mycoplasma gallisepticum thyA-nrdFEI locus

Sequencing of Mycoplasma gallisepticum genome fragment containing thymidylate synthase and ribonucleotide reductase gene clusters reveals both its unusual organization and gene content. Sequence analysis indicates the presence of a gene whose product can be considered as a fusion of two full size proteins: the N-terminal part shows significant similarity to mycoplasmal dihydrofolate reductases, while the C-terminal part of the polypeptide chain shows significant similarity to eukaryotic deoxycytidylate deaminase. Phylogenetic analysis has suggested that the C-terminal part of the M. gallisepticum fusion gene and eukaryotic deoxycytidylate deaminase genes are xenologous. No chromosomal regions encoding peptides similar to the C-terminal part of this fusion protein were found in completely sequenced genomes of Mycoplasma genitalium and Mycoplasma pneumoniae. Genes for ribonucleoside diphosphate reductase alpha chain (nrdE), NrdI protein (nrdI), and ribonucleoside diphosphate reductase beta chain (nrdF) have an opposite direction of transcription with respect to genes for thymidylate synthase (thyA), and dihydrofolate reductase-deoxycytidylate deaminase fusion protein.

Regulatory potential, phyletic distribution and evolution of ancient, intracellular small-molecule-binding domains

Central cellular functions such as metabolism, solute transport and signal transduction are regulated, in part, via binding of small molecules by specialized domains. Using sensitive methods for sequence profile analysis and protein structure comparison, we exhaustively surveyed the protein sets from completely sequenced genomes for all occurrences of 21 intracellular small-molecule-binding domains (SMBDs) that are represented in at least two of the three major divisions of life (bacteria, archaea and eukaryotes). These included previously characterized domains such as PAS, GAF, ACT and ferredoxins, as well as three newly predicted SMBDs, namely the 4-vinyl reductase (4VR) domain, the NIFX domain and the 3-histidines (3H) domain. Although there are only a limited number of different superfamilies of these ancient SMBDs, they are present in numerous distinct proteins combined with various enzymatic, transport and signal-transducing domains. Most of the SMBDs show considerable evolutionary mobility and are involved in the generation of many lineage-specific domain architectures. Frequent re-invention of analogous architectures involving functionally related, but not homologous, domains was detected, such as, fusion of different SMBDs to several types of DNA-binding domains to form diverse transcription regulators in prokaryotes and eukaryotes. This is suggestive of similar selective forces affecting the diverse SMBDs and resulting in the formation of multidomain proteins that fit a limited number of functional stereotypes. Using the "guilt by association approach", the identification of SMBDs allowed prediction of functions and mode of regulation for a variety of previously uncharacterized proteins.

Lyme disease in central Europe

Lyme borreliosis is a fascinating disease, the aetiopathology of which is not yet completely known. Different subspecies of Borrelia burgdorferi sensu lato are responsible for the variable clinical course of the disease. Some new cutaneous (alopecia) and ocular (photophobia and retinal vasculitis) manisfestations have been described and the largest prospective study on erythema migrans during pregnancy was published during the last year. Optimal therapy of Lyme borreliosis is still lacking, but doxycycline, amoxicillin, penicillin, and ceftriaxone are recommended most frequently.

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.

Two independent transcriptional units control the complex and simultaneous expression of the bmp paralogous chromosomal gene family in Borrelia burgdorferi

The chromosomal paralogous gene family 36 encodes for four lipoproteins with high amino acid homology that are expressed in vivo in humans and animals and are immunogenic. Transcriptional analysis of the bmp gene cluster indicated that all four genes of this cluster are expressed in vitro and constitute two transcriptional units with a complex pattern of transcription, including alternative monocistronic and polycistronic messages. One unit consists of bmpD, whose transcription is coupled to the transcription of the ribosomal protein genes, rpsG and rpsL. The second unit includes bmpC, bmpA and bmpB. The simultaneous expression of the four bmp genes in Borrelia burgdorferi suggests that their gene products may have either different or complementary functions. Primer extension experiments identified promoters for bmpD, bmpC and bmpA, but not for bmpB. The concentration of gene-specific mRNA paralleled its promoter homology to the Escherichia coli sigma70 promoter. The linkage of bmpD expression to rpsL and rpsG suggests that the expression of this gene may be controlled by growth-related global regulation mechanisms in B. burgdorferi. These results indicate that the bmp family constitutes a good model for the investigation of complex regulation of chromosomal gene expression in this bacteria.

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.

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.