# Genome Sequence of the Saprophyte Leptospira biflexa Provides Insights into the Evolution of Leptospira and the Pathogenesis of Leptospirosis

Mathieu Picardeau, Dieter M. Bulach, Christiane Bouchier, Richard L.

Zuerner, Nora Zidane, Peter J. Wilson, Sophie Creno, Elizabeth S. Kuczek, Simona Bommezzadri, John C. Davis, et al.

► To cite this version:

Mathieu Picardeau, Dieter M. Bulach, Christiane Bouchier, Richard L. Zuerner, Nora Zidane, et al. Genome Sequence of the Saprophyte Leptospira biflexa Provides Insights into the Evolution of Leptospira and the Pathogenesis of Leptospirosis. PLoS ONE, 2008, 3 (2), pp.e1607. 10.1371/journal.pone.0001607. hal-03020737

HAL Id: hal-03020737

https://hal.science/hal-03020737

Submitted on 24 Nov 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
# Genome Sequence of the Saprophyte *Leptospira biflexa* Provides Insights into the Evolution of *Leptospira* and the Pathogenesis of Leptospirosis

Mathieu Picardeau¹,<fnref n="9" />, Dieter M. Bulach²,³,⁴,⁵, Christiane Bouchier⁵, Richard L. Zuerner⁶, Nora Zidane⁵, Peter J. Wilson⁷, Sophie Creno⁵, Elizabeth S. Kuczek⁷, Simona Bommezzadri¹, John C. Davis⁷, Annette McGrath⁷, Matthew J. Johnson⁷, Caroline Boursaux-Eude⁸, Torsten Seemann², Zoé Rouy<fnref n="9" />, Ross L. Coppel²,³, Julian I. Rood²,³,⁴, Aurélie Lajus<fnref n="9" />, John K. Davies²,³,⁴, Claudine Médigue<fnref n="9" />,¹⁰, Ben Adler²,³,⁴*

¹Unité de Biologie des Spirochètes, Institut Pasteur, Paris, France, ²Victorian Bioinformatics Consortium, Monash University, Clayton, Victoria, Australia, ³Australian Research Council Centre of Excellence in Structural and Functional Microbial Genomics, Department of Microbiology, Monash University, Clayton, Victoria, Australia, ⁴Australian Bacterial Pathogenesis Program, Department of Microbiology, Monash University, Clayton, Victoria, Australia, ⁵Plate-forme Génomique, Institut Pasteur, Paris, France, ⁶Bacterial Diseases of Livestock Research Unit, National Animal Disease Center (NADC), Agricultural Research Service (ARS), United States Department of Agriculture (USDA), Ames, Iowa, United States of America, ⁷Australian Genome Research Facility, Gehrmann Laboratories, University of Queensland, St. Lucia, Queensland, Australia, ⁸Plate-forme Intégration et Analyse génomique, Institut Pasteur, Paris, France, ⁹Commissariat à l'Energie Atomique (CEA), Direction des Sciences du Vivant, Laboratoire de Génomique Comparative, Institut de Génomique, Genoscope, Evry, France, ¹⁰Centre National de la Recherche Scientifique (CNRS) UMR8030, Génomique Métabolique, Evry, France

## Abstract

*Leptospira biflexa* is a free-living saprophytic spirochete present in aquatic environments. We determined the genome sequence of *L. biflexa*, making it the first saprophytic *Leptospira* to be sequenced. The *L. biflexa* genome has 3,590 protein-coding genes distributed across three circular replicons: the major 3,604 chromosome, a smaller 278-kb replicon that also carries essential genes, and a third 74-kb replicon. Comparative sequence analysis provides evidence that *L. biflexa* is an excellent model for the study of *Leptospira* evolution; we conclude that 2052 genes (61%) represent a progenitor genome that existed before divergence of pathogenic and saprophytic *Leptospira* species. Comparisons of the *L. biflexa* genome with two pathogenic *Leptospira* species reveal several major findings. Nearly one-third of the *L. biflexa* genes are absent in pathogenic *Leptospira*. We suggest that once incorporated into the *L. biflexa* genome, laterally transferred DNA undergoes minimal rearrangement due to physical restrictions imposed by high gene density and limited presence of transposable elements. In contrast, the genomes of pathogenic *Leptospira* species undergo frequent rearrangements, often involving recombination between insertion sequences. Identification of genes common to the two pathogenic species, *L. borgpetersenii* and *L. interrogans*, but absent in *L. biflexa*, is consistent with a role for these genes in pathogenesis. Differences in environmental sensing capacities of *L. biflexa*, *L. borgpetersenii*, and *L. interrogans* suggest a model which postulates that loss of signal transduction functions in *L. borgpetersenii* has impaired its survival outside a mammalian host, whereas *L. interrogans* has retained environmental sensory functions that facilitate disease transmission through water.

**Citation:** Picardeau M, Bulach DM, Bouchier C, Zuerner RL, Zidane N, et al (2008) Genome Sequence of the Saprophyte *Leptospira biflexa* Provides Insights into the Evolution of *Leptospira* and the Pathogenesis of Leptospirosis. PLoS ONE 3(2): e1607. doi:10.1371/journal.pone.0001607

**Editor:** Dana Davis, University of Minnesota, United States of America

**Received** November 16, 2007; **Accepted** January 17, 2008; **Published** February 13, 2008

This is an open-access article distributed under the terms of the Creative Commons Public Domain declaration which stipulates that, once placed in the public domain, this work may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose.

**Funding:** This work was supported by the following: National Health and Medical Research Council, Australia; Australian Research Council; U.S. Department of Agriculture; Institut Pasteur, Paris, France; French Ministry of Research ANR Jeunes Chercheurs (no. 05-JCJC-0105-01); Consortium National de Recherche en Génomique (no. RNG-20040057); ACI IMPBio 2004 (MicroScope project). Funders had no role in the planning, performance or analysis of the work.

**Competing Interests:** The authors have declared that no competing interests exist.

## Introduction

The genus *Leptospira* contains pathogenic and saprophytic species that differ in their capacity for survival in a vast array of environments that range from soil and water [1] to the tissues of mammalian hosts during acute and chronic infection [2]. Typically, long term colonization by pathogenic *Leptospira* of the proximal renal tubules of mammalian maintenance host species provides a persistent source of infection; thus, pathogenic *Leptospira* is shed in the urine of chronically infected animals, facilitating transmission to naïve hosts [2]. Leptospirosis is of considerable

importance to international public health, with more than half a million cases reported annually due largely to environmental exposure to pathogenic *Leptospira* species, with mortality rates of up to 25% in some outbreaks. In addition, leptospirosis in production animals results in a significant economic burden worldwide [2].

Recent applications of molecular taxonomy techniques to this genus reveal extensive genetic diversity within *Leptospira*, with more than 16 pathogenic and saprophytic species recognized [3,4]. A significant challenge in the future will be to more precisely correlate these genetic differences with biological differences. The relationship between leptospiral genome content, pathogenesis

<figure><img src="image_2.png" /></figure>

<footnote>*E-mail: Ben.Adler@med.monash.edu.au</footnote>

<footnote n="9">⁹These authors contributed equally to this work.</footnote>
and the ability to survive in diverse environmental niches is a particularly important area of investigation, highlighted by our recent findings. A process of genome erosion and loss of gene function in *L. borgpetersenii* serovar Hardjo [5] limits its viability outside a mammalian host and likely impairs disease transmission through water, a route commonly used by *L. interrogans* to infect new hosts.

To gain insight into the genetic potential of *Leptospira* and to help identify genes that contribute to long-term survival in surface water, we determined the genome sequences of two *L. biflexa* strains, the first saprophytic *Leptospira* species to be characterized by genome analysis. Comparison of these data to genomes of the pathogenic species *L. borgpetersenii* [5], and *L. interrogans* [6,7] provides an opportunity to identify features that are unique to pathogenic and saprophytic species, thereby providing new experimental directions and novel perspectives on leptospiral evolution, environmental persistence and the causation of disease.

## Results and Discussion

*Leptospira* genomes vary in replicon content and genetic organization

The genome of *Leptospira biflexa* serovar Patoc strain Patoc1 (Ames strain) consists of three replicons with a total of 3,956,086 base pairs (bp) (Figure 1). The two larger replicons share extensive similarity to the two chromosomes that comprise the genomes of *L. borgpetersenii* and *L. interrogans* and are therefore referred to as chromosome I (CI; 3,603,977 bp; GC% 38.89) and chromosome II (CII; 277,995 bp; GC% 39.27). *L. biflexa* possesses a third circular replicon (74,114 bp; GC% 37.47), that we designate p74, not found in the previously sequenced pathogenic *Leptospira* species. (The complete genomic sequences of *L. biflexa* serovar Patoc strain Patoc1, strains Paris and Ames, have been deposited in GenBank under the Accession Numbers: CP000777, CP000778, CP000779, CP000786, CP000787 and CP000788).

A total of 3,590 protein-coding genes (CDSs) was identified in *L. biflexa* (Table 1). Most of these genes are located on CI, including two *rfi* genes, two *rrl* genes and two *rrs* genes, coding for 5S, 23S and 16S rRNA molecules respectively. These rRNA genes are not linked to

each other, a feature unusual among most bacteria but common among *Leptospira*. Similar to the slower-growing pathogenic species *L. interrogans* and *L. borgpetersenii*, each of which has a relatively low number (37) of transfer RNA (tRNA) genes, the faster-growing *L. biflexa* surprisingly has only 35 tRNA genes. This finding indicates that the growth rate of *Leptospira* is not restricted by the low number of tRNA and rRNA genes as previously suggested [7], but rather may be due to other differing metabolic capacities amongst *Leptospira* spp. As reported for other *Leptospira* spp. [8], essential genes such as *gltB* (glutamate synthase) and *asd* (aspartate semialdehyde dehydrogenase) are carried on CII in *L. biflexa*.

It is unclear if replicon III (p74) functions as a chromosome, i.e. carries genes essential for survival, or is an extrachromosomal element or plasmid. Thirteen genes on p74 have orthologs located on CI in pathogenic *Leptospira* (Table 2). For example, *recBCD* are located on CI in *L. interrogans* and *L. borgpetersenii*, but are located on p74 in *L. biflexa*. Mutation of these housekeeping genes in other bacterial species can affect transformation competence and viability [9], suggesting that p74 is essential for the survival of *L. biflexa*.

GC skew analysis (Figure 1) suggests that CI, CII and p74 are theta-type replicons that replicate bi-directionally from a unique origin. Each *L. biflexa* replicon encodes its own partition proteins from origin-proximal genes that may recognize and interact with a replication-specific binding site. The replication origin of CI resembles typical circular chromosomes from other bacteria, i.e. it is adjacent to *dnaA*, which encodes the initiator protein, and other genes such as *dnaN*, *recF* and *gyrAB*. In contrast, the replication origins of CII, replicon p74 and the previously reported leptophage LE1 [10] resemble phage and plasmid replicons; they contain both a partition operon and a downstream putative *rep* gene. For p74, the predicted replication origin based on GC skew analysis includes a partition operon (pLEPBI0001/LBF_5000 and pLEPBI0002/LBF_5001) and an adjacent gene, pLEPBI0003/LBF_5003, the product of which shares 63% similarity with the leptophage LE1 Rep protein. Regardless of the similarity between the LE1 Rep protein and the product of pLEPBI003, p74 and the LE1 prophage can co-exist (data not shown). Furthermore, we demonstrated that this region of p74 directs autonomous

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Circular maps of the three *L. biflexa* replicons. (1) the coordinates in bp beginning at 0 = oriC; (2) dark pink: genes unique to *L. biflexa*, not found in *L. interrogans* serovar Copenhageni and *L. borgpetersenii* serovar Hardjobovis (identity >40% over 80% of the length of the smallest protein). (3) dark purple: genes found in *L. biflexa*, *L. interrogans* and *L. borgpetersenii* (identity >40% over 80% of the length of the smallest protein). (4) red: genes found in *L. biflexa* and *L. borgpetersenii*, but not in *L. interrogans* (identity >40% over 80% of the length of the smallest protein). (5) brown: genes found in *L. biflexa* and *L. interrogans*, but not in *L. borgpetersenii* (identity >40% over 80% of the length of the smallest protein). (6) blue: genes found in *L. biflexa* and other sequenced spirochetes (*Borrelia afzelii* PKo, *Borrelia burgdorferi*, *Borrelia garinii*, *Treponema denticola* and *Treponema pallidum*) (identity >40% over 80% of the length of the smallest protein). (7) The innermost ring shows GC skew; positive skew is shown in grey, and negative skew is shown in black.
doi:10.1371/journal.pone.0001607.g001</figcaption></figure>

<figure><img src="image_4.png" /></figure>
<table><caption>**Table 1.** Summary of genome features of pathogenic and saprophytic *Leptospira*</caption><thead><tr><th rowspan="2">Features</th><th colspan="2">¥L. borgpetersenii</th><th colspan="2">¥L. interrogans</th><th colspan="2">¥L. biflexa</th><th rowspan="2">P74</th><th rowspan="2">LE-1 prophage<sup>b</sup></th></tr><tr><th>CI</th><th>CII</th><th>CI</th><th>CII</th><th>CI</th><th>CII</th></tr></thead><tbody><tr><td>Size (bp)</td><td>3,614,456</td><td>317,335</td><td>4,277,185</td><td>350,181</td><td>3,603,977</td><td>277,995</td><td>74,116</td><td>73,623</td></tr><tr><td>G+C content (%)</td><td>41.0</td><td>41.2</td><td>35.1</td><td>35.0</td><td>38.9</td><td>39.3</td><td>37.5</td><td>38.5</td></tr><tr><td>Protein-coding percentage</td><td>80</td><td>80</td><td>74.9</td><td>75.5</td><td>92.3</td><td>93.3</td><td>90.9</td><td>93.4</td></tr><tr><td><strong>Protein coding sequences</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>CDS<sup>a</sup></td><td>2,607</td><td>237</td><td>3,105</td><td>274</td><td>3,268</td><td>266</td><td>56</td><td>82</td></tr><tr><td>With assigned function</td><td>1,644</td><td>135</td><td>1,817</td><td>159</td><td>2,042</td><td>141</td><td>31</td><td>19</td></tr><tr><td>Conserved hypothetical</td><td>373</td><td>32</td><td>484</td><td>34</td><td>464</td><td>43</td><td>5</td><td>2</td></tr><tr><td>Unique hypothetical</td><td>590</td><td>70</td><td>804</td><td>81</td><td>762</td><td>82</td><td>20</td><td>61</td></tr><tr><td><strong>Transposases</strong></td><td>215</td><td>26</td><td>26</td><td>0</td><td>8</td><td>1</td><td>1</td><td>0</td></tr><tr><td><strong>Pseudogenes</strong></td><td>340</td><td>28</td><td>38</td><td>3</td><td>32</td><td>1</td><td>0</td><td>0</td></tr><tr><td><strong>Transfer RNA genes</strong></td><td>37</td><td>0</td><td>37</td><td>0</td><td>35</td><td>0</td><td>0</td><td>0</td></tr><tr><td><strong>Ribosomal RNA genes</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>23S</td><td>2</td><td>0</td><td>2</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td></tr><tr><td>16S</td><td>2</td><td>0</td><td>2</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td></tr><tr><td>5S</td><td>1</td><td>0</td><td>1</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td></tr></tbody></table>

¥L. borgpetersenii serovar Hardjo strain L550, L. interrogans serovar Copenhageni strain Fiocruz, L. biflexa serovar Patoc strain Ames

aexcluding transposases and pseudogenes

b[11,38]

doi:10.1371/journal.pone.0001607.t001

replication in *L. biflexa* (Figure S1). This will facilitate the construction of new shuttle vectors for the genetics of *Leptospira* and their use in co-transformation experiments with the *L. biflexa-E. coli* shuttle vector derived from LE1 [11].

We suggest that the CII and p74 replicons evolved from a stabilized circular intermediate of a progenitor phage related to LE1 [10]. Intermediates in this process may be related to the LaiGI-1 Genomic Island [12], which exists both as an element integrated into CI of *L. interrogans* serovar Lai or as an autonomously replicating plasmid. Stabilization of these intermediates could have occurred through illegitimate recombination resulting in incorporation of essential genes from CI into the

smaller replicons. Evidence for this model is provided by the presence of genes located on the CI replicon in pathogenic *Leptospira* on either CII or p74 in *L. biflexa*, as noted above.

The genus *Leptospira* is renowned for the stability of its agglutinating antigens during *in vitro* culture, with examples of strains maintaining serovar identity during more than 80 years of propagation [2], implying considerable genomic stability in the absence of selective pressure for antigenic change. Furthermore, more than 2% of the *L. interrogans* serovar Copenhageni genome is dedicated to LPS biosynthesis. In *L. biflexa* about 1.4% of the genome encodes LPS biosynthesis functions with essentially all relevant CDS encoded on the same strand, and these genes are

<table><caption>**Table 2.** CDS from Replicon III (p74) that have an ortholog in Chromosome I in other *Leptospira*</caption><thead><tr><th>Stop</th><th>start</th><th>locus_tag</th><th>ortholog_tag</th><th>product</th></tr></thead><tbody><tr><td>4267</td><td>3794</td><td>LBF_5005</td><td>SPN2759</td><td>Conserved hypothetical protein</td></tr><tr><td>5304</td><td>4993</td><td>LBF_5007</td><td>SPN2285</td><td>Conserved hypothetical protein</td></tr><tr><td>7485</td><td>5608</td><td>LBF_5009</td><td>SPN2142</td><td>Serine phosphatase RsbU, regulator of sigma subunit</td></tr><tr><td>10908</td><td>11909</td><td>LBF_5013</td><td>SPN2858</td><td>ABC-type Fe3+-siderophore transport system, permease component</td></tr><tr><td>16746</td><td>18548</td><td>LBF_5018</td><td>SPN2289</td><td>Exodeoxyribonuclease V, alpha subunit</td></tr><tr><td>18545</td><td>22168</td><td>LBF_5019</td><td>SPN2290</td><td>Exodeoxyribonuclease V, beta subunit</td></tr><tr><td>22171</td><td>25470</td><td>LBF_5020</td><td>SPN2291</td><td>Exodeoxyribonuclease V, gamma subunit</td></tr><tr><td>41270</td><td>42064</td><td>LBF_5030</td><td>SPN0228</td><td>Bacteriophage-related protein*</td></tr><tr><td>53226</td><td>52999</td><td>LBF_5037</td><td>SPN1718</td><td>Conserved hypothetical protein</td></tr><tr><td>60761</td><td>60456</td><td>LBF_5044</td><td>SPN3221</td><td>Antitoxin of toxin-antitoxin stability system</td></tr><tr><td>61204</td><td>60767</td><td>LBF_5045</td><td>SPN3222</td><td>Hypothetical protein</td></tr><tr><td>62047</td><td>63093</td><td>LBF_5047</td><td>SPN1129</td><td>Homoserine kinase</td></tr><tr><td>63112</td><td>64008</td><td>LBF_5048</td><td>SPN1151</td><td>GGDEF domain receiver component of a two-component response regulator</td></tr></tbody></table>

*Ortholog found on Chromosome II or Chromosome I in other *Leptospira*.

doi:10.1371/journal.pone.0001607.t002

<figure><img src="image_5.png" /></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Venn diagram showing numbers of unique and shared genes amongst *L. interrogans*, *L. borgpetersenii* and *L. biflexa*. Orthologous CDS were identified in a pair-wise fashion using Whole-Genome Reciprocal Best-Hit BLAST Analysis [37]. Manual curation ensured a one to one relationship for orthologous CDS, particularly in situations where sets of paralogous CDS existed and in addition evaluated the nature of the relationship between CDS with reciprocal best-hits but low expect values. This analysis was performed using the *L. interrogans* serovar Copenhageni strain Fiocruz, *L. borgpetersenii* serovar Hardjo strain L550 and *L. biflexa* serovar Patoc strain Ames genome sequences.
doi:10.1371/journal.pone.0001607.g002</figcaption></figure>

identical in the Ames and Paris strains. In contrast, non-LPS encoding regions in the two *L. biflexa* strains maintained separately for 17 years show some evidence of minor accumulation of changes. Comparison of the Paris and Ames strains of *L. biflexa* reveals four indels, three of which are insertions of ISLbi1 elements into the coding regions of genes in the Ames strain (LBF_0259, LBF_2295 and LBF_2512). The first two genes encode proteins of unknown function, and the third encodes a protein with a predicted role in lipid metabolism; the significance of these insertions is not known. In the Paris strain, the insertion element is restricted to one copy each on CII and p74; these copies are maintained in the Ames strain. The fourth difference occurs in CII in which an additional 250 bp are found in an intergenic region in the Ames strain. Aside from these differences, the genomes of the Paris and Ames strains are virtually identical; in this article the term *L. biflexa* will refer to the Ames strain unless stated otherwise.

Approximately two thirds of the genes in *L. biflexa* have orthologous genes in the pathogens *L. interrogans* and *L. borgpetersenii* (Figure 2), consistent with a common origin for leptospiral saprophytes and pathogens. The genes conserved across the three species are distributed in the two chromosomes of all *Leptospira* species and strains (Table 3). The high sequence identity of the small chromosome *parA* gene product (>87%) between all

*Leptospira* spp. suggests a shared ancestry for this replicon. Moreover, there is no difference in inter-chromosomal DNA parameters such as GC%, codon preferences and gene density within the three *L. biflexa* replicons, despite the relatively wide range of these parameters in pathogenic *Leptospira*. By contrast, codon preferences of genes of the leptospiral bacteriophage LE1 and its host, *L. biflexa* [10] are different (data not shown). These data suggest a long-standing relationship between the three *L. biflexa* replicons and perhaps a more recent acquisition of LE1.

## The evolution of the *Leptospira* species

Phylogenetic analysis based on comparison of 16S rRNA sequences indicates that members of the family *Leptospiraceae* form the deepest branch in spirochete evolution, with divergence of saprophytic and pathogenic *Leptospira* likely being the result of a single event [13]. Further diversification of species within each of these two evolutionary branches of *Leptospira* is supported by multilocus sequence analysis and DNA renaturation kinetics [3,14,15].

We used comparative genome analysis to help identify key features showing patterns of variation consistent with the action of selective evolutionary pressure within *Leptospira*. The *L. biflexa* and *L. borgpetersenii* genomes are similar in size (Table 1), but the gene density in *L. biflexa* is much higher, probably as a result of IS-mediated genome erosion in *L. borgpetersenii*. In contrast, the *L. interrogans* genome is larger, probably reflecting the added genetic information required for survival both within mammalian hosts and aquatic environments, whereas *L. biflexa* and *L. borgpetersenii* are restricted to aquatic and mammalian host environments, respectively. The dearth of IS-elements in the *L. biflexa* genome (five IS-elements) is in stark contrast to their abundance in the *L. interrogans* (36 and 69 IS-elements in the Fiocruz L1-130 and Lai 56601 strains, respectively) and *L. borgpetersenii* (167 IS-elements) genomes. The presence of large numbers of IS-elements is an indicator of genome plasticity in *Leptospira* species. Taken together, our results suggest that the *L. biflexa* gene order is more likely to have a closer relationship to the progenitor genome for the genus *Leptospira*. There is low syntenity between the sequenced pathogenic *Leptospira* (Figure 3) [5], despite the short evolutionary distance separating them.

Gene and functional redundancy is also more common in the pathogens in comparison to the saprophytes (Table 4). The pathogens have more paralogs (excluding transposases), with 203 (5.01%) and 438 (10.52%) paralogs for *L. interrogans* serovar Copenhageni Fiocruz L1-130 and *L. borgpetersenii* serovar Hardjo L550, respectively, compared to 62 paralogs (1.65%) in the saprophyte *L. biflexa*. We identified 43, 67, and 53 loci with gene duplication in the genomes of *L. biflexa*, *L. interrogans*, and *L. borgpetersenii*, respectively. The genome of *L. biflexa* therefore appears to be more stable than those of *L. interrogans* and *L. borgpetersenii*. Moreover, the *L. biflexa* genome (92%) has a greater

<table><caption>**Table 3.** Distribution of the orthologs over the two chromosomes of *Leptospira* spp.</caption><thead><tr><th></th><th><i>L. biflexa</i> strain Patoc1</th><th><i>L. interrogans</i> strain Fiocruz L1-139</th><th><i>L. interrogans</i> strain Lai 56601</th><th><i>L. borgpetersenii</i> strain L550</th><th><i>L. borgpetersenii</i> strain JB197</th></tr></thead><tbody><tr><td>No of CDS shared between C1 replicons (1)</td><td>1411 (41.58%)</td><td>1429 (39.13%)</td><td>1429 (37.41%)</td><td>1482 (38.73%)</td><td>1448(38.20%)</td></tr><tr><td>No of CDS shared between CII replicons (2)</td><td>80 (28.46%)</td><td>83 (27.66%)</td><td>83 (26.17%)</td><td>82 (24.47%)</td><td>82 (25.07%)</td></tr></tbody></table>

<figure><img src="image_7.png" /></figure>

<footnote>(1) Number of CDS (orthologs) found in the large chromosome (CII) of one leptospire that are also found in the large chromosomes of the other four leptospiral large chromosomes. (2) Number of CDS (orthologs) found in the small chromosome (CII) of one leptospire that are also found in the small chromosomes of the other four leptospiral large chromosomes. Putative orthologous relations between two genomes are defined as gene couples satisfying the bi-directional best hit (BBH) criterion or a blastP alignment threshold, a minimum of 40% sequence identity on 80% of the length of the smallest protein. CDS, coding regions.
doi:10.1371/journal.pone.0001607.t003</footnote>
A

<figure><img src="image_8.png" /></figure>

B

<figure><img src="image_9.png" /></figure>

C

<figure><img src="image_10.png" /></figure>

**Figure 3. Synteny plot between the five Leptospira genomes.** The line plots were obtained using synteny results between the large CI(A) or small CII(B) chromosomes of *L. biflexa* serovar Patoc strain Patoc1, *L. interrogans* serovar Lai strain 56601, *L. interrogans* serovar Copenhageni strain Fiocruz L1-130, *L. borgpetersenii* serovar Hardjo strain L550, and *L. borgpetersenii* serovar Hardjo strain JB197. A line plot (C) compares synteny between *L. borgpetersenii* serovar Hardjo strain JB197 and *L. interrogans* serovar Copenhageni strain Fiocruz L1-130. Comparative analysis was performed using the MaGe interface [35] in the SpiroScope database (https://www.genoscope.cns.fr/agc/mage). The minimum size of the synteny groups is set to five genes. In green: synteny groups are organized on the same strand; in red: synteny groups are organized on opposite strands. doi:10.1371/journal.pone.0001607.g003

<figure><img src="image_11.png" /></figure>
<table><caption>**Table 4.** Comparative genomics of *Leptospira* spp. Putative orthologous relations between two genomes are defined as gene couples satisfying the bi-directional best hit (BBH) criterion or a blastP alignment threshold, a minimum of 40% sequence identity on 80% of the length of the smallest protein. Putative paralogous relations between two genomes are defined as gene couples satisfying the bi-directional best hit (BBH) criterion or a blastP alignment threshold, a minimum of 60% sequence identity on 80% of the length of the smallest protein.</caption><thead><tr><th>Paralogs/orthologs</th><th>L. biflexa strain<br>Patoc 1</th><th>L. interrogans strain<br>Fiocruz L1-130</th><th>L. borgpetersenii strain<br>L550</th><th>L. borgpetersenii strain<br>JB197</th><th>L. interrogans strain<br>Lai 56601</th></tr></thead><tbody><tr><td>L. biflexa strain Patoc 1</td><td>62 (1.5%)</td><td>1650 (44.04%)</td><td>1635 (43.64%)</td><td>1631 (43.53%)</td><td>1652 (44.10%)</td></tr><tr><td>L. interrogans strain Fiocruz L1-130</td><td>1633 (40.35%)</td><td>203 (5.01%)</td><td>2913 (71.97%)</td><td>2907 (71.83%)</td><td>3745 (92.53%)</td></tr><tr><td>L. borgpetersenii strain L550</td><td>1674 (40.23%)</td><td>3084 (74.11%)</td><td>438 (10.52%)</td><td>3936 (94.59%)</td><td>3077 (73.94%)</td></tr><tr><td>L. borgpetersenii strain JB197</td><td>1636 (39.73%)</td><td>3052 (74.13%)</td><td>3922 (95.26%)</td><td>354 (8.59%)</td><td>3044 (73.93%)</td></tr><tr><td>L. interrogans strain Lai 56601</td><td>1638 (39.60%)</td><td>3778 (91.34%)</td><td>2942 (71.13%)</td><td>2934 (70.93%)</td><td>204 (4.93%)</td></tr></tbody></table>

doi:10.1371/journal.pone.0001607.t004

gene density than the *L. borgpetersenii* (80%) and *L. interrogans* genomes (75%) (Table 1). Greater gene density can contribute to a relatively stable gene order. For example, the spirochete *Borrelia* has a gene density of more than 92% and the *B. garinii* and *B. burgdorferi* genomes are essentially collinear [16].

## Genetic determinants involved in survival in the environment

The spectrum of ecological niches occupied by diverse *Leptospira* species raises questions as to how the capacity to survive in these diverse environments has evolved within different species in the genus *Leptospira*. We propose that their common progenitor had a genome more like *L. biflexa*. Subsequently, acquisition of genes that enabled *Leptospira* to infect mammals would have expanded the range of environments that it could successfully occupy. Loss of environmental survival genes then would lead to dependence on a mammalian host and eventually return the genome to a smaller size (i.e. *L. borgpetersenii*), consistent with our model for IS-mediated genome degradation [5].

The extensive repertoire of genes (Table 5) encoding proteins involved in signal transduction in *L. biflexa* (287 CDS) compared with *L. interrogans* (214 CDS) and *L. borgpetersenii* (167 CDS) is consistent with an enhanced metabolic capability in *L. biflexa* reflected by its environmental habitat and likely contributes to its enhanced growth rate relative to the pathogens. An analogous situation is seen in *Mycobacterium*, where *M. marinum* which occupies both animal and aquatic environments has many more genes encoding environmental sensing and metabolic proteins than the closely related, obligate human pathogen *M. tuberculosis* [17,18].

The genomes of *Leptospira* spp. contain a number of genes involved in the production of exopolysaccharides. These genes, such as those encoding glycosyltransferases, alginate biosynthesis, and lipopolysaccharide transport systems, may contribute to colonization of both biotic and abiotic surfaces. When investigating the formation of biofilms on solid surfaces, we have observed the production of a strong biofilm by *L. biflexa* and *L. interrogans* (Unpublished results). The formation of such a biofilm is consistent with the life of saprophytic species in water; it may facilitate *L. biflexa* occupying particular environmental niches. The presence of biofilms may also play an important role in chronic carriage of the pathogen *L. interrogans* in animal reservoirs by facilitating colonization of the renal tubules. Interestingly, genes involved in alginate biosynthesis are present in both *L. biflexa* (11 genes) and *L. interrogans* (8 genes), but are absent in *L. borgpetersenii*, a finding consistent with the reduced environmental survival of *L. borgpetersenii* [5].

## Comparative genomics of pathogenic and saprophytic *Leptospira*

The molecular mechanisms of pathogenesis in leptospirosis remain almost entirely unknown. The set of leptospiral, pathogen-specific genes, defined as those with no orthologous gene in *L. biflexa*, is likely to be enriched for genes that play a role in pathogenesis. Significantly, the majority (893) of these 1,431 genes has no known function (Table S1abc), suggesting the presence of pathogenic mechanisms unique to *Leptospira*. Among the remaining 538 genes, there are several genes encoding response regulators and environment sensing proteins; these genes are likely to represent adaptations that enable survival in environments not encountered by the saprophytic strain. There is an expansion of genes encoding leucine-rich repeat (LLR) proteins from one gene in *L. biflexa* to eight and 18 genes in *L. borgpetersenii* and *L. interrogans*, respectively. Although these LLR proteins have no obvious function, in *Treponema denticola* the LRR protein LrrA appears to have roles in attachment and to, and penetration of, host tissues [19]; the diversity of LRR proteins in pathogenic *Leptospira* may be important for successful infection of a wide variety of mammalian host species.

The regulation of transcription differs significantly between the saprophytic and pathogenic strains as indicated by the presence of more than 50 saprophyte-specific, and more than 20-pathogen specific, transcription regulators, whereas there are only 27 transcription regulators that are common to all of the sequenced *Leptospira* species.

Because the bacterial surface is an interface between the pathogen and the host, any differences in cell-surface proteins might reflect variation in pathogenesis mechanisms among *Leptospira* spp., which contain a relatively low repertoire of trans-membrane proteins in their outer membranes [20]. However, the outer membranes contain a predominance of lipoproteins, that may be either surface-exposed or located in the periplasm. Predicted lipoproteins are prominent in both saprophytic (164 predicted lipoproteins) and pathogenic *Leptospira* (184 and 130 predicted lipoproteins in *L. interrogans* and *L. borgpetersenii* respectively). Despite these similarities, there are significant differences, most notably the absence in *L. biflexa* of an ortholog of the major outer membrane lipoprotein, LipL32 [20]. Moreover, 89 *L. biflexa* lipoprotein genes have no orthologs in other *Leptospira* species, and more than 90 lipoproteins from the pathogenic species have no orthologs in the *L. biflexa* genome. In addition to LipL32, other characterized lipoproteins that have no orthologs in *L. biflexa* include LipL41, LipL36 and several LipL45 related proteins [20]. However, LipL21, which was reported as a pathogen-specific

<figure><img src="image_12.png" /></figure>
<table><caption>**Table 5.** Distribution of general protein functions between leptospiral species based on the COG function classification scheme.</caption><thead><tr><th>*COG Function Classification</th><th>¥L. biflexa</th><th>¥L. borgpetersenii</th><th>¥L. interrogans</th></tr></thead><tbody><tr><td colspan="4"><strong>INFORMATION STORAGE AND PROCESSING</strong></td></tr><tr><td>J Translation, ribosomal structure and biogenesis</td><td>154</td><td>174</td><td>153</td></tr><tr><td>K Transcription</td><td>166</td><td>104</td><td>109</td></tr><tr><td>L Replication, recombination and repair</td><td>94</td><td>91</td><td>102</td></tr><tr><td>B Chromatin structure and dynamics</td><td>2</td><td>2</td><td>2</td></tr><tr><td colspan="4"><strong>CELLULAR PROCESSES AND SIGNALING</strong></td></tr><tr><td>D Cell cycle control, cell division, chromosome partitioning</td><td>21</td><td>22</td><td>22</td></tr><tr><td>V Defense mechanisms</td><td>39</td><td>32</td><td>37</td></tr><tr><td>T Signal transduction mechanisms</td><td>287</td><td>167</td><td>214</td></tr><tr><td>M Cell wall/membrane/envelope biogenesis</td><td>230</td><td>199</td><td>218</td></tr><tr><td>N Cell motility</td><td>93</td><td>84</td><td>89</td></tr><tr><td>U Intracellular trafficking, secretion, and vesicular transport</td><td>71</td><td>73</td><td>71</td></tr><tr><td>O Posttranslational modification, protein turnover, chaperones</td><td>105</td><td>96</td><td>100</td></tr><tr><td colspan="4"><strong>METABOLISM</strong></td></tr><tr><td>C Energy production and conversion</td><td>132</td><td>115</td><td>119</td></tr><tr><td>G Carbohydrate transport and metabolism</td><td>91</td><td>76</td><td>91</td></tr><tr><td>E Amino acid transport and metabolism</td><td>163</td><td>136</td><td>150</td></tr><tr><td>F Nucleotide transport and metabolism</td><td>46</td><td>52</td><td>52</td></tr><tr><td>H Coenzyme transport and metabolism</td><td>119</td><td>112</td><td>120</td></tr><tr><td>I Lipid transport and metabolism</td><td>101</td><td>83</td><td>99</td></tr><tr><td>P Inorganic ion transport and metabolism</td><td>120</td><td>72</td><td>88</td></tr><tr><td>Q Secondary metabolites biosynthesis, transport and catabolism</td><td>35</td><td>23</td><td>27</td></tr><tr><td colspan="4"><strong>POORLY CHARACTERIZED</strong></td></tr><tr><td>R General function prediction only</td><td>311</td><td>237</td><td>294</td></tr><tr><td>S Function unknown</td><td>174</td><td>157</td><td>192</td></tr><tr><td>CDS Not Classified (not related to any COG)</td><td>1,266</td><td>931</td><td>1,245</td></tr><tr><td>Total CDS (count excludes transposases and pseudogenes)</td><td>3,590</td><td>2,843</td><td>3,378</td></tr></tbody></table>

lipoprotein [21] based on antibody reactivity, has an ortholog in *L. biflexa* with 50% similarity. This level of similarity is consistent with a different function for the LipL21 orthologs in the two *Leptospira* species.

Several putative virulence factors previously identified in pathogenic *Leptospira* spp. are not present in *L. biflexa*, including the Lig surface proteins containing immunoglobulin-like repeats predicted to play a role in the adhesion to host tissues [22]. Similarly, LfhA, a putative factor H binding protein [23] that has also been shown to bind the extracellular matrix protein laminin [24], is shared among pathogenic *Leptospira*, but is lacking in *L. biflexa*. Although the genome of *L. biflexa* contains putative hemolysins [25], its genome is devoid of genes encoding enzymes capable of degrading tissues, such as the range of sphingomyelinases found in pathogenic species [26,27,28]. The role of sphingomyelinases in the pathogenesis of leptospirosis has been controversial; are they key virulence factors or do they merely play a role in nutrient acquisition? Their absence in *L. biflexa* strongly supports their involvement in survival within mammalian hosts. Interestingly, the membrane protein, Loa22, the only protein to date that has been shown genetically to be required for virulence in *L. interrogans* [29], has a *L. biflexa* ortholog with 73% similarity.

Its role in either pathogenic or saprophytic species is unknown, but its presence in the saprophytic species suggests that it is involved in survival rather than being a direct virulence factor and is consistent with the common progenitor hypothesis.

*L. interrogans* is the most frequently reported agent of human leptospirosis. The disease is also generally more severe with *L. interrogans* than with *L. borgpetersenii* [2]. On this basis we propose that the subset of *L. interrogans* genes that have no orthologs in either *L. biflexa* or *L. borgpetersenii* may contain virulence factors that are responsible for the more severe form of leptospirosis. Other subsets that may be enriched for genes involved in particular biological functions include those genes that have orthologs in *L. biflexa* and *L. interrogans* and not *L. borgpetersenii* which may contain genes involved in survival outside the animal host. While the loss of many genes from the *L. borgpetersenii* genome has occurred through genome reduction [5], the presence of 265 unique genes in *L. borgpetersenii* (Figure 2; Table S1abc) indicates that these genomes have also gained some additional genes during the course of their evolution. The lateral acquisition of genetic material is often associated with IS-elements. While there is an association between IS-elements and genes that are isolate- or species-specific, in particular *in L. interrogans*, no mechanism has been determined.

<figure><img src="image_13.png" /></figure>

<footnote>*Each COG assignment has been manually curated to ensure consistent classification across orthologous proteins. A feature of the COG scheme is that some COGs have multiple functional classifications.</footnote>

<footnote>¥L. borgpetersenii serovar Hardjo strain L550, L. interrogans serovar Copenhageni strain Fiocruz, L. biflexa serovar Patoc strain Ames</footnote>

<footnote>doi:10.1371/journal.pone.0001607.t005</footnote>
Other possible mechanisms for horizontal gene acquisition may include the involvement of bacteriophage, as integrated prophages or bacteriophage remnants are present.

## The Core Leptospiral Genes

Saprophytic and pathogenic *Leptospira* species belong to two distinct phylogenetic groups, leading us to conclude that the 2,052 genes shared by both groups constitute the core genome of this genus (Figure 2). As expected, many of the functional categories that are involved in essential housekeeping functions, such as DNA and RNA metabolism, protein processing and secretion, cell structure, cellular processes, and energetic and intermediary metabolism, are represented in the core gene set.

The presence of orthologous genes in all the sequenced leptospiral species is an indicator that these genes were acquired prior to the radiation of the genus, and as such, is a strong indicator that these genes have not been laterally acquired. The substantial proportion of the total leptospiral genes that are in this category is perhaps an indicator that lateral transfer of genes into the genus *Leptospira* is a minor contributor to the overall genetic composition of the genus and an indicator that the genus has undergone an extended period of 'genetic isolation'. This notion is supported by the fact that approximately 20% of the core leptospiral genes are unique to the genus. However, this does not preclude horizontal acquisition of some genes or gene clusters.

## Implications for genetic studies

Despite the development of basic tools such as transposon mutagenesis for the pathogenic species of *Leptospira*, targeted gene inactivation is not yet possible. Therefore, there is likely to be continued interest in the use of *L. biflexa* as a model bacterium for genetic analysis [30]. Knowledge of the distribution of orthologous genes in *L. biflexa* will be an important resource for the elucidation of function for genes common to pathogenic and saprophytic species.

## Materials and Methods

### Sequencing and annotation of the genome of *L. biflexa*

The strain, *L. biflexa* serovar Patoc strain Patoc1, was initially isolated from stream water [31], maintained in the collection of the National Reference Center of *Leptospira* (Institut Pasteur, Paris, France) and designated the Paris strain. A second strain was derived from the same source but kept in the culture collection at the National Animal Disease Center (NADC), Ames, IA since 1990 and is designated the Ames strain. Each strain was colony purified before growth for genomic DNA isolation. For each isolate, genomic DNA was randomly sheared by nebulization (HydroShear, GeneMachines) and the ends repaired enzymatically. Small fragments (~1.5–4 kb) were ligated either to a derivative of plasmid pGEM7Zf+ (Promega) (Paris strain) or pSMART-HC and pSMART-LC (Lucigen Corp.) (Ames strain). Large (35–45 kb) DNA fragments were ligated to fosmid pCC1FOS (Epicentre, Madison, WI) (Paris strain). Intermediate sized DNA fragments (7–12 kb) were prepared from partially BamHI-digested Ames strain DNA and ligated into pZERO-1 (Invitrogen). Plasmid DNA preparation was performed with the TempliPhi DNA sequencing template amplification kit (GE Healthcare -Biosciences) or SprintPrep plasmid preparation kits (Agencourt Bioscience). Fosmid DNA purification was performed with the Montage BAC Miniprep96 Kit (Millipore) [25]. Sequencing reactions were performed from both ends of DNA templates using ABI PRISM BigDye Terminator cycle sequencing ready reactions kits and run on a 3700 or a 3730 xl Genetic Analyzer (Applied Biosystems) at the Genomics platform (Pasteur Genopole

Île-de-France), the Australian Genome Research Facility, Brisbane, Australia or at the Genomics Facility at the NADC. Base-calls from sequence data were made using Phred [32]. Sequences not meeting our production quality criteria (at least 100 bases called with a quality over 20) were discarded. The traces were assembled using Phrap and Consed [33].

Whole genome shotgun sequencing was performed until approximately 6× genome coverage was achieved. Autofinish [34] was used to design primers for improving regions of low quality sequence and for primer walking along templates that spanned the gaps between contigs. Several strategies were used to orientate contigs to enable directed PCR-based approaches to span gaps between contigs. These strategies included Blast-based approaches that identified contig ends with hits to the same gene or to genes within the same locus, provided that there was conservation of gene order in both the *L. interrogans* and the *L. borgpetersenii* genomes. Repeated sequences longer than 600 bases, such as IS-elements and the *rrl* and *rn* genes, were identified and curated manually to ensure that at least two templates spanned these regions, thus confirming the assembly of these regions. At this stage, single, circularised contigs representing CII and p74 were identified. Combinational PCR was used to close the gaps between the final contigs making up CI. Outward-directed primers were designed for each of the contig ends; the primer sequences were subsequently checked and confirmed to be unique on the genome. The combinational PCR process required approximately 600 PCR reactions pairing each of the primers. In each instance where a gap was spanned, the size of the PCR product was less than 2 kb. The sequence of these PCR products was determined and added to the assembly to enable the closing of CI. PCR with independent primers was used to confirm the joins determined by the combinational PCR. In addition, for the Paris strain five fosmid clones were completely sequenced by transposon-assisted sequencing (Finnzymes, TGS II kit). Seven fosmid clones were also selected by Southern hybridization and sequenced in order to verify the 23S RNA and 16S RNA regions and 4 repeat regions. Validation of the final assembly was achieved by comparison of *in silico* digestion patterns with macrorestriction patterns obtained by PFGE with *Not*I and *Asc*I. For each genome, the error rate was less than 1 error per 10,000 bp.

The complete genome sequence was obtained from 58,663 and 40,260 sequences for the Paris and Ames strains respectively (giving >8× and 7× coverage). For the Paris strain, coding sequences (CDSs) likely to encode proteins were predicted with the AMIGene software [35]. Annotation was performed as described previously [25] using the MaGe annotation platform [36]. All the data were stored in SpiroScope, a relational database which is publicly available (http://www.genoscope.cns.fr/agc/mage). The Ames strain was annotated as described previously for the *L. borgpetersenii* serovar Hardjo genome using the Wasabi interactive platform [5].

## Supporting Information

### Figure S1

Found at: doi:10.1371/journal.pone.0001607.s001 (0.05 MB DOC)

### Table S1 Supplementary Table 1 a, b and c

Found at: doi:10.1371/journal.pone.0001607.s002 (0.46 MB XLS)

## Acknowledgments

We thank David Alt, Richard Hornsby, A. Magnier and A. Blanchard for technical support.

<figure><img src="image_14.png" /></figure>
## Author Contributions

Conceived and designed the experiments: BA CB MP RZ DB PW EK AM CB-E. Performed the experiments: CM ZR NZ SC SB JD MJ AL.

## References

1. Trueba G, Zapata S, Madrid K, Cullen P, Haake D (2004) Cell aggregation: a mechanism of pathogenic *Leptospira* to survive in fresh water. International Microbiology 7: 35–40.

2. Faine S, Adler B, Bolin C, Perolat P (1999) *Leptospira* and Leptospirosis. Melbourne, Australia: MediSci.

3. Brenner DJ, Kaufmann AF, Sulzer KR, Steigerwalt AG, Rogers FC, et al. (1999) Further determination of DNA relatedness between serogroups and serovars in the family *Leptospiraceae* with a proposal for *Leptospira alexanderi* sp. nov. and four new *Leptospira* genospecies. International Journal of Systematic Bacteriology 49: 839–858.

4. Levett PN (2001) Leptospirosis. Clinical Microbiology Reviews 14: 296–326.

5. Bulach DM, Zuerner RL, Wilson P, Seemann T, McGrath A, et al. (2006) Genome reduction in *Leptospira borgpetersenii* reflects limited transmission potential. Proceedings of the National Academy of Sciences of the United States of America 103: 14560–14565.

6. Nascimento ALTO, Ko AI, Martins EAL, Monteiro-Vitorello CB, Ho PL, et al. (2004) Comparative genomics of two *Leptospira interrogans* serovars reveals novel insights into Physiology and Pathogenesis. Journal of Bacteriology 186: 2164–2172.

7. Ren SX, Fu G, Jiang XG, Zeng R, Miao YG, et al. (2003) Unique physiological and pathogenic features of *Leptospira interrogans* revealed by whole-genome sequencing. Nature 422: 888–893.

8. Zuerner RL, Herrmann JL, Saint Girons I (1993) Comparison of genetic maps for two *Leptospira interrogans* serovars provides evidence for two chromosomes and intraspecies heterogeneity. Journal of Bacteriology 175: 5445–5451.

9. Kickstein E, Harms K, Wackernagel W (2007) Deletions of *recBCD* or *recD* influence genetic transformation differently and are lethal together with a *recJ* deletion in *Acinetobacter baylyi*. Microbiology 153: 2259–2270.

10. Bourhy P, Frangeul L, Couve E, Glaser P, Saint Girons I, et al. (2005) Complete nucleotide sequence of the LE1 prophage from the spirochete *Leptospira biflexa* and characterization of Its replication and partition functions. Journal of Bacteriology 187: 3931–3940.

11. Saint Girons I, Bourhy P, Ottone C, Picardeau M, Yelton D, et al. (2000) The LE1 bacteriophage replicates as a plasmid within *Leptospira biflexa*: construction of an *L. biflexa-Escherichia coli* shuttle vector. Journal of Bacteriology 182: 5700–5705.

12. Bourhy P, Salaun L, Lajus A, Medigue C, Boursaux-Eude C, et al. (2007) A Genomic island of the pathogen *Leptospira interrogans* serovar Lai can excise from its chromosome. Infection and Immunity 75: 677–683.

13. Paster BJ, Dewhirst FE, Weisburg WG, Tordoff LA, Fraser GJ, et al. (1991) Phylogenetic analysis of the spirochetes. Journal of Bacteriology 173: 6101–6109.

14. Haake DA, Suchard MA, Kelley MM, Dundoo M, Alt DP, et al. (2004) Molecular evolution and mosaicism of leptospiral outer membrane proteins involves horizontal DNA transfer. Journal of Bacteriology 186: 2818–2828.

15. Ramadass P, Jarvis BD, Corner RJ, Penny D, Marshall RB (1992) Genetic characterization of pathogenic *Leptospira* species by DNA hybridization. International Journal of Systematic Bacteriology 42: 215–219.

16. Glockner G, Lehmann R, Romualdi A, Pradella S, Schulte-Spechtel U, et al. (2004) Comparative analysis of the *Borrelia garinii* genome. Nucleic Acids Research 32: 6038–6046.

17. Stamm LM, Brown EJ (2004) *Mycobacterium marinum*: the generalization and specialization of a pathogenic mycobacterium. Microbes and Infection 6: 1418–1428.

18. Stinear TP, Seemann T, Pidot S, Frigui W, Reysset G, et al. (2007) Reductive evolution and niche adaptation inferred from the genome of *Mycobacterium ulcerans*, the causative agent of Buruli ulcer. Genome Research 17: 192–200.

Analyzed the data: RC BA JR JD MP RZ DB AM TS. Wrote the paper: BA JR MP RZ DB.

19. Ikegami A, Honma K, Sharma A, Kuramitsu H (2004) Multiple functions of the leucine-rich repeat protein LrrA of *Treponema denticola*. Infection and Immunity 72: 4619–4627.

20. Cullen PA, Haake DA, Adler B (2004) Outer membrane proteins of pathogenic spirochetes. FEMS Microbiology Reviews 28: 291–318.

21. Cullen PA, Haake DA, Bulach DM, Zuerner RL, Adler B (2003) LipL21 Is a novel surface-exposed lipoprotein of pathogenic *Leptospira* species. Infection and Immunity 71: 2414–2421.

22. Choy HA, Kelley MM, Chen TL, Moller AK, Matsunaga J, et al. (2007) Physiological osmotic induction of *Leptospira interrogans* adhesion: LigA and LigB bind extracellular matrix proteins and fibrinogen. Infection and Immunity 75: 2441–2450.

23. Verma A, Hellwage J, Artiushin S, Zipfel PF, Kraiczy P, et al. (2006) LfhA, a novel Factor H-binding protein of *Leptospira interrogans*. Infection and Immunity 74: 2659–2666.

24. Barbosa AS, Abreu PAE, Neves FO, Atzingen MV, Watanabe MM, et al. (2006) A newly identified leptospiral adhesin mediates attachment to laminin. Infection and Immunity 74: 6356–6364.

25. Louvel H, Bommezzadri S, Zidane N, Boursaux-Eude C, Creno S, et al. (2006) Comparative and functional genomic analyses of iron transport and regulation in *Leptospira* spp. Journal of Bacteriology 188: 7893–7904.

26. Segers RP, van Gestel JA, van Eys GJ, van der Zeijst BA, Gaastra W (1992) Presence of putative sphingomyelinase genes among members of the family *Leptospiraceae*. Infection and Immunity 60: 1707–1710.

27. Segers RP, van der Drift A, de Nijs A, Corcione P, van der Zeijst BA, et al. (1990) Molecular analysis of a sphingomyelinase C gene from *Leptospira interrogans* serovar hardjo. Infection and Immunity 58: 2177–2185.

28. del Real G, Segers RP, van der Zeijst BA, Gaastra W (1989) Cloning of a hemolysin gene from *Leptospira interrogans* serovar hardjo. Infection and Immunity 57: 2588–2590.

29. Ristow P, Bourhy P, McBride F, via Weykamp da C, Figueira C, et al. (2007) The OmpA-like protein Loa22 Is essential for leptospiral virulence. PLoS Pathogens 3: e97.

30. Louvel H, Picardeau M (2007) Genetic manipulation of *Leptospira biflexa*. In: Coico R, Kowalik TF, Quarles JM, Stevenson B, Taylor R, eds. Current Protocols in Microbiology. Hoboken: Wiley.

31. Babudieri B (1961) Studio serologico del gruppo Semaranga-Patoc di Leptospira biflexa.; Cagliari-Sassari. pp 408–414.

32. Ewing B, Hillier L, Wendl MC, Green P (1998) Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Research 8: 175–185.

33. Gordon D, Abajian C, Green P (1998) Consed: a graphical tool for sequence finishing. Genome Research 8: 195–202.

34. Gordon D, Desmarais C, Green P (2001) Automated finishing with Autofinish. Genome Research 11: 614–625.

35. Bocs S, Crueveiller S, Vallenet D, Nuel G, Medigue C (2003) AMIGene: Annotation of Microbial Genes. Nucleic Acids Research 31: 3723–3726.

36. Vallenet D, Labarre I, Rouy Z, Barbe V, Bocs S, et al. (2006) MaGe - a microbial genome annotation system supported by synteny result. Nucleic Acids Research 34: 53–65.

37. Fuchsman CA, Rocab G (2006) Whole-genome reciprocal BLAST analysis reveals that *Planctomycetes* do not share an unusually large number of genes with *Eukarya* and *Archaea*. Applied and Environmental Microbiology 72: 6841–6844.

38. Saint Girons I, Margarita D, Amouriaux P, Baranton G (1990) First isolation of bacteriophages for a spirochaete: potential genetic tools for *Leptospira*. Research in Microbiology 141: 1131–1138.

<figure><img src="image_15.png" /></figure>