<title>Introduction</title> <p>Leptospirosis is one of the most widespread zoonotic diseases in the world and is caused by pathogenic spirochetes within the genus <italic>Leptospira</italic>. Spirochetes belong to an ancient branch of eubacteria, with <italic>Leptospira</italic> representing its deepest division ##REF##1917844##[1]##. <italic>Leptospira</italic> are genetically diverse bacteria. Genetic classification of this genus is based on DNA homology and divides pathogenic <italic>Leptospira</italic> into seven main species: <italic>L. interrogans</italic>, <italic>L. borgpetersenii</italic>, <italic>L. weilii</italic>, <italic>L. noguchii</italic>, <italic>L. santarosai</italic>, <italic>L. kirschneri and L. alexanderi</italic>
##UREF##0##[2]##–##UREF##1##[4]##. In addition, there are currently eleven recognized species with a saprophytic or intermediate pathogenic status, including the saprophytic species <italic>L. biflexa</italic> and <italic>L. meyeri</italic>, and <italic>L. fainei and L. inadai</italic> exemplifying species with an intermediate status ##REF##8409911##[5]##–##REF##9734039##[9]##. Whole genome sequencing of <italic>L. interrogans</italic> serovars Lai and Copenhageni and two strains of <italic>L. borgpetersenii</italic> serovar Hardjo has revealed the occurrence of frequent gene rearrangements and fragmentation, perhaps indicating a rapid adaptation to new environments by pathogenic <italic>Leptospira</italic>
##REF##16973745##[10]##–##REF##12712204##[12]##. It has been proposed that genome reduction detected in <italic>L. borgpetersenii</italic> reflects lower environmental survivability corresponding to limited potential for indirect transmission ##REF##16973745##[10]##, in contrast to <italic>L. interrogans</italic>, a species that frequently passes through surface water between mammalian hosts ##UREF##2##[13]##.</p> <p>We previously characterized the <italic>S10-spc-α</italic> ribosomal protein cluster of <italic>L. interrogans</italic> serovar Lai ##REF##10620683##[14]##. The cluster consists of 17.5 kb comprising 32 genes that, with the exception of <italic>fus</italic>, <italic>tuf</italic>, <italic>secY</italic>, <italic>adk</italic> and <italic>infA</italic>, code for ribosomal proteins. The <italic>secY</italic> gene codes for preprotein translocase, a gene that has diagnostic value and potential for resolving taxonomic questions in <italic>Leptospira</italic>
##REF##8409911##[5]##, ##REF##10620683##[14]##. Genetic organization of ribosomal proteins is highly conserved and a prototypical S10 locus may predate divergence of Archaea and Bacteria ##REF##9395406##[15]##. However, translocation of several genes throughout the <italic>S10-spc-α</italic> locus differentiates Gram-positive from Gram-negative bacteria ##REF##9395406##[15]##. The genetic organization of the <italic>L. interrogans S10-spc-α</italic> locus is unique, as it contains all genes found in the <italic>Escherichia coli</italic> locus, and all genes except <italic>map</italic> that are found in the <italic>Bacillus subtilis</italic> locus ##REF##10620683##[14]##. The <italic>L. interrogans S10-spc-α</italic> locus is not typical of other spirochetes; several genes found in the <italic>S10-spc-α</italic> locus of <italic>L. interrogans</italic> are translocated to different portions of the <italic>Borrelia burgdorferi</italic> and <italic>Treponema pallidum</italic> genomes ##REF##10620683##[14]##. Considering the high plasticity of the <italic>Leptospira</italic> genome ##REF##16973745##[10]##, ##REF##15028702##[11]##, ##REF##7690025##[16]##, it is unclear if genetic organization of the <italic>S10-spc-α</italic> locus is conserved amongst <italic>Leptospira</italic>, or if the genetic organization shared among <italic>Borrelia</italic> and <italic>Treponema</italic> may occur among some <italic>Leptospira</italic> species, and predate divergence of <italic>Leptospira</italic> from other spirochete genera.</p> <p>In this study, we examined genetic organization and content of the <italic>S10-spc-α</italic> locus in <italic>Leptospira</italic>, and report that this locus is highly conserved throughout the genus. These data suggest that maintenance of the <italic>S10-spc-α</italic> operon structure is essential regardless of the extent of other rearrangements that have occurred during <italic>Leptospira</italic> evolution. Comparative sequence analysis of four segments of the <italic>S10-spc-α</italic> locus provides new information on phylogenetic relationships between pathogenic <italic>Leptospira</italic>.</p>
<title>Materials and Methods</title> <title>Bacterial strains and media</title> <p>
<italic>Leptospira</italic> strains used in this study were from the reference collections of the WHO/FAO/OIE Collaborating Center for Reference and Research on Leptospirosis at KIT Biomedical Research, Amsterdam, The Netherlands, and the USDA Leptospirosis Reference Center at the National Centers for Animal Health, Ames, USA (##SUPPL##1##Table S1##). Bacteria were propagated at 30°C in EMJH liquid media as described by Ellinghausen and McCullough ##REF##14266934##[29]## as modified by Johnson and Harris ##REF##6027998##[30]##.</p> <title>DNA extraction</title> <p>
<italic>Leptospira</italic> were grown to late log phase, harvested by centrifugation, and genomic DNA was extracted using a QIAamp DNA mini kit (Qiagen, Germany) following the manufacturer's instructions. DNA concentration was determined using a Nano-Drop-1000 spectrophotometer (ThermoFisher Scientific, USA) and by visual comparison with Smart Ladder SF (Eurogentec S.A., Belgium) after agarose gel electrophoresis in 1.5% agarose gels, stained with ethidium bromide according to standard procedures ##UREF##5##[31]##.</p> <title>PCR analysis</title> <p>Adjacent and overlapping fragments from the whole <italic>S10-spc-α</italic> locus were amplified by PCR from various <italic>Leptospira</italic> strains using primers listed in ##TAB##0##Table 1## and ##SUPPL##2##S2##. Several primers were designed by cross-species alignment of available <italic>L. interrogans</italic> and <italic>L. borgpetersenii</italic> genome sequences ##REF##16973745##[10]##–##REF##12712204##[12]## and access to the <italic>L. biflexa</italic> Patoc I genome sequence before publication (D. Bulach ##REF##18270594##[17]##). In addition, an iterative approach was used to develop primers useful for sequencing <italic>secY</italic> by identifying conserved regions suitable for amplification of adjacent variable regions across divergent species for which the genome sequences are yet unavailable. Primer sets were designed to produce a series of overlapping amplification products to ensure the presence and correct location of genes in the locus.</p> <p>PCR amplifications were done on a PTC-100 Peltier Thermal Cycler (MJ Research, USA) using the following program: denaturation for 5 min. at 94°C, followed by 34 cycles consisting of annealing, 1 min at 52°C, primer extension, 2 min at 72°C, denaturation, 1 min at 94°C. PCR products were separated by agarose gel electrophoresis and visualized as described above.</p> <title>Sequencing</title> <p>PCR amplification products were purified by QIAquick PCR purification kits (Qiagen Corp.) prior to DNA sequencing. Nucleotide sequences were determined by dye termination reactions separated on ABI Prism 310 and ABI 3700 (Applied Biosystems, USA) DNA sequencers. Sequencing was done on both complementary and forward strands and repeated at least twice to obtain reliable sequence data. Sequence data were edited using Sequencher (Gene Codes Corp., USA).</p> <title>Phylogenetic analysis: Binary Analysis of PCR data</title> <p>The presence (1) or absence (0) of correctly amplified fragments within the <italic>S10-spc-α</italic> locus, for each of the analyzed species, was codified in a discrete binary 40-character matrix covering a complete set of 24 taxa, representing eleven <italic>Leptospira</italic> species (##SUPPL##1##Table S1##). The characters were weighted proportionally to fragment size and assumed sequence homology for fragments with identical estimated size. The data matrix was analyzed under parsimony criteria using the branch and bound algorithm; support for branches in the unrooted tree was estimated by bootstrapping (100 replicates) with the program PAUP* v. 4.0b10 ##UREF##6##[32]##. The inferred phylogenetic relationships are based on both gene organization and sequence variation within the complete <italic>S10-spc-α</italic> locus. Phylogenetic signals contained in this data set were evaluated by <italic>g</italic>
<sub>1</sub> estimation (g<sub>1</sub> = −0.947). The negative skew of the distribution of three lengths, under parsimony criterion, is originated from trees with low scores based in highly informative data ##REF##1624764##[33]##.</p> <title>Phylogenetic analysis: Comparative Sequence Analysis</title> <p>Sequence data from four loci within the <italic>S10-spc-α</italic> locus were obtained to conduct a distance and parsimony-based phylogenetic analysis of pathogenic <italic>Leptospira</italic> using MEGA4 and PAUP* v. 4.0b10, respectively. Nucleotide diversity and diverse sequence parameters were obtained with MEGA4 ##REF##17488738##[34]## and DNASP ##REF##14668244##[35]##. The hypothesis that all mutations are selectively neutral was tested using Tajima's D test ##REF##2513255##[36]## implemented in DNASP. The confidence limits of <italic>D</italic> (two-tailed test) was obtained assuming that <italic>D</italic> follows the beta distribution and the confidence limits given in equation 47 and ##TAB##1##Table 2##, respectively in Tajima, 1989 ##REF##2513255##[36]##. Confidence intervals were also determined for Tajima's <italic>D</italic> by computer simulations using the coalescent algorithm. In distance analysis, midpoint rooted trees were obtained by the neighbor-joining method with Tamura-Nei distances ##REF##3447015##[37]##, ##REF##8336541##[38]## and the cluster support was estimated by bootstrapping with 2000 replicates ##REF##28561359##[39]##. The gaps were ignored only when they are included in the two sequences compared, using the pairwise-deletion option. In parsimony analysis, a branch-and-bound search was used with 2000 bootstraps. The homogeneity of the four partitioned data sets was evaluated using the incongruence-length difference test ##UREF##7##[40]## implemented in PAUP* v. 4.0b10. An initial tree was inferred from a data set that concatenated all the available sequenced fragments, i.e. data from the <italic>rplE</italic>, <italic>rpsN</italic>, <italic>rpsH</italic> (primer pair 301–300), <italic>rplB</italic>, <italic>rpsS</italic> (primer pair 621–625), <italic>rplC</italic>, <italic>rplD</italic> (primer pair 624–650) and <italic>secY</italic> (primer pairs G1–G2) (##TAB##0##Table 1##). This data set represents 1684 bp from each of nine <italic>Leptospira</italic> species and 19 representative strains, including <italic>L. biflexa</italic> strain Patoc I as an outgroup. Concurrently, sequence data from each of the individual fragments used in the concatenated set were analyzed separately using identical analysis methodology, to search for topological local incongruence responsible for the low support of particular nodes in the initial tree. In addition, two <italic>secY</italic> fragments (spanning the G1–G2 and SecY II–IV primer sets, respectively) were used for the reconstruction of phylogenetic relationships between 131 strains, using <italic>L. biflexa</italic> strain Patoc I as an outgroup. The <italic>secY</italic> sequences were stripped to a standard size of 1289 bp, whereas the sequences derived from G1–G2 were significantly shorter: 245 bp. Using this latter fragment, a comparison of monophyletic clustering and resolution of <italic>Leptospira</italic> species respect to the standard 16S rDNA data was done.</p> <title>ACT Alignment</title> <p>Alignment of <italic>S10-spc-α</italic> locus sequences from <italic>L. interrogans</italic>, <italic>L. borgpetersenii</italic>, and <italic>L. biflexa</italic> was done using BLASTN ##REF##2231712##[41]## with settings adjusted to identify regions having ≥80% sequence identity. Data were visualized using the Artemis Comparison Tool (ACT) ##REF##15976072##[42]##.</p>
<title>Results</title> <title>Amplification of the S10-spc-α locus of L. interrogans</title> <p>Correctly sized fragments as deduced from the positions of the primer pairs on the locus (##TAB##0##Table 1##) were obtained from all six <italic>L. interrogans</italic> strains (Lai, M20, RGA, Hond Utrecht IV, Pomona and Hardjoprajitno) for each of the 40 primer pairs tested. These data indicate that the <italic>S10-spc-α</italic> locus is well conserved in <italic>L. interrogans</italic> (##SUPPL##3##Table S3##). Remarkably, the amplification pattern of <italic>L. inadai</italic> serovar Malaya strain H6 was identical to that of <italic>L. interrogans</italic>, a finding that we note below indicates this strain was incorrectly classified previously as <italic>L. inadai</italic>.</p> <title>Comparative PCR analysis of the S10-spc-α locus in pathogenic Leptospira</title> <p>Amplification patterns of different <italic>L. borgpetersenii</italic> and <italic>L. kirschneri</italic> strains shared a high level of identity (one and two differences, respectively). However, marked strain differences were found within the species <italic>L. santarosai</italic> (8), <italic>L. noguchii</italic> (9), <italic>L. weilii</italic> (15) and <italic>L. alexanderi</italic> (14). Predictably, because genetic relatedness is used to differentiate <italic>Leptospira</italic> species, the amplification profiles varied depending on the species from which the template DNA was isolated (##SUPPL##3##Table S3##). These data show that strains composing these species likely have higher sequence variation within the <italic>S10-spc-α</italic> locus than that seen in <italic>L. interrogans</italic>. To confirm that failed PCR amplifications were due to sequence variation at or near the primer annealing sites, and not a disruption of gene synteny, a series of additional primers were designed that directed amplification from conserved sequences in adjacent genes through the regions in question. Amplification using these additional primer sets confirmed that all genes initially identified in the <italic>L. interrogans S10-spc-α</italic> locus were present throughout the same locus of all pathogenic <italic>Leptospira</italic> species. This conserved organization extends as far as <italic>fus</italic>, encoding EF-G at the 5′ end of the locus, through <italic>rpsD</italic> at the 3′ end of the locus. Thus the genetic organization of the <italic>S10-spc-α</italic> locus is conserved in all pathogenic <italic>Leptospira</italic> spp. with no signs of disruptions or translocations of genes within the locus.</p> <title>Comparative PCR analysis of the S10-spc-α locus of non-pathogenic Leptospira</title> <p>Attempts to perform PCR analysis of DNA from <italic>Leptospira</italic> species with saprophytic or intermediate (i.e. questionable) pathogenic status frequently failed to generate products or yielded anomalous sized amplicons. These data imply a marked divergence in the <italic>S10-spc-α</italic> sequence content from pathogenic <italic>Leptospira</italic> (##SUPPL##3##Table S3##). Interestingly, the amplification profile of <italic>L. meyeri</italic> strain ICF is consistent with a pathogenic status whereas the profile of <italic>L. meyeri</italic> strain Veldrat Semarang 173 is more similar to those of the saprophytic and intermediate species <italic>L. biflexa</italic>, <italic>L. fainei</italic>, and <italic>L. inadai</italic>.</p> <p>To determine if the genetic composition of this segment of the genome is different between saprophytic and pathogenic <italic>Leptospira</italic>, the corresponding regions of the <italic>L. biflexa</italic>, <italic>L. interrogans</italic>, and <italic>L. borgpetersenii</italic> genomes (GenBank accession numbers for <italic>L. interrogans</italic>, <italic>L. borgpetersenii</italic> and <italic>L. biflexa</italic> are AE016823, CP000348, CP000786, respectively) were aligned with BLAST and the results visualized by ACT (##SUPPL##0##Fig. S1##). These data show that saprophytic and pathogenic <italic>Leptospira</italic> have the same organization in the <italic>S10-spc-α</italic> locus, and the lack of successful PCR amplification is likely due to extensive sequence drift within the genus.</p> <title>Phylogenetic analysis from binary data</title> <p>The parsimony criterion was used to infer phylogenetic relationships within <italic>Leptospira</italic> from binary data. The most parsimonious tree generated from these data shows two prominent well-supported clades: 1) a basal clade, with bootstrapping value of 81%, that includes <italic>L. fainei</italic> and <italic>L. inadai</italic>, two species with intermediate pathogenic status, and the saprophytic <italic>L. meyeri</italic> strain Veldrat Semarang 173; and 2) a sister clade, supported with a 68% bootstrap value, that contains pathogenic <italic>Leptospira</italic> species (##FIG##0##Fig. 1##). Within the pathogenic clade, relationships among <italic>L. alexanderi</italic>, <italic>L. santarosai</italic>, and <italic>L. weilii</italic> species, are poorly resolved. In contrast, <italic>L. interrogans</italic>, <italic>L. kirschneri</italic>, <italic>L. borgpetersenii</italic>, and <italic>L. noguchii</italic> are clustered in a well-supported clade where <italic>L. kirschneri</italic> and <italic>L. interrogans</italic> appear as closely related siblings.</p> <p>Surprisingly, there are two exceptions to the predicted distribution of strains. <italic>L. inadai</italic> strain H6 clusters with <italic>L. interrogans</italic>, and <italic>L. meyeri</italic> strain ICF branches within pathogenic species suggesting a pathogenic status for these two strains. As we note in the Discussion section, we believe strain H6 is incorrectly classified as <italic>L. inadai</italic>.</p> <title>Multilocus sequence typing (MLST)</title> <p>Phylogenetic analysis was done on four conserved loci within the <italic>S10-spc-α</italic> locus and compared to available <italic>Leptospira rrs</italic> sequence data (GenBank accession numbers EU365895-EU365966). DNA amplification of target sequences from the intermediate strains <italic>L. fainei</italic> strain BUT 6, <italic>L. inadai</italic> strain 10, and <italic>L. meyeri</italic> strain Veldrat Semarang 173 was not successful. Therefore, these strains were not included in the analysis. The sequences for the saprophytic strain Patoc I were deduced from its genome sequence ##REF##18270594##[17]##. None of the analyzed sequences are significantly deviated from neutral expectations (<italic>P</italic>&gt;0.1). The shortest G1–G2 fragment (245 bp) showed the highest nucleotide diversity, π value of 0.14, whereas in the 300–301 fragment π was 0.09. Congruently, the mean divergence values (D) for pairwise comparisons ranged from 0.103 to 0.171 for the 300–301 and G1–G2 fragment, respectively. The lowest phylogenetic signal was obtained for 300–301 sequences. In contrast, the 621–625 and 624–650 fragments showed a phylogenetic signal slightly higher than the combined data set (##TAB##1##Table 2##).</p> <p>Distance and parsimony analysis yielded identical or similar topologies and bootstrapping values were comparable for the concordant nodes, although they were generally lower in parsimony trees. Alternative branching patterns in parsimony trees (with bootstrap value&lt;50%) occurred in nodes showing the lowest bootstrap support in distance topologies.</p> <p>In the composite tree (##FIG##1##Fig. 2E##), pathogenic strains were separated into two well-supported clades that are similar, but not identical to clades resolved in the binary tree. One clade consists of the sister sub-clades containing <italic>L. interrogans</italic> and <italic>L. noguchii</italic>, with <italic>L. kirschneri</italic> located in a basal position. This clade is consistently recovered in all topologies (##FIG##1##Fig. 2##), with the exception of the tree based on the 621–625 fragment (##FIG##1##Fig. 2B##) and parsimony topology generated from G1–G2 sequences (##FIG##1##Fig. 2D##), where <italic>L. kirschneri</italic> and <italic>L. noguchii</italic> swap their positions. The close relationships of these species are also apparent through comparative analysis using 16S rDNA sequence data (##FIG##1##Fig. 2F##) and independent binary data (##FIG##0##Fig. 1##). The second clade of the composite tree contains <italic>L. borgpetersenii</italic>, <italic>L. santarosai</italic>, <italic>L, alexanderi</italic>, <italic>L. weilii</italic> and <italic>L. meyeri</italic> strain ICF. Although, the branching pattern within this clade has lower support, the sibling relationship between <italic>L. alexanderi</italic> and <italic>L. weilii</italic> is well conserved. The relative positions of <italic>L. borgpetersenii</italic>, <italic>L. meyeri</italic> and <italic>L. santarosai</italic> are uncertain and vary depending on the data set and method of analysis (##FIG##1##Fig. 2##). In the tree inferred from the G1–G2 locus, the Celledoni and Sarmin strains of <italic>L. weilii</italic> are located in separate clades.</p> <p>The repeated findings that placed <italic>L. inadai</italic> strain H6 within the <italic>L. interrogans</italic> cluster, suggested that this strain is probably misclassified and belongs to <italic>L. interrogans</italic>. To rule out that an incorrect strain was used in our study, we repeated the sequence analysis with an H6 strain originating from the CDC collection used to establish the current taxonomic description of <italic>Leptospira</italic>
##UREF##0##[2]##. Results with the CDC H6 reference strain were identical with results obtained with our strain excluding an error in our collection.</p> <p>The <italic>S10-spc-α</italic> locus encodes ribosomal proteins that interact with rRNA, therefore ribosomal protein and rRNA sequences are expected to have parallel phylogenies. Because <italic>rrs</italic> is a well-accepted target for phylogenetic analysis we constructed a phylogenetic tree from available <italic>rrs</italic> sequence data. The <italic>rrs</italic> based phylogenetic tree is similar to the locus-deduced tree, showing close relationships between the species <italic>L. interrogans</italic>, <italic>L. kirschneri</italic> and <italic>L. noguchii</italic>. Both the clade support and genetic divergence among other <italic>Leptospira</italic> species based on <italic>rrs</italic> sequence data was lower than for <italic>S10-spc-α</italic> locus data alone, a finding consistent with a slower rate of sequence drift in rRNA than ribosomal protein genes.</p> <title>Phylogeny of secY versus its G1–G2 domain</title> <p>The 20-mer primers G1 and G2 amplify a 285 bp fragment of <italic>secY</italic>, and these primers were developed previously as a diagnostic PCR for the detection of <italic>Leptospira</italic> DNA ##REF##8409911##[5]##. A 245 bp fragment flanked by the G1–G2 primers has been shown previously to be a useful tool for discriminating between species ##UREF##3##[18]##–##REF##15863283##[21]##. This study provides an opportunity to broaden the evaluation of the G1–G2 domain by comparing the discriminative value of this domain with the majority of the <italic>secY</italic> sequence. Sequences for <italic>secY</italic> were obtained from 131 <italic>Leptospira</italic> strains (GenBank accession numbers EU357938–EU358070). The phylogenetic tree produced from <italic>secY</italic> sequence data was compared to a tree derived from the extracted sequences of the 245 bp fragment flanked by primers G1 and G2 (##FIG##2##Fig. 3##). These two trees are similar, resolve <italic>Leptospira</italic> species, and discriminate between strains. With few exceptions, all strains clustered with other members of the same species as determined by DNA-DNA hybridization analysis ##UREF##1##[4]##. Because of a limitation presented by the original G2/G2 primer pair, it does not amplify DNA from <italic>L. kirschneri</italic>; two new primers were designed (SecYII and SecYIV) that flank the G1 and G2 annealing sites. These primers amplify <italic>secY</italic> sequences from all pathogenic strains (data not shown).</p>
<title>Discussion</title> <p>Whole genome sequence analyses of different <italic>Leptospira</italic> species reveal extensive plasticity, including rearrangements, duplications, and disruptions of otherwise conserved segments of the genome ##REF##16973745##[10]##–##REF##12712204##[12]##, ##REF##18270594##[17]##. Previously, we demonstrated that <italic>L. interrogans</italic> strain Lai contained a large ribosomal protein locus spanning the <italic>S10</italic>, <italic>spc</italic>, and <italic>α</italic> loci identified in widely divergent eubacterial genera ##REF##10620683##[14]##. Notably, this entire locus is transcribed from either of two promoters upstream of <italic>fus</italic>, the first gene in the operon, and comprises one of the longest known prokaryotic transcripts ##REF##10620683##[14]##. In the present study, we show that genetic content and organization of the <italic>S10-spc-α</italic> locus is well conserved across the genus <italic>Leptospira</italic>, a finding that is somewhat remarkable given the extent of rearrangements that have disrupted synteny during <italic>Leptospira</italic> evolution. The conserved <italic>S10-spc-α</italic> organization includes the presence of the 5′ <italic>fus</italic> gene coding for elongation factor EF-G and the genes <italic>adk</italic> (adenylate kinase), <italic>infA</italic> (IF1), and <italic>rpsD</italic> (S4) located at the 3′ end of the locus, genes that are dispersed in the <italic>B. burgdorferi</italic> and <italic>T. pallidum</italic> genomes ##REF##10620683##[14]##. The genetic organization of the <italic>Leptospira S10-spc-α</italic> locus is unique among spirochetes ##REF##10620683##[14]##, and the data presented in this work support phylogenetic evidence that suggests <italic>Leptospira</italic> are one of the oldest branches in spirochete evolution. Conservation of the <italic>Leptospira S10-spc-α</italic> locus is in stark contrast to the unique organization of rRNA genes, where the <italic>rrs</italic>, <italic>rrl</italic>, and <italic>rrn</italic> genes are not closely linked to each other, but are dispersed throughout the larger of two chromosomes comprising the <italic>Leptospira</italic> genome ##REF##7690025##[16]##, ##REF##1447129##[22]##. Despite a lack of synteny for the ribosomal RNA genes, rRNA genes show limited sequence divergence. Generally, rRNA sequence conservation is a consequence of low tolerance to change due to structural constraints within the ribosome and a requirement to maintain specific binding sites for ribosomal proteins ##REF##3309345##[23]##, ##REF##6265904##[24]##.</p> <p>PCR analysis of the <italic>S10-spc-α</italic> locus showed a number of regions that were more consistently amplified than other regions (##SUPPL##3##Table S3##), suggesting that either this locus has undergone rearrangements or that sequence drift affected the efficiency of primers to faithfully bind template from diverse species. Alignment of genomic sequences spanning the <italic>S10-spc-α</italic> locus showed that the genetic organization of this locus is conserved among pathogenic and saprophytic <italic>Leptospira</italic> (##SUPPL##0##Fig. S1##). Therefore, variable success in amplifying regions of the <italic>S10-spc-α</italic> locus from diverse <italic>Leptospira</italic> species is likely due to sequence drift; <italic>Leptospira</italic> species have substantial differences in sequence composition as shown by DNA∶DNA hybridization analysis ##UREF##0##[2]##. Additionally, the PCR primers were designed primarily from the available genomic sequences of two pathogenic <italic>Leptospira</italic> serovars, and our results may be biased due to the divergence between pathogenic and saprophytic species. The binary PCR data positioned <italic>Leptospira</italic> species into two clades; one clade contained only pathogenic species, while the other contained both saprophytic species and species with intermediate pathogenic potential. One important aspect of our findings is confirmation that <italic>L. fainei</italic>, <italic>L. inadai</italic>, and <italic>L. meyeri</italic>, known to present a group of <italic>Leptospira</italic> with intermediate pathogenic potential, form a distinct cluster separate from true pathogenic species, suggesting the presence of three distinct lines of evolution within this genus.</p> <p>We selected four loci within the <italic>S10-spc-α</italic> locus that were consistently amplified from <italic>Leptospira</italic> species in initial studies to perform phylogeny studies. Phylogenetic trees deduced from the separate loci as well as from the concatenated sequence were similar, and resulted in trees each having two clades, results similar to those obtained from the binary PCR data. The clades contained branches that, with few exceptions, reflected species designations based on <italic>rrs</italic> sequence analysis ##REF##15090524##[25]##, MLST analysis ##REF##17121682##[26]##, multilocus enzyme electrophoresis (MLEE) ##UREF##4##[27]##, and DNA homology data ##UREF##0##[2]##.</p> <p>Three anomalies were found during comparison of the binary and sequence-based phylogenetic trees. First, the two strains of species <italic>L. meyeri</italic> were separated into different branches. Strain ICF was positioned in the pathogenic clusters whereas Veldrat Semarang 173 appeared in the saprophytic/intermediate pathogen cluster. This is consistent with previous reports that ICF is a pathogenic strain and Veldrat Semarang 173 is a saprophytic one ##REF##8409911##[5]##, ##REF##17020799##[28]##. The findings of this work imply that <italic>L. meyeri</italic> is composed of strains with different pathogenic potential. A second anomaly detected in this work affects the classification of strain H6. Strain H6 was designated a member of <italic>L. inadai</italic> based on DNA∶DNA reassociation analysis ##UREF##0##[2]##, but MLEE data contradicted this finding ##UREF##4##[27]##. To exclude the possibility that the discrepancy in previous studies, and in our work, was due to contamination, we analyzed strain H6 from both our collections and the reference collection at CDC used to develop the current species designations using DNA hybridization data ##UREF##0##[2]##, and found both strains had identical sequences to <italic>L. interrogans</italic>. Consequently, we recommend that strain H6 be reclassified as <italic>L. interrogans</italic>. The third anomaly involves <italic>L. weilii</italic> strains Celledoni and Sarmin. These two strains are separated into separate clades in the G1–G2 sequence-based tree (##FIG##1##Fig. 2E##). However, these two strains share the same clade based on analyses using binary data or sequence data. We believe that gene duplication and recombination events might have facilitated horizontal transfer of all or part of <italic>secY</italic> (corresponding to the G1–G2 region). It should be noted that genes contained in the <italic>S10-spc-α</italic> locus are duplicated in <italic>L. borgpetersenii</italic> strain L550, but are found as unique copy genes in all other sequenced <italic>Leptospira</italic> genomes, including <italic>L. borgpetersenii</italic> strain JB197 ##REF##16973745##[10]##. Thus, duplication of this locus and subsequent DNA acquisition via horizontal genetic transfer could facilitate stable integration of divergent <italic>secY</italic> genes.</p> <p>The <italic>S10-spc-α</italic> locus includes the <italic>secY</italic> gene encoding preprotein translocase. Primer pair G1–G2 is positioned within this gene and directs amplification of a 285 bp fragment from all pathogenic species except <italic>L. kirschneri</italic>
##REF##8409911##[5]##. Although it has been suggested that this small fragment has a high discriminating power making it useful for a quick speciation ##UREF##3##[18]##–##REF##15863283##[21]##, data supporting that contention is fragmentary. One goal of this study was to determine if analysis of the G1–G2 region provided sufficient information for <italic>Leptospira</italic> spp. discrimination. Data generated in the present study is more comprehensive than previous reports; phylogenetic trees based on the G1–G2 segment are in accordance with <italic>rrs</italic> based trees, showing that analysis of this small fragment can be used to identify species.</p> <p>Genetic analysis of the <italic>S10-spc-α</italic> locus contributes to a better understanding of <italic>Leptospira</italic> evolution. Trees generated from analysis of sequence data generated here provide analysis of more conserved loci than those studied previously ##REF##15090524##[25]##, and may be more useful in comparing evolution of the genus. A conserved, yet distinct genetic organization of this locus provides additional support for the early divergence of <italic>Leptospira</italic> from other spirochetes. Finally, from a practical standpoint, we demonstrate that analysis of a 245 bp segment of <italic>secY</italic> is suitable for rapid identification of <italic>Leptospira</italic> species.</p>
<p>Conceived and designed the experiments: BV AAA RZ JQ RAH. Performed the experiments: BV AAA RZ. Analyzed the data: NA DB JQ. Wrote the paper: BV AAA RZ NA JQ RAH.</p> <p>
<italic>S10-spc-α</italic> is a 17.5 kb cluster of 32 genes encoding ribosomal proteins. This locus has an unusual composition and organization in <italic>Leptospira interrogans</italic>. We demonstrate the highly conserved nature of this region among diverse <italic>Leptospira</italic> and show its utility as a phylogenetically informative region. Comparative analyses were performed by PCR using primer sets covering the whole locus. Correctly sized fragments were obtained by PCR from all <italic>L. interrogans</italic> strains tested for each primer set indicating that this locus is well conserved in this species. Few differences were detected in amplification profiles between different pathogenic species, indicating that the <italic>S10-spc-α</italic> locus is conserved among pathogenic <italic>Leptospira</italic>. In contrast, PCR analysis of this locus using DNA from saprophytic <italic>Leptospira</italic> species and species with an intermediate pathogenic capacity generated varied results. Sequence alignment of the <italic>S10-spc-α</italic> locus from two pathogenic species, <italic>L. interrogans</italic> and <italic>L. borgpetersenii</italic>, with the corresponding locus from the saprophyte <italic>L. biflexa</italic> serovar Patoc showed that genetic organization of this locus is well conserved within <italic>Leptospira</italic>. Multilocus sequence typing (MLST) of four conserved regions resulted in the construction of well-defined phylogenetic trees that help resolve questions about the interrelationships of pathogenic <italic>Leptospira</italic>. Based on the results of <italic>secY</italic> sequence analysis, we found that reliable species identification of pathogenic <italic>Leptospira</italic> is possible by comparative analysis of a 245 bp region commonly used as a target for diagnostic PCR for leptospirosis. Comparative analysis of <italic>Leptospira</italic> strains revealed that strain H6 previously classified as <italic>L. inadai</italic> actually belongs to the pathogenic species <italic>L. interrogans</italic> and that <italic>L. meyeri</italic> strain ICF phylogenetically co-localized with the pathogenic clusters. These findings demonstrate that the <italic>S10-spc-α</italic> locus is highly conserved throughout the genus and may be more useful in comparing evolution of the genus than loci studied previously.</p>
<title>Supporting Information</title>
<p>We thank Marc Suchard and David Haake for their continuous support and critical reading of the manuscript, David Alt and Richard Hornsby for expert technical support, and Mark Wilson for providing <italic>Leptospira</italic> cultures.</p>
<fig id="pone-0002752-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002752.g001</object-id><label>Figure 1</label><caption><title>Consensus tree based on PCR amplification data.</title><p>Majority-rule consensus tree elaborated under the parsimony criterion and based on binary data (absence/presence) coded from amplification patterns in the <italic>S10-spc-α</italic> locus for different <italic>Leptospira</italic> species. Numbers on nodes are bootstrap support after 100 replicates. Only bootstrap values above or equal to 50% are shown. Species included in the sequence analysis are coded in color. <italic>L. biflexa</italic> was used as the outgroup. CI = 0.346.</p></caption></fig> <fig id="pone-0002752-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002752.g002</object-id><label>Figure 2</label><caption><title>Phylogenetic trees based on Tamura-Nei distances and elaborated using the Neighbor-Joining method.</title><p>Distances were calculated from the 300–301 (A), 621–625 (B), 624–650 (C) and G1–G2 (D) sequence fragments within the <italic>S10-spc-α</italic> locus of pathogenic species of <italic>Leptospira</italic>. The total evidence was combined and analyzed under identical conditions (E). In addition, data available from 16S rDNA (<italic>rrs</italic>) sequences were used to obtain an alternative hypothesis for the relationships of diverse <italic>Leptospira</italic> strains (F). Dotted lines show alternative branching patterns, with bootstrapping values ≥50%, obtained in the consensus majority rule tree obtained by parsimony criterion. Numbers above branches represent the percentage of bootstrapping results (2000 replicates). Trees are drawn to scale as indicated by the bar depicted below each tree; bars represent the estimated distance in units of the number of base substitutions per site. The scale the 16S rRNA-based tree is expanded relative to other loci. <italic>L. biflexa</italic> was used as the outgroup.</p></caption></fig> <fig id="pone-0002752-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002752.g003</object-id><label>Figure 3</label><caption><title>Circular phylogenetic trees based in Tamura-Nei distances and elaborated using Neighbor-Joining method.</title><p>Distances were calculated from G1–G2 (A) restricted sequences or the <italic>secY</italic> sequences (B), and are based on analysis of 131 strains of pathogenic species of <italic>Leptospira</italic>. Numbers above branches represent the percentage of bootstrapping results (2000 replicates). Only bootstrap values above or equal to 50% are shown. <italic>L. biflexa</italic> was used as the outgroup. Dots indicate strains with divergent positions compared to those from DNA-DNA reassociation analysis ##UREF##0##[2]##.</p></caption></fig>
<table-wrap id="pone-0002752-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002752.t001</object-id><label>Table 1</label><caption><title>Primer pairs and positions in the <italic>S10-spc-α</italic> locus of <italic>L. interrogans</italic> serovar Lai.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Fragment</td><td align="left" rowspan="1" colspan="1">Primer pair</td><td align="left" rowspan="1" colspan="1">Position</td><td align="left" rowspan="1" colspan="1">Genes</td><td align="left" rowspan="1" colspan="1">Fragment</td><td align="left" rowspan="1" colspan="1">Primer pair</td><td align="left" rowspan="1" colspan="1">Position</td><td align="left" rowspan="1" colspan="1">Genes</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">1.</td><td align="left" rowspan="1" colspan="1">737-745</td><td align="left" rowspan="1" colspan="1">843-1435</td><td align="left" rowspan="1" colspan="1">urp</td><td align="left" rowspan="1" colspan="1">35.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">301-258</td><td align="left" rowspan="1" colspan="1">11601-12581</td><td align="left" rowspan="1" colspan="1">rplE-rplF</td></tr><tr><td align="left" rowspan="1" colspan="1">2.</td><td align="left" rowspan="1" colspan="1">740-744</td><td align="left" rowspan="1" colspan="1">1305-1873</td><td align="left" rowspan="1" colspan="1">urp</td><td align="left" rowspan="1" colspan="1">36.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">301-191</td><td align="left" rowspan="1" colspan="1">11601-12948</td><td align="left" rowspan="1" colspan="1">rplE-rplF</td></tr><tr><td align="left" rowspan="1" colspan="1">3.</td><td align="left" rowspan="1" colspan="1">748-751</td><td align="left" rowspan="1" colspan="1">1759-2493</td><td align="left" rowspan="1" colspan="1">fus</td><td align="left" rowspan="1" colspan="1">37.</td><td align="left" rowspan="1" colspan="1">314-191</td><td align="left" rowspan="1" colspan="1">12348-12948</td><td align="left" rowspan="1" colspan="1">rpsH, rplF</td></tr><tr><td align="left" rowspan="1" colspan="1">4. <xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">752-751</td><td align="left" rowspan="1" colspan="1">2269-2493</td><td align="left" rowspan="1" colspan="1">fus</td><td align="left" rowspan="1" colspan="1">38.</td><td align="left" rowspan="1" colspan="1">314-428c</td><td align="left" rowspan="1" colspan="1">12348-14047</td><td align="left" rowspan="1" colspan="1">rpsH-rpmD</td></tr><tr><td align="left" rowspan="1" colspan="1">5. <xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">752-729</td><td align="left" rowspan="1" colspan="1">2269-2832</td><td align="left" rowspan="1" colspan="1">fus</td><td align="left" rowspan="1" colspan="1">39.</td><td align="left" rowspan="1" colspan="1">314-430c</td><td align="left" rowspan="1" colspan="1">12348-14372</td><td align="left" rowspan="1" colspan="1">rpsH-rplO</td></tr><tr><td align="left" rowspan="1" colspan="1">6. <xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">735-729</td><td align="left" rowspan="1" colspan="1">2406-2832</td><td align="left" rowspan="1" colspan="1">fus</td><td align="left" rowspan="1" colspan="1">40.</td><td align="left" rowspan="1" colspan="1">802-R1c</td><td align="left" rowspan="1" colspan="1">12735-13445</td><td align="left" rowspan="1" colspan="1">rplF, rplR</td></tr><tr><td align="left" rowspan="1" colspan="1">7. <xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">735-743</td><td align="left" rowspan="1" colspan="1">2406-3304</td><td align="left" rowspan="1" colspan="1">fus</td><td align="left" rowspan="1" colspan="1">41.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">802-428c</td><td align="left" rowspan="1" colspan="1">12735-14047</td><td align="left" rowspan="1" colspan="1">rplF-rpmD</td></tr><tr><td align="left" rowspan="1" colspan="1">8.</td><td align="left" rowspan="1" colspan="1">735-667</td><td align="left" rowspan="1" colspan="1">2406-4394</td><td align="left" rowspan="1" colspan="1">fus, tuf</td><td align="left" rowspan="1" colspan="1">42.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">802-430c</td><td align="left" rowspan="1" colspan="1">12735-14372</td><td align="left" rowspan="1" colspan="1">rplF-rplO</td></tr><tr><td align="left" rowspan="1" colspan="1">9. <xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">743c-706</td><td align="left" rowspan="1" colspan="1">3304-3814</td><td align="left" rowspan="1" colspan="1">fus</td><td align="left" rowspan="1" colspan="1">43.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">191c-428c</td><td align="left" rowspan="1" colspan="1">12948-14047</td><td align="left" rowspan="1" colspan="1">rplF-rpmD</td></tr><tr><td align="left" rowspan="1" colspan="1">10.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">743c-667</td><td align="left" rowspan="1" colspan="1">3304-4394</td><td align="left" rowspan="1" colspan="1">fus, tuf</td><td align="left" rowspan="1" colspan="1">44.</td><td align="left" rowspan="1" colspan="1">191c-430c</td><td align="left" rowspan="1" colspan="1">12948-14372</td><td align="left" rowspan="1" colspan="1">rplF-rplO</td></tr><tr><td align="left" rowspan="1" colspan="1">11.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">800-660</td><td align="left" rowspan="1" colspan="1">3683-4327</td><td align="left" rowspan="1" colspan="1">fus, tuf</td><td align="left" rowspan="1" colspan="1">45.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">R1-428c</td><td align="left" rowspan="1" colspan="1">13427-14047</td><td align="left" rowspan="1" colspan="1">rplR-rpmD</td></tr><tr><td align="left" rowspan="1" colspan="1">12.</td><td align="left" rowspan="1" colspan="1">800-667</td><td align="left" rowspan="1" colspan="1">3683-4394</td><td align="left" rowspan="1" colspan="1">fus, tuf</td><td align="left" rowspan="1" colspan="1">46.</td><td align="left" rowspan="1" colspan="1">R1-430c</td><td align="left" rowspan="1" colspan="1">13427-14372</td><td align="left" rowspan="1" colspan="1">rplR-rplO</td></tr><tr><td align="left" rowspan="1" colspan="1">13.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">657-654</td><td align="left" rowspan="1" colspan="1">4350-5255</td><td align="left" rowspan="1" colspan="1">tuf</td><td align="left" rowspan="1" colspan="1">47.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">428-430c</td><td align="left" rowspan="1" colspan="1">14047-14372</td><td align="left" rowspan="1" colspan="1">rpmD, rplO</td></tr><tr><td align="left" rowspan="1" colspan="1">14.</td><td align="left" rowspan="1" colspan="1">657-624c</td><td align="left" rowspan="1" colspan="1">4350-5976</td><td align="left" rowspan="1" colspan="1">tuf-rplC</td><td align="left" rowspan="1" colspan="1">48.</td><td align="left" rowspan="1" colspan="1">428-G2c</td><td align="left" rowspan="1" colspan="1">14047-15468</td><td align="left" rowspan="1" colspan="1">rpmD-secY</td></tr><tr><td align="left" rowspan="1" colspan="1">15.</td><td align="left" rowspan="1" colspan="1">659-648</td><td align="left" rowspan="1" colspan="1">4438-5465</td><td align="left" rowspan="1" colspan="1">tuf</td><td align="left" rowspan="1" colspan="1">49.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">430-G2c</td><td align="left" rowspan="1" colspan="1">14372-15468</td><td align="left" rowspan="1" colspan="1">rplO, secY</td></tr><tr><td align="left" rowspan="1" colspan="1">16.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">732-624c</td><td align="left" rowspan="1" colspan="1">5240-5976</td><td align="left" rowspan="1" colspan="1">tuf-rplC</td><td align="left" rowspan="1" colspan="1">50</td><td align="left" rowspan="1" colspan="1">634-635</td><td align="left" rowspan="1" colspan="1">14643-16387</td><td align="left" rowspan="1" colspan="1">
<italic>rplO-adk</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">17.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">647-618</td><td align="left" rowspan="1" colspan="1">5297-5806</td><td align="left" rowspan="1" colspan="1">tuf, rpsJ</td><td align="left" rowspan="1" colspan="1">51.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">443-G2c</td><td align="left" rowspan="1" colspan="1">15276-15468</td><td align="left" rowspan="1" colspan="1">
<italic>secY</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">18.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">647-624c</td><td align="left" rowspan="1" colspan="1">5297-5976</td><td align="left" rowspan="1" colspan="1">tuf-rplC</td><td align="left" rowspan="1" colspan="1">52.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">443-G1</td><td align="left" rowspan="1" colspan="1">15276-15752</td><td align="left" rowspan="1" colspan="1">
<italic>secY</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">19.</td><td align="left" rowspan="1" colspan="1">624-650</td><td align="left" rowspan="1" colspan="1">5976-6790</td><td align="left" rowspan="1" colspan="1">rplC, rplD</td><td align="left" rowspan="1" colspan="1">53.<xref ref-type="table-fn" rid="nt102">**</xref>
</td><td align="left" rowspan="1" colspan="1">SecYII-SecYIV</td><td align="left" rowspan="1" colspan="1">15289-15946</td><td align="left" rowspan="1" colspan="1">
<italic>secY</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">20.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">624-644</td><td align="left" rowspan="1" colspan="1">5976-7151</td><td align="left" rowspan="1" colspan="1">rplC-rplW</td><td align="left" rowspan="1" colspan="1">54.</td><td align="left" rowspan="1" colspan="1">G2-G1</td><td align="left" rowspan="1" colspan="1">15468-15752</td><td align="left" rowspan="1" colspan="1">
<italic>secY</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">21.</td><td align="left" rowspan="1" colspan="1">624-621c</td><td align="left" rowspan="1" colspan="1">5976-7847</td><td align="left" rowspan="1" colspan="1">rplC-rplB</td><td align="left" rowspan="1" colspan="1">55.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">G2-444</td><td align="left" rowspan="1" colspan="1">15468-15970</td><td align="left" rowspan="1" colspan="1">
<italic>secY</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">22.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">651-644</td><td align="left" rowspan="1" colspan="1">6883-7151</td><td align="left" rowspan="1" colspan="1">rplD, rplW</td><td align="left" rowspan="1" colspan="1">56.</td><td align="left" rowspan="1" colspan="1">G2-429</td><td align="left" rowspan="1" colspan="1">15468-16353</td><td align="left" rowspan="1" colspan="1">
<italic>secY-adK</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">23.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">643-621c</td><td align="left" rowspan="1" colspan="1">7138-7847</td><td align="left" rowspan="1" colspan="1">rplW, rplB</td><td align="left" rowspan="1" colspan="1">57.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">G2-400</td><td align="left" rowspan="1" colspan="1">15468-16640</td><td align="left" rowspan="1" colspan="1">
<italic>secY-infA</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">24.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">622-621c</td><td align="left" rowspan="1" colspan="1">7689-7847</td><td align="left" rowspan="1" colspan="1">rplB</td><td align="left" rowspan="1" colspan="1">58.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">260-458c</td><td align="left" rowspan="1" colspan="1">16616-18104</td><td align="left" rowspan="1" colspan="1">
<italic>infA-rpsD</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">25.</td><td align="left" rowspan="1" colspan="1">621-625</td><td align="left" rowspan="1" colspan="1">7847-8504</td><td align="left" rowspan="1" colspan="1">rplB, rpsS</td><td align="left" rowspan="1" colspan="1">59.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">458-507</td><td align="left" rowspan="1" colspan="1">18104-18696</td><td align="left" rowspan="1" colspan="1">
<italic>rpsD, rpoA</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">26.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">621-605c</td><td align="left" rowspan="1" colspan="1">7847-9082</td><td align="left" rowspan="1" colspan="1">rplB-rpsC</td><td align="left" rowspan="1" colspan="1">60.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">458-504</td><td align="left" rowspan="1" colspan="1">18104-19376</td><td align="left" rowspan="1" colspan="1">
<italic>rpsD, rpoA</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">27.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">605-460</td><td align="left" rowspan="1" colspan="1">9082-10196</td><td align="left" rowspan="1" colspan="1">rpsC-rpmC</td><td align="left" rowspan="1" colspan="1">61.</td><td align="left" rowspan="1" colspan="1">450-479</td><td align="left" rowspan="1" colspan="1">18163-19264</td><td align="left" rowspan="1" colspan="1">
<italic>rpsD, rpoA</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">28.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">801-803c</td><td align="left" rowspan="1" colspan="1">10105-10965</td><td align="left" rowspan="1" colspan="1">rpmC-rplX</td><td align="left" rowspan="1" colspan="1">62.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">477-504</td><td align="left" rowspan="1" colspan="1">18584-19376</td><td align="left" rowspan="1" colspan="1">
<italic>rpoA</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">29.</td><td align="left" rowspan="1" colspan="1">801-301c</td><td align="left" rowspan="1" colspan="1">10105-11601</td><td align="left" rowspan="1" colspan="1">rpmC-rplE</td><td align="left" rowspan="1" colspan="1">63.</td><td align="left" rowspan="1" colspan="1">477-501c</td><td align="left" rowspan="1" colspan="1">18584-19791</td><td align="left" rowspan="1" colspan="1">
<italic>rpoA, rplQ</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">30.</td><td align="left" rowspan="1" colspan="1">801-300</td><td align="left" rowspan="1" colspan="1">10105-12110</td><td align="left" rowspan="1" colspan="1">rpmC-rpsH</td><td align="left" rowspan="1" colspan="1">64.</td><td align="left" rowspan="1" colspan="1">503-480</td><td align="left" rowspan="1" colspan="1">18862-19621</td><td align="left" rowspan="1" colspan="1">
<italic>rpoA, rplQ</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">31.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">310-309</td><td align="left" rowspan="1" colspan="1">10167-10672</td><td align="left" rowspan="1" colspan="1">rpmC-rplN</td><td align="left" rowspan="1" colspan="1">65.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">478-501c</td><td align="left" rowspan="1" colspan="1">19371-19791</td><td align="left" rowspan="1" colspan="1">
<italic>rpoA, rplQ</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">32.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">310-277</td><td align="left" rowspan="1" colspan="1">10167-11107</td><td align="left" rowspan="1" colspan="1">rpmC-rplX</td><td align="left" rowspan="1" colspan="1">66.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">478-502</td><td align="left" rowspan="1" colspan="1">19371-20341</td><td align="left" rowspan="1" colspan="1">
<italic>rpoA, rplQ,</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">33.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">310-301c</td><td align="left" rowspan="1" colspan="1">10167-11601</td><td align="left" rowspan="1" colspan="1">rpmC-rplE</td><td align="left" rowspan="1" colspan="1">67.</td><td align="left" rowspan="1" colspan="1">501-502</td><td align="left" rowspan="1" colspan="1">19791-20341</td><td align="left" rowspan="1" colspan="1">
<italic>rplQ</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">34.<xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">301-300</td><td align="left" rowspan="1" colspan="1">11601-12110</td><td align="left" rowspan="1" colspan="1">rplE-rpsH</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0002752-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002752.t002</object-id><label>Table 2</label><caption><title>Nucleotide diversity, divergent estimations and parameters estimated from the sequences of 4 fragments in diverse <italic>Leptospira</italic> species.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Fragment</td><td align="left" rowspan="1" colspan="1">Sites</td><td align="left" rowspan="1" colspan="1">Polymorphic Sites</td><td align="left" rowspan="1" colspan="1">Mean D</td><td align="left" rowspan="1" colspan="1">π</td><td align="left" rowspan="1" colspan="1">Θ per site</td><td align="left" rowspan="1" colspan="1">Tajima's D</td><td align="left" rowspan="1" colspan="1">
<italic>P</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>-g1</italic>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">300–301</td><td align="left" rowspan="1" colspan="1">469</td><td align="left" rowspan="1" colspan="1">179</td><td align="left" rowspan="1" colspan="1">0.1030</td><td align="left" rowspan="1" colspan="1">0.0891</td><td align="left" rowspan="1" colspan="1">0.1011</td><td align="left" rowspan="1" colspan="1">−1.2016</td><td align="left" rowspan="1" colspan="1">&gt;0.1</td><td align="left" rowspan="1" colspan="1">−0.5670</td></tr><tr><td align="left" rowspan="1" colspan="1">621–625</td><td align="left" rowspan="1" colspan="1">479</td><td align="left" rowspan="1" colspan="1">176</td><td align="left" rowspan="1" colspan="1">0.1190</td><td align="left" rowspan="1" colspan="1">0.1034</td><td align="left" rowspan="1" colspan="1">0.1306</td><td align="left" rowspan="1" colspan="1">−0.8672</td><td align="left" rowspan="1" colspan="1">&gt;0.1</td><td align="left" rowspan="1" colspan="1">−0.9404</td></tr><tr><td align="left" rowspan="1" colspan="1">624–650</td><td align="left" rowspan="1" colspan="1">491</td><td align="left" rowspan="1" colspan="1">226</td><td align="left" rowspan="1" colspan="1">0.1440</td><td align="left" rowspan="1" colspan="1">0.1205</td><td align="left" rowspan="1" colspan="1">0.1623</td><td align="left" rowspan="1" colspan="1">−1.0739</td><td align="left" rowspan="1" colspan="1">&gt;0.1</td><td align="left" rowspan="1" colspan="1">−1.0876</td></tr><tr><td align="left" rowspan="1" colspan="1">G1–G2</td><td align="left" rowspan="1" colspan="1">245</td><td align="left" rowspan="1" colspan="1">91</td><td align="left" rowspan="1" colspan="1">0.1710</td><td align="left" rowspan="1" colspan="1">0.1434</td><td align="left" rowspan="1" colspan="1">0.1381</td><td align="left" rowspan="1" colspan="1">0.1606</td><td align="left" rowspan="1" colspan="1">&gt;0.1</td><td align="left" rowspan="1" colspan="1">−0.5914</td></tr></tbody></table></alternatives></table-wrap>
<supplementary-material content-type="local-data" id="pone.0002752.s001"><label>Figure S1</label><caption><p>Alignment of the L. biflexa, L. interrogans and L. borgpetersenii S10-spc-a genome sequences. Regions of greater than 80% sequence identity are shown as blue (between L. interrogans and L. borgpetersenii) and red (L. biflexa and L. interrogans). White regions indicate segments where sequence identity drops below 80%. Regions of similarity were determined using Blastn under default settings except the -m 8 output option was used. The display was generated using ACT. Note that the orientation of these sequences shown in the figure is consistent with the genomic sequence data in GenBank and are inverted relative to the direction of transcription. GenBank accession numbers for the genomes of L. interrogans, L. borgpetersenii and L. biflexa are AE016823, CP000348, CP000786, respectively.</p><p>(2.24 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002752.s002"><label>Table S1</label><caption><p>Leptospira strains used for the S10-spc-α locus study.</p><p>(0.25 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002752.s003"><label>Table S2</label><caption><p>All primers used for the S10-spc-α locus analysis.</p><p>(0.09 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002752.s004"><label>Table S3</label><caption><p>Amplification through the S10-spc-α operon of Leptospira spp. Positive and negative PCR scores for amplification reactions along the locus from various strains.</p><p>(0.05 MB DOC)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><label>*</label><p>Fragments used in the phylogenetic analysis from the binary data.</p></fn><fn id="nt102"><label>**</label><p>Primer pair used to produced G1–G2 sequences from all pathogenic species.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>This study was supported by EC grant ICA4-CT-2001-10086, and by the U.S. Department of Agriculture.</p></fn></fn-group>
<graphic id="pone-0002752-t001-1" xlink:href="pone.0002752.t001"/> <graphic xlink:href="pone.0002752.g001"/> <graphic id="pone-0002752-t002-2" xlink:href="pone.0002752.t002"/> <graphic xlink:href="pone.0002752.g002"/> <graphic xlink:href="pone.0002752.g003"/>
<media xlink:href="pone.0002752.s001.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002752.s002.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002752.s003.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002752.s004.doc"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["2"], "element-citation": ["\n"], "surname": ["Brenner", "Kaufmann", "Sulzer", "Steigerwalt", "Rogers"], "given-names": ["DJ", "AF", "KR", "AG", "FC"], "year": ["1999"], "article-title": ["Further determination of DNA relatedness between serogroups and serovars in the family Leptospiraceae with a proposal for "], "italic": ["Leptospira alexanderi sp. nov.", "Leptospira"], "source": ["Int J Syst Bacteriol 49 Pt"], "volume": ["2"], "fpage": ["839"], "lpage": ["858"]}, {"label": ["4"], "element-citation": ["\n"], "surname": ["Yasuda", "Steigerwalt", "Sulzer", "Kaufmann", "Rogers"], "given-names": ["P", "A", "K", "A", "F"], "year": ["1987"], "article-title": ["Deoxyribonucleic acid relatedness between serogroups and serovars in the family Leptospiraceae with proposals for seven new "], "italic": ["Leptospira"], "source": ["Int J Syst Bacteriol"], "volume": ["37"], "fpage": ["407"], "lpage": ["415"]}, {"label": ["13"], "element-citation": ["\n"], "surname": ["Faine", "Adler", "Bolin", "Perolat"], "given-names": ["S", "B", "C", "P"], "year": ["1999"], "article-title": ["\n"], "italic": ["Leptospira"], "publisher-loc": ["Melbourne, Australia"], "publisher-name": ["MediSci"], "fpage": ["272"]}, {"label": ["18"], "element-citation": ["\n"], "surname": ["De los A Valverede", "Raminez", "Montes de Octa", "Goris", "Ahmed"], "given-names": ["MA", "JM", "LG", "MGA", "N"], "year": ["2008"], "article-title": ["Arenal, a new "], "italic": ["Leptospira"], "source": ["Infection, Genetics and Evolution"], "comment": ["In press"]}, {"label": ["27"], "element-citation": ["\n"], "surname": ["Letocart", "Boerlin", "Boerlin-Petzold", "Goudet", "Baranton"], "given-names": ["M", "P", "F", "J", "G"], "year": ["1999"], "article-title": ["Genetic structure of the genus Leptospira by multilocus enzyme electrophoresis."], "source": ["Int J Syst Bacteriol"], "volume": ["49"], "fpage": ["231"], "lpage": ["238"]}, {"label": ["31"], "element-citation": ["\n"], "surname": ["Maniatis", "Fritsch", "Sambrook"], "given-names": ["T", "EF", "J"], "year": ["1983"], "article-title": ["Molecular Cloning: A Laboratory Manual."], "publisher-loc": ["Cold Spring Harbor, NY"], "publisher-name": ["Cold Spring Harbor Laboratory"], "fpage": ["545"]}, {"label": ["32"], "element-citation": ["\n"], "surname": ["Swofford"], "given-names": ["DL"], "year": ["1998"], "article-title": ["PAUP*: Phylogenetic Analysis Using Parsimony (and Other Methods) 4.0 Beta."], "publisher-loc": ["Sunderland, MA"], "publisher-name": ["Sinauer Associates, Inc"]}, {"label": ["40"], "element-citation": ["\n"], "surname": ["Farris", "K\u00e4llersj\u00f6", "Kluge", "Bult"], "given-names": ["JS", "M", "AG", "C"], "year": ["1995"], "article-title": ["Constructing a significance test for incongruence."], "source": ["Syst Biol"], "volume": ["44"], "fpage": ["570"], "lpage": ["572"]}]
PLoS One. 2008 Jul 16; 3(7):e2752