# What Makes a Bacterial Species Pathogenic?: Comparative Genomic Analysis of the Genus Leptospira.

Derrick E Fouts, Michael A Matthias, Haritha Adhikarla, Ben Adler, Luciane

Amorim-Santos, Douglas E Berg, Dieter Bulach, Alejandro Buschiazzo,
Yung-Fu Chang, Renee L Galloway, et al.

► To cite this version:

Derrick E Fouts, Michael A Matthias, Haritha Adhikarla, Ben Adler, Luciane Amorim-Santos, et al. What Makes a Bacterial Species Pathogenic?: Comparative Genomic Analysis of the Genus Leptospira.. PLoS Neglected Tropical Diseases, 2016, 10 (2), pp.e0004403. 10.1371/journal.pntd.0004403. pasteur-01436457

HAL Id: pasteur-01436457

https://pasteur.hal.science/pasteur-01436457v1

Submitted on 16 Apr 2019

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>
<figure><img src="image_2.png" /></figure>

CrossMark
click for updates

RESEARCH ARTICLE

# What Makes a Bacterial Species Pathogenic?: Comparative Genomic Analysis of the Genus *Leptospira*

Derrick E. Fouts¹*, Michael A. Matthias², Haritha Adhikarla³, Ben Adler⁴, Luciane Amorim-Santos³,⁵, Douglas E. Berg², Dieter Bulach⁶, Alejandro Buschiazzo⁷,⁸, Yung-Fu Chang⁹, Renee L. Galloway¹⁰, David A. Haake¹¹,¹², Daniel H. Haft¹, Rudy Hartskeerl¹³, Albert I. Ko³,⁵, Paul N. Levett¹⁴, James Matsunaga¹¹,¹², Ariel E. Mechaly⁷, Jonathan M. Monk¹⁵, Ana L. T. Nascimento¹⁶,¹⁷, Karen E. Nelson¹, Bernhard Palsson¹⁵, Sharon J. Peacock¹⁸, Mathieu Picardeau¹⁹, Jessica N. Ricaldi²⁰, Janjira Thaipandungpanit²¹, Elsio A. Wunder, Jr.³,⁵, X. Frank Yang²², Jun-Jie Zhang²², Joseph M. Vinetz²,²⁰,²³*

1 J. Craig Venter Institute, Rockville, Maryland, United States of America, 2 Division of Infectious Diseases, Department of Medicine, University of California San Diego School of Medicine, La Jolla, California, United States of America, 3 Department of Epidemiology of Microbial Diseases, Yale School of Public Health, New Haven, Connecticut, United States of America, 4 Australian Research Council Centre of Excellence in Structural and Functional Microbial Genomics, Department of Microbiology, Monash University, Clayton, Australia, 5 Centro de Pesquisas Gonçalo Moniz, Fundação Oswaldo Cruz/MS, Salvador, Bahia, Brazil, 6 Victorian Bioinformatics Consortium, Monash University, Clayton, Victoria, Australia, 7 Institut Pasteur de Montevideo, Laboratory of Molecular and Structural Microbiology, Montevideo, Uruguay, 8 Institut Pasteur, Department of Structural Biology and Chemistry, Paris, France, 9 Department of Population Medicine & Diagnostic Sciences, College of Veterinary Medicine, Cornell University, Ithaca, New York, United States of America, 10 Centers for Disease Control and Prevention (DHHS, CDC, OID, NCEZID, DHCP, BSPB), Atlanta, Georgia, United States of America, 11 VA Greater Los Angeles Healthcare System, Los Angeles, California, United States of America, 12 David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America, 13 WHO/FAO/OIE and National Collaborating Centre for Reference and Research on Leptospirosis, KIT Biomedical Research, Royal Tropical Institute (KIT), Amsterdam, The Netherlands, 14 Government of Saskatchewan, Disease Control Laboratory Regina, Canada, 15 Department of Bioengineering, University of California, San Diego, La Jolla, California, United States of America, 16 Centro de Biotecnologia, Instituto Butantan, São Paulo, SP, Brazil, 17 Programa Interunidades em Biotecnologia, Instituto de Ciências Biomédicas, USP, São Paulo, SP, Brazil, 18 Department of Medicine, University of Cambridge, Cambridge, United Kingdom, 19 Institut Pasteur, Biology of Spirochetes Unit, National Reference Centre and WHO Collaborating Center for Leptospirosis, Paris, France, 20 Instituto de Medicina Tropical Alexander von Humboldt; Facultad de Medicina Alberto Hurtado, Universidd Peruana Cayetano Heredia, Lima, Peru, 21 Faculty of Tropical Medicine, Mahidol University, Bangkok, Thailand, 22 Department of Microbiology and Immunology, Indiana University School of Medicine, Indianapolis, Indiana, United States of America, 23 Instituto de Medicina "Alexander von Humboldt," Universidad Peruana Cayetano Heredia, Lima, Peru

* Current address: National Library of Medicine, National Institutes of Health, Bethesda, Maryland, United States of America

* dfouts@jcvi.org (DEF); jvinetz@ucsd.edu (JMV)

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

OPEN ACCESS

**Citation:** Fouts DE, Matthias MA, Adhikarla H, Adler B, Amorim-Santos L, Berg DE, et al. (2016) What Makes a Bacterial Species Pathogenic?: Comparative Genomic Analysis of the Genus *Leptospira*. PLoS Negl Trop Dis 10(2): e0004403. doi:10.1371/journal.pntd.0004403

**Editor:** Pamela L. C. Small, University of Tennessee, UNITED STATES

**Received:** June 19, 2015

**Accepted:** January 3, 2016

**Published:** February 18, 2016

**Copyright:** This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.

**Data Availability Statement:** All primary sequence data are available on GenBank and the accession numbers are located in Table 1 of the manuscript.

**Funding:** This project has been funded in whole or part with federal funds from the National Institute of Allergy and Infectious Diseases, National Institutes of Health, Department of Health and Human Services under Contract Number HHSN272200900007C. This work was also supported in part by the following U.S. Public Health Service grants: U19AI115658 (JMV), R01AI108276 (JMV), D43TW007120 (JMV), K24AI068903 (JMV), R21AI115273 (MAM),

## Abstract

Leptospirosis, caused by spirochetes of the genus *Leptospira*, is a globally widespread, neglected and emerging zoonotic disease. While whole genome analysis of individual pathogenic, intermediately pathogenic and saprophytic *Leptospira* species has been reported, comprehensive cross-species genomic comparison of all known species of infectious and non-infectious *Leptospira*, with the goal of identifying genes related to pathogenesis and
R01AI052473 (AIK), U01AI088752 (AIK), R25TW009338 (AIK), R01TW009504 (AIK), and R01AI121207 (AIK). In addition, support to the A. Buschiazzo team was provided in part by grants FSA_1_2013_1_12557 and ALI_1_2014_1_4982 from ANII (Uruguay). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

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

mammalian host adaptation, remains a key gap in the field. Infectious *Leptospira*, comprised of pathogenic and intermediate pathogenic *Leptospira*, evolutionarily diverged from non-infectious, saprophytic *Leptospira*, as demonstrated by the following computational biology analyses: 1) the definitive taxonomy and evolutionary relatedness among all known *Leptospira* species; 2) genomically-predicted metabolic reconstructions that indicate novel adaptation of infectious *Leptospira* to mammals, including sialic acid biosynthesis, pathogen-specific porphyrin metabolism and the first-time demonstration of cobalamin (B12) autotrophy as a bacterial virulence factor; 3) CRISPR/Cas systems demonstrated only to be present in pathogenic *Leptospira*, suggesting a potential mechanism for this clade's refractoriness to gene targeting; 4) finding *Leptospira* pathogen-specific specialized protein secretion systems; 5) novel virulence-related genes/gene families such as the Virulence Modifying (VM) (PF07598 paralogs) proteins and pathogen-specific adhesins; 6) discovery of novel, pathogen-specific protein modification and secretion mechanisms including unique lipoprotein signal peptide motifs, Sec-independent twin arginine protein secretion motifs, and the absence of certain canonical signal recognition particle proteins from all *Leptospira*; and 7) and demonstration of infectious *Leptospira*-specific signal-responsive gene expression, motility and chemotaxis systems. By identifying large scale changes in infectious (pathogenic and intermediate pathogenic) vs. non-infectious *Leptospira*, this work provides new insights into the evolution of a genus of bacterial pathogens. This work will be a comprehensive roadmap for understanding leptospirosis pathogenesis. More generally, it provides new insights into mechanisms by which bacterial pathogens adapt to mammalian hosts.

## Author Summary

Leptospirosis is an emerging and re-emerging globally important zoonotic infectious disease caused by spirochetes of the genus *Leptospira*. This genus is complex, with members that cause lethal human disease, yet mechanisms that underlie pathogenesis remain obscure. *Leptospira* species are divided into those that are infectious for mammals, and those that are non-infectious environmental saprophytes. Based on biological characteristics and molecular phylogeny, infectious *Leptospira* are further divided into pathogenic and intermediate pathogenic members. The pan-genus genomic analysis of 20 *Leptospira* species reported here shows the evolutionary relationship of the different *Leptospira* clades, and various genetic factors related to virulence and pathogenesis. Infectious *Leptospira* show key adaptations to mammals, for example sialic acid biosynthesis, pathogen-specific porphyrin metabolism, and the observation that pathogenic *Leptospira* are vitamin B12 autotrophs, able to synthesize it from a simple amino acid precursor, L-glutamine. A large novel protein family of unknown function—the Virulence Modifying proteins—is found uniquely in pathogenic *Leptospira*. Similarly, the CRISPR/Cas system was only found in pathogenic *Leptospira*. A comparative genomic analysis of a complex bacterial genus allowed us to identify large-scale changes that provides new insights into general processes by which bacteria evolve to become pathogenic.
Introduction

Leptospirosis is a globally widespread zoonotic disease with important health consequences for humans and domesticated animals [1, 2]. Infectious *Leptospira* have significant affinity for specific mammals but vary in how strictly they adapt to specific hosts [1]. Rodent reservoirs (e.g., reservoir hosts (rats, mice) do not exhibit disease but have long-term renal colonization and excrete organisms in the urine, which is key to leptospiral ecology and its life cycle. Infected livestock (e.g., cattle, pigs) and companion animals (e.g., dogs) may suffer fetal loss and acute kidney, liver and lung injury in response to infection. Infected humans variably exhibit clinical manifestations including asymptomatic infection [3] with or without long-term renal carriage [4], undifferentiated fever, renal failure, jaundice, hemorrhage (especially the severe pulmonary hemorrhage syndrome), meningitis, shock and death.

Past taxonomy divided the *Leptospira* genus into a single pathogenic and a single saprophytic species denoted as *L. interrogans* and *L. biflexa*, respectively, which, in turn, were divided into more than 250 serovars based on the cross-agglutinin absorption (CAAT) assay [1, 5]. In the 1990s, DNA hybridization (DDH) identified 17 ‘genomospecies’ [6], which also distinguished DDH from serovar. DDH complemented by molecular methods and experimental studies have since confirmed the existence of at least 22 species [7–12], and grouping of species as infectious (sometimes referred to as group I and group II pathogens, corresponding to “pathogenic” and “intermediately pathogenic”, respectively) and non-infectious (“saprophytic”) [13]. Technical challenges in performing DDH has led to the development of many different molecular approaches to species identification [14, 15]. The International Committee on Systematics of Prokaryotes, Subcommittee on the Taxonomy of *Leptospiraceae* recently agreed that genome sequence comparison should replace DDH for species definition [16]. Such methods include sequence-based phylogeny and calculation of *in silico* genomic similarities between isolates by using draft genomes [17–20].

Leptospiral typing is important for carrying out outbreak investigations and in identifying likely mammalian host reservoir sources of infection. Two commonly used molecular methods performed are pulsed-field gel electrophoresis and multilocus sequencing typing (MLST). MLST has the advantage that it reflects the underlying population genetic structure, is reproducible, is robustly supported by experimental data, and even can be used directly to identify infecting *Leptospira* in clinical samples [21–27]. Genome sequencing, which has become widely available, together with automated tools that assign MLST sequence types directly from sequence data, has demonstrated an important potential for typing, with the expectation that automated analysis tools will become sufficiently user-friendly for rapid and efficient whole genome analysis and comparison, including phylogenetic analysis based on the identification of single nucleotide polymorphisms (SNPs) in the core genome.

The goal of the *Leptospira* Genome Project, initiated in 2011, has been to obtain and compare whole genome information for all known *Leptospira* species. Among the goals of this analysis are the following: i) identifying *Leptospira* pathogenesis mechanisms that might explain heterogeneity in clinical manifestations of leptospirosis; ii) understanding the relationship of genomic content and context to pathogenesis; iii) determining the definitive evolutionary relationship of *Leptospira* towards understanding how infectious *Leptospira* diverged from saprophytes; and iv) identifying common antigens for improving diagnosis and vaccine development. Prior to this project, there were 9 known pathogenic *Leptospira* species, 5 intermediate *Leptospira* species, and 6 saprophytic *Leptospira* species [28], for which whole genome sequence analysis was available for two pathogenic species (two serovars of *L. interrogans*, Lai [29] and Copenhageni [30], two serovars of *L. borgpetersenii* [31]), one intermediate pathogen *L. licerasiae* [32], and one saprophyte species, *L. biflexa* [33]. Since the advent of the present
large-scale project the whole genome sequence of another pathogen, *L. santarosai* serovar Sher-
man has been reported [34] but without comparative analysis.

The present study reports a systematic comparative genome analysis of the 20 *Leptospira* species known when this project began. These species comprise the pathogenic, intermediate pathogenic and saprophytic clades, defined by 16S rRNA gene sequence [2, 35, 36] and complemented by DNA-DNA hybridization [7, 11, 32]. This analysis focuses on the main genomic features and content distinguishing infectious from non-infectious *Leptospira*, and on how specific genes and gene families have ramified in the pathogenic and intermediate pathogenic *Leptospira*.

**Methods**

*Leptospira* isolates

A globally representative collection of the 20 *Leptospira* species known at the advent of this
project was analyzed here, and provided by members of the leptospirosis research community
*(Table 1).*

DNA preparation from isolates

A standard operating procedure was established for all contributing laboratories to follow in preparing DNA for whole genome sequencing. *Leptospira* were considered to be like Gram-negative bacteria for the purpose of DNA extraction because of the presence of lipopolysaccharide and a thin peptidoglycan cell wall. Either ~10¹² bacterial cells or 30 mL of the densest possible culture of *Leptospira* (in EMJH medium) were centrifuged and the pellet resuspended in 180 µl Buffer ATL (all buffer abbreviations are according to the manufacturer and related specifically to components of the kit) (Qiagen Tissue Kit, Valencia, CA, USA), and then purified according to the protocol “Purification of Total DNA from Animal Tissues (Spin-Column Protocol),” including the use of proteinase K, thorough vortexing throughout the procedure, and incubated at 56°C until the cells were completely lysed, according to the manufacturer’s instructions. Lysis was usually complete in 1–3 hr. Before adding Buffer AL, 10 µl of RNAs cocktail was added (mixture of two highly purified ribonucleases, RNase A (500 U/ml) and RNase T1 (20,000 U/ml); Ambion, Life Technologies, Carlsbad, CA), and incubated for 30–60 min at 37°C. After vortexing for 15 s, Buffer AL was added to the sample, which was again mixed thoroughly by vortexing, followed by adding ethanol (96–100%), and mixed again by vortexing. It was considered essential that the sample, Buffer AL, and ethanol were mixed immediately and thoroughly by vortexing or pipetting to yield a homogeneous solution. After this point samples were handled with large bore, genomic DNA-compatible tips. The mixture from step 3 (including any precipitate) was pipetted into the DNeasy Mini spin column placed in a 2 ml collection tube (provided), centrifuged at 6000 x g for 1 min and the flow-through discarded. The DNeasy Mini spin column was placed in a new 2 ml collection tube, and washed with 500 µl Buffer AW1. The DNeasy Mini spin column was placed in a new 2 ml collection tube, washed with 500 µl Buffer AW2, and centrifuge for 3 min at 20,000 x g (14,000 rpm) to dry the DNeasy membrane; remaining ethanol was considered to interfere with sequencing reactions and processing. The DNeasy Mini spin column was eluted twice by adding 2 x 100 µl Buffer AE directly onto the DNeasy membrane, incubating at room temperature for 1 min, and then centrifuged for 1 min at 6000 x g (8000 rpm). DNA was shipped on dry ice to JCVI for sequencing. Quality control included certification of intact, high molecular weight DNA and required 15–30 µg for fragment libraries, complementarily documented by agarose gel image containing a DNA Mass Ladder, OD260/280 determination, and an estimated DNA concentration from a fluorometric assay (SYBR Green, Quant-IT PicoGreen dsDNA Assay Kits).
<table><caption>Table 1. Leptospiral species metadata and genome statistics.</caption><thead><tr><th>Species</th><th>Serovar</th><th>Strain</th><th>Lifestyle</th><th>Source</th><th>Country of Origin</th><th>Contributing Lab</th><th>Bioproject</th><th>WGS Master or GenBank Accession</th><th>BEI/ATCC Size (Mbp)</th><th># Contigs</th><th>#genes</th><th>#proteins</th><th>Contig N50</th><th>#Scaffolds</th><th>Average Coverage</th><th>G+C %</th><th>Locus Tag</th><th>Finishing status<sup>a</sup></th></tr></thead><tbody><tr><td>alexcanderi</td><td>Marmhaa 3</td><td>L 60<sup>t</sup></td><td>pathogenic</td><td>Homo sapiens</td><td>China</td><td>Hartskeerl</td><td>PRJNA74069</td><td>AHMTO00000000.2</td><td>NR-22256</td><td>4.22</td><td>65</td><td>4582</td><td>4541</td><td>156.869</td><td>14</td><td>58.1x</td><td>40.2%</td><td>LEP1GSC062 IHQD</td></tr><tr><td>alstoni</td><td>Pingchang</td><td>80-412</td><td>pathogenic</td><td>Frog</td><td>China</td><td>Galloway</td><td>PRJNA77154</td><td>AHD00000000.2</td><td>NR-35981</td><td>4.44</td><td>67</td><td>4392</td><td>4380</td><td>152.980</td><td>33</td><td>57.6x</td><td>42.5%</td><td>LEP1GSC193 IHQD</td></tr><tr><td>borgpetersenii</td><td>Javanica</td><td>UI 09931</td><td>pathogenic</td><td>Homo sapiens</td><td>Thailand</td><td>Peacock</td><td>PRJNA74113</td><td>AHNP00000000.2</td><td>NR-20171</td><td>3.89</td><td>15</td><td>4053</td><td>4010</td><td>468.398</td><td>8</td><td>60.4x</td><td>40.1%</td><td>LEP1GSC103 IHQD</td></tr><tr><td>internogens*</td><td>Copenhageni</td><td>Fiocruz L1-130</td><td>pathogenic</td><td>Homo sapiens</td><td>Brazil</td><td></td><td>ATCC</td><td></td><td>AE016823/4</td><td>4.63</td><td>2</td><td>3762</td><td>3667</td><td></td><td>2</td><td></td><td>35.0%</td><td>LIC</td><td>Complete</td></tr><tr><td>kirschneri</td><td>Cyptopteri</td><td>3522 C<sup>t</sup></td><td>pathogenic</td><td>Chiroptera (bat)</td><td>Indonesia</td><td>Hartskeerl</td><td>PRJNA74057</td><td>AHMNO00000000.2</td><td>NR-22255</td><td>4.41</td><td>24</td><td>4029</td><td>3986</td><td>415.138</td><td>13</td><td>58.0x</td><td>35.9%</td><td>LEP1GSC049 IHQD</td></tr><tr><td>kimetyi</td><td>Malaysia</td><td>Bejo-Isog<sup>t</sup></td><td>pathogenic</td><td>environmental</td><td>Malaysia</td><td>Hartskeerl</td><td>PRJNA74061</td><td>AHMPO00000000.2</td><td>NR-22254</td><td>4.42</td><td>4</td><td>4282</td><td>4238</td><td>3,808.637</td><td>4</td><td>57.2x</td><td>44.8%</td><td>LEP1GSC052 IHQD</td></tr><tr><td>noguchii</td><td>Panama</td><td>CZ 21<sup>t</sup></td><td>pathogenic</td><td>Dolichis marsupialis (opossum)</td><td>Panama</td><td>Hartskeerl</td><td>PRJNA67231</td><td>AKWV00000000.2</td><td>NR-22283</td><td>4.71</td><td>35</td><td>4565</td><td>4520</td><td>280.299</td><td>12</td><td>58.2x</td><td>35.5%</td><td>LEP1GSC059 IHQD</td></tr><tr><td>santarosaii</td><td>Shermani</td><td>1342K<sup>t</sup></td><td>pathogenic</td><td>Homo sapiens</td><td>Panama</td><td>Hartskeerl</td><td>PRJNA178172</td><td>AHB00000000.2</td><td>Not Deposited</td><td>3.99</td><td>68</td><td>4049</td><td>4002</td><td>108.025</td><td>40</td><td>62.8x</td><td>41.8%</td><td>LEP1GSC048 IHQD</td></tr><tr><td>weili</td><td>undetermined</td><td>LNT 1234</td><td>pathogenic</td><td>Homo sapiens</td><td>Thailand</td><td>Peacock</td><td>PRJNA74087</td><td>AHNC00000000.2</td><td>NR-20183</td><td>4.26</td><td>83</td><td>4478</td><td>4436</td><td>172.307</td><td>24</td><td>59.7x</td><td>40.8%</td><td>LEP1GSC086 IHQD</td></tr><tr><td>broomii</td><td>Hurstbridge</td><td>539<sup>t</sup></td><td>intermediate</td><td>Homo sapiens</td><td>Denmark</td><td>Hartskeerl</td><td>PRJNA74059</td><td>AHMNO00000000.2</td><td>NR-22253</td><td>4.40</td><td>12</td><td>4249</td><td>4205</td><td>2,457.338</td><td>6</td><td>54.0x</td><td>43.0%</td><td>LEP1GSC050 IHQD</td></tr><tr><td>fairei</td><td>Hurstbridge</td><td>BUT 6<sup>t</sup></td><td>intermediate</td><td>Sus scrofa (wild boar)</td><td>Australia</td><td>Hartskeerl</td><td>PRJNA167230</td><td>AKWZ00000000.2</td><td>NR-22252</td><td>4.29</td><td>13</td><td>4157</td><td>4113</td><td>579.278</td><td>5</td><td>59.1x</td><td>43.5%</td><td>LEP1GSC058 IHQD</td></tr><tr><td>inadai</td><td>Lyme</td><td>10<sup>t</sup></td><td>intermediate</td><td>Homo sapiens</td><td>USA</td><td>Hartskeerl</td><td>PRJNA74055</td><td>AHMMP00000000.2</td><td>NR-22251</td><td>4.46</td><td>27</td><td>4314</td><td>4264</td><td>810.442</td><td>14</td><td>55.9x</td><td>44.6%</td><td>LEP1GSC047 IHQD</td></tr><tr><td>licerasiae</td><td>Varilai</td><td>VAR 010<sup>t</sup></td><td>intermediate</td><td>Homo sapiens</td><td>Peru</td><td>Vinetz</td><td>PRJNA74167</td><td>AHOO00000000.2</td><td>ATCC</td><td>4.21</td><td>14</td><td>3974</td><td>3931</td><td>553.148</td><td>4</td><td>58.3x</td><td>41.1%</td><td>LEP1GSC185 IHQD</td></tr><tr><td>wolffii</td><td>Khorat</td><td>Khorat-H2<sup>t</sup></td><td>intermediate</td><td>Homo sapiens</td><td>Thailand</td><td>Hartskeerl</td><td>PRJNA167232</td><td>AKWX00000000.2</td><td>NR-22250</td><td>4.40</td><td>23</td><td>4252</td><td>4206</td><td>916.783</td><td>10</td><td>58.5x</td><td>45.6%</td><td>LEP1GSC061 IHQD</td></tr><tr><td>biflexa<sup>t</sup></td><td>Patoc</td><td>Patoc I (Paris)</td><td>saprophytic</td><td>stream water</td><td>Italy</td><td></td><td>PRJNA56893</td><td>CP0007867</td><td>ATCC</td><td>3.95</td><td>2</td><td>3774</td><td>3725</td><td></td><td>3</td><td>6x</td><td>38.9%</td><td>LEPBI Complete</td></tr><tr><td>meyeri</td><td>Hardjo</td><td>Went 5</td><td>saprophytic</td><td>unknown</td><td>Canada</td><td>Galloway</td><td>PRJNA167225</td><td>AKXE00000000</td><td>NR-29052</td><td>4.19</td><td>7</td><td>4011</td><td>3969</td><td>2,478.877</td><td>3</td><td>61.1x</td><td>38.0%</td><td>LEP1GSC017 IHQD</td></tr><tr><td>terpstraee</td><td>Hualin</td><td>LT 11-33<sup>t</sup></td><td>saprophytic</td><td>unknown</td><td>China</td><td>Hartskeerl</td><td>PRJNA178717</td><td>AOGW00000000.2</td><td>ATCC</td><td>4.09</td><td>23</td><td>3932</td><td>3889</td><td>372.285</td><td>10</td><td>67.1x</td><td>38.2%</td><td>LEP1GSC203 IHQD</td></tr><tr><td>vanthieli</td><td>Holland</td><td>WaZ</td><td>saprophytic</td><td>water</td><td>Netherlands</td><td>Galloway</td><td>PRJNA177160</td><td>AOGY00000000.2</td><td>ATCC</td><td>4.23</td><td>89</td><td>4251</td><td>4205</td><td>164.554</td><td>28</td><td>58.1x</td><td>38.9%</td><td>LEP1GSC199 IHQD</td></tr><tr><td>wolbachii</td><td>Codice</td><td>CCD</td><td>saprophytic</td><td>unknown</td><td>USA</td><td>Galloway</td><td>PRJNA177156</td><td>AOGZ00000000.2</td><td>700522</td><td>4.08</td><td>25</td><td>3956</td><td>3911</td><td>2,540.234</td><td>9</td><td>70.5x</td><td>39.2%</td><td>LEP1GSC195 IHQD</td></tr><tr><td>yanagavae</td><td>Sao Paulo</td><td>Sao Paulo</td><td>saprophytic</td><td>water</td><td>Brazil</td><td>Hartskeerl</td><td>PRJNA178716</td><td>AOGX00000000.2</td><td>ATCC</td><td>4.06</td><td>47</td><td>4011</td><td>3967</td><td>330.086</td><td>19</td><td>60.4x</td><td>38.2%</td><td>LEP1GSC202 IHQD</td></tr></tbody></table>

§ Improved High-Quality Draft (IHQD); High-Quality Draft (HQD)

* J. Bacteriol. 186 (7), 2164–2172 (2004)

†PLoS ONE 3 (2), E1607 (2008)

Ttype strain

doi:10.1371/journal.pntd.0004403.t001
## Genome sequencing, draft assembly and annotation

The genomes of 17 *Leptospira* species (the whole genome sequences of the remaining 3 species studied in the present analysis, *L. interrogans* serovar Copenhageni strain L1-130 and *L. biflexa* serovar Patoc Strain Patoc I, and *L. licerasiae* serovar Varillal strain VAR010, already were published [33, 37, 38]) were sequenced at JCVI by a combination of Illumina HiSeq (2x100 bp) and 454 FLX Titanium. Briefly, paired-end libraries were constructed for each sequencing technology from random nebulized genomic DNA in the 300–800 bp (Illumina) and 2–3 kb (454) size ranges. Sequence reads were generated with a target average read depth of ~ 20–30 fold (454) and ~60-fold (Illumina) coverage. Sequences for all 18 strains were assembled using the Celera Assembler version 6.1 [38], and ordered using NUCmer [39] to align the contigs to the best-matching closed *Leptospira* reference genome. All 18 new genome sequences underwent manual gap closure to elevate the genome status to improved high-quality draft (Table 1). Contigs were annotated for protein- and RNA-encoding features using the JCVI automated annotation pipeline essentially as described [40–43] except HMMs were run using HMMER3 [44].

## Phylogenetic analysis

16S rRNA trees were generated by first creating a multiple sequence alignment to the bacterial 16S rRNA reference alignment using Ribosomal Database Project release 10 (RDP-X) [45]. The aligned FASTA sequences were downloaded and trimmed to remove gapped columns using Belvu (v2.31) [46]. Based on the alignment, a bootstrapped Maximum-likelihood tree was subsequently inferred using phylopFasta, an in-house wrapper script [47] for the Phylip program [48, 49].

SecY trees were created by first aligning *secY* nucleotide sequences using Clustal Omega [50] with 100 combined guide-tree/HMM iterations. The multiple sequence alignment was trimmed to remove gapped columns and a bootstrapped Maximum-likelihood tree was inferred as was done for the 16S rRNA trees.

The nucleotide sequences of 7 MLST housekeeping genes were extracted from the 20 genomes, and sequence types (STs) assigned using the MLST website (http://leptospira.mlst.net/) [51]. A multiple sequence alignment of the concatenated sequences of 7 MLST loci was performed using ClustalW implemented in MEGA version 5 [52]. A maximum likelihood tree was re-constructed using an algorithm implemented in PhyML version 3.0.1 [53]. The model of sequence evolution used was the generalized time-reversible (GTR) model with gamma-distributed rate variation. The CLC Main Work Bench version 7.0 was used to edit and display the tree (Qiagen, USA).

Universal protein marker trees were constructed using a set of 39 proteins that are universally conserved among bacteria and produce monophyletic phylogenies, suggesting that they undergo minimal horizontal transfer (S1 Table) [54–56]. Protein sequences were aligned using ClustalW (v1.83) [57] using default settings. The alignment was trimmed to remove gapped columns using *trimAl* (v1.2r59) using–nogaps option and– fasta output option [58]. Aligned and trimmed predicted amino acid sequences of each species were concatenated as described previously [40], in the following order: AspS, FusA, GyrB, InfB, LepA, LeuS, PyrG, RplA, RplB, RplC, RplE, RplF, RplK, RplM, RplN, RplO, RplP, RplR, RplV RpoA, RpoB, RpsB, RpsC, RpsD, RpsE, RpsG, RpsH, RpsI, RpsK, RpsL, RpsM, RpsO, RpsQ, SecY, SerS, TopA, TsaD, Tuf, and YchF. The resulting alignment of 11241 amino acids was used to generate a Maximum-likelihood tree from 100 bootstrapped using raxmlFasta, an in-house wrapper script for the raxmlHPC (v7.0.4) [59].

A pan-genome tree was constructed using the mean of the BLASTP Score Ratio (BSR) as described previously [60]. The PanOCT output file
*100_pairwise_BSR_distance_matrix_phylip.txt* was used as input for Neighbor [48, 49] to build an unrooted UPGMA Neighbor-Joining tree. This PanOCT output file is a Phylip-style distance matrix derived from the pairwise mean BSR of core proteins present in 100% of genomes where a single value is presented for each pair of genomes in the pan-genome.

## *In silico* DNA-DNA hybridization

Genome relatedness among *Leptospira* strains was determined pairwise from fully or partially sequenced genomes using the genome-to-genome distances (GGD) calculator (S2 Table) [61]. This analysis was complemented by *in silico* DNA-DNA hybridization as previously reported [11].

## Pan-genome analysis

Clusters of orthologous proteins were generated using version ver3_13 of PanOCT [62]. Since PanOCT does not place paralogs into its ortholog clusters, but does produce a paralogs.txt file that specifies which clusters are paralogs, an in-house PERL script, paralogs_matchtable.pl, was created to merge paralogous clusters. This approach was necessary because analysis of core and novel genes has historically been defined for clusters containing all paralogs [63–68]. The R script, compute_pangenome.R, from Park et al. [67] and paralog_matchtable.pl were used to construct the pan-genome, core and novel genes plot. We initially chose not to compute permutations in genome order for the reasons described in [69]. As a consequence of a lack of permutations, compute_pangenome.R was modified to load in a defined genome order of addition.

Orthologous protein content was compared and illustrated in a Venn diagram that was constructed using output from an in-house PERL script, create_meta_groupings.pl that uses output from PanOCT and a file that describes how genomes are to be grouped. Genomes were grouped by whether they are infectious (group I or group II), non-infectious (saprophytic), or an outgroup (*Leptonema illini*). Since there were multiple genomes per group (except for *Lep-tonema*), clusters were counted if there was a majority (50%), all-but-one, or all protein members from a particular group or groups. Clusters not matching these criteria were not counted.

## Metabolic reconstructions

Four draft metabolic network reconstructions were created for representative species chosen from pathogenic, intermediate and saprophytic clades including the following: *L. interrogans* (serovar Copenhageni strain L1-130), *L. licerasiae* (serovar Varillal strain VAR010), *L. biflexa* (serovar Patoc strain Patoc 1 (Ames)) and *L. kmetyi* (serovar Malaysia strain Bejo-Iso9). *L. kmetyi* was included in these comparisons because in addition to having recently been reported to infect humans in the Caribbean islands [70, 71], initial genome examination suggested that *L. kmetyi* could belong to a transitional group between the group I and group II pathogens and distinct from the other group I pathogens, here represented by *L. interrogans*. The reconstructions were built using the ModelSeed framework [72].

The COBRApy toolbox [73] was used to perform Flux Balance Analysis (FBA) [74] simulations and constraint-based analyses using the gurobi linear programming solver [75]. The constraint-based model consists of an *S* matrix composed of distinct metabolites and reactions including exchange and biomass reactions (S3 Table). Each of the reactions has an upper and lower bound on the flux it can carry. Reversible reactions have an upper bound of 1,000 mmol gD*W*⁻¹h⁻¹ and a lower bound of -1,000 mmol gD*W*⁻¹h⁻¹, making them practically unconstrained, while irreversible reactions have a lower bound of zero. By default, the biomass reaction was set as the objective to be maximized. The exchange reactions that allow for extracellular metabolites to pass in and out of the system were defined such that a positive flux indicates flow out. The GapFind MILP algorithm [76] encoded in the COBRApy Toolbox was
performed for the models unable to grow in minimal medium (biomass objective function equal to zero), to find exchange reactions allowing for *in vitro* growth, indicating strain-specific auxotrophies.

An *in silico* minimal medium was constructed that supported growth for all of the *Leptospira* models, consisting of trace elements (magnesium, manganese zinc, sulfate, calcium, copper, phosphate, cobalt, chlorine, potassium, ferrous and ferric iron and ammonia), water, oxygen, heme, CO₂, the vitamins thiamine (B1), folate and menaquinone, glycerol and the fatty acids lauric acid, stearic acid, and decanoate and aminoethanol, meso-2,6, diaminopimelate, as well a variety of amino acids (glutamate, aspartate, tyrosine, phenylalanine, asparagine (S3 Table).

## Genomic location and genetic organization of vitamin B12 biosynthesis-related gene clusters in *Leptospira*

PanOCT data were used to identify the cob I/III and cob II gene clusters in infectious and non-infectious *Leptospira* (when present), which encode proteins predicted to participate in B12 transport or synthesis. To determine the genomic locations of the *btuB* and cob II and I/III clusters, a custom *Postgresql* database was created using the annotated genomes of 20 species. Orthologs were identified with blastp (-F "m S" -s T) and conserved genomic neighborhood using the Prokaryotic Sequence homology Analysis Tool (PSAT) [77].

## Results

### General genomic descriptions

The genomes of 17 *Leptospira* spp. isolates were newly determined for this study, representing 8 pathogenic, 4 intermediate and 5 saprophytic clades and used in pan-genomic comparative analyses along with the previously reported genomes of *L. interrogans* serovar Copenhageni strain Fiocruz L1-130 [30], *L. biflexa* serovar Patoc strain Patoc I (Paris) [33] and *L. licerasiae* serovar Varillal strain VAR010T [32]. Thirteen of these isolates were sequenced to a genome finishing status [78, 79] of "Improved High-Quality Draft" (IHQD) and 5 to a status a "High-Quality Draft" (HQD) (Table 1). To achieve a genome finishing status of IHQD, manual finishing was conducted consisting of contig sequence extension, sequence gap closure, and PCR to link physical ends. On average, the genomes assembled into 36 contigs [range 4 (*L. kmetyi*) to 89 (*L. vanthielii*)], with an average genome size of 4.26 Mbp in length [range 3.89 Mbp (*L. borgpetersenii* to 4.71 Mbp (*L. noguchii*)) at an average of 59.7-fold sequence coverage. The average G+C% was 40.7% [range 35.5% (*L. noguchii*) to 45.6% (*L. wolfii*)]. These genomes were predicted to encode an average of 4,197 protein-coding sequences per genome [range 3,932 (*L. terpstrae*) to 4,582 (*L. alexanderi*)].

### Phylogenetic analysis of the *Leptospira* Genus to determine evolutionary relatedness

Twenty genome sequences (17 new, 3 previously published) of isolates representing 20 of the 22 known *Leptospira* spp. (Table 1) were used to determine phylogenetic distances between species; two recently reported species (*L. idonii* [80] and *L. mayottensis* [81]) are not included here. Phylogenetic relationships among all *Leptospira* species were analyzed in five independent ways (Fig 1): A, a core set of 39 concatenated genes coding for housekeeping proteins (universal markers); B, a pan-genus set of 1350 proteins; C, 3) multilocus sequence typing (MLST)[51]; D, 16S rRNA (highly conserved); and E, *secY* (highly variable)[70]. *Leptonema illini* was used as the outgroup for all analyses. Each approach yielded different nodes and
<figure><img src="image_4.png" /></figure>

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

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

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

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

**Fig 1. Phylogenetic analysis of Leptospira species.** Consensus maximum-likelihood trees are depicted using multiple alignments of 16S rRNA (A), secY (B), MLST (C) and 39 concatenated protein data sets (D). The numbers along the branches denote percent occurrence of nodes among 100 bootstrap replicates. A pan-genome tree was generated based on the mean of the BLASTP score ratio of core 1135 proteins (E). The scale bar represents the number of nucleotide (A-C), amino acid (D & E) substitutions.

doi:10.1371/journal.pntd.0004403.g001

branches of the species, except for 16S rRNA sequences, for which deduced phylogeny did not discriminate between *L. meyeri* and *L. yanagawae*. Additionally, these five approaches revealed three clades correctly clustering members of the nine pathogenic (group I) and five intermediate pathogenic (group II) and six non-pathogenic (saprophytic) species. Only the trees based
on *secY*, the universal markers and the pan-genome clearly separated the closely related pathogenic species *L. interrogans*, *L. kirschneri* and *L. noguchii* from the other 6 pathogenic species. As expected, phylogenetic positions shifted to some extent between the single locus-based analyses and became more consistent using the multi-locus approaches. The following pairs of species showed close relationships: the pathogenic species (except in the *secY* tree) *L. interrogans* and *L. kirschneri*, *L. alexanderi* and *L. weilii*; the group II species, *L. inadai* and *L. broomii*; and the non-pathogenic (except 16S and *secY* trees) *L. wolbachii* and *L. vanthielii*.

Furthermore, the genome relatedness between pairs of representative strains of each of the 20 species from fully or partially sequenced genomes confirmed the genetic relatedness among *Leptospira* species as established by DDH (Fig 1/ S2 Table).

## Pan-genome analysis

The pan-genome is defined as a core set of genes shared by all isolates plus a variable set of genes shared by a subset of isolates, and strain-specific or novel genes. Based on these 20 representative genomes and raw PanOCT output, the size of the core- and pan-genome was determined to contain 1,764 and 17,477 genes, respectively; however, with paralogs collapsed, the size of the core- and pan-genome was determined to be 1,592 and 13,822 genes, respectively (S1A Fig). The number of species-specific or novel genes ranged from 233 to 892 (S1A Fig) for each new genome added. After the addition of the third genome (*L. noguchii*), the size of the core gene set plateaued, while the pan-genome continued to rise (S1A Fig).

To determine whether the *Leptospira* pan-genome is open or closed (as defined below, the number of new genes identified (i.e., unique or strain-specific genes) for each genome added was determined and fit to a power law function (n = κN⁻⁰) as described previously [63]. Conceptually, a pan-genome is closed when sequencing the genomes of additional isolates fails to increase gene number (i.e., the entire gene repertoire has been discovered) [82]. The exponent (α) indicates whether the pan-genome is open (α ≤ 1) or closed (α > 1) [83]. Using this equation, the pan-genome of *Leptospira* was inferred to be open (α = 0.49 ± 0.02) (S1B Fig). From an exponential decay function, the number of new genes predicted for each genome (species) added was extrapolated and calculated to be 409 ± 12 on average (S1B Fig).

The distribution of protein clusters representing gene families among the three groups (pathogenic, intermediate, saprophytic) is depicted in a Venn diagram (Fig 2A). Because there were multiple genomes per group (except for the *Leptonema* outgroup), clusters were counted if there was a majority (50%), all-but-one, or all protein cluster members from a particular group or groups. Focusing on the majority criteria, pathogens and intermediates had nearly equal numbers of group-specific genes (416 and 424, respectively), and the highest number of shared genes between two groups (369). Binary comparisons of pathogens and intermediates with saprophytes revealed just 52, and 78 genes shared, respectively. When comparing only *Leptospira*-specific genes, the core genome was comprised of 737 genes, with the majority of genes being shared with *Leptonema illini*. Closer examination of species-specific genes showed that pathogenic *Leptospira* have more species-specific genes on average (637±129) than do intermediates (418±126) or saprophytes (321±90). *L. noguchii* sv. Panama str. CZ 214T had the greatest number of species-specific genes among species compared in this study. To understand the function of genes shared among infectious *Leptospira*, the distribution of protein functions was examined for clusters shared among infectious and non-infectious *Leptospira* (Fig 2B). The only functional category dominated by pathogenic *Leptospira* was “mobile and extrachromosomal elements.” The functional categories that stood out most among genes shared between pathogens and intermediates was “biosynthesis of cofactors, prosthetic groups, and carriers” and “fatty acid and phospholipid metabolism.” Saprophyte-specific genes
A.

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

B.

<figure><img src="image_10.png" /><figcaption>**Fig 2.** Pan-genomic comparisons of 20 *Leptospira* species. **Panel A:** Orthologous protein clusters were binned, counted and placed into a Venn diagram by whether clusters contained proteins from genomes in each of the three *Leptospira* groups: pathogenic (A), intermediate (B), saprophytic (C) and the *Leptonema* outgroup (D). Clusters were counted if there was a majority (50%), all-but-one, or all protein members from a particular group or groups (separated by colons). Singleton clusters, representing species-specific or strain-specific genes are noted in circles surrounding the Venn diagram. Clusters not matching any of these criteria or containing at least one protein from another group were considered as ambiguous groupings. The Venn diagram is not to scale. **Panel B:** Protein clusters unique to pathogenic, intermediate, and saprophytic groups or shared only between pathogenic and intermediate groups were counted by main functional role categories. See key for group colors.</figcaption></figure>

doi:10.1371/journal.pntd.0004403.g002
dominated 10 of the 16 functional role categories, many of which were involved in central
intermediary and energy metabolism, gene regulation, signal transduction, protein fate, cell
envelope, and transport functions.

## Protein secretion systems

All *Leptospira* clades are predicted to have type II protein secretion systems but do not appear to contain a type III secretion system. Lipoproteins in *Leptospira* have particular pathogenetic significance because of their potential as vaccine targets and virulence factors involved in host-pathogen interactions [84–93].

### An unusual sec system in *Leptospira*

*Leptospira* contains the *Sec* system for signal-peptide-containing proteins and signal peptidase to remove the signal peptide at the time of secretion. Genes encoding the signal recognition particle (SRP) protein Ffh and receptor protein FtsY were not found in any *Leptospira* genome, nor was the SRP structural RNA. Generally, the lack of SRP and its receptor is unusual in bacteria, although the system is missing in the genus *Dehalococcoides* and also, apparently, in the uncultivated marine lineage SAR86 [94]. The narrow, elongated spiral shape limits the distance a ribosome can be from the *Leptospira* plasma membrane and may obviate the need for translation arrest by SRP. However, in looking for novel features near *Sec* system genes in *Leptospira* showed a novel gene inserted between the normally consecutive genes for *Sec* system proteins *SecY* and *YajC*, encoding a non-globular protein with an N-terminal signal peptide and a transmembrane segment towards the C-terminus, with the majority of residues in between consisting of low complexity, poorly conserved sequence especially rich in *Lys*, *Glu*, and *Asn*. No homologs to this low-complexity protein occur outside the *Leptospira* genus. We postulate that this novel gene could be involved in protein secretion.

### Unusual sec-independent (twin-arginine) translocation system

Twin-arginine translocation (TAT) in prokaryotes allows completion of complete protein folding prior to Sec-independent secretion through the plasma membrane. Except in the halo-philes, where high salt outside the cell explains the need for folding prior to export, TAT substrates tend to be redox cofactor-binding proteins [95]. These proteins fold and bind their cofactors before crossing the membrane. The *tatA* and *tatC* gene-encoded components of the translocase are evident in *Leptospira* genomes, but the twin-arginine signal itself proved elusive. The TIGRFAMs collection [96] hidden Markov model (HMM) TIGR01409 finds no sequence scoring near the trusted cutoff in any species of *Leptospira*, nor in *Leptonema illini*. However, alignments of full-length homologs in *Leptospira* to recognizable TAT translocation substrates from other lineages could be extended into the N-terminal signal region. Such alignments often showed a *Lys-Arg* dipeptide in the *Leptospira* sequence aligned to the *Arg-Arg* motif of recognizable TAT signal sequences. This observation triggered a review of all candidate families of TAT translocation substrates in *Leptospira*, and produced iterative refinement of the lineage-specific TAT signal, and a catalog of TAT substrates.

Eleven protein families were confirmed as TAT substrates by multiple criteria, including strong conservation of the putative TAT signal within the protein family, alignment to non-spirochete homologs that extended N-terminally into the TAT signal region, and strong sequence similarity of the putative TAT signal motif, usually RKxFL, across the different *Lep*-tospira putative TAT translocation substrate families. A continuous 18-residue stretch from each protein in each of these was used to construct a seed alignment of *Leptospira* TAT signal sequences, including the modified Twin-Arg motif and the adjacent hydrophobic region.
<figure><img src="image_11.png" /><figcaption>**Fig 3.** TAT signal sequence in *Leptospira sp.* The X-axis shows position in an ungapped alignment. The Y-axis shows information content, measured in bits.</figcaption></figure>

doi:10.1371/journal.pntd.0004403.g003

Comparative analysis predicted a conserved TAT signal sequence in all *Leptospira* [97] (Fig 3). Eleven protein families comprising the defining alignment and for two additional families are strong candidates for TAT-dependent, Sec-independent translocation (S4 Table). Only one of the 13 families, the PhoX alkaline phosphatase family, was observed to be largely restricted to pathogenic species of *Leptospira* [98].

Inspecting the family of LIC_10874 (a 4Fe-4S dicluster domain protein family of LIC_10874) within and outside the genus *Leptospira* demonstrated conservation of the putative TAT signal in both, and the substitution in *Leptospira* of the second Arg by Lys, as in other families. This family is notable, however, because in multiple species from phylogenetically distant clades, translation start sites can be assigned with high confidence, and the TAT signal begins rather far (some 50 residues) from that start. Member sequences in this family all share a well-conserved prefix domain, ~50 amino acids in length, between the start of translation, and star of the recognizable TAT signal.

An unexpected feature in the *Leptospira* TAT system cassette is a probable serine phosphatase encoded next to *tatC* (family TIGR04400), which either overlaps it or is present within five base pairs in all 20 leptospiral species examined. It is not known whether this putative phosphatase is involved in Sec-independent translocation *per se*, rather than in its regulation, or in some unrelated process.

## Lipoprotein secretion

Cleavage of pro-lipoproteins by the type II signal peptidase occurs within a short lipobox sequence, which includes the invariant cysteine that is targeted for covalent modification with lipids. In *E. coli* lipoproteins, the -1 position immediately preceding the peptidase cleavage site is highly conserved, being occupied by the small nonpolar amino acids Ala or Gly in the vast majority of cases [99]. A previous analysis noted somewhat larger residues, such as Asn, Ser, and Cys, at the -1 position of experimentally-verified spirochete lipoproteins [100].
Examination of *L. interrogans* lipoprotein orthologs in saprophytic and intermediate species revealed a number of unexpected amino acids at the -1 position (S5 Table). For example, the bulky amino acid Tyr was found in the -1 position in LipL21 in all intermediate *Leptospira* spp., whereas those from pathogenic and saprophytic species possess typical -1 residues Ala and Ser, respectively. Conversely, some lipoproteins of saprophytic or intermediate species with expected amino acids at the -1 position were found to have orthologs in *L. interrogans* with variant amino acids at the -1 position (S5 Table). For example, the outer membrane lipoprotein Loa22 of saprophytic species has the allowed residue Asn at -1 while its orthologs in pathogenic and intermediate species has Leu or Phe at the -1 position. Similarly, LIC11088 orthologs in most intermediate and saprophytic species possess permitted residues at -1, whereas the pathogens have Gln or charged residues. Thus, the availability of genome sequences from across the genus *Leptospira* has confirmed the much higher flexibility in the leptospiral lipobox and is anticipated to lead to redefinition of the pan-leptospiral lipobox to accommodate increased amino acid flexibility at the -1 position.

The substrate specificities of the first two enzymes in lipoprotein biogenesis, prolipoprotein diacylglyceryl transferase (Lgt) and signal peptidase II (Lsp), are likely to be influenced by amino acids at the -1 position relative to the lipoprotein cleavage site [101]. For this reason, these enzymes would be expected to possess novel structural properties that allow recognition of an expanded set of residues at the -1 position of the lipobox. Consistent with this notion, Lgt orthologs of all 20 *Leptospira* strains lack the signature sequence that defines most Lgt proteins (Prosite accession PS01311) [102]. Interestingly, *Leptonema illini*, the bacterium most closely related to *Leptospira* among sequenced organisms, harbors two Lgt paralogs: one quite close in sequence to leptospiral Lgt orthologs and another with signature sequence similar to Lgt of all other organisms. This arrangement suggests duplication of *lgt* in an ancestor common to *Lep*-tospira and *Leptonema*, with subsequent loss of the latter and functional divergence of the for-mer to accommodate bulkier -1 lipobox residues. Similarly, Lsp of *Leptospira* species possesses an extra 22 or 24 residues at a position corresponding to a location within the second periplas-mic loop of *E. coli* Lsp [103], which is missing from the Lsp sequence of other bacteria with well-characterized lipoproteins, including the spirochetes *B. burgdorferi* and *T. pallidum*. The sequence features of leptospiral Lgt and Lsp suggest the presence of novel structural features at the active sites of these enzymes consistent with variability at the -1 position of the leptospiral lipobox.

## Metabolic reconstructions

*In silico*, genomically-based metabolic network reconstructions were created for four represen-tative *Leptospira* species occupying different clades: *L. interrogans* and *L. kmetyi* (pathogen; group I), *L. licerasiae* (intermediate pathogen; group II) and *L. biflexa* (non-pathogenic). *L. kmetyi* was chosen for analysis because preliminary genomic inspection suggested unusual fea-tures of this species with regard to pathogenesis-related genes (*vide infra*). The base *in silico* media and default computational bounds (S3 Table) represent every compound allowed to enter the system for cellular uptake to allow all models to produce biomass. Removal of some of these compounds leads to species-specific growth (a unique model-predicted auxotrophy). The use of a steady-state assumption does not allow the flux balance analysis to take into account specific concentrations of a given metabolite but, rather, the predicted rate of uptake, secretion or transformation. The default uptake bounds for each metabolite are provided (S3 Table). Negative bounds represent entry of the metabolite into the extra-cellular compartment where they can then be consumed by the model. The bounds are only constraints on the maxi-mum rate of consumption of a given compound. The actual rate of consumption is predicted
<figure><img src="image_12.png" /><figcaption>**Fig 4.** Core and pan-metabolic capabilities of the *Leptospira* genus. The core and pan metabolic content was determined for genome-scale metabolic models of 4 different *Leptospira* species. A) Core content, illustrated by the intersection of the Venn diagram, shared with all species. The pan content consists of all content in any model and includes the core content. The Venn diagram is not to scale. B) Classification of reactions in the core and pan reactomes by metabolic subsystem.</figcaption></figure>

doi:10.1371/journal.pntd.0004403.g004

by the model. Units are in mmol gDW-1 hr-1. Analogously to core and pan genomes, the reaction content of each model can be used to construct core and pan metabolic networks. The core network consists of those reactions occurring in all representatives of *Leptospira*, while the pan network consists of all reactions that can potentially occur in any individual *Leptospira* species (Fig 4A).

Major differences in the metabolic networks between the four species arose in amino acid metabolism, biosynthesis of cofactors and vitamins and carbohydrate metabolism (Fig 4B). These large groups were further divided into specific pathways, specifically those of porphyrin metabolism, folate metabolism, starch and sucrose metabolism, as well as phenylalanine and tyrosine metabolism. A large difference was observed in porphyrin metabolism, specifically for the biosynthesis of cobalamin (vitamin B12) in *L. interrogans* (see below).

The conversion of static metabolic network reconstructions into computable mathematical models allows computation of phenotypes based on the content of each reconstruction [104]. Thus, the four strain-specific reconstructed networks were converted into genome-scale metabolic models that allow for the computational/simulation prediction of phenotype. This set of genomic scale models (GEMs) allows for a meaningful interpretation of the content of each reconstruction and allows for the prediction of each strain's different metabolic capabilities [105]. Because reactions belonging to the amino acid metabolism subsystem made up the majority of reactions in the pan-reactome, it was hypothesized that these capabilities may reflect functional differences between different *Leptospira* species. To test this hypothesis, different minimal media formulations were created *in silico* and used to test each model's growth capabilities. The models predicted all of the *Leptospira* tested to be auxotrophic for the amino acids aspartate, histidine and asparagine as well as vitamins B₁ (thiamin) and vitamin K2 (menaquinone).

Beyond the auxotrophies predicted to be shared by all 4 *Leptospira* models, potential species-specific auxotrophies for other vitamins and amino acids were also identified. All of the strain-specific models were predicted to be auxotrophic for phenylalanine except for *L. interrogans*, which was predicted to have the enzyme prephenate dehydrogenase encoded for by *novF* that converts chorismate to prephanate, a precursor to tyrosine and phenylalanine. The model for *L. interrogans* lacks prephanate oxidoreductase, which would predict inability to convert prephanate to tyrosine. Only *L. kmetyi* was found to have the enzymatic machinery capable of synthesizing tyrosine from phenylalanine. Among these representative species of the three
clades, only *L. interrogans* was predicted to be an L-glutamate auxotroph due to the lack of L-glutamate oxidoreductase.

Additional major differences were observed between the pathogenic and the non-pathogen *Leptospira*. A major difference in the lysine biosynthesis pathway was observed for the models of the pathogens, *L. kmetyi* and *L. interrogans*. Only these models possessed the *dapABCDE* genes required to convert L-aspartate 4-semialdehyde to LL-2,6, diaminopimelate required for peptidoglycan and lysine biosynthesis. Therefore both *L. licerasiae* and *L. biflexa* were predicted to be LL-2,6, diaminopimelate auxotrophs. Furthermore, only the pathogens *L. interrogans* and *L. kmetyi* possessed a full folate (vitamin B9) biosynthesis pathway using as precursor guanosine 5'-triphosphate. *L. biflexa* and *L. licerasiae* could produce vitamin B₂ (riboflavin), but lack the reactions to convert it to folate including dihydroneopterin aldolase encoded for by *folB*; therefore, the models for *L. biflexa* and *L. licerasiae* were folate auxotrophs while the models for *L. kmetyi* and *L. interrogans* were not.

## Vitamin B12 biosynthesis

The vitamin B12 biosynthesis genes in infectious *Leptospira* are grouped into two clusters: cob I/III and cob II (Table 2). Though the exact number of reactions for each pathway in *Leptospira* remains unknown, in *Salmonella enterica* Typhimurium cob I comprises genes for the biosynthesis of adenosylcobinamide, cob II genes for the synthesis of the lower axial ligand 5,6-dimethylbenzimidazole (DMB) and a third cluster cob III, the nucleotide loop that joins DMB to the corrin ring to complete B12 biosynthesis. In infectious *Leptospira*, cob II is a five-gene cluster that includes three genes *cobTSC* that participate in the synthesis of DMB (and two genes that may or may not participate in B12 biosynthesis), suggesting that the first cluster encodes enzymes for the synthesis of adenosylcobinamide guanine diphosphate. The first 12 genes encode enzymes that participate in the synthesis of the corrin ring (cob I) whereas the last five, enzymes for the addition of the nucleotide loop (cob III). Intriguingly, cob III of the infectious species includes a gene *cbiZ* encoding an enzyme that participates in an alternative cobinamide salvage pathway first described in the archeon, strain Göl [106].

Cob I/III gene clusters in the sequenced pathogenic *Leptospira* vary in length, from 16 in *L. santarosai* Shermani 1342K<sup>T</sup> to 19 in *L. alexanderi* Manhoa 3L60<sup>T</sup> and *L. borgpetersenii* Javanica UI 0993. Presumably owing to repeated *in vitro* sub-culture in medium containing B12, several genes have been inactivated or deleted in the strains tested. For example, the Javanica UI 0993 cobI/III cluster contains a gene fragment resulting from a premature stop codon in a gene encoding a histidine phosphatase superfamily branch 1 (*hps_1*) protein present in all other pathogenic *Leptospira* including other *L. borgpetersenii* strains (Hardjo L550 and JB197); and, cob I/III in *L. alexanderi* 3L 60T, contains a disrupted cobyrinic acid a,c-diamide synthase, inactivated by a frame shift mutation. A gene encoding a flavodoxin reductase present in cobI/III of other infectious has been deleted in *L. santarosai* 1342KT and *L. alstoni* 80-412. The cob I/III clusters of *L. kmetyi* (group I) and all group II pathogens contain three genes, two encoding a putative cobalt transporter (*cbtBA*) and a gene encoding an additional *hps_1* protein (Table 2).

Cob II also varies in length among all infectious (pathogenic and intermediate pathogenic) *Leptospira*, from three genes in *L. broomii* Hurstbridge 5399 and *L. inadai* Lyme 10<sup>T</sup> to seven in *L. noguchii* Panama CZ214<sup>T</sup>. The genes comprising the cob II cluster in non-pathogenic *Leptospira* are found in two discrete clusters in non-pathogenic species (e.g., LEPBI_I2857 and LEPBI_I2858, LEPBI_I2938-LEPBI_I2940) in *L. biflexa* (Table 2), suggesting that homologs in pathogenic *Leptospira* were acquired *en bloc* after the divergence of pathogenic and non-pathogenic *Leptospira*.
Leptospiral glycobiology: structure and diversity of rfb/O-antigen loci, lipid A, and sialic acid biosynthesis-encoding regions

**General features of Leptospira rfb loci.** Lipopolysaccharide (LPS) has long been a major focus of leptospiral microbiology not because of its (low potency) endotoxigenic activities (see below in the lipid A section) and because, notably, leptospiral LPS is the basis for serovar identification and vaccine development [36, 107–109]. Because of the importance of LPS in leptospiral biology, we carried out a comprehensive analysis of the genomic locations, structures and neighborhood of leptospiral rfb loci, also known as the O-antigen loci, in 20 species of *Leptospira*.

Using previously described leptospiral O-antigen gene clusters as a guide [110–112], we identified and schematized all clades of *Leptospira rfb* loci (Fig 5; S2 Fig depicts the *rfb* locus in genomic context). Of the genomes representing 20 *Leptospira* species, 17 known serovars were compared. The O-antigen biosynthesis gene clusters were located in three different genomic locations and ranged in size from 3,768 bp (*L. wolffii* sv. Khorat) to 121,402 bp (*L. alexanderi* sv. Manhao3). This region in *L. wolffii* sv. Khorat is now the smallest predicted leptospiral *rfb* biosynthesis cluster, consisting of just 4 genes, replacing the locus of *L. licerasiae* sv. Varillal [32]. All pathogenic leptospiral species, and the intermediates *L. inadai*, *L. broomii* and *L. fainei* have their *rfb* loci located in the same genomic position, sandwiched between a copper-binding protein on the left and the ribosomal protein S6 on the right (S2 Fig). The same protein-encoding genes (viz., *MarR* and *DASS*) define the start and end of the O-antigen cluster, respectively. Of the *rfb* loci of pathogenic *Leptospira*, serovars Manhao 3, Javanica, and Pingchang were most similar in size and gene content (Fig 5). Notably, *L. broomii* [113] (and *L. fainei* [114–116], both of which are serovar Hurstbridge, had nearly identical *rfb* gene clusters, predicting that *L. broomii* would also be serovar Hurstbridge, and confirmed by serology [113]. The presence of a specific serovar in different species has been previously observed in isolates of both *L. interrogans* and *L. borgpetersenii* serovar Hardjo, which have highly similar gene content in their O-antigen biosynthetic loci [117]. The *rfb* gene cluster of saprophytic *Leptospira* is downstream of the gene encoding ribosomal protein S6, lacks DASS, and is smaller (median 60,710 bp vs. ~99,520 bp) than pathogenic *Leptospira rfb* gene clusters (S2 Fig). O-antigen gene loci in serovars Varillal and Khorat are located in a third location, between *murC* and *purK*, consistent with a novel branching in the phylogenetic tree.

Consistently, the downstream flanking genes in the *rfb* loci are far more conserved than the upstream genes (Fig 5; S2 Fig). This finding is especially true for the pathogenic serovars and between three of the five intermediate serovars represented. For pathogenic serovars and serovar Lyme, a conserved block of genes is involved in O polysaccharide processing via the Wzy-dependent pathway. This export system was also identified in two saprophytic serovars (e.g., Holland and Codice). Overall, 12 of the representative 20 species genomes encoded the Wzy-dependent system and one genome (*L. licerasiae* sv. Varillal only encoded a putative flippase (Wzx) with no identifiable Wzy ortholog. There were no orthologs of Wzz, the O-antigen chain length determinant, in the 20 genomes studied. Also conserved in the 3-prime region of only those serovars with the Wzy-dependent pathway is a gene encoding a protein with homology to *E. coli WcaJ/S. enterica WbaP*, which are members of the PHPT family of polyisoprenyl-phosphate hexose-1-phosphate transferases that function to transfer glycosyl-1-phosphate to a lipid undecaprenol carrier, initiating formation of the O-unit in O-antigen assembly. In *L. borgpetersenii* serovar Hardjobovis, this protein, encoded by *orfH13*, is an UND-pp-galactosyltransferase [118].

The other major pathway of O-antigen polysaccharide biosynthesis is the *Wzm/Wzt*-encoded or ABC-transporter dependent pathway [119, 120]. Six of the 20 representative genomes encoded orthologs of *Wzm* and *Wzt*. Two of the genomes were from intermediate
<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /><figcaption>**Fig 5. Structure of Leptospira rfb locus gene clusters.** The rfb region and beginning and ending CDSs (blue) 9 of pathogenic (A), 5 intermediate (B), and 6 saprophytic (C) representative Leptospira species were compared. rfb region CDSs are labeled by locus identifier and colored by functional role categories as noted in the boxed key. Gene symbols, when present, are noted above their respective genes. BLASTP matches between CDSs are colored by protein percent identity (see key).</figcaption></figure>

doi:10.1371/journal.pntd.0004403.g005
and four were from saprophytic groups, representing 5 known serovars (e.g. Hurstbridge, Sao-
paulo, Hualin, Hardjo type Went, and Patoc). The genome analysis did not provide a clear
indication of the export system used by serovar Khorat. One other known O-antigen biosyn-
thesis pathway the synthase-dependent pathway [119]. BLASTP searches of WbbE and WbbF
from the only known example of this pathway, from the plasmid-encoded O:54 antigen of S.
*enterica* serovar Borreze [121], failed to identify any homologs in the representative 20 *Leptos-
pira* genomes. It is possible that serovar Khorat uses a novel mechanism for O -antigen
biosynthesis.

A dTDP-rhamnose biosynthesis gene cluster, encoding *rfbABCD* was found in the con-
served 3-prime end of the predicted O-antigen biosynthetic gene clusters of only pathogenic
*Leptospira* spp. serovar Saopaulo, found in the saprophytic species *L. yanagawae*, encoded
homologs of all four of these genes, but in the order *rfbCBAD*, where the genes *rfbABC* appear
to have been inverted (Fig 5C). The genes *rfbAB* and *rfbC*, were found in a different location
with *rfbC* separated by several genes in serovars Hualin, Hardjo, and Patoc. These same isolates
lacked an *rfbD* homolog. Only *rfbAB* homologs were identified in serovar Hurstbridge.

**L. licerasiae-type surface polysaccharide cassettes.** We previously reported that *L. licera-siae*, which is antigenically unique, lacks the type of extremely large O-antigen biosynthesis region found in *L. interrogans* and nearly all other *Leptospira* [122]. Instead, the one serovar of *L. licerasiae*, Varillal, has a six-gene cluster with three glycosyltransferase genes between two normally adjacent, convergently transcribed genes: the *murC* gene involved in cell wall biosyn-thesis and *purK* gene of purine biosynthesis. *Leptospira wolffii* had a similar genomic *rfb* locus, again with a six-gene cluster positioned between *murC* and *purK* (Fig 5B); antigenic relatedness to *L. licerasiae* serovar Varillal remains to be confirmed experimentally. The first glycosyltrans-ferase in this cassette, LEP1GSC185_2122 (GenBank EIE02925) in *L. licerasiae* and LEP1GSC061_3728 (GenBank EPG64090) in *L. wolffii*, are highly conserved and would be a useful marker for this extremely small O-antigen gene cluster. No other protein in the replace-ment six-gene cassette is conserved across the different variants. Genes in these regions have no close homologs in any other *Leptospira*, in the O-antigen region or anywhere else, support-ing the notion that these cassettes provide unique carbohydrate chemistry and serology, and is not simply an unusual gene neighborhood for otherwise common leptospiral enzymes.

**Lipid A biosynthesis.** The lipid A of leptospiral LPS is not as potent an endotoxin as lipid A moieties of other bacteria such as the Enterobacteriaceae or *Neisseria* spp.; the mechanistic explanation for this observation is that *L. interrogans* lipid A has different acyl chains and novel phosphorylation on the position of the lipid A that abrogate endotoxinogenicity [123]. The lipid A biosynthetic pathway of *L. interrogans* serovar Lai involves 13 enzymes, encoded by genes *lpxA*, *lpxC*, *lpxD1*, *lpxD2*, *lpxB1*, *lpxB2*, *lpxK*, *kdtA*, *kdsB1*, *kdsB2*, *lnt*, *kdsA* (also found as *waaA*) and *htrB*. The presence and homology of amino acid (aa) sequences of these enzymes was compared between 21 different species and/or serovars of *Leptospira* spp classi-fied in three different groups: pathogenic (PT, 10 species), intermediate (IM, 5 species) and non-pathogenic or commensal (NP, 6 species). Most proteins were found in all *Leptospira* spe-cies (S6 Table). However *lpxB2*, was found only in 4 pathogenic species/serovar and 1 non-pathogenic, *lpxD2* was not found in intermediate species/serovar and *htrB* was only present in 1 pathogenic and 1 non-pathogenic species/serovar. The *kdsB1* and *kdsB2* were only found in two species/serovar (*L. interrogans* sv. Lai and *L. inadai* sv. Lyme), all other species/serovar had only one *kdsB* that showed a higher level of similarity with *kdsB2* from *L. interrogans* sv. Lai than with *kdsB1*. Although we found that some genomes lack one or two lipid A biosynthetic genes (e.g. *lpxD2* and *kdsB2*), the computation analysis is still consistent with functional bio-synthetic pathways still being present in all species, because, for the genomes lacking one of the duplicate genes, the remaining ones (e. g. *lpxD1* and *kdsB1*) may be able to complement the
function of the lipid A biosynthetic pathway. Another possibility is that the genes are present
in the genomes, but we missed the genes because of gene diverge or gaps in genome sequence
obtained. Finally, the variable presence of lipid A biosynthesis genes may relate to some as yet
undiscovered structural differences in lipid A moieties among *Leptospira*.

The predicted protein sequences of individual lipid A biosynthesis pathway were nearly
identical among *Leptospira* as predicted using an identity matrix (S7A Table). The homology
between two sequences is expressed within the range of 0 to 1 (identical or completely homolo-
gous). The results presented hereby are expressed as the mean of homology values within each
group, compared to the pathogenic species group. The *lpxA* amino acid sequence was found in
all species, although the average similarity within pathogens was 0.928, while the homology of
intermediates and saprophytes was 0.694 and 0.581, respectively, when compared to pathogen
sequences. This analysis was carried out for each amino acid sequence (S7B Table).

**Sialic acids as post-translational modifications restricted to pathogenic Leptospira.**

Previous studies have demonstrated that pathogenic *Leptospira* endogenously synthesize Neu5Ac, the most common sialic acid, and that an observed gene fusion event suggested that *L. interrogans* uses a Neu5Ac biosynthetic pathway that is more similar to that of animals than to other bacteria. Lectin-based affinity purification of NulO-modified molecules, followed by mass spectrometric identification suggested post-translational modification of surface lipoproteins, including the putative virulence factor Loa22 [124, 125]. In the genomes analyzed for this study, 3 of the 9 pathogens had the complete cluster of genes involved in the production of sialic acids; 3 more lacked 1 gene in the cluster (Table 3). *L. weilii* contains only 2 genes from the cluster (*spsE* and *rfbB3*) and *L. kirschneri* and *L. noguchii* have only the *spsE* gene

All genomes, except *L. licerasiae* and *L. wolffii* have a N-acetylneuraminic (sialic) acid synthetase (*spsE*) gene (NP_711790.1). Phylogenetic analysis of this protein shows 2 distinct groups (S3 Fig). The first group contains the proteins from pathogens that contain the whole cluster. These proteins are related to the synthetases involved in the production of legionaminic acids. The second group contains the proteins from the intermediate species, the saprophytes and the pathogens *L. kirschneri* and *L. noguchii*. This group of synthetases is related to those producing pseudaminic acids.

The lack of a second sialic acid synthetase (NP_711794.1) in *L. kirschneri* and *L. noguchii* differentiates these pathogens from *L. interrogans*, which does contain this gene. These synthetases contain a phosphatase domain in addition to the NeuB domain, which suggests an animal-like Neu5Ac biosynthetic pathway. The pathogen *L. weilii* lacks NP_711794.1 but saprophyte *L. vanthielli* contains a similar synthetase but one that is missing the N-terminal transferase domain present in leptospiral pathogens. Finally, a UDP-N-acetylglucosamine diphosphorylase (NP_714003.1) was found in all leptospiral genomes studied. This gene is not located within the sialic acid gene cluster, and is also annotated as a MobA-like NTP transferase domain, therefore its role in sialic acid biosynthesis is unclear.

The sialic acid biosynthetic genes in leptospiral pathogens have some notable characteristics. *L. alexanderi* lacks O-acyltransferase (*neuD*) and this species and *L. borgpetersenii* have a truncated version of a nucleoside-diphosphate-sugar epimerase (NP_711787.1) (S3 Fig; Table 3). Only *L. santarosai* has a N-acetylneuraminic (sialic) acid synthetase with a phosphoglycerate dehydrogenase domain. Notably, none of the intermediate or saprophyte species contain the metabolic machinery to synthesize sialic acids, confirming previous suggestions [124].

Leptospiral mobile elements: phage and CRISPR-Cas systems

**Phage.** Bacteriophages are abundant biological entities that have significant effects on bacterial evolution. Some estimates suggest that there are approximately ten-fold more phages
than bacteria [126, 127]. However, our current knowledge of phages infecting *Leptospira* spp. is limited. Three distinct *L. biflexa* phages were isolated from sewage water in Paris. Morphological analysis by electron microscopy revealed that these three phages belong to the *Myoviridae* family and seem to be morphologically similar with polyhedral heads and contractile tails [128]. One of these phages, the 74-kb LE1 prophage, was shown to replicate as a double-stranded circular replicon in *L. biflexa* [129]. The genome of LE1 has a GC content of 36%, similar to that of *Leptospira* spp., and most of the 79 predicted ORFs display no similarity to known ORFs, but 21 ORFs appeared to be organized in clusters that might encode head and tail structural proteins and immunity repressor proteins [130].

Next generation sequencing and refinement of computational methods have allowed comparative genome analysis to discover new prophage and genomic islands [131]. A few phage related genomic islands have thus been characterized in *L. interrogans* and *L. licerasiae* [122, 132, 133] and it was previously shown that one of these genomic islands can excise from the *L. interrogans* chromosome [133].

To determine the distribution of prophages within the *Leptospira* genus, Phage_Finder [60] was run under both strict (-S) and non-strict modes to identify predicted prophage regions. Phage_Finder predicted a total of 14 major prophage regions across the 20 genomes, most of which were found to be shared between the *Leptospira* species (Table 4). Among the prophage sequences, the LE1-like prophage is found in many genomes, suggesting that double-stranded DNA tailed phages, which are the most frequently observed phages in bacteria [134], are common phages infecting *Leptospira*. The presence of numerous phage-associated sequences in the genome of pathogens and intermediates, in comparison to the saprophytes, suggests that phages have played an important part in the evolution of these lineages, as has been experimentally shown in *L. biflexa* [135].

Further experimental studies of *Leptospira* phages, which would include both electron microscopic visualization and production of phage *in vitro*, will be important to determine whether recombinant *Leptospira* phage might be useful for genetic manipulation studies of different *Leptospira* species, particularly pathogens and intermediates.

**CRISPR/Cas systems.** Three described types of CRISPR/Cas systems are common in *Leptospira* genomes and only found in infectious members of the genus (Table 5): the *E. coli* (type I-E), DVULG (type I-C), and MYXAN systems [96, 136, 137]. A single sequenced genome, *L. inadai* serovar Lyme str. 10, has the recently described PreFran type, which has been found in *Prevotella* and *Francisella* [96]. Four of the 20 representative strains have components of two CRISPR/Cas types, suggesting that for some isolates there is redundancy in CRISPR/Cas machinery. Surprisingly, half of the 20 representative *Leptospira* strains contained predicted CRISPR repeats, which were more common in pathogens and intermediates than in saprophytes, which lacked CRISPR systems. In none of the six saprophytes examined—*L. biflexa*, *L. meyeri*, *L. wolbachii*, *L. vanthielii*, *L. yanagawae*, and *L. terpstrae*—were CRISPR/Cas systems detected, suggesting that these species rely on some other mechanism for escaping phage/plasmid attack. In these saprophytes, we were also unable to detect sequences encoding prophage, while CRISPR systems and prophage occurred together in several, but not all, representative pathogenic and intermediate strains.

When present, between one and six CRISPR repeat arrays were detected, containing between three and 25 spacer sequences (Table 5). Since CRISPR spacer sequences in other organisms are known to target phage sequences for destruction, we wondered if any of the 239 predicted spacer sequences targeted any of the known *Leptospira* spp. phage or predicted prophages. A database containing the nucleotide sequences of the 19 predicted prophages from this study plus the LE1 phage and the prophage from Qin et al. [32, 132] was constructed and used to search all 239 predicted spacer sequences using BLAST+ 2.2.30 [138]. Upon filtering
<table><caption>Table 5. CRISPR/Cas Systems Identified in Leptospira Species Representatives.</caption><thead><tr><th rowspan="2">Species</th><th rowspan="2">Lifestyle</th><th rowspan="2">Contig</th><th rowspan="2">Direct Repeat Consensus</th><th rowspan="2"># Spacers</th><th rowspan="2">Location</th><th colspan="5">CRISPR System Genome Properties Evidence*</th></tr><tr><th>Dvulg (type I-C)</th><th>Ecoli (type I-E)</th><th>Myxan</th><th>PreFran</th></tr></thead><tbody><tr><td rowspan="3">alexanderi</td><td rowspan="3">pathogenic</td><td>AHMT02000065</td><td>CGGTTCAACCCACGCATGTSGAGAATAG</td><td>16</td><td>6749..7777</td><td>0</td><td>1</td><td>0</td><td>0</td></tr><tr><td>AHMT02000052</td><td>CGGTTCAACCCACGTGTGGGAAAG</td><td>25</td><td>65616..67218</td><td></td><td></td><td></td><td></td></tr><tr><td>AHMT02000009</td><td>... AACCCACGTGTGGGAAAG</td><td>6</td><td>11327..11714</td><td></td><td></td><td></td><td></td></tr><tr><td>alstoni</td><td>pathogenic</td><td>AOHD02000029</td><td>GTGCTCAACGCCTAACGGCATCAAAGGTATGTTCAG</td><td>14</td><td>329..1366</td><td>0</td><td>0.125</td><td>0.444</td><td>0</td></tr><tr><td rowspan="2">borgpetersenii</td><td rowspan="2">pathogenic</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>0</td><td>0.875</td><td>0</td><td>0</td></tr><tr><td>AE016823.1</td><td>GTGCTCAACGCCTAACGGCATCAAAGTTATATTTCAG</td><td>3</td><td>1138348..1134101</td><td>1</td><td>0</td><td>0.889</td><td>0</td></tr><tr><td rowspan="2">interrogans</td><td rowspan="2">pathogenic</td><td>AE016823.1</td><td>TTCCTAAGAAATAGGGAAATTTAAAAA</td><td>4</td><td>1451041..1451345</td><td></td><td></td><td></td><td></td></tr><tr><td>AHMN02000020</td><td>GTGCTCAACGCCTTACGGCATCAAAGTTAAGTTTCAG</td><td>5</td><td>13791..14191</td><td>1</td><td>0</td><td>0.556</td><td>0</td></tr><tr><td rowspan="2">kirschneri</td><td rowspan="2">pathogenic</td><td>AHMN02000010</td><td>ATTTCTTTCTCTTTAAGAGGAAGTGGATTGAAAC</td><td>3</td><td>88470..88720</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>AKWY02000022</td><td>GTTTCAATCCACTCTTCTCTAARAAAGAGAGATTAAAC</td><td>12</td><td>56305..57198</td><td>1</td><td>0</td><td>0</td><td>0</td></tr><tr><td rowspan="4">knetyi</td><td rowspan="4">pathogenic</td><td>AKWY02000014</td><td>AKWY02000014</td><td>GTTTCTTTCTCTTTGAGGAAAGTGGATTGAAAC</td><td>5</td><td>7603..7996</td><td></td><td></td><td></td></tr><tr><td>AKWY02000003</td><td>AKWY02000003</td><td>GTTTCTTTCTCTTTGAGGAAAGTGGATTGAAAC</td><td>10</td><td>55925..56687</td><td></td><td></td><td></td></tr><tr><td>AKWY02000021</td><td>AKWY02000021</td><td>GTTTCAATCCACTTCCTTTTAAACCAAGGAAGAAC</td><td>7</td><td>48415..48952</td><td></td><td></td><td></td></tr><tr><td>AKWY02000021</td><td>AKWY02000021</td><td>GTTTCTTTCTCTTTACTTTAGAGGAAGTGGATTGAAAC</td><td>8</td><td>315674..316283</td><td></td><td></td><td></td></tr><tr><td rowspan="4">noguchii</td><td rowspan="4">pathogenic</td><td>AOHB02000026</td><td>GTGTTCAACCCACGCATGTGGGAAAG</td><td>5</td><td>19990..20321</td><td>0</td><td>1</td><td>0.889</td><td>0</td></tr><tr><td>AOHB02000050</td><td>CGGTTCAACCCACGCATGTGGGAAAG</td><td>9</td><td>16484..17061</td><td></td><td></td><td></td><td></td></tr><tr><td>AOHB02000050</td><td>GTGCTCAACGCCTTACGGCATCAAAGTTATATTTCAG</td><td>3</td><td>50240..50492</td><td></td><td></td><td></td><td></td></tr><tr><td>AOHB02000050</td><td>..TGCTCAACGCCTTACGGCATCAAAGTTATATTTCAG</td><td>3</td><td>52981..53232</td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">weili</td><td rowspan="2">pathogenic</td><td>AOHB02000003</td><td>CTTTTCCCCACATGCGTGGGGTTGAACCG</td><td>3</td><td>56764..56975</td><td></td><td></td><td></td><td></td></tr><tr><td>AHNC02000050</td><td>CTTTTCCCCACATGCGTGGGGTTGAACCG</td><td>7</td><td>80411..80865</td><td>0</td><td>0.875</td><td>0</td><td>0</td></tr><tr><td rowspan="2">broomii</td><td rowspan="2">intermediate</td><td>AHNC02000082</td><td>CTTTTCCCCACACACGTGGGGTTGAACCG</td><td>16</td><td>222247..223251</td><td></td><td></td><td></td><td></td></tr><tr><td>AHMO02000008</td><td>CTGAAACTAACTTTGATGCCGAAAGCGTGTTGAGCAT</td><td>11</td><td>908793..909615</td><td>0</td><td>0</td><td>1</td><td>0</td></tr><tr><td rowspan="2">fainei</td><td rowspan="2">intermediate</td><td>AKWZ02000010</td><td>CGGTTCATCCCCACGTGTGGGAAATAG</td><td>5</td><td>1694300..1694633</td><td>0</td><td>1</td><td>0</td><td>0</td></tr><tr><td>AKWZ02000010</td><td>CGGTTCATCCCCACGTGTGGGAAATAG</td><td>6</td><td>1695946..1696339</td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">inadai</td><td rowspan="2">intermediate</td><td>AKWZ02000004</td><td>CTATTCCCCACACGTGGGGATGAACCG</td><td>18</td><td>200759..201886</td><td></td><td></td><td></td><td></td></tr><tr><td>AHMM02000017</td><td>GTTTCAATTCAGATGGTTCGATTAAAG</td><td>22</td><td>566571..568061</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><td>licerasiae</td><td>intermediate</td><td>NA</td><td>ATCTACAAAAGTAGAAATTCTTCCTCTCTTTAGAG</td><td>9</td><td>625117..625713</td><td></td><td></td><td></td><td></td></tr><tr><td>wolffii</td><td>intermediate</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>biflexa</td><td>saprophytic</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>meyerii</td><td>saprophytic</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>terpstraes</td><td>saprophytic</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>vanthieli</td><td>saprophytic</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>wolbachii</td><td>saprophytic</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>yanagawae</td><td>saprophytic</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>0</td><td>0</td><td>0</td><td>0</td></tr></tbody></table>

†Predicted using CRISPRFinder (Grissa et al., 2007).

* Evidence value ranges from 0 (absent) to 1 (present; all steps of the genome property were found)

NA = Not applicable because no "confirmed" CRISPR repeats were identified

doi:10.1371/journal.pntd.0004403.t005
L. noguchii_AKWY02000022_spacer10: 56913-AAGATGACCTTGCGGGACATCGCCTTTCTTATTCT-56948

L. santarosai_AOHB02000043PFPR01: 15582-AAGATGACCTTGCGGGACATCGCCTTTCTcTATTCT-15547

L. noguchii_AKWY02000021_R3_spacer1: 423886-AAGATGAGAACGGAAAACCAAAATACTCTTTCAA-423919

L. santarosai_AOHB02000043PFPR01: 22939-AAGATGAGAACGGAAAACCAAAATACTCTTTCAA-22972

L. noguchii_AKWY02000021_R3_spacer1: 423886-AAGATGAgAACGGAAAACCaAAATACTCTTTCAA-423919

L. weilii_AHNC02000035PFPR01: 49925-AAGATGAcAACGGAAAACCgAAATACTCcTTCAA-49892

L. weilii_AHNC02000082_spacer2: 222338-ATACAGCGATTACATTCGGATTCAGTATTTTT-222369

L. weilii_AHNC02000035PFPR01: 50179-ATACAGCGATTACATTCGGATTCAGTATTTTT-50148

**Fig 6.** CRISPR Spacer Sequences that Recognize *Leptospira* Predicted Prophages. The CRISPR sequences are shown, which correspond to specific prophage accession numbers as listed in Table 4.

doi:10.1371/journal.pntd.0004403.g006

the data for matches spanning the entire spacer sequence, with 3 or fewer mismatches and with a bitscore of 30+ revealed three spacer sequences matching two predicted prophage sequences (Fig 6). Two different *L. noguchii* spacers matched the same predicted *L. santarosai* prophage (Table 5). One of the same *L. noguchii* spacers also matched a predicted *L. weilii* prophage and was also recognized by an *L. weilii* spacer (Fig 6).

## Virulence and Survival Mechanisms

**Adhesion to Extracellular Matrix (ECM).** The presence of genes encoding putative adhesive proteins through the 20 sequenced species was analyzed by BLAST and comparative genome analysis (S8 Table). The widespread distribution of these genes within the *Leptospira* genus suggests that their functions arose independent of mammalian adaptation, but any potential role in adaptation to an environmental lifestyle remains speculative. These predicted adhesion-related proteins were generally distributed among the 20 species except for the predicted adhesin-encoding gene *LenB*, identified only in pathogenic *L. interrogans* serovars Lai and Copenhagen strains; in the saprophyte *L. meyeri*; and *Lsa27*, *Lsa21*, *LipL53* present in two, three and five pathogen species, respectively. Three predicted adhesin-encoding genes were restricted to pathogenic species: *Lsa30*, *Lsa44* and *Mfn6*. *Lsa30* was present in all species, *Lsa44* was absent in *L. interrogans* serovar Lai and *L. weilii*, and *Mfn6* was absent in *L. weilii* and *L. alexanderi*. The genes encoding *Lsa23*, *Lsa26*, *Lsa33*, *Lsa45*, *Lsa66*, *LipL32* and *Mfn1* were found in all infectious species (pathogens and intermediates) but absent in saprophytes. The *Len* protein family members were variably distributed. *Lsa24/LenA* was identified in all sequenced strains, while *LenB*, *LenC*, *LenD*, *LenE* and *LenF* were found in both pathogens and saprophytes; *LenD* was also found in the intermediate *L. wolffii*. *Lsa23*, *Lsa26*, *Lsa33*, *Lsa45*, *Lsa66*, *LipL32* and *Mfn1* were found in all pathogenic and intermediate but are absent in saprophyte strains. The *Len* protein family showed a random distribution among the genome
species: *Lsa24/LenA* was identified in all sequenced strains, while *LenB*, *LenC*, *LenD*, *LenE* and *LenF* were found in both pathogenic and saprophytic *Leptospira*, and *LenD* was also found in *L. wolffii*, an intermediate. The adhesins *Lsa20*, *Lsa25*, *Lsa36*, *Lsa63*, *TlyC*, *OmpL1*, *OmpL37*, *OmpL47*, *Mfn7* and *rLIC12976* were identified in all pathogenic, intermediate and saprophyte species; *Mfn9* was found in all except in *L. santarosai*. Also listed are plasminogen- and complement regulator-binding proteins (S8 Table) that have functions related to the predicted proteins listed above.

**Complement evasion and ECM degradation via metalloproteases.** Two leptospiral proteases have been suggested as virulence factors: thermolysin and collagenase. Thermolysins are members of the M4 metalloprotease family that can be identified bioinformatically by the presence of two N-terminal propeptide (FTP, and PepSY) and two C-terminal protease domains (Peptidase_M4 and Peptidase_M4_C) [139]. Using Pfam HMMs targeting these four domains (e.g., PF07504, PF03413, PF012868 and PF01447), we identified LIC13322 and four additional predicted thermolysin orthologs only in pathogenic *Leptospira* spp.: LIC10715, LIC13320, LIC13321, and LEP1GSC059_0182 (S9 Table), primarily among *L. interrogans*, *L. kirschneri* and *L. noguchii*. LEP1GSC059_0182 was found only in one species, *L. noguchii*. No thermolysin ortholog was found in intermediate or saprophytic species.

Collagenase has been suggested to be a virulence factor in *Leptospira* based on observed *in vivo* expression, detection of specific anti-collagenase antibodies induced by infection, and the effects of *ColA* mutagenesis and complementation on traversal of cell monolayers and outcome of experimental animal infection [140, 141]. Comparative genomic analysis identified two collagenase genes, restricted to pathogenic *Leptospira* spp.: orthologs of one (LIC_12760) were found in all pathogens except *L. kmetyi*; an additional paralog, EMN46521, was restricted to *L. weilii* and *L. alexanderi*, and based on nearly identical size and closely related amino acid sequences, likely arose by gene duplication (S9 Table). The implications of this latter finding for pathogenesis are unclear.

**Resistance to oxidative stress.** Three enzyme systems have conventionally been associated with the ability of pathogenic bacteria to defend against host-derived oxidative stress-related mediators such as hydrogen peroxide and superoxide radicals: catalases, peroxidases and superoxide dismutase. While catalases generate water and oxygen from H₂O₂, peroxidases generate water and an oxygen radical.

*KatA* and another but uncharacterized predicted catalase ortholog (LEP1GSC062_4039) were only found in pathogenic *Leptospira*, suggesting that these enzymes play an important role for *Leptospira* living within the mammalian host (Table 6). Conversely, superoxide dismutase was not found in pathogenic *Leptospira* but only in saprophytic *Leptospira*, suggesting either that pathogenic *Leptospira* are not exposed to oxygen radicals in the environment, or, more likely, that this clade of *Leptospira* has developed alternative ways to detoxify oxygen radicals.

**Immunodominant proteins of *Leptospira*.** Previously published protein microarray analysis demonstrated the presence of immunodominant proteins of *L. interrogans* serovar Cophenhageni using sera from confirmed leptospirosis cases in Bahia state, Brazil. The top 24 immunogenic hits from this analysis were analyzed throughout the genus *Leptospira*, focusing on the presence of orthologs and their amino acid similarities (S10A Table). Only 1 of these genes (a methyltransferase, NC_005823.1) was restricted to *L. interrogans*, and only 1 (Lig A) was found in only *L. interrogans* and *kirschneri*, with the caveat that a 56% homologous LigA domain was found in *L. alstoni*, a leptospire with unclear disease potential. Orthologs of 20 of these 24 hits were detected in all *Leptospira* species but with variable amino acid similarities, suggesting that species- (and perhaps serovar-) specific protein microarrays might be necessary for accurate assessment of immune responses induced by different *Leptospira* in humans.
<table><caption>Table 6. Proteins involved in the resistance to oxidative stress.</caption><thead><tr><th>Cluster</th><th>Representative</th><th>GenBank Accession</th><th>Function</th><th>L. interrogans sv. Copenhageni st. Floucazi.1-190</th><th>L. kirschneri sv. Cyprinoperli st. 3522 CT</th><th>L. noglechi sv. Panama st. 214T</th><th>L. alstoni sv. Pringotang st. 80-4-2</th><th>L. weillav. undetermined st. LNT 1254</th><th>L. alexanderi sv. Manhao st. 3. L 607</th><th>L. borgpetersenii sv. Jovenica st. U10691</th><th>L. santarosai sv. Shermani st. 1342KT</th><th>L. ametyi sv. Malaysia at. Bepo-ho9T</th><th>L. fairei sv. Hurstbridge st. BUT 6T</th><th>L. broomia sv. Hurstbridge st. 5399T</th><th>L. wolffii sv. Khorat at. Khonak-H2T</th><th>L. licerassae sv. Varillii at. VAR 010</th><th>L. woldi sv. Lyne st. 107</th><th>L. woldi sv. Collete at. CDC</th><th>L. yanagawae sv. Ssao at. Paulo</th><th>L. belfza sv. Ssao at. Paulo</th><th>L. vanthelli sv. Holland at. Waz</th><th>L. terestra sv. Hulatin at. LT 1-137</th><th>L. mayeri sv. Hulatin at. HT 1-137</th></tr><tr><th></th><th></th><th></th><th></th><th>P</th><th>P</th><th>P</th><th>P</th><th>P</th><th>P</th><th>P</th><th>P</th><th>I</th><th>I</th><th>I</th><th>I</th><th>I</th><th>S</th><th>S</th><th>S</th><th>S</th><th>S</th></tr></thead><tbody><tr><td>LEP1/GSC202_0889</td><td>EQA63414</td><td>caliolase (KAA)</td><td></td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>LEP1/GSC195_0257</td><td>EQA67809</td><td>caliolase/cytochrome (KAG)</td><td></td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>LEP1/GSC202_4039</td><td>EQA67078</td><td>caliolase</td><td></td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>LEP1/GSC206_1809</td><td>EMN42918</td><td>di-kerne/cytochrome C peroxidase</td><td></td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>LEP1/GSC206_3491</td><td>EMN4421</td><td>di-kerne/cytochrome C peroxidase</td><td></td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>LEP1/GSC205_4289</td><td>EQA52716</td><td>di-kerne/cytochrome C peroxidase</td><td></td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>LEP1/GSC203_3912</td><td>EMY61934</td><td>di-kerne/cytochrome C peroxidase</td><td></td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>LEP1/GSC202_2886</td><td>EQA62707</td><td>glutathione peroxidase</td><td></td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>LEP1/GSC202_4239</td><td>EQA53474</td><td>glutathione peroxidase</td><td></td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>LEP1/GSC205_3788</td><td>EQA72367</td><td>glutathione peroxidase</td><td></td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>LEP1/GSC195_3278</td><td>EQO97317</td><td>superoxide dismutase</td><td></td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>

P = Pathogenic species; I = Intermediate species; S = Saprophytic species

doi:10.1371/journal.pntd.0004403.t006
Further, to validate such arrays with well-defined sera from leptospirosis cases, identification of the infecting leptospire will be necessary.

The extensively-studied leptospiral immunoglobulin-like (Lig) protein family is comprised of three proteins, LigA, LigB and LigC [142, 143], which have bacterial immunoglobulin-like (Big) repeat domains, a motif found in virulence factors of other bacterial pathogens [144, 145]. The three genes encoding Lig proteins were believed to be pathogen-specific [142, 146-148]. Comparative analysis of 20 genomes confirmed that *ligA* and *ligB* were present exclusively in pathogenic *Leptospira*: *ligB* was identified in all pathogenic species, *ligA* was found in three of the nine pathogenic species, *L. alstoni*, *L. kirschneri* and *L. interrogans* (S10B Table). While *ligC* was found in the five intermediate species and five of nine pathogenic species (S10B Table), none of the *lig* genes were identified in genomes of saprophytes. *ligC* was previously identified as a pseudogene from sequence analysis of a limited number of strains [142, 147, 148].

The unique structure of LigA, LigB and LigC proteins, which includes a large number of tandem Big domains, is conserved across species for which *lig* genes were found. *ligB* and *ligC* encode molecules which are comprised of a lipobox sequence, 12 tandem Big2 type domains and a C-terminal non-Big domain (S10B Table), whereas *ligA* encodes a protein with 13 tandem Big domains which lacks a C-terminal non-Big domain. Of note, further prospection of the genome sequence and ortholog families identified four additional genes that encode Big2 and Big3_4 domain-type containing proteins (S10B Table). These genes are different from the conventional *lig* gene family in that they encode proteins with a small (1-2) number of Big domains. Interestingly, a gene (LIC13050) encoding a protein with two Big3_4 domains was found in all *Leptospira* species, including the saprophytes.

**PF07598 paralogous gene family.** Previous work identified a group I-specific family of proteins corresponding to Pfam model PF07598 [149] that was expressed *in vivo* in a hamster model of acute leptospirosis [149], and expanded in strains, e.g., Copenhageni and Lai, that commonly cause severe disease, suggesting that these proteins contribute to *Leptospira* virulence. These prior studies focused on finding PF07598 orthologs to the *L. interrogans* Lai attenuated strain in the 20 representative *Leptospira* spp., but did not look for strain-specific homologs that match PF07598 HMM. To identify novel PF07598 family members, we identified clusters of protein orthologs from our pan-genome run that matched PF07598 above trusted cut-offs; this analysis also included matches within the genomes of previously sequenced strains *L. borgpetersenii* Hardjo L550 and JB197 (and previously annotated as conserved hypothetical proteins). At least 26 distinct orthologs ranging in length from 47 (LEP1GSC049_1303 unique to *L. kirschneri*) amino acids to 651 (LEP1GSC193_2756; *L. alstoni* 80-412) were identified (Table 7; only homologs longer than 200 amino acids shown). As previously reported, *L. santarosai* 1342KT contains two distinct homologs, while *L. kirschneri* 3522 CT at least 15 and *L. noguchii* CZ214 at least 14 including 5 (LEP1GSC059_0232, LEP1GSC059_3018, LEP1GSC059_3019, LEP1GSC059_3599 and LEP1GSC059_3600) without an apparent ortholog in any of the other strains tested (Table 7). *L. borgpetersenii* Javanica (4 total) and the previously sequenced genomes of two Hardjo strains (3) contain two orthologs in common, while Javanica has two distinct copies (LEP1GSC103_4030 and LEP1GSC103_0672) not present in Hardjo and both Hardjo strains share an ortholog (LBJ_1339 and LBL_1564) not present in Javanica. In addition, *L. interrogans* L1-130 contains an ortholog (LIC_10639) shared with *L. noguchii* CZ214T and *L. kirschneri* 3522 CT not present in *L. interrogans* 56601 consistent with the hypothesis that serovar Lai has lost this ortholog.

To better understand the evolution of this paralogous gene family, a phylogenetic tree of all PF07598 members detected in infectious *Leptospira* was constructed using homologs longer than 200 amino acids shown (Fig 7A). A complex web of lineage specific gene duplications and loss was revealed. For example, as highlighted in (Fig 7A), successive gene duplications and
<figure><img src="image_15.png" /><figcaption>**Fig 7. Phylogenetic Relationship of PF07598 Paralogous Family in Leptospira.** (A) Unrooted bootstrapped phylogenetic tree; (*) Gene duplication event; (**) gene duplication event; (***) gene deletion. (B) Principal components analysis was used to arrange PF07598 family members. Color legend indicates the PF07598 family members from specific serovars depicted as diamonds. Arrowheads indicate *L. noguchii*-specific orthologs. Only PF07598 family members longer than 200 amino acids are included in the analysis. Clusters (A, B and C) were defined by K-means clustering with Kendall rank correlation.</figcaption></figure>

doi:10.1371/journal.pntd.0004403.g007

subsequent gene loss have led to four distinct clusters containing proteins in *L. interrogans*, *L. kirschneri* and *L. noguchii* PF07598 family members. The initial duplication event led to divergence of the first orthologous group comprising LEP1GSC059_0224 (*L. noguchii*) and LIC_12985 and LA0591 (*L. interrogans*) and apparent loss of the corresponding ortholog in *L. kirschneri* (*). The second, an *L. interrogans*-specific event (**), led to the divergence of two orthologous groups comprised of LA_0589, LIC_12986 and LEP1GSC049_3370, and LA_3388, LIC_10778 and LEP1GSC049_0186, respectively. *L. interrogans* Lai 56601 has seemingly lost an ortholog belonging to a group containing LIC10639 and LEP1GSC049_1381
comprised of proteins belonging to *L. interrogans* Copenhageni and *L. kirschneri*, respectively; and *L. kirschneri* 3522CT, *L. alexanderi* and both *L. borgpetersenii* Hardjo strains seem to have lost an ortholog present in the other pathogenic strains (***). This pattern of species and sero-var specific gene duplication and gene deletion occurs throughout the tree. K-means clustering with Kendall rank correlation grouped the *L. interrogans* Lai orthologs into three clusters comprising family members with >90% inclusion probability: LA_1400 and LA_1402 (cluster A); LA_0589, LA_0591, LA_0835 and LA_3388 (cluster B); and LA_0620, LA_0769, LA_0934, LA_2628, LA_3271 and LA_3490 (cluster C) (Fig 7B).

**Motility and chemotaxis.** Since motility is required for pathogenesis [150, 151], it is plausible that there are differences in motility and chemotaxis gene content that distinguish infectious from non-infectious species. We identified a total of 76 CDSs encoding proteins involved in leptospiral motility and chemotaxis, using the annotated genome of the *Leptospira interrogans* serovar Copenhageni strain Fiocruz L1-130 as a reference. We established the amino acid sequence identity of CDSs in 20 *Leptospira* genomes based on their respective orthologs in the strain Fiocruz L1-130 genome (S4 Fig; S11 Table). Among these, 37 CDSs were predicted to encode proteins in the basal body assembly and export apparatus; 7 CDSs were predicted to encode proteins in the flagellar hook assembly; 7 CDSs were predicted to encode proteins involved in the filament assembly; and 25 CDSs were predicted to encode proteins in chemotaxis (S4 Fig; S11 Table).

Proteins involved in motility were highly conserved among all the 20 *Leptospira* species according to BLAST analysis and PanOCT ortholog clusters (S4 Fig). The filament is the portion of the flagella which demonstrated the highest amino acid sequence identity, with a mean of 97.9%, 86.4 and 72.4% amino acid sequence identity among pathogenic, intermediate and saprophytic species, respectively (S11 Table). The ORFs that encode the flagellar hook proteins also demonstrate high amino acid sequence conservation, with an average sequence identity of 86.0% and 61.7% in pathogenic and saprophytic species, respectively (S11 Table). FliK, a bi-functional protein involved in determining hook length and modulating export-pathway specificity at the hook-filament checkpoint [152, 153], was the only protein that showed a low level of identity among the three species groups, including within the pathogenic species (69.7%, S4 Fig and S11 Table).

Although CDSs encoding basal body proteins showed the lower amino-acid identity among motility genes, the average identity was high, ranging from 93, 73, and 60% within pathogenic, intermediate and saprophytic species (S4 Fig and S11 Table). In this category, three proteins showed 50% identity or lower when comparing pathogenic species with intermediates and saprophytes. The protein FlgA is involved in the P-ring formation, whereas the FliO and FliJ are involved in the export apparatus. In addition, five CDSs showed an amino acid sequence identity below 50% between pathogenic and saprophytic species, which were CDSs encoding proteins FlgH and FlgL involved in the L- and P-ring formation, respectively, proteins FlgN and FlhX involved in the export apparatus, and the FliG1 protein, which is involved in the motor switch. The P- and L-ring form the outer cylinder and acts as a bushing for the central rod [153, 154] and is believed to participate only passively in the motor mechanism, while the FliG1 protein is believed to be partly responsible for the asymmetrical rotation of the flagella [155].

CDSs encoding chemotaxis proteins are highly conserved among pathogenic species (87% amino acid sequence identity). In contrast, ORFs encoding such proteins are less conserved when comparing pathogenic species with intermediates and saprophytes species groups and have lower amino acid sequence identity (48 and 43%, respectively, S4 Fig and S11 Table). More than 70% of the orthologs of chemotaxis proteins within intermediate and saprophyte species had less than 50% amino acid sequence identity when compared to pathogenic species. Among these proteins, the majority were methyl-accepting protein (MCP) homologs, but
include also chemotaxis regulators like *cheA*, *cheR*, *cheB* and *cheY* (S4 Fig and S11 Table).

Whereas orthologs associated with construction of the flagellar filament are conserved across pathogenic, intermediate and saprophytic species, ORFs encoding two chemotaxis proteins, one MCP homolog and cheR1, in pathogenic species had no orthologs in intermediate and saprophytic species. Furthermore, one MCP homolog protein in pathogenic and intermediate species had no orthologs in saprophyte species, suggesting a degree of divergence with respect to chemotaxis between pathogenic, intermediate and saprophyte *Leptospira* species.

## Gene regulation and sensory transduction

**Alternative sigma (σ) factors.** σ factors are a class of proteins constituting essential dissociable subunits of prokaryotic RNA polymerase. σ factors provide promoter recognition specificity to the polymerase and contribute to DNA strand separation. All bacterial species have a housekeeping σ-factor (σ⁷⁰) responsible for transcription from the majority of promoters. Most bacteria encode additional alternative σ-factors that redirect RNAP to distinct sets of promoters, which can contribute both directly and indirectly to environmental adaptation and bacterial virulence. In addition to a housekeeping sigma factor σ⁷⁰ (LIC11701, RpoD), all *Leptospira* species have an alternative sigma factors σ⁵⁴ (LIC11545, RpoN) involved in nitrogen and many cellular and environmental regulations, σF involved in flagella gene expression (LIC11380 (FliA), 5–11 extracytoplasmic function (ECF) σ factors (σE) involved in regulation of membrane and periplasmic stress, and more than 30 anti-σ regulators (S12 Table).

Leptospiral species differ in σ-factors. First, pathogenic *Leptospira* have two activators (enhancer-binding protein, EBP) for σ⁵⁴, whereas saprophytic *Leptospira* species has only one. σ⁵⁴ is a unique sigma factor that is phylogenetically different from other σ actors. It recognizes a unique –24/–12 promoter sequence (instead of -35/-10 sequence for σ⁷⁰) and its activation always requires an activator, EBP. Signals feed into EBP and activate σ⁵⁴-dependent genes. Each EBP-σ⁵⁴ pairs responds to different signals and activates a set of genes. Our analyses show that while all pathogenic *Leptospira* appear to have two activators (herein named as *Leptospira enhancer-binding protein A* and *B* that can be denominated EBP-A and EBP-B, saprophytic *Leptospira* have only one EBP (EBP-A) (S12 Table). Although the upstream signals and downstream targets remain to be elucidated, we speculate that LepA-σ⁵⁴ modulates a group of genes involved in environmental survival for both pathogenic and saprophytic *Leptospira*, whereas LepB-σ⁵⁴ is important for pathogenic *Leptospira* species to adapt to host environment. Second, pathogenic and saprophytic leptospiral species differ in ECF σ factors. Pathogenic and intermediate pathogenic *Leptospira* have 9 to 10 ECFs, and saprophytic species often have 5 ECFs. One ECF (LIC10599) is only found in highly pathogenic *Leptospira*, while 2 ECFs are only associated with saprophytic *Leptospira* (S5 Fig; S12 Table). Lastly, *Leptospira* have more than 30 regulators predicted as anti-σfactors, anti-anti-σfactors, and regulators of anti-anti-σfactors. Although their functions remain unclear, some of these regulators may modulate ECF functions as observed in *B. subtilis*. Nevertheless, there are some obvious differences in their distributions among *Leptospira* species (S12 Table). It is conceivable that the ECFs and regulators of σ factors present only in saprophytic *Leptospira* are involved in responding to environmental stress, whereas the ECFs and regulators of σ factors present only in pathogenic *Leptospira* are likely important for *Leptospira's* life cycle in mammalian hosts.

**Two Component Systems (TCS).** TCSs are the predominant molecular switches controlling signaling events in bacteria. Typically, TCSs consist of a sensor histidine kinase (HK) and an effector response regulator (RR). A single polypeptide merging both components results in hybrid histidine kinases (HHKs). HKs and RRs are usually found adjacent to each other in the genome. Orphan TCS proteins are unpaired HKs/RRs, which work with their cognate partners
<table><caption>Table 8. Summary of two component systems identified in 20 *Leptospira* genomes.</caption><thead><tr><th>Species</th><th>Non-Orphan HK</th><th>Non-Orphan RR</th><th>Orphan HK</th><th>Orphan RR</th><th>Hybrid HK</th><th>Hpt</th><th>Total</th><th>TCS to genome size <sup>a</sup></th></tr></thead><tbody><tr><td><strong>Pathogenic</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td><em>L. interrogans</em></td><td>12</td><td>12</td><td>13</td><td>22</td><td>12</td><td>3</td><td><strong>74</strong></td><td>15.99</td></tr><tr><td><em>L. kirschneri</em></td><td>9</td><td>9</td><td>12</td><td>24</td><td>12</td><td>3</td><td><strong>69</strong></td><td>15.64</td></tr><tr><td><em>L. noguchii</em></td><td>12</td><td>12</td><td>12</td><td>23</td><td>12</td><td>3</td><td><strong>74</strong></td><td>15.70</td></tr><tr><td><em>L. alstoni</em></td><td>15</td><td>15</td><td>14</td><td>24</td><td>14</td><td>4</td><td><strong>86</strong></td><td>19.38</td></tr><tr><td><em>L. weilli</em></td><td>13</td><td>13</td><td>11</td><td>19</td><td>11</td><td>3</td><td><strong>70</strong></td><td>16.41</td></tr><tr><td><em>L. alexanderi</em></td><td>11</td><td>11</td><td>12</td><td>20</td><td>11</td><td>3</td><td><strong>68</strong></td><td>16.09</td></tr><tr><td><em>L. borgpetersenii</em></td><td>12</td><td>12</td><td>12</td><td>21</td><td>11</td><td>3</td><td><strong>71</strong></td><td>18.22</td></tr><tr><td><em>L. santarosai</em></td><td>13</td><td>13</td><td>11</td><td>20</td><td>11</td><td>3</td><td><strong>71</strong></td><td>17.80</td></tr><tr><td><em>L. kmetyi</em></td><td>18</td><td>18</td><td>18</td><td>25</td><td>19</td><td>3</td><td><strong>101</strong></td><td>22.85</td></tr><tr><td><strong>Intermediate</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td><em>L. fainei</em></td><td>19</td><td>19</td><td>15</td><td>27</td><td>18</td><td>4</td><td><strong>102</strong></td><td>23.79</td></tr><tr><td><em>L. broomii</em></td><td>21</td><td>21</td><td>14</td><td>26</td><td>21</td><td>4</td><td><strong>107</strong></td><td>24.34</td></tr><tr><td><em>L. wolffii</em></td><td>18</td><td>18</td><td>19</td><td>30</td><td>16</td><td>4</td><td><strong>105</strong></td><td>23.85</td></tr><tr><td><em>L. licerasiae</em></td><td>15</td><td>15</td><td>14</td><td>36</td><td>16</td><td>5</td><td><strong>101</strong></td><td>23.98</td></tr><tr><td><em>L. inadai</em></td><td>17</td><td>17</td><td>14</td><td>28</td><td>17</td><td>5</td><td><strong>98</strong></td><td>21.98</td></tr><tr><td><strong>Saprophytic</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td><em>L. wolbachii</em></td><td>17</td><td>17</td><td>13</td><td>33</td><td>20</td><td>4</td><td><strong>104</strong></td><td>25.47</td></tr><tr><td><em>L. yanagawae</em></td><td>14</td><td>14</td><td>16</td><td>32</td><td>18</td><td>4</td><td><strong>98</strong></td><td>24.16</td></tr><tr><td><em>L. biflexa</em></td><td>14</td><td>14</td><td>14</td><td>34</td><td>19</td><td>4</td><td><strong>99</strong></td><td>25.05</td></tr><tr><td><em>L. vanthielii</em></td><td>15</td><td>15</td><td>12</td><td>33</td><td>19</td><td>4</td><td><strong>98</strong></td><td>23.15</td></tr><tr><td><em>L. terpstrae</em></td><td>16</td><td>16</td><td>15</td><td>34</td><td>22</td><td>4</td><td><strong>107</strong></td><td>26.14</td></tr><tr><td><em>L. meyeri</em></td><td>16</td><td>16</td><td>13</td><td>36</td><td>21</td><td>5</td><td><strong>107</strong></td><td>25.54</td></tr></tbody></table>

doi:10.1371/journal.pntd.0004403.t008

that are far apart in the genome. In addition, multistep phosphorelays may include intermediate histidine phosphotransferase proteins (Hpt), adding further complexity to TCS networks. Pathogenic, intermediate and saprophytic *Leptospira* species encode an unusually large and diverse set of TCSs, including orphan HKs and RRs, HHKs and Hpts in addition to classical paired HK/RR systems (Table 8 and S13 Table). Of note, more than 60% of the TCS genes found in *Leptospira* genomes encode non-classical orphan HK, orphan RR, HHK and Hpt proteins (Table 8). Overall, pathogenic species had the lowest average number of TCS genes (76), while saprophytic ones had the largest (102) (Table 8). Genome size-normalized TCS data revealed that pathogenic *Leptospira* species have roughly 35% less TCS genes in comparison to intermediate and saprophytic species (S6 Fig and Table 8). Additionally, pathogenic species had a proportionally lower number of strain-specific TCS genes compared to intermediate and saprophytic species (S7 Fig and S14 Table). We also identified a core set of 16 TCS genes shared among all the *Leptospira* genomes being analyzed (S7 Fig). Half of these TCS genes, conserved among all *Leptospira* species, were orphan HK/RRs (S15 Table). Taking into account their high conservation throughout the species, irrespective of saprophytic or pathogenic mode of lives, this core set of TCSs probably regulates pivotal cellular pathways in *Leptospira*.

Although there were 15 TCS genes conserved among pathogen and intermediate species, we did not identify genes that were shared between saprophyte:pathogen or saprophyte:intermediate species (S7 Fig). This finding is in agreement with the previous observation that

<footnote><sup>a</sup> Number of TCSs / Total genome size (Mb) respective to each *Leptospira* species.</footnote>

<footnote>Abbreviations: HK, histidine kinase; RR, response regulator; Hpt, histidine phosphotransferase; TCSs, two component systems; Mb, megabase.</footnote>
intermediate *Leptospira spp.* are more closely related to pathogens rather than to saprophytes, and that the gene order is more conserved in pathogenic and intermediate species [32]. Among the 15 TCS genes found in both pathogenic and intermediate species, 5 were orphan RRs, 4 were orphan HKs, 2 were HHKs, and only 1 was a classical HK:RR pair (S13 Table).

It is worth noting that *L. kmetyi* and *L. alstoni* were outliers with respect to being pathogenic species that harbored the largest number of species-specific TCS genes (27 and 17, respectively) and the largest overall number of TCSs within this species group (S7 Fig and S14 Table). In contrast, *L. alexanderi* did not contain any species-specific TCS genes (S7 Fig and S14 Table). The seven TCS genes present in all pathogenic species (S16 Table) may represent moieties involved in common signaling pathways and may play a role in host colonization and pathogenic mechanisms.

## Discussion

Here we used comparative whole genome analysis to answer the overall question, “what makes a bacterial genus pathogenic?” This analysis delineated the definitive phylogenetic relationship among 20 species of *Leptospira*, and demonstrated that infectious species and clades of *Leptospira* contain unique genes that are not found in non-infectious *Leptospira* (summarized in Table 9).

In a general sense, the comprehensiveness of this analysis is fundamentally important for understanding large-scale evolutionary mechanisms by which saprophytic bacteria acquire genes to enable infectiousness and pathogenicity. More specifically, considering how complex the genus *Leptospira* is—among the most complex genera of pathogenic bacteria—our analyses indicate that many genetic events over evolutionary time have given rise to pathogenic *Leptospira* of diverse biological properties. The genus *Leptospira* contains non-infectious environmental saprophytes and those members infectious to mammals. Infectious *Leptospira* are subdivided further into phylogenetically separated groups: pathogens (group I) and

<table><caption>**Table 9. Summary of Some Key Genomic Differences in *Leptospira* Species that Suggest Role in Evolution from Saprophyte to Infectious Pathogen.**</caption><thead><tr><td>Name of feature</td><td>General function (Table/Fig)</td><td>Comments</td></tr></thead><tbody><tr><td>Rfb locus</td><td>Serovar-specific polysaccharide biosynthesis on lipopolysaccharide (Fig 5; S2 Fig)</td><td>Rfb loci more complex in pathogenic Leptospira than in intermediates and saprophytes, possibly reflecting adaptation to mammalian environments and pathogen-host cell interactions</td></tr><tr><td>CRISPR/Cas systems</td><td>Defense against exogenous nucleic acids (Table 5; Fig 6)</td><td>Found only in pathogenic and intermediate Leptospira</td></tr><tr><td>Proteolytic enzymes (thermolysin, collagenase)</td><td>Degradation of host proteins (eg., complement components, interstitial proteins including collagen) (S9 Table)</td><td>Immune evasion, invasion of mammalian host tissues</td></tr><tr><td>Sialic acid biosynthetic genes</td><td>Predicted to modify proteins with neuraminic acid (Table 3; S3 Fig)</td><td>Complete loci found only in pathogenic Leptospira; sialic acid-modified proteins predicted to be involved in pathogen-host interactions</td></tr><tr><td>Catalase/peroxidase/superoxide dismutase</td><td>Detoxification of oxidative radicals (Table 6)</td><td>Catalase only found in pathogenic Leptospira; superoxide dismutase absent in pathogenic Leptospira</td></tr><tr><td>Extracellular matrix (ECM)-binding proteins</td><td>Adhesion (S8 Table)</td><td>Some members found only in pathogenic and intermediates; some found only in pathogenic; some found in all clades. Potentially intermediate Leptospira adhesion to host cells (in case of infectious Leptospira); may be involved movement of Leptospira in and through tissue interstitium; may generally promoting invasion/colonization processes</td></tr><tr><td>Cobalamin biosynthesis</td><td>Production of vitamin B12 (Table 2)</td><td>Infectious Leptospira predicted to be autotrophic for B12 biosynthesis, predicted to be important during in vivo mammalian infection</td></tr><tr><td>PF07598 family</td><td>Unknown (Table 7)</td><td>Upregulated in vivo; found only in pathogenic Leptospira</td></tr></tbody></table>

doi:10.1371/journal.pntd.0004403.t009
intermediate pathogens (group II). Previous systems (pathogenomic) approaches began to
yield insights into the large scale set of genes that enable infectious *Leptospira* to adhere, invade,
colonize, persist, evade the immune system and cause disease in mammalian reservoirs hosts as
well as accidental hosts [32, 149, 156]. By identifying novel gene families, differences in patho-
gen clade-related gene content, and key potential metabolic differences among infectious *Lep-
tospira* species as well as *Leptospira* species that have contrasting potential for causing severe
disease, for example, the present work contributes towards the next generation of leptospirosis
pathogenesis experimental studies. The data and analyses resulting from this *Leptospira*
Genome Project will contribute to new research directions in diagnostics, vaccine and thera-
apeutics development to prevent and ameliorate leptospirosis, with One Health relevance for
the health of humans and animals of veterinary importance alike.

This comparative analysis of the genus *Leptospira* assessed phylogenetic relationships
among species in several independent ways, including single locus, multilocus [51] and whole
genome approaches. All approaches robustly confirmed the separation of the 20 *Leptospira*
species into three clades: pathogens, intermediate pathogens and saprophytes; infectious *Lep-
tospira* include members of the pathogen and intermediate pathogen clades. Whole genome
analysis produced consistent dendrograms, similar but not identical to those obtained by
multi-locus sequence typing (MLST) [51]; MLST has been most useful for characterizing *Lep-
tospira* isolates [22, 51, 70, 157, 158] but also has been used to identify *Leptospira* strains
directly in clinical samples. Speciation based on Bayesian analysis of 16s rDNA (rrs) gene
sequences has become generally accepted [35] for *Leptospira*, especially differentiation of path-
ogen from non-pathogen (potential contaminant). However, the *Leptospira* 16S rRNA gene is
highly conserved so that species cannot be further subdivided. For example, 16S rDNA
sequence-deduced phylogeny could not distinguish *L. meyeri* from *L. yanagawae*. However,
the data provided here and elsewhere demonstrate that the gold standard for future definitive
taxonomical definition of any *Leptospira* isolate will be based on whole genome sequence-
based *in silico* DNA-DNA hybridization [11]. Genome-to-genome distances (GGD) analysis
confirmed *in silico* DNA-DNA hybridization results. As previously determined by the classical
*in vitro*-performed DNA-DNA hybridization studies, each of the representative strains repre-
sents a distinct species (estimated hybridization between pair of strains <70%) (S2 Table). For
example, *L. interrogans* serovar Icterohaemorrhagiae strain M20 is phylogenetically related to
*L. kirschneri* strain 3522C (estimated hybridization 42.30% ± 2.53) and *L. noguchii* strain
CZ214 (estimated hybridization 37.80% ± 2.49), while other pathogenic, intermediate and sap-
rophytic species are distantly related to *L. interrogans* (S2 Table). One serovar may belong to
more than one species (Fig 5) so that serovar does not have precise taxonomic implications,
although the present analysis is limited by not exploring within-species O-antigen loci in the
genomic data from the 320 isolates for which genomic sequence information was generated.
Finally, serovars belonging to the same species had GGD values higher than 70% DDH similar-
ity. *In silico* DDH values therefore accurately reflect whole genome relatedness and may be
used for the purpose of species delineation [159, 160], thus replacing the classical DNA-DNA
hybridization technique which, we argue, is now obsolete. The use of genome sequences also
provides reusable data and reproducible results. The GC content and a set of core genes
(including ribosomal genes) can also be extracted from genome sequences to verify that the
data are phylogenetically consistent.

## Lipopolysaccharide *rfb* biosynthetic loci in 20 *Leptospira* species

A remarkable feature of *Leptospira*—one that has often dominated the study of *Leptospira* over
the past century—is the serologically-determined variety of *Leptospira* serovars. Serovar
identification based on serological agglutination tests previously was the basis of *Leptospira* taxonomy, and continues to be important for understanding leptospirosis epidemiology related to human and animal infection and mammalian reservoir hosts. The present comparative genome analysis of 20 *Leptospira* species, as well as hundreds of additional *Leptospira* whole genome sequences not yet completely analyzed, demonstrates genetic determinants of the anti-genic carbohydrates (putatively O-antigens). Therefore molecular analysis of genes in the *rfb* O-antigen locus rather than serological tests will be the future basis of serovar identification.

Leptospiral LPS has long been known to be non-endotoxigenic, the basis for which is thought to be due to unusual modifications of the core Lipid A component of LPS resulting in altered Toll-like receptor-mediated innate immune responses [161-164]. Leptospiral serovar, determined by LPS antigenicity, seems to have some association with mammalian host prediction—for example, serovar Copenhageni with *Rattus* spp., serovar Canicola with dogs [1, 2, 165]—but a causal, mechanistic role here remains to be experimentally demonstrated. Indirect evidence from other microbial systems suggests the speculative hypothesis that the diversity of leptospiral LPS may be driven by environmental ecology-mediated selection pressures [166, 167], such as has been reported in *Salmonella* spp. which has many (>2000) serotypes). The genus *Leptospira* has been reported to contain more than 300 serovars [1, 2, 165], the basis for which remains essentially unexplored.

## Metabolic reconstructions

Here we report the first genomically-predicted metabolic network analysis [168, 169] of *Leptospira*, comparing members of the pathogen, intermediate pathogen and saprophyte clades. These large-scale reconstructions allow classification of the conserved metabolic capabilities (core metabolic network) and the unique metabolic capabilities (pan metabolic network). These reconstructions can be further converted to metabolic models of metabolism to probe metabolic capabilities computationally.

The most striking differences between infectious/pathogenic *Leptospira* and non-pathogens arose in porphyrin and vitamin biosynthetic capabilities. *L. interrogans* was shown to have a nearly complete vitamin B12 biosynthetic pathway that enables *de novo* B12 synthesis from an L-glutamate precursor, while *L. biflexa* completely lacked this pathway. Only pathogenic *Leptospira*—*L. interrogans* and *L. kmetyi*—were predicted to have a full folate biosynthetic pathway. These differences in biosynthetic capabilities may allow such pathogens to survive in nutrient-limited niches within the mammalian host. These observations are consistent with previous observations that found that *L. interrogans* serovar Canicola can grow *in vitro* in the absence of B12 but not B₁ [170], but in contrast with observations of others who concluded that pathogenic *Leptospira* could grow in the absence of B₁ but not B12 [171].

An open question in *Leptospira* biology is why *L. interrogans* grows more slowly than do intermediate pathogens and saprophytes, such as *L. licerasiae* and *L. biflexa*, which grow rapidly in defined EMJH medium [38]. The metabolic network model of *L. interrogans* was shown to lack L-glutamate oxidoreductase, an enzyme involved in recruiting ammonia as a nitrogen source [172, 173]), predicting a lower growth yield compared to the other *Leptospira* models in our *in-silico* minimal media analysis. The model of *L. biflexa* predicted the greatest yield with this reaction because this *Leptospira* contains the L-aspartate ammonia-lyase reaction allowing it to convert L-aspartate into fumarate and ammonia, in addition to using this component solely for biomass generation. These observations hint at one possible solution to the question of different growth rates, but model-guided experimentation is required to validate this prediction. Predictions made using these metabolic networks depend on an accurate reactome [174], and must be validated experimentally. If the models reported here are further curated and
experimentally validated they would be the first such mode for a pathogenic spirochete. Such an approach will yield fundamental insights into diverse metabolic capabilities of this diverse phylum, including *Borrelia* spp., which cause Lyme borreliosis and relapsing fever, and *Treponema* spp., which cause syphilis, yaws, periodontitis and other diseases.

## Vitamin autotrophy

Vitamin B12 (cobalamin) is the largest and most complex of natural organometallic cofactors and coenzymes, its *de novo* synthesis requiring ~30 energetically costly enzymatic steps [175]. Mammals have evolved highly complex, regulated mechanisms to absorb, transport and store cobalamin [176], reminiscent of the baroque processes involved in iron uptake, transport and storage in humans [177] [178]. We found that the genetic machinery for B12 autotrophy is found in infectious but not saprophytic *Leptospira*, leading us to speculate that such autotrophy allows *Leptospira* to infect mammals in the face of B12 sequestration by the host.

All *Leptospira* survive in the external environment, but our analysis predicts that only pathogenic strains make cobalamin *de novo* from L-glutamate, suggesting that this process is critical *in vivo*. Important in this context, cobalamin absorption and utilization in mammals is mediated by an elaborate set of carrier proteins, receptors and transporters [176] that generally are presumed necessary to process and protect this very large molecule. Considering the mechanistic details of cobalamin handling in mammals and the B12 autotrophy of infectious *Leptospira*, we hypothesize that mammalian B12 systems deprive invasive microbes of cobalamin, akin to the role of iron absorption, transport and sequestration known to sequester iron from pathogens, which have evolved siderophore mechanisms to acquire iron *in vivo* in mammalian hosts. Comparative analysis of *de novo* cobalamin biosynthesis in *Leptospira* predicts that infectious *Leptospira* are autotrophic for synthesizing this compound while saprophytes are auxotrophic, and suggests lines of experimentation to explore further the details cobalamin biosynthesis in *Leptospira*. The significance of the absence of complete cobalamin biosynthetic pathways in some group I *Leptospira* remains unclear.

## Detoxification of reactive oxygen species suggests resistance to host defense and differences in ecological niche

Previous comparative biochemical studies of spirochetes demonstrated catalase activity only in pathogenic *Leptospira* (all of which were previously classified inclusively as *L. interrogans*), and superoxide dismutase activity only in the saprophytic *L. biflexa*; peroxidase activity was present in both clades [179]. An important finding in our comparative genome analysis was that the *Leptospira* catalases, *KatA* and the putative catalase ortholog (LEP1GSC062_4039), were only found in pathogenic *Leptospira*, while the single leptospiral superoxide dismutase gene *sod* was found only in saprophytic *Leptospira*. Catalase has classically been associated with resistance to phagocyte-produced oxidative burst-mediated killing of pathogens (*viz*. hydrogen peroxide), typically intracellularly after phagocytosis. The presence of catalases only in pathogenic *Leptospira* suggests the testable hypothesis that this enzyme class may be involved in intracellular resistance to intracellular host cell killing, following on the published observations of pathogenic *Leptospira* within phagolysosomes [180, 181]; whether pathogenic *Leptospira* survive and proliferate in this subcellular compartment has not been conclusively demonstrated. Conversely, the absence of *sod* in pathogenic *Leptospira* suggests that this clade occupies an environmental niche not exposed to oxidative radicals. Presumably the observation that *sod* is restricted to saprophytes suggests that this clade is exposed to a different context in which oxidative radicals are found in the environment.
## Adhesion

Extracellular matrix (ECM)-binding proteins potentially mediate *Leptospira* adhesion to mammalian host cells, movement of *Leptospira* in and through tissue interstitium, and generally promoting invasion/colonization processes [37, 182–195]. A diverse array of ECM-binding proteins has been identified suggesting a redundancy of adhesion molecules that are probably part of invasion strategies of *Leptospira*. Indeed, many putative adhesins are multifunctional as they bind plasminogen and generate plasmin [185, 196, 197], increasing proteolytic processes associated with infection, or they could participate in immune evasion strategies by interacting with complement regulators (see below). Moreover, several of these proteins are recognized by human leptospirosis serum samples indicating their expression during infection. A caveat is that almost all these putative adhesins have been identified based on binding studies with recombinant proteins. Although site-directed mutagenesis of pathogenic *Leptospira* spp remains difficult, relatively straightforward methods have been developed for functional analysis of putative adhesins genes through gain-of-function studies in *L. biflexa* [92] [198] [199]. Many putative adhesins are present in saprophytic *Leptospira* but the presence of DNA sequences does not mean that the proteins are expressed. Although non-pathogenic *Leptospira* species may encode the genes for putative adhesins, transcripts of *ompL1* and the protein itself, for example, have not been detected in *L. biflexa* serovar Patoc (190). In any event, adhesion to environmental biotic or abiotic structures may be part of the biology of saprophytic *Leptospira* but such a concept remains speculative at this time [200].

Many adhesins are multipurpose proteins as they bind PLG and produce PLA. The generation of fully active PLA, aside from its other functions, also contributes to the enhanced degradation of complement components. Indeed, it has been shown that in *Leptospira*, PLA decreases C3b and human IgG deposition, most probably through their degradation, thereby hampering opsonization, restricting complement antibacterial functions [201]. Another mechanism of complement evasion is through the acquisition of host regulators of complement activation. Surface microbial proteins that bind to complement inhibitors and activate them permit pathogens to inhibit the complement response on the bacteria [202]. Binding of *Leptospira* to factor H (FH), factor H-like protein (FHL-1) and C4 binding protein (C4BP) has been reported [203, 204] and several complement regulators-binding proteins have been identified [205–211]. Lsa23 is an example of multifunctional protein capable of binding ECM, PLG/PLA and complement regulators and should play a role in leptospiral virulence [211].

## Immune evasion via proteolysis of complement

The protein family includes several metalloprotease members that are considered virulence factors in several pathogens [212]. Aureolysin, a zinc-dependent metalloprotease of *S. aureus*, acts in synergy with host regulators to inactivate C3 thus potentially inactivating host immune response [213]. Indeed, in the case of *Leptospira*, thermolysins were only found in pathogenic strains. Recently, experimental evidence of the sequence LIC13322 encoding a putative thermolysin on the direct degradation of complement factors has been reported, suggesting its role in immune evasion by pathogenic leptospiral strains [214]. Thus, it seems that pathogenic *Leptospira*, like other successful pathogens, utilize at least two strategies to circumvent the complement system: acquisition of host complement inhibitors and degradation of complement components, either thorough PLG/PLA generation or by the presence of bacterial proteases.

## Sialic acids

We show here a distinctive presence of the sialic acid cluster in most pathogenic *Leptospira* species and notably absent from intermediately pathogenic and saprophytic *Leptospira* species.
This fact suggests a role for sialic acids as virulence determinants. The data shown are consis-
tent with previously published data in which *L. interrogans* and *L. alexanderi* produced di-acet-
ylated nonulosonic acids and *L. santarosai* did not produce identifiable nonulosonic acid
species, and this observation could be associated with some species missing particular enzymes
of the pathway [215].

The pathogens that contain the whole sialic acid biosynthesis cluster have N-acetylneurami-
nic acid synthetases predicted to produce legionaminic acids, while other species have a cluster
predicted to produce pseudaminic acid. Pseudaminic acid has been shown to be required for
flagella biogenesis in *Campylobacter* spp. and *Helicobacter* spp. and function as virulence fac-
tors [216]. Leptospiral flagella are located between the inner membrane (IM) and outer mem-
brane (OM) and drive motility [217]. The presence of a pseudaminic acid pathway in all
species suggests glycosylation could have a similar role in flagella biogenesis.

Legionaminic acid has been associated with virulence in *Legionella pneumophila* and *Cam-
pylobacter coli* where its cell surface location seems to be involved in adhesion, cell-cell interac-
tion and immune evasion [218]. We hypothesize that similar mechanisms occur in pathogenic
leptospires.

PF07598 paralogous gene family

A novel gene family—first identified in a pathogenomic screen of *L. interrogans* serovar Lai
[219]—was found to have orthologs present in pathogenic *Leptospira* but not in intermediate
or saprophytic *Leptospira*. In *L. interrogans* serovar Lai, these genes were reported to be upregu-
lated *in vivo* [219] and here we report that the numbers of paralogs varies among the patho-
gens, with the *L. interrogans*, *L. kirschneri* and *L. noguchi* having the most. These observations
suggest that the PF07598 genes contribute to leptospiral virulence but the mechanism(s) by
which they do so remain to be elucidated; no functional annotations for the PF07598 gene fam-
ily is yet possible. Experimental studies of this gene family will likely provide insight into lepto-
spirosis pathogenesis.

Motility and chemotaxis

Although motility is essential for pathogenesis of *Leptospira* [150, 151], all *Leptospira* spiro-
chetes, including those belonging to intermediate and saprophytic species, are motile. Consis-
tent with this observation, we found that pathogenic, intermediate and saprophytic species of
*Leptospira* have all the genes necessary to assemble a functional flagellar apparatus. Further-
more, the flagella-encoding genes are highly conserved within the genus, indicating that patho-
genic and non-pathogenic *Leptospira* do not differ significantly with respect to their flagella
apparatus and structure.

In contrast to what we found with regard to motility, there was high diversity with respect to predicted amino acid sequence identity among genes encoding chemotaxis proteins from different species of *Leptospira*. Furthermore, not all chemotaxis proteins are present in all the species, which corresponds to different chemotactic behaviors observed in pathogenic and saprophytic *Leptospira* [220]. The majority of the diversity among the chemotaxis proteins was observed in MCPs, trans-membrane sensor proteins that trigger the intracellular signal transduction in bacterial chemotaxis [221], and located at the cell poles near the basal body and flagellar motor as other chemotaxis proteins [217, 222]. This finding, together with the observation that there is a higher diversity among genes associated with the basal body and that asymmetrical periplasmic flagellar rotation occurs with the interaction of basal body and chemotaxis proteins, suggests that the sensing and chemotactic response regulated by this proteins may impact their survival in specific environments, including their ability to infect a mammalian host.
## Gene regulation

Transcriptional regulation is the most common form of regulation in bacteria, often via a specific transcriptional regulator (activator or repressor) or an alternative sigma factor. In this study, we found that *Leptospira* have only three types of alternative sigma ($\sigma$) factors ($\sigma^5^4$, $\sigma^F$, $\sigma^E$), which is fewer than *E. coli* which have genome sizes similar to those of *Leptospira*. Our analysis revealed that there is a distinct difference in the $\sigma^5^4$ regulatory networks between pathogenic and saprophytic *Leptospira* species. All pathogenic species have two $\sigma^5^4$ regulatory networks, LepA-$\sigma^5^4$ and LepB-$\sigma^5^4$, while saprophytic *Leptospira* only have LepA-$\sigma^5^4$. Historically, $\sigma^5^4$ is known to be involved in nitrogen assimilation, and now is well recognized to regulate diverse functions in response to various stimuli [223]. In addition, $\sigma^5^4$ has also been shown to be essential for infection in some pathogenic bacteria. For example, in another spirochetal pathogen, *B. burgdorferi*, the EBP activator, Rrp2, and $\sigma^4$, controls production of RpoS which in turn, governs expression of many virulence factors important for mammalian infection such as OspC [224, 225]. Thus, it is logical to postulate that LepB-$\sigma^5^4$, is involved in survival in natural environments for both pathogenic and saprophytic *Leptospira*, while LepA-$\sigma^5^4$ plays a role in survival in the host for pathogenic species. This hypothesis merits experimental testing.

Extracytoplasmic function $\sigma$ factors (ECF $\sigma$ factors) are the most diverse alternative $\sigma$ factors found in many bacteria [226, 227]. Many bacteria contain multiple ECF $\sigma$ factors. For example, *Pseudomonas aeruginosa* has more than 19 ECF $\sigma$ factors [226]. Based on sequence analysis, ECF $\sigma$ factors have been grouped into over 40 classes [226]. Our analyses showed that *Leptospira* have 5–10 ECF $\sigma$ factors, and pathogenic *Leptospira* have 5 more ECF $\sigma$ factors than saprophytic species (S11 Table), which is consistent with the more complex life cycle of pathogenic species than of saprophytic species. All *Leptospira* spp. have one copy of ECF31 and ECF43 with unknown functions. Pathogenic *Leptospira* have additional 5 unclassified ECF $\sigma$ factors. Saprophytic *Leptospira* have one copy of ECF41 and ECF42 that are not found in pathogenic species. Although functions of ECF41 and ECF42-type remain unknown, one report showed that one of the ECF41 $\sigma$ factors, *SigJ* in *Mycobacterium tuberculosis*, is involved in resistance to hydrogen peroxide [228]. It is unclear whether ECF41 in saprophytic *Leptospira* (LEPBI_I1070) has a similar function as *SigJ*, and if so, how it contributes to the survival of saprophytic *Leptospira* in the environment.

The activity of ECF $\sigma$ factors is often regulated by an anti-$\sigma$ factor, a transmembrane protein that binds and inhibits the activity of ECF $\sigma$ [229]. Cleavage of anti-$\sigma$ factor by proteases leads to release and activate $\sigma^E$. Extracellular signals regulate this intramembrane proteolysis often via an anti-anti-$\sigma$ factor (or called anti-$\sigma$ antagonist). Our analyses revealed that both pathogenic and saprophytic *Leptospira* have more than 30 $\sigma^E$ regulators. Among them, 17 are only found in pathogenic/intermediate *Leptospira*, while 19 are found solely in saprophytic *Leptospira* (S11 Table). These differences likely reflect the variety of signals sensed by pathogenic and saprophytic *Leptospira*.

In addition to alternative $\sigma$ factors, both pathogenic and saprophytic *Leptospira* species have many putative transcriptional regulators, far more than is found in other pathogenic spirochetes such as *Borrelia burgdorferi* and *Treponema pallidum* [230, 231]. Our initial analyses of transcriptional regulators among *Leptospira* species did not yield a distinct pattern of correlation with pathogenicity. Further *in silico* and experimental analyses to confirm the prediction and more importantly, to determine their regulatory role in *Leptospira*, is needed.

## Sensory transduction

*Leptospira* species have a high number of two-component sensory systems (TCSs) (70–100) compared to *Borrelia*, *Treponema* and *Bradyspira* (6–20). The number of TCS genes found in a
particular *Leptospira* species strongly correlated with the likely diversity of ecological niches that the species encounters, a phenomenon observed in other bacteria [203]. The lower number of TCSs found in pathogenic species may be linked to the process of host adaptation, whereas the larger numbers of unique TCSs in intermediates and, even more in saprophytes, may be instru- mental for sensing and adapting to a more diverse range of environmental conditions. Regard- less, almost all of the pathogenic species encode more than 70 TCS genes, indicating that the *Leptospira* pathogen requires a highly complex network of signaling processes for its life cycle. Interestingly, the proportion of TCS genes encoding orphan HK/RR and HHK proteins are higher in *Leptospira* (> 60%) than in other bacteria where orphan TCS proteins are unusual. These findings suggest that branched signaling pathways may be relevant in this genus and could confer added physiological advantages to *Leptospira* under specific circumstances.

One limitation in the present analytical approach is the difficulty in defining a robust and confident automatic method to segregate orthologous clusters among all *Leptospira* strains, and especially those that relate to TCS function. As in most *in silico* analyses, further biochemical experiments are needed to confirm the role of the various TCS categories identified in this study. Another limitation of this cross-species comparative analysis is that differences in serovars/ strains with *Leptospira* species were not studied; such analyses will be a future priority given the strength of the approach and the depth of existing data. This *in silico* approach will also not iden- tify novel virulence factors nor mechanisms of pathogenesis based on sequence analysis alone.

In summary, the large-scale comparative genomic analysis of 20 *Leptospira* species has pro- vided broad insights into how infectious members of this genus acquired the genes necessary to acquire pathogenicity and virulence, placing these species within a definitive phylogeny. Novel, *Leptospira* species-specific genes and gene families were identified. Genomically-based metabolic reconstruction predictions predict novel adaptation of infectious *Leptospira* to mam- mals (summarized in Table 9), including sialic acid biosynthesis, pathogen-specific porphyrin metabolism and the first-time demonstration of riboswitch-regulated cobalamin (B12) autot- rophy as a bacterial virulence factor. Only pathogenic *Leptospira* contain CRISPR/Cas systems, suggesting not only a potential mechanism for this clade’s refractoriness to gene targeting but also possible novel means to be able to genetically modify pathogenic *Leptospira*. Whether restriction modification systems might contribute to gene targeting has yet to be analyzed in detail, but the publicly available whole genome data sets provided in support of the present work will be contribute to carrying out such analyses. A novel virulence-related genes/gene family epitomized by the PF07598 group of paralogs suggests adaptation and diversification of this protein family within the pathogenic clade. Identifying large scale changes in infectious (pathogenic and intermediate pathogenic) as compared to non-infectious *Leptospira* has yield large-scale, novel insights into the evolution of a bacterial pathogen, provides the basis for new directions in leptospirosis pathogenesis research. It also makes novel genomic and pathoge- nomic contributions to the field of bacterial pathogenesis, which is of general interest.

**Supporting Information**

S1 Fig. **Pan-genome, core and novel genes of the 20 sequenced *Leptospira* species.** The blue and red lines denote the pan-genome and core genes as genomes are added in the order noted on along the x-axis (A). The bars indicate the number of novel gene families discovered for each genome added. The color of the bars illustrate the three main groupings of *Leptospira*: pathogenic (red), intermediate (blue), and saprophytic (green). The number of novel genes discovered with the addition of each new genome (B) was estimated using a pan-genome model based on the original model presented by Tetellin et al. [63]. Purple circles are the median of each distribution (grey circles). Power law (red lines) and exponential (blue lines) regressions
were plotted to determine (α), and tg(θ), respectively. The exponent (α) indicates whether the pan-genome is open (α ≤ 1) or closed (α > 1) [83] and tg(θ) denotes the average extrapolated number of strain-specific/novel genes.

(PDF)

**S2 Fig. Flanking Genes Surrounding the *Leptospira rfb* locus gene clusters.** The *rfb* region and flanking CDSs (blue) 9 of pathogenic (A), 5 intermediate (B), and 6 saprophytic (C) representative *Leptospira* species were compared. *rfb* region CDSs are labeled by locus identifier and colored by functional role categories as noted in the boxed key. Gene symbols, when present, are noted above their respective genes. BLASTP matches between CDSs are colored by protein percent identity (see key).

(PDF)

**S3 Fig. Phylogenetic analysis of leptospiral N-acetylneuraminic (Sialic) Acid Synthetase (NeuB) protein sequences.** Maximum-likelihood tree shows pathogens (red lines), intermediates (green lines) and saprophytes (blue lines). Numbers denote node support. A red box highlights those proteins that are part of a complete sialic acid cluster.

(PDF)

**S4 Fig. Heat map of ORFs encoding 51 motility and 25 chemotaxis proteins identified in analysis of the 20 *Leptospira* genomes.** ORFs are identified according to their *L. interrogans* serovar Copenhageni strain Fiocruz L1-130 number. The heat map shows the degree of amino acid sequence identity of ORFs with their respective orthologs in the *L. interrogans* strain Fiocruz L1-130 genome.

(PDF)

**S5 Fig. Comparison of ECF Sigma (σ) Factors Among *Leptospira*.** Venn diagram showing distribution of ECF σfactors unique or shared among the pathogenic (*L. interrogans* L1-130), intermediate pathogenic (*L. kmetyi*) and saprophytic (*L. biflexa*) species. The number and locus ID of ECF σ factors that are unique or shared among these *Leptospira* species are labeled in each sector of the diagram.

(PDF)

**S6 Fig. Normalized number of Leptospiral two component systems by genome size.** The number of TCS genes was normalized per Mbp genome (y-axis) of representative Leptospiral species (x-axis). See key for shading of pathogenic, intermediate and saprophyte genomes.

(PDF)

**S7 Fig. Venn diagram showing the distribution of TCS genes among *Leptospira* species.** The ratios depicted inside each one of the major groupings, correspond to the number of TCS ortholog genes present in the [majority:all-but-one:all] species of that particular group. True cut-off values for these Figs correspond to the presence of the gene in 50% (majority), 90% (all but one) or 100% (all) of the particular group of species. Sequence clusters that do not match the indicated cut-off value or those from unexpected groupings are included in the “ambiguous grouping” set. Singleton clusters, representing species-specific genes are noted in circles surrounding the Venn diagram.

(PDF)

**S1 Table. Universal protein markers.**

(PDF)

**S2 Table. Estimates of genome relatedness of *Leptospira* species.**

(XLSX)
**S3 Table.** Metabolites and reactions including exchange and biomass reactions used for metabolic reconstructions.

(XLSX)

**S4 Table.** Candidate substrates for the leptospiral TAT protein secretion system.

(XLSX)

**S5 Table.** Examination of the -1 position of leptospiral lipobox sequences.

(XLSX)

**S6 Table.** Identification of Lipid A biosynthesis proteins in *Leptospira*.

(XLSX)

**S7 Table.** A. Identity matrix comparison of lipid A biosynthesis pathway genes across the genus *Leptospira*. B. Homology comparison of amino acid sequences of enzymes involved in the synthesis of cell wall lipid A from *Leptospira* species.

(XLSX)

**S8 Table.** Leptospiral proteins involved in adhesion to extracellular matrix, plasminogen binding and complement evasion.

(XLSX)

**S9 Table.** Proteases with a potential role in host-pathogen interactions.

(XLSX)

**S10 Table.** A. Amino acid sequence identity comparison of proteins orthologous to *Leptospira interrogans* serovar Copenhageni immunodominant proteins. B. Amino acid identity of leptospiral proteins to Lig proteins and domains.

(XLSX)

**S11 Table.** Comparison of leptospiral proteins involved in motility and chemotaxis.

(XLSX)

**S12 Table.** Sigma factors and accessory proteins involved in gene regulation in *Leptospira* species.

(XLSX)

**S13 Table.** Inventory of two component systems in *Leptospira* species.

(XLSX)

**S14 Table.** Identification of strain-specific two component systems in *Leptospira* species.

(XLSX)

**S15 Table.** Core three component system proteins conserved among *Leptospira* species.

(XLSX)

**S16 Table.** Two component systems specific to pathogenic *Leptospira* species.

(XLSX)

## Acknowledgments

We are grateful to the J. Craig Venter Institute sequencing, bioinformatics and IT departments for supporting the infrastructure required to determine the genome sequences, annotation and pan-genome and other analyses carried out in this project.
Author Contributions

Conceived and designed the experiments: DEF MAM JMV. Performed the experiments: DAF MAM. Analyzed the data: DEF MAM HA BA LAS DEB DB AB YFC RLG DAH DHH RH AIK PNL JM AEM JMM ALTN KEN BP SJP MP JNR JT EAW XFY JJZ JMV. Contributed reagents/materials/analysis tools: DEF MAM HA BA LAS DEB DB AB YFC RLG DAH DHH RH AIK PNL JM AEM JMM ALTN KEN BP SJP MP JNR JT EAW XFY JJZ JMV. Wrote the paper: DEF MAM HA BA LAS DEB DB AB YFC RLG DAH DHH RH AIK PNL JM AEM JMM ALTN KEN BP SJP MP JNR JT EAW XFY JJZ JMV.

References

1. Levett PN. Leptospirosis. Clin Microbiol Rev. 2001; 14(2):296-326. PMID: 11292640

2. Bharti AR, Nally JE, Ricaldi JN, Matthias MA, Diaz MM, Lovett MA, et al. Leptospirosis: A zoonotic disease of global importance. Lancet Infect Dis. 2003; 3:757-71. PMID: 14652202

3. Ashford DA, Kaiser RM, Spiegel RA, Perkins BA, Weyant RS, Bragg SL, et al. Asymptomatic infection and risk factors for leptospirosis in Nicaragua. Am J Trop Med Hyg. 2000; 63(5-6):249-54. PMID: 11421372.

4. Ganoza CA, Matthias MA, Saito M, Cespedes M, Gotuzzo E, Vinetz JM. Asymptomatic renal colonization of humans in the peruvian Amazon by Leptospira. PLoS Negl Trop Dis. 2010; 4(2):e612. Epub 2010/02/27. doi: 10.1371/journal.pntd.0000612 PMID: 20186328; PubMed Central PMCID: PMC2826405.

5. Dikken H, Kmety E. Serological typing methods of leptospires. In: B T., NJ R., editors. Methods in Microbiology. 11. London: Academic Press; 1978. p. 259-307.

6. Brenner DJ, Kaufmann AF, Sulzer KR, Steigerwalt AG, Rogers FC, Weyant RS. 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. 1999; 49 Pt 2:839-58. Epub 1999/05/13. PMID: 10319510.

7. Brenner DJ, Kaufmann AF, Sulzer KR, Steigerwalt AG, Rogers FC, Weyant RS. 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. Int J Syst Bacteriol. 1999; 49:839-58. PMID: 10319510

8. Ramadass P, Jarvis BDW, Corner RJ, Penny D, Marshall RB. Genetic characterization of pathogenic Leptospira species by DNA hybridization. Int J Syst Bacteriol. 1992; 42:215-9. PMID: 1581182

9. Ramadass P, Jarvis BDW, Corner RJ, Cinco M, Marshall RB. DNA relatedness among strains of Leptospira biflexa. Int J Syst Bacteriol 1990; 40:231-5. PMID: 2397191

10. Yasuda PH, Steigerwalt AG, Sulzer KR, Kaufmann AF, Rogers FC, Brenner DJ. Deoxyribonucleic acid relatedness between serogroups and serovars in the family *Leptospiraceae* with proposals for seven new *Leptospira* species. Int J Syst Bacteriol. 1987; 37:407-15.

11. Bourhy P, Collet L, Brisse S, Picardeau M. Leptospira mayottensis sp. nov., a pathogenic species of the genus Leptospira isolated from humans. Int J Syst Evol Microbiol. 2014; 64(Pt 12):4061-7. doi: 10.1099/ijs.0.066597-0 PMID: 25249563.

12. Saito M, Villanueva SY, Kawamura Y, Iida K, Tomida J, Kanemaru T, et al. Leptospira idonii sp. nov., isolated from environmental water. International journal of systematic and evolutionary microbiology. 2013; 63(Pt 7):2457-62. doi: 10.1099/ijs.0.047233-0 PMID: 23203626.

13. Ko AI, Goarant C, Picardeau M. Leptospira: the dawn of the molecular genetics era for an emerging zoonotic pathogen. Nat Rev Microbiol. 2009; 7(10):736-47. Epub 2009/09/17. nrmicro2208 [pii] doi: 10.1038/nrmicro2208 PMID: 19756012.

14. Ahmed A, Grobusch MP, Klatser PR, Hartskeerl RA. Molecular approaches in the detection and characterization of Leptospira. J Bacteriol Parasitol. 2011;S5-002.

15. Cerqueira GM, Picardeau M. A century of Leptospira strain typing. Infect Genet Evol. 2009; 9:760-8. doi: 10.1016/j.meegid.2009.06.009 PMID: 19540362

16. Levett PN, Smythe L. International Committee on Systematics of Prokaryotes Subcommittee on the taxonomy of Leptospiraceae. Minutes of the closed meeting, 9 October 2013, Fukuoka, Japan. Int J Syst Evol Micr. 2014; 64:in press.

17. Auch AF, vonJan M, Klenk HP G M. Digital DNA-DNA hybridization for microbial species delineation by means of genome-to-genome sequence comparison. Stand Genomic Sci. 2010; 2:117-34. doi: 10.4056/sigs.531120 PMID: 21304684
18. Goris J, Konstantinidis KT, Klappenbach JA, Coenye T, Vandamme P, Tiedje JM. DNA-DNA hybridization values and their relationship to whole-genome sequence similarities. Int J Syst Evol Micr. 2007; 57:81–91.

19. Wolf YI, Rogozin I, Grishin N, Tatusov R, Koonin E. Genome trees constructed using five different approaches suggest new major bacterial clades. BMC Evol Biol. 2001; 1:8. PMID: 11734060

20. Auch AF, Klenk HP, Göker MS. Standard operating procedure for calculating genome-to-genome distances based on high-scoring segment pairs. Stand Genomic Sci. 2010; 2:142–8. doi: 10.4056/sigs.541628 PMID: 21304686

21. Boonsilp S, Thaipadungpanit J, Amornchai P, Wuthiekanun V, Bailey MS, Holden MT, et al. A single multilocus sequence typing (MLST) scheme for seven pathogenic Leptospira species. PLoS neglected tropical diseases. 2013; 7(1):e1954. Epub 2013/01/30. doi: 10.1371/journal.pntd.0001954 PMID: 23359622; PubMed Central PMCID: PMC3554523.

22. Ahmed A, Ferreira AS, Hartskeerl RA. Multilocus sequence typing (MLST): markers for the traceability of pathogenic Leptospira strains. Methods in molecular biology. 2015; 1247:349–59. doi: 10.1007/978-1-4939-2004-4_25 PMID: 25399108.

23. Nalam K, Ahmed A, Devi SM, Francalacci P, Baig M, Sechi LA, et al. Genetic affinities within a large global collection of pathogenic Leptospira: implications for strain identification and molecular epidemiology. PLoS One. 2010; 5(8):e12637. doi: 10.1371/journal.pone.0012637 PMID: 20805987; PubMed Central PMCID: PMC2929200.

24. Romero EC, Blanco RM, Galloway RL. Analysis of multilocus sequence typing for identification of Leptospira isolates in Brazil. J Clin Microbiol. 2011; 49(11):3940–2. doi: 10.1128/JCM.01119-11 PMID: 21880969; PubMed Central PMCID: PMC3209071.

25. Ahmed A, Thaipadungpanit J, Boonsilp S, Wuthiekanun V, Nalam K, Spratt BG, et al. Comparison of two multilocus sequence based genotyping schemes for Leptospira species. PLoS Negl Trop Dis. 2011; 5(11):e1374. doi: 10.1371/journal.pntd.0001374 PMID: 22087342; PubMed Central PMCID: PMC3210738.

26. Agampodi SB, Moreno AC, Vinetz JM, Matthias MA. Utility and limitations of direct multi-locus sequence typing on qPCR-positive blood to determine infecting Leptospira strain. The American journal of tropical medicine and hygiene. 2013; 88(1):184–5. doi: 10.4269/ajtmh.2012.12-0526 PMID: 23208890; PubMed Central PMCID: PMCPMC3541733.

27. Chiani Y, Jacob P, Varni V, Landolt N, Schmeling MF, Pujato N, et al. Isolation and clinical sample typing of human leptospirosis cases in Argentina. Infect Genet Evol. 2016; 37:245–51. doi: 10.1016/j. meegid.2015.11.033 PMID: 26658064.

28. Lehmann JS, Matthias MA, Vinetz JM, Fouts DE. Leptospiral pathogenomics. Pathogens. 2014; 3 (2):280–308. doi: 10.3390/pathogens3020280 PMID: 25437801; PubMed Central PMCID: PMC4243447.

29. Ren SX, Fu G, Jiang XG, Zeng R, Miao YG, Xu H, et al. Unique physiological and pathogenic features of Leptospira interrogans revealed by whole-genome sequencing. Nature. 2003; 422(6934):888–93. doi: 10.1038/nature01597 PMID: 12712204.

30. Nascimento AL, Ko AI, Martins EA, Monteiro-Vitarello CB, Ho PL, Haake DA, et al. Comparative genomics of two *Leptospira interrogans* serovars reveals novel insights into physiology and pathogenesis. J Bacteriol. 2004; 186(7):2164–72. Epub 2004/03/19. PMID: 15028702; PubMed Central PMCID: PMC374407.

31. Bulach DM, Zuerner RL, Wilson P, Seemann T, McGrath A, Cullen PA, et al. Genome reduction in *Leptospira borgpetersenii* reflects limited transmission potential. Proc Natl Acad Sci U S A. 2006; 103 (39):14560–5. Epub 2006/09/16. doi: 10.1073/pnas.0603979103 PMID: 16973745; PubMed Central PMCID: PMC1599999.

32. Ricaldi JN, Fouts DE, Selengut JD, Harkins DM, Patra KP, Moreno A, et al. Whole Genome Analysis of *Leptospira licerasiae* Provides Insight into Leptospiral Evolution and Pathogenicity. PLoS Negl Trop Dis. 2012; 6(10):e1853. doi: 10.1371/journal.pntd.0001853 PMID: 23145189; PubMed Central PMCID: PMC3493377.

33. Picardeau M, Bulach DM, Bouchier C, Zuerner RL, Zidane N, Wilson PJ, 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):e1607. Epub 2008/02/14. doi: 10.1371/journal.pone.0001607 PMID: 18270594; PubMed Central PMCID: PMC2229662.

34. Chou LF, Chen YT, Lu CW, Ko YC, Tang CY, Pan MJ, et al. Sequence of Leptospira santarosai serovar Shermani genome and prediction of virulence-associated genes. Gene. 2012; 511(2):364–70. doi: 10.1016/j.gene.2012.09.074 PMID: 23041083.

35. Matthias MA, Diaz MM, Campos KJ, Calderon M, Willig MR, Pacheco V, et al. Diversity of bat-associated Leptospira in the Peruvian Amazon inferred by bayesian phylogenetic analysis of 16S ribosomal
DNA sequences. Am J Trop Med Hyg. 2005; 73(5):964-74. PMID: 16282313; PubMed Central PMCID: PMC2270400.

36. Levett PN. Systematics of leptospiraceae. Current topics in microbiology and immunology. 2015; 387:11-20. doi: 10.1007/978-3-662-45059-8_2 PMID: 25388130.

37. Nascimento AL, Ko AI, Martins EA, Monteiro-Vitorello CB, Ho PL, Haake DA, et al. Comparative genomics of two Leptospira interrogans serovars reveals novel insights into physiology and pathogenesis. J Bacteriol. 2004; 186(7):2164-72. Epub 2004/03/19. PMID: 15028702; PubMed Central PMCID: PMC374407.

38. Matthias MA, Ricaldi JN, Cespedes M, Diaz MM, Galloway RL, Saito M, et al. Human leptospirosis caused by a new, antigenically unique Leptospira associated with a Rattus species reservoir in the Peruvian Amazon. PLoS Negl Trop Dis. 2008; 2(4):e213. Epub 2008/04/03. doi: 10.1371/journal.pntd.0000213 PMID: 18382606; PubMed Central PMCID: PMC2271056.

39. Delcher AL, Phillippy A, Carlton J, Salzberg SL. Fast algorithms for large-scale genome alignment and comparison. Nucleic acids research. 2002; 30(11):2478-83. Epub 2002/05/30. PMID: 12034836; PubMed Central PMCID: PMC117189.

40. Fouts DE, Mongodin EF, Mandrell RE, Miller WG, Rasko DA, Ravel J, et al. Major Structural Differences and Novel Potential Virulence Mechanisms from the Genomes of Multiple *Campylobacter* Species. PLoS Biol. 2005; 3(1):e15. PMID: 15660156; PubMed Central PMCID: PMC539331.

41. Fouts DE, Tyler HL, DeBoy RT, Daugherty S, Ren Q, Badger JH, et al. Complete Genome Sequence of the N2-Fixing Broad Host Range Endophyte *Klebsiella pneumoniae* 342 and Virulence Predictions Verified in Mice. PLoS Genet. 2008; 4(7):e1000141. PMCID: PMC2453333. PMID: 18654632; PubMed Central PMCID: PMC2453333. doi: 10.1371/journal.pgen.1000141

42. Davidsen T, Beck E, Ganapathy A, Montgomery R, Zafar N, Yang Q, et al. The comprehensive microbial resource. Nucleic acids research. 2010; 38 part 1:1–10. Epub 2009/11/07. doi: 10.1093/nar/gkp912 PMID: 19892825; PubMed Central PMCID: PMC2808947.

43. Chen Y, Stine OC, Badger JH, Gil AI, Nair GB, Nishibuchi M, et al. Comparative Genomic Analysis of *Vibrio parahaemolyticus*: Serotype Conversion and Virulence. BMC Genomics. 2011; 12:294. Epub 2011/06/08. doi: 10.1186/1471-2164-12-294 PMID: 21645368; PubMed Central PMCID: PMC3130711.

44. Finn RD, Clements J, Eddy SR. HMMER web server: interactive sequence similarity searching. Nucleic acids research. 2011; 39(Web Server issue):W29-37. Epub 2011/05/20. doi: 10.1093/nar/gkr367 PMID: 21593126; PubMed Central PMCID: PMC3125773.

45. Cole JR, Wang Q, Cardenas E, Fish J, Chai B, Farris RJ, et al. The Ribosomal Database Project: improved alignments and new tools for rRNA analysis. Nucleic acids research. 2009; 37 part 1:1-10. Epub 2008/11/14. doi: 10.1093/nar/gkn879 PMID: 19004872; PubMed Central PMCID: PMC2686447.

46. Sonnhammer EL, Hollich V. Scoredist: a simple and robust protein sequence distance estimator. BMC Bioinformatics. 2005; 6:108. Epub 2005/04/29. doi: 10.1186/1471-2105-6-108 PMID: 15857510; PubMed Central PMCID: PMC1131889.

47. phylipFasta—Wrapper for the Phylip Package Written in Ruby. Available from: https://github.com/jhbadger/phyloFasta.

48. Felsenstein J. PHYLIP—Phylogeny Inference Package (Version 3.2). Cladistics. 1989; 5:164-6.

49. Felsenstein J. PHYLIP (Phylogeny Inference Package) version 3.69. Distributed by the author, Department of Genome Sciences, University of Washington, Seattle. 2009.

50. Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011; 7:539. Epub 2011/10/13. doi: 10.1038/msb.2011.75 PMID: 21988835; PubMed Central PMCID: PMC3261699.

51. Boonsilp S, Thaipadungpanit J, Amornchai P, Wuthiekanun V, Bailey MS, Holden MT, et al. A single multilocus sequence typing (MLST) scheme for seven pathogenic *Leptospira* species. PLoS Negl Trop Dis. 2013; 7(1):e1954. Epub 2013/01/30. doi: 10.1371/journal.pntd.0001954 PMID: 23359622; PubMed Central PMCID: PMC3554523.

52. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Molecular biology and evolution. 2011; 28(10):2731-9. Epub 2011/05/07. doi: 10.1093/molbev/msr121 PMID: 21546353; PubMed Central PMCID: PMC3203626.

53. Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Systematic biology. 2003; 52(5):696-704. Epub 2003/10/08. PMID: 14530136.
54. Brown JR, Douady CJ, Italia MJ, Marshall WE, Stanhope MJ. Universal trees based on large combined protein sequence data sets. Nat Genet. 2001; 28(3):281-5. Epub 2001/06/30. doi: 10.1038/90129 PMID: 11431701.

55. Santos SR, Ochman H. Identification and phylogenetic sorting of bacterial lineages with universally conserved genes and proteins. Environ Microbiol. 2004; 6(7):754-9. Epub 2004/06/10. doi: 10.1111/j.1462-2920.2004.00617.x PMID: 15186354.

56. Ciccarelli FD, Doerks T, von Mering C, Creevey CJ, Snel B, Bork P. Toward automatic reconstruction of a highly resolved tree of life. Science. 2006; 311(5765):1283-7. Epub 2006/03/04. doi: 10.1126/science.1123061 PMID: 16513982.

57. Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic acids research. 1994; 22(22):4673-80. Epub 1994/11/11. PMID: 7984417; PubMed Central PMCID: PMC308517.

58. Capella-Gutierrez S, Silla-Martinez JM, Gabaldon T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009; 25(15):1972-3. Epub 2009/06/10. doi: 10.1093/bioinformatics/btp348 PMID: 19505945; PubMed Central PMCID: PMC2712344.

59. Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006; 22(21):2688-90. Epub 2006/08/25. doi: 10.1093/bioinformatics/btl446 PMID: 16928733.

60. Fouts DE. Phage_Finder: Automated Identification and Classification of Prophage Regions in Complete Bacterial Genome Sequences. Nucleic Acids Res. 2006; 34(20):5839-51. PMID: 17062630; PubMed Central PMCID: PMC1635311.

61. Meier-Kolthoff JP, Auch AF, Klenk HP, Göker M. Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinformatics. 2012; (14):60.

62. Fouts DE, Brinkac L, Beck E, Inman J, Sutton G. PanOCT: Automated Clustering of Orthologs Using Conserved Gene Neighborhood for Pan-Genomic Analysis of Bacterial Strains and Closely Related Species. Nucleic Acids Res. 2012; 40(22):e172. Epub 2012/08/21. doi: 10.1093/nar/gks757 PMID: 22904089; PubMed Central PMCID: PMC3526259.

63. Tettelin H, Massignani V, Cieslewicz MJ, Donati C, Medini D, Ward NL, et al. Genome analysis of multiple pathogenic isolates of *Streptococcus agalactiae*: implications for the microbial "pan-genome". Proc Natl Acad Sci U S A. 2005; 102(39):13950-5. Epub 2005/09/21. doi: 10.1073/pnas.0506758102 PMID: 16172379; PubMed Central PMCID: PMC1216834.

64. Hogg JS, Hu FZ, Janto B, Boissy R, Hayes J, Keefe R, et al. Characterization and modeling of the *Haemophilus influenzae* core and supragenomes based on the complete genomic sequences of Rd and 12 clinical nontypeable strains. Genome Biol. 2007; 8(6):R103. Epub 2007/06/07. doi: 10.1186/gb-2007-8-6-r103 PMID: 17550610; PubMed Central PMCID: PMC2394751.

65. Rasko DA, Rosovitz MJ, Myers GS, Mongodin EF, Fricke WF, Gajer P, et al. The pangenome structure of *Escherichia coli*: comparative genomic analysis of *E. coli* commensal and pathogenic isolates. J Bacteriol. 2008; 190(20):6881-93. Epub 2008/08/05. doi: 10.1128/JB.00619-08 PMID: 18676672; PubMed Central PMCID: PMC2566221.

66. Davie JJ, Earl J, de Vries SP, Ahmed A, Hu FZ, Bootsma HJ, et al. Comparative analysis and supra-genome modeling of twelve *Moraxella catarrhalis* clinical isolates. BMC Genomics. 2011; 12:70. Epub 2011/01/29. doi: 10.1186/1471-2164-12-70 PMID: 21269504; PubMed Central PMCID: PMC3045334.

67. Park J, Zhang Y, Buboltz AM, Zhang X, Schuster SC, Ahuja U, et al. Comparative genomics of the classical *Bordetella* subspecies: the evolution and exchange of virulence-associated diversity amongst closely related pathogens. BMC Genomics. 2012; 13:545. Epub 2012/10/12. doi: 10.1186/1471-2164-13-545 PMID: 23051057; PubMed Central PMCID: PMC3533505.

68. Gordienko EN, Kazanov MD, Gelfand MS. Evolution of pan-genomes of *Escherichia coli*, *Shigella* spp., and *Salmonella enterica*. J Bacteriol. 2013; 195(12):2786-92. Epub 2013/04/16. doi: 10.1128/JB.02285-12 PMID: 23585535; PubMed Central PMCID: PMC3697250.

69. Jacobsen A, Hendriksen RS, Aaresturp FM, Ussery DW, Friis C. The *Salmonella enterica* pan-genome. Microb Ecol. 2011; 62(3):487-504. Epub 2011/06/07. doi: 10.1007/s00248-011-9880-1 PMID: 21643699; PubMed Central PMCID: PMC3175032.

70. Bourhy P, Herrmann Storck C, Theodose R, Olive C, Nicolas M, Hochedez P, et al. Serovar diversity of pathogenic Leptospira circulating in the French West Indies. PLoS Negl Trop Dis. 2013; 7(3): e2114. doi: 10.1371/journal.pntd.0002114 PMID: 23516654; PubMed Central PMCID: PMC3597474.

71. Hochedez P, Escher M, Decoussy H, Pasgrimaud L, Martinez R, Rosine J, et al. Outbreak of leptospirosis among canyoning participants, Martinique, 2011. Euro Surveill. 2013; 18(18):20472. PMID: 23725775.
72. Henry CS, DeJongh M, Best AA, Frybarger PM, Linsay B, Stevens RL. High-throughput generation, optimization and analysis of genome-scale metabolic models. Nat Biotechnol. 2010; 28(9):977-82. Epub 2010/08/31. doi: 10.1038/nbt.1672 PMID: 20802497.

73. Ebrahim A, Lerman JA, Palsson BO, Hyduke DR. COBRApy: COnstraints-Based Reconstruction and Analysis for Python. BMC Syst Biol. 2013; 7:74. Epub 2013/08/10. doi: 10.1186/1752-0509-7-74 PMID: 23927696; PubMed Central PMCID: PMC3751080.

74. Orth JD, Thiele I, Palsson BO. What is flux balance analysis? Nat Biotechnol. 2010; 28(3):245-8. Epub 2010/03/10. doi: 10.1038/nbt.1614 PMID: 20212490; PubMed Central PMCID: PMC3108565.

75. Gurobi Optimization I. Gurobi Optimizer Reference Manual 2014 [2014 Nov 17]. Available from: http://www.gurobi.com.

76. Kumar VS, Maranas CD. GrowMatch: an automated method for reconciling in silico/in vivo growth predictions. PLoS Comput Biol. 2009; 5(3):e1000308. Epub 2009/03/14. doi: 10.1371/journal.pcbi.1000308 PMID: 19282964; PubMed Central PMCID: PMC2645679.

77. Fong C, Rohmer L, Radey M, Wasnick M, Brittnacher MJ. PSAT: a web tool to compare genomic neighborhoods of multiple prokaryotic genomes. BMC Bioinformatics. 2008; 9:170. Epub 2008/03/28. doi: 10.1186/1471-2105-9-170 PMID: 18366802; PubMed Central PMCID: PMC2358893.

78. Chain PS, Grafham DV, Fulton RS, Fitzgerald MG, Hostetler J, Muzny D, et al. Genomics. Genome project standards in a new era of sequencing. Science. 2009; 326(5950):236-7. Epub 2009/10/10. doi: 10.1126/science.1180614 PMID: 19815760; PubMed Central PMCID: PMC3854948.

79. Nelson KE, Weinstock GM, Highlander SK, Worley KC, Creasy HH, Wortman JR, et al. A catalog of reference genomes from the human microbiome. Science. 2010; 328(5981):994-9. Epub 2010/05/22. doi: 10.1126/science.1183605 PMID: 20489017; PubMed Central PMCID: PMC2940224.

80. Saito M, Villanueva SY, Kawamura Y, Iida K, Tomida J, Kanemaru T, et al. Leptospira idonii sp. nov., isolated from environmental water. Int J Syst Evol Microbiol. 2013; 63(Pt 7):2457-62. Epub 2012/12/04. doi: 10.1099/ijs.0.047233-0 PMID: 23203626.

81. Bourhy P, Collet L, Brisse S, Picardeau M. Leptospira mayottensis sp. nov., a pathogenic species of the genus Leptospira isolated from humans. Int J Syst Evol Microbiol. 2014; 64(Pt 12):4061-7. Epub 2014/09/25. doi: 10.1099/ijs.0.066597-0 PMID: 25249563.

82. Vernikos G, Medini D, Riley DR, Tettelin H. Ten years of pan-genome analyses. Current opinion in microbiology. 2014; 23C:148-54. Epub 2014/12/09. doi: 10.1016/j.mib.2014.11.016 PMID: 25483351.

83. Tettelin H, Riley D, Cattuto C, Medini D. Comparative genomics: the bacterial pan-genome. Current opinion in microbiology. 2008; 11(5):472-7. Epub 2008/12/18. PMID: 19086349.

84. Nally JE, Whitelegge JP, Aguilera R, Pereira MM, Blanco DR, Lovett MA. Purification and proteomic analysis of outer membrane vesicles from a clinical isolate of Leptospira interrogans serovar Copen-hageni. Proteomics. 2005; 5(1):144-52. PMID: 15672460.

85. Haake DA, Martinich C, Summers TA, Shang ES, Pruetz JD, McCoy AM, et al. Characterization of leptospiral outer membrane lipoprotein LipL36: downregulation associated with late-log-phase growth and mammalian infection. Infection and immunity. 1998; 66(4):1579-87. PMID: 9529084

86. Barnett JK, Barnett D, Bolin CA, Summers TA, Wagar EA, Cheville NF, et al. Expression and distribution of leptospiral outer membrane components during renal infection of hamsters. Infection and immunity. 1999; 67(2):853-61. PMID: 9916100

87. Haake DA, Mazel MK, McCoy AM, Milward F, Chao G, Matsunaga J, et al. Leptospiral outer membrane proteins OmpL1 and LipL41 exhibit synergistic immunoprotection. Infection and immunity. 1999; 67(12):6572-82. PMID: 10569777

88. Haake DA, Chao G, Zuerner RL, Barnett JK, Barnett D, Mazel M, et al. The leptospiral major outer membrane protein LipL32 is a lipoprotein expressed during mammalian infection. Infection and immunity. 2000; 68(4):2276-85. PMID: 10722630.

89. Guerreiro H, Croda J, Flannery B, Mazel M, Matsunaga J, Galvao Reis M, et al. Leptospiral proteins recognized during the humoral immune response to leptospirosis in humans. Infection and immunity. 2001; 69(8):4958-68. PMID: 11447174

90. Matsunaga J, Sanchez Y, Xu X, Haake DA. Osmolarity, a key environmental signal controlling expression of leptospiral proteins LigA and LigB and the extracellular release of LigA. Infection and immunity. 2005; 73(1):70-8. PMID: 15618142.

91. Matsunaga J, Lo M, Bulach DM, Zuerner RL, Adler B, Haake DA. Response of Leptospira interrogans to physiologic osmolarity: relevance in signaling the environment-to-host transition. Infection and immunity. 2007; 75(6):2864-74. Epub 2007/03/21. IAI.01619-06 [pii] doi: 10.1128/IAI.01619-06 PMID: 17371863; PubMed Central PMCID: PMC1932867.
92. Pinne M, Matsunaga J, Haake DA. Leptospiral outer membrane protein microarray, a novel approach to identification of host ligand-binding proteins. J Bacteriol. 2012; 194(22):6074-87. doi: 10.1128/JB.01119-12 PMID: 22961849; PubMed Central PMCID: PMC3486348.

93. Ristow P, Bourhy P, da Cruz McBride FW, Figueira CP, Huerre M, Ave P, et al. The OmpA-like protein Loa22 is essential for leptospiral virulence. PLoS Pathog. 2007; 3(7):e97. PMID: 17630832.

94. Dupont CL, Rusch DB, Yooseph S, Lombardo MJ, Richter RA, Valas R, et al. Genomic insights to SAR86, an abundant and uncultivated marine bacterial lineage. Isme J. 2012; 6(6):1186-99. Epub 2011/12/16. doi: 10.1038/ismej.2011.189 PMID: 22170421; PubMed Central PMCID: PMC3358033.

95. Rose RW, Bruser T, Kissinger JC, Pohlschroder M. Adaptation of protein secretion to extremely high-salt conditions by extensive use of the twin-arginine translocation pathway. Molecular microbiology. 2002; 45(4):943-50. Epub 2002/08/16. PMID: 12180915.

96. Haft DH, Selengut JD, Richter RA, Harkins D, Basu MK, Beck E. TIGRFAMs and Genome Properties in 2013. Nucleic acids research. 2013; 41 (Database issue):D387-95. Epub 2012/12/01. doi: 10.1093/nar/gks1234 PMID: 23197656; PubMed Central PMCID: PMC3531188.

97. Crooks GE, Hon G, Chandonia JM, Brenner SE. WebLogo: a sequence logo generator. Genome Res. 2004; 14(6):1188-90. Epub 2004/06/03. doi: 10.1101/gr.849004 PMID: 15173120; PubMed Central PMCID: PMC419797.

98. Drozd M, Gangaiah D, Liu Z, Rajashekara G. Contribution of TAT system translocated PhoX to *Cam-pylobacter jejuni* phosphate metabolism and resilience to environmental stresses. PLoS One. 2011; 6 (10):e26336. Epub 2011/10/27. doi: 10.1371/journal.pone.0026336 PMID: 22028859; PubMed Central PMCID: PMC3197622.

99. Gonnet P, Rudd KE, Lisacek F. Fine-tuning the prediction of sequences cleaved by signal peptidase II: a curated set of proven and predicted lipoproteins of *Escherichia coli* K-12. Proteomics. 2004; 4 (6):1597-613. Epub 2004/06/03. doi: 10.1002/pmic.200300749 PMID: 15174130.

100. Setubal JC, Reis M, Matsunaga J, Haake DA. Lipoprotein computational prediction in spirochaetal genomes. Microbiology. 2006; 152(Pt 1):113-21. Epub 2005/12/31. doi: 10.1099/mic.0.28317-0 PMID: 16385121; PubMed Central PMCID: PMC2667199.

101. Paetzel M, Dalbey RE, Strynadka NC. Crystal structure of a bacterial signal peptidase apoenzyme: implications for signal peptide binding and the Ser-Lys dyad mechanism. J Biol Chem. 2002; 277 (11):9512-9. Epub 2001/12/14. doi: 10.1074/jbc.M110983200 PMID: 11741964.

102. Sigrist CJ, de Castro E, Cerutti L, Cuche BA, Hulo N, Bridge A, et al. New and continuing developments at PROSITE. Nucleic acids research. 2013; 41 (Database issue):D344-7. Epub 2012/11/20. doi: 10.1093/nar/gks1067 PMID: 23161676; PubMed Central PMCID: PMC3531220.

103. Munoa FJ, Miller KW, Beers R, Graham M, Wu HC. Membrane topology of *Escherichia coli* prolipoprotein signal peptidase (signal peptidase II). J Biol Chem. 1991; 266(26):17667-72. Epub 1991/09/15. PMID: 1894646.

104. Bordbar A, Monk JM, King ZA, Palsson BO. Constraint-based models predict metabolic and associated cellular functions. Nat Rev Genet. 2014; 15(2):107-20. Epub 2014/01/17. doi: 10.1038/nrg3643 PMID: 24430943.

105. Monk JM, Charusanti P, Aziz RK, Lerman JA, Premyodhin N, Orth JD, et al. Genome-scale metabolic reconstructions of multiple *Escherichia coli* strains highlight strain-specific adaptations to nutritional environments. Proc Natl Acad Sci U S A. 2013; 110(50):20338-43. Epub 2013/11/28. doi: 10.1073/pnas.1307797110 PMID: 24277855; PubMed Central PMCID: PMC3864276.

106. Woodson JD, Escalante-Semerena JC. CbiZ, an amidohydrolase enzyme required for salvaging the coenzyme B12 precursor cobinamide in archaea. Proc Natl Acad Sci U S A. 2004; 101(10):3591-6. doi: 10.1073/pnas.0305939101 PMID: 14990804; PubMed Central PMCID: PMC373507.

107. Adler B. History of leptospirosis and leptospira. Current topics in microbiology and immunology. 2015; 387:1-9. doi: 10.1007/978-3-662-45059-8_1 PMID: 25388129.

108. Adler B. Vaccines against leptospirosis. Current topics in microbiology and immunology. 2015; 387:251-72. doi: 10.1007/978-3-662-45059-8_10 PMID: 25388138.

109. Faine M, Adler B, Bolin C, Perolat P. *Leptospira* and Leptospirosis. Melbourne: MedScience; 1999.

110. Bulach DM, Kalambaheti T, de la Pena-Moctezuma A, Adler B. Lipopolysaccharide biosynthesis in Leptospira. J Mol Microbiol Biotechnol. 2000; 2(4):375-80. PMID: 11075908

111. de la Pena-Moctezuma A, Bulach DM, Adler B. Genetic differences among the LPS biosynthetic loci of serovars of Leptospira interrogans and Leptospira borgpetersenii. FEMS Immunol Med Microbiol. 2001; 31(1):73-81. Epub 2001/07/31. S0928-8244(01)00245-0 [pii]. PMID: 11476985.

112. de la Pena-Moctezuma A, Bulach DM, Kalambaheti T, Adler B. Comparative analysis of the LPS biosynthetic loci of the genetic subtypes of serovar Hardjo: Leptospira interrogans subtype Hardjoprajitno
and Leptospira borgpetersenii subtype Hardjobovis. FEMS microbiology letters. 1999; 177(2):319-
26. PMID: 10474199

113. Levett PN, Morey RE, Galloway RL, Steigerwalt AG. Leptospira broomii sp. nov., isolated from humans with leptospirosis. Int J Syst Evol Microbiol. 2006; 56(Pt 3):671-3. Epub 2006/03/04. 56/3/671 [pii] doi: 10.1099/ijs.0.63783-0 PMID: 16514048.

114. Perolat P, Chappel RJ, Adler B, Baranton G, Bulach DM, Billinghurst ML, et al. *Leptospira fainei* sp. nov., isolated from pigs in Australia. Int J Syst Bacteriol. 1998; 48 Pt 3:851-8. PMID: 9734039.

115. Petersen AM, Boye K, Blom J, Schlichting P, Krogfelt KA. First isolation of Leptospira fainei serovar Hurstbridge from two human patients with Weil's syndrome. J Med Microbiol. 2001; 50(1):96-100. PMID: 11192512

116. Chappel RJ, Khalik DA, Adler B, Bulach DM, Faine S, Perolat P, et al. Serological titres to Leptospira fainei serovar hurstbridge in human sera in Australia. Epidemiol Infect. 1998; 121(2):473-5. PMID: 9825801.

117. de la Pena-Moctezuma A, Bulach DM, Kalambaheti T, Adler B. Comparative analysis of the LPS bio-synthetic loci of the genetic subtypes of serovar Hardjo: *Leptospira interrogans* subtype Hardjoprajitno and *Leptospira borgpetersenii* subtype Hardjobovis. FEMS microbiology letters. 1999; 177(2):319-26. Epub 1999/09/04. PMID: 10474199.

118. Bulach DM, Kalambaheti T, de la Pena-Moctezuma A, Adler B. Functional analysis of genes in the rfb locus of *Leptospira borgpetersenii* serovar Hardjo subtype Hardjobovis. Infection and immunity. 2000; 68(7):3793-8. Epub 2000/06/17. PMID: 10858186; PubMed Central PMCID: PMC101650.

119. Raetz CR, Whitfield C. Lipopolysaccharide endotoxins. Annual review of biochemistry. 2002; 71:635-700. Epub 2002/06/05. doi: 10.1146/annurev.biochem.71.110601.135414 PMID: 12045108; PubMed Central PMCID: PMC2569852.

120. Cuthbertson L, Mainprize IL, Naismith JH, Whitfield C. Pivotal roles of the outer membrane polysaccharide export and polysaccharide copolymerase protein families in export of extracellular polysaccharides in gram-negative bacteria. Microbiology and molecular biology reviews: MMBR. 2009; 73 (1):155-77. Epub 2009/03/05. doi: 10.1128/MMBR.00024-08 PMID: 19258536; PubMed Central PMCID: PMC2650888.

121. Keenleyside WJ, Perry M, Maclean L, Poppe C, Whitfield C. A plasmid-encoded rfbO:54 gene cluster is required for biosynthesis of the O:54 antigen in Salmonella enterica serovar Borreze. Molecular microbiology. 1994; 11(3):437-48. Epub 1994/02/01. PMID: 7512186.

122. Ricaldi JN, Fouts DE, Selengut JD, Harkins DM, Moreno A, Lehmann JS, et al. Whole genome analysis of *Leptospira licerasiae* provides insight into Leptospiral evolution and pathogenicity. PLoS Negl Trop Dis. 2012; 6:e1853. doi: 10.1371/journal.pntd.0001853 PMID: 23145189

123. Que-Gewirth NL, Ribeiro AA, Kalb SR, Cotter RJ, Bulach DM, Adler B, et al. A methylated phosphate group and four amide-linked acyl chains in *Leptospira interrogans* lipid A. The membrane anchor of an unusual lipopolysaccharide that activates TLR2. J Biol Chem. 2004; 279(24):25420-9. Epub 2004/03/27. doi: 10.1074/jbc.M400598200 PMID: 15044492; PubMed Central PMCID: PMC2556802.

124. Ricaldi JN, Matthias MA, Vinetz JM, Lewis AL. Expression of sialic acids and other nonulosonic acids in *Leptospira*. BMC Microbiol. 2012; 12:161. Epub 2012/08/03. doi: 10.1186/1471-2180-12-161 PMID: 22853805; PubMed Central PMCID: PMC3438082.

125. Ristow P, Bourhy P, da Cruz McBride FW, Figueira CP, Huerre M, Ave P, et al. The OmpA-like protein Loa22 is essential for leptospiral virulence. PLoS pathogens. 2007; 3(7):e97. Epub 2007/07/17. doi: 10.1371/journal.ppat.0030097 PMID: 17630832; PubMed Central PMCID: PMC1914066.

126. Brussow H, Hendrix RW. Phage genomics: small is beautiful. Cell. 2002; 108:13-6. PMID: 11792317

127. Rohwer F. Global phage diversity. Cell. 2003; 113:141. PMID: 12705861

128. Saint Girons I, Margarita D, Amouriaux P, Baranton G. First isolation of bacteriophages for a spirochaete: potential genetic tools for *Leptospira*. Res Microbiol. 1990; 141:1131-8. PMID: 2092364

129. Saint Girons I, Bourhy P, Ottone C, Picardeau M, Yelton D, Hendrix RW, et al. The LE1 bacteriophage replicates as a plasmid within *Leptospira biflexa*: construction of an *L. biflexa-Escherichia coli* shuttle vector. J Bacteriol. 2000; 182:5700-5. PMID: 11004167

130. Bourhy P, Frangeul L, Couve E, Glaser P, Saint Girons I, Picardeau M. Complete nucleotide sequence of the LE1 prophage from the spirochete *Leptospira biflexa* and characterization of its replication and partition functions. J Bacteriol. 2005; 187:3931-40. PMID: 15937155

131. Dobrindt U, Hochhut B, Hentschel U, Hacker J. Genomic islands in pathogenic and environmental microorganisms. Nat Rev Microbiol. 2004; 2:414-24. PMID: 15100694

132. Qin JH, Zhang Q, Zhang ZM, Zhong Y, Yang Y, Hu BY, et al. Identification of a novel prophage-like gene cluster actively expressed in both virulent and avirulent strains of Leptospira interrogans serovar Lai. Infect Immun. 2008; 76:2411-9. doi: 10.1128/IAI.01730-07 PMID: 18362131
133. Bourhy P, Salaün L, Lajus A, Médigue C, Boursaux-Eude C, Picardeau M. A genomic island of the pathogen Leptospira interrogans serovar Lai can excise from its chromosome. Infect Immun. 2007; 75:677-83. PMID: 17118975

134. Ackermann H-W. 5500 Phages examined in the electron microscope. Arch Virol. 2006; 152:227-43. PMID: 17051420

135. Saint Girons I, Bourhy P, Ottone C, Picardeau M, Yelton D, Hendrix RW, et al. The LE1 bacteriophage replicates as a plasmid within leptospira biflexa: construction of an L. biflexa-escherichia coli shuttle vecto. J Bacteriol. 2000; 182(20):5700-5. PMID: 11004167

136. Haft DH, Selengut J, Mongodin EF, Nelson KE. A guild of 45 CRISPR-associated (Cas) protein families and multiple CRISPR/Cas subtypes exist in prokaryotic genomes. PLoS Comput Biol. 2005; 1(6): e60. Epub 2005/11/18. doi: 10.1371/journal.pcbi.0010060 PMID: 16292354; PubMed Central PMCID: PMC1282333.

137. Makarova KS, Haft DH, Barrangou R, Brouns SJ, Charpentier E, Horvath P, et al. Evolution and classification of the CRISPR-Cas systems. Nature reviews Microbiology. 2011; 9(6):467-77. Epub 2011/05/10. doi: 10.1038/nrmicro2577 PMID: 21552286; PubMed Central PMCID: PMC3380444.

138. Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009; 10:421. Epub 2009/12/17. doi: 10.1186/1471-2105-10-421 PMID: 20003500; PubMed Central PMCID: PMC2803857.

139. Nickerson NN, Joag V, McGavin MJ. Rapid autocatalytic activation of the M4 metalloprotease aureo-lysin is controlled by a conserved N-terminal fungalysin-thermolysin-propeptide domain. Mol Microbiol. 2008; 69(6):1530-43. doi: 10.1111/j.1365-2958.2008.06384.x PMID: 18673454.

140. Kassegne K, Hu W, Ojcius DM, Sun D, Ge Y, Zhao J, et al. Identification of collagenase as a critical virulence factor for invasiveness and transmission of pathogenic Leptospira species. J Infect Dis. 2014; 209(7):1105-15. doi: 10.1093/infdis/jit659 PMID: 24277745.

141. Janwitthayanan W, Keelawat S, Payungporn S, Lowanitchapat A, Suwancharoen D, Poovorawan Y, et al. In vivo gene expression and immunoreactivity of Leptospira collagenase. Microbiol Res. 2013; 168(5):268-72. doi: 10.1016/j.micres.2012.12.005 PMID: 23305770.

142. Matsunaga J, Barocchi MA, Croda J, Young TA, Sanchez Y, Siqueira I, et al. Pathogenic *Leptospira* species express surface-exposed proteins belonging to the bacterial immunoglobulin superfamily. Molecular microbiology. 2003; 49(4):929-45. Epub 2003/08/02. PMID: 12890019; PubMed Central PMCID: PMC1237129.

143. Chang YF, Chen CS, Palaniappan RU, He H, McDonough SP, Barr SC, et al. Immunogenicity of the recombinant leptospiral putative outer membrane proteins as vaccine candidates. Vaccine. 2007; 25 (48):8190-7. Epub 2007/10/16. doi: 10.1016/j.vaccine.2007.09.020 PMID: 17936448.

144. Hamburger ZA, Brown MS, Isberg RR, Bjorkman PJ. Crystal structure of invasin: a bacterial integrin-binding protein. Science. 1999; 286(5438):291-5. Epub 1999/10/09. PMID: 10514372.

145. Luo Y, Frey EA, Pfuetzner RA, Creagh AL, Knoechel DG, Haynes CA, et al. Crystal structure of enteropathogenic *Escherichia coli* intimin-receptor complex. Nature. 2000; 405(6790):1073-7. Epub 2000/07/13. doi: 10.1038/35016618 PMID: 10890451.

146. Palaniappan RU, Chang YF, Jusuf SS, Artiushin S, Timoney JF, McDonough SP, et al. Cloning and molecular characterization of an immunogenic LigA protein of *Leptospira interrogans*. Infection and immunity. 2002; 70(11):5924-30. Epub 2002/10/16. PMID: 12379666; PubMed Central PMCID: PMC130282.

147. McBride AJ, Cerqueira GM, Suchard MA, Moreira AN, Zuerner RL, Reis MG, et al. Genetic diversity of the Leptospiral immunoglobulin-like (Lig) genes in pathogenic Leptospira spp. Infection, genetics and evolution: journal of molecular epidemiology and evolutionary genetics in infectious diseases. 2009; 9 (2):196-205. Epub 2008/11/26. doi: 10.1016/j.meegid.2008.10.012 PMID: 19028604; PubMed Central PMCID: PMC2812920.

148. Cerqueira GM, McBride AJ, Picardeau M, Ribeiro SG, Moreira AN, Morel V, et al. Distribution of the leptospiral immunoglobulin-like (lig) genes in pathogenic *Leptospira* species and application of *ligB* to typing leptospiral isolates. J Med Microbiol. 2009; 58(Pt 9):1173-81. Epub 2009/06/17. doi: 10.1099/jmm.0.009175-0 PMID: 19528180; PubMed Central PMCID: PMC2887549.

149. Lehmann JS, Fouts DE, Haft DH, Cannella AP, Ricaldi JN, Brinkac L, et al. Pathogenomic Inference of Virulence-Associated Genes in *Leptospira interrogans*. PLoS neglected tropical diseases. 2013; 7 (10): e2468. doi: 10.1371/journal.pntd.0002468 PMID: 24098822

150. Lambert A, Picardeau M, Haake DA, Sermswan RW, Srikram A, Adler B, et al. FlaA proteins in *Lep-tospira interrogans* are essential for motility and virulence but are not required for formation of the fla-gellum sheath. Infection and immunity. 2012; 80(6):2019-25. Epub 2012/03/28. doi: 10.1128/IAI.00131-12 PMID: 22451522; PubMed Central PMCID: PMC3370569.