Molecular characterization of Treponema pallidum subsp. pallidum in Switzerland and France with a new multilocus sequence typing scheme

PLOS ONE, Jul 2018

Syphilis is an important public health problem and an increasing incidence has been noted in recent years. Characterization of strain diversity through molecular data plays a critical role in the epidemiological understanding of this re-emergence. We here propose a new high-resolution multilocus sequence typing (MLST) scheme for Treponema pallidum subsp. pallidum (TPA). We analyzed 30 complete and draft TPA genomes obtained directly from clinical samples or from rabbit propagated strains to identify suitable typing loci and tested the new scheme on 120 clinical samples collected in Switzerland and France. Our analyses yielded three loci with high discriminatory power: TP0136, TP0548, and TP0705. Together with analysis of the 23S rRNA gene mutations for macrolide resistance, we propose these loci as MLST for TPA. Among clinical samples, 23 allelic profiles as well as a high percentage (80% samples) of macrolide resistance were revealed. The new MLST has higher discriminatory power compared to previous typing schemes, enabling distinction of TPA from other treponemal bacteria, distinction between the two main TPA clades (Nichols and SS14), and differentiation of strains within these clades.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

http://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0200773&type=printable

Molecular characterization of Treponema pallidum subsp. pallidum in Switzerland and France with a new multilocus sequence typing scheme

July Molecular characterization of Treponema pallidum subsp. pallidum in Switzerland and France with a new multilocus sequence typing scheme Linda Grillova 0 1 Tanika Bawa 1 Lenka Mikalova 0 1 Angèle Gayet-Ageron 1 Kay Nieselt 1 Michal Strouhal 0 1 Patrice Sednaoui 1 Tristan Ferry 1 Matthias Cavassini 1 Stephan Lautenschlager 1 Fabrizio Dutly 1 Marta Pla-DÂõaz 1 Michael KruÈ tzen 1 3 Fernando Gonza lez-Candelas 1 Homayoun C. Bagheri 1 David SÏ majs 0 1 Natasha Arora 1 Philipp P. Bosshard 1 2 0 Department of Biology, Masaryk University , Brno , Czech Republic , 2 Department of Fundamental Neuroscience, University of Geneva , Geneva , Switzerland , 3 Division of Clinical Epidemiology, University of Geneva Hospitals and Faculty of Medicine , Geneva, Switzerland, 4 Center for Bioinformatics , University of TuÈbingen , TuÈ bingen, Germany, 5 Medical Biology, Institute Alfred Fournier , Paris, France, 6 Hospices civils de Lyon, Lyon , France , 7 Department of Infectious Diseases, Centre Hospitalier Universitaire Vaudois, Lausanne, Switzerland, 8 Department of Dermatology , Triemlispital, Zurich , Switzerland , 9 IMD Institut fuÈ r medizinische & molekulare Diagnostik AG, Zurich, Switzerland, 10 Unidad Mixta Infecci oÂn y Salud Pu blica FISABIO/Universidad de Valencia, CIBER in Epidemiology and Public Health , Valencia , Spain 1 Editor: Baochuan Lin, Defense Threat Reduction Agency , UNITED STATES 2 Department of Dermatology, University Hospital Zurich , Zurich , Switzerland , 16 Faculty of Medicine, University of Zurich , Zurich , Switzerland 3 Department of Anthropology, University of Zurich , Zurich , Switzerland , 12 Repsol Technology Center , Madrid , Spain , 13 Institute for Evolutionary Biology and Environmental Studies, University of Zurich , Zurich , Switzerland , 14 Zurich Institute of Forensic Medicine, University of Zurich , Zurich , Switzerland Syphilis is an important public health problem and an increasing incidence has been noted in recent years. Characterization of strain diversity through molecular data plays a critical role in the epidemiological understanding of this re-emergence. We here propose a new high-resolution multilocus sequence typing (MLST) scheme for Treponema pallidum subsp. pallidum (TPA). We analyzed 30 complete and draft TPA genomes obtained directly from clinical samples or from rabbit propagated strains to identify suitable typing loci and tested the new scheme on 120 clinical samples collected in Switzerland and France. Our analyses yielded three loci with high discriminatory power: TP0136, TP0548, and TP0705. Together with analysis of the 23S rRNA gene mutations for macrolide resistance, we propose these loci as MLST for TPA. Among clinical samples, 23 allelic profiles as well as a high percentage (80% samples) of macrolide resistance were revealed. The new MLST has higher discriminatory power compared to previous typing schemes, enabling distinction of TPA from other treponemal bacteria, distinction between the two main TPA clades (Nichols and - Data Availability Statement: All relevant data are within the paper and its Supporting Information files. Funding: This work was supported by the University of Zurich, and the Grant Agency of the Czech Republic [GA17-25455S] to DS and [J1725589Y] to MS, and by the Ministry of Health of the Czech Republic [17-31333A] to DS. Repsol Technology Center provided support in the form of salary for an author [HB], but did not have any SS14), and differentiation of strains within these clades. additional role in the study choice, design, data collection and analysis, decision to publish, or preparation of the manuscript. The specific roles of authors are articulated in the `author contributions' section. Introduction Treponema pallidum subsp. pallidum (TPA), the bacterium responsible for syphilis, causes more than 5.6 million new syphilis cases annually and an increase in number of cases has been noted worldwide in the last few years [ 1 ]. The recent global outbreaks of syphilis infections are driven by multiple factors including suboptimal diagnostics in the early phase of infection, patients not adopting safer sex rules, mutations conferring macrolide resistance, and the lack of a syphilis vaccine [ 2 ]. However, a comprehensive epidemiological understanding of the underlying factors and the patterns characterizing the re-emergence is still missing. Molecular typing is required in order to discriminate genetic variants and investigate their potential association with phenotypic features such as pathogenicity and antibiotic resistance. Furthermore, molecular typing provides population structure data that can be used to map the geographical distribution of strain types, draw associations between strains and specific demographic groups, and examine transmission dynamics. In some cases, molecular typing can also be a valuable diagnostic tool, as illustrated by recent cases of syphilis-bejel confusion that were resolved through genetic data [ 3, 4 ]. The molecular characterization of TPA is limited by the lack of a continuous cultivation system under in vitro conditions. Therefore, research efforts have been primarily directed at strains propagated in laboratory animals, leading to the publication of a few TPA whole genomes since 1998 [5±10]. The data from these genomes as well as that of individual genes from clinical samples have been used in the development of molecular typing schemes (Table 1). Due to the low level of variation across TPA strains (below 0.03%), these schemes generally target the most variable loci, in order to gain sufficient discriminatory power. These loci include genes predicted to code for outer membrane proteins (TP0136, TP0548) or virulence factors (tpr genes), genes comprising tandem repeats (arp and rpsA), as well as the genes associated with resistance to macrolides (23S rRNA genes) [15±17]. The most commonly used schemes are the CDCT and ECDCT with over 3,000 clinical isolates typed worldwide (reviewed in [ 18, 19 ]). However, these typing schemes were recently questioned due to possible intra-strain variability of the arp and tpr loci. These two loci were found to reveal subtype discrepancies in 11 out of 18 parallel samples taken from sero-positive individuals [20]. An alternative typing scheme, the SBMT, has been used to type over 170 clinical samples in Argentina, Belgium and the Czech Republic [12, 19, 21±23]. One of the limitations of all these typing schemes is that they generally have low discriminatory power for strains of the SS14-clade [ 19, 22 ], the more dominant of the two major clades identified RFLPb of tpr genes, determination of repeats in arp gene Sanger sequencing RFLPb of tpr genes, determination of repeats in arp gene, Sanger sequencing of TP0548 RFLPb of tpr genes, determination of repeats in arp and rpsA genes TP0136, TP0548, TP0705, Sanger sequencing 23S rRNA Year of introduction 1998 2006 2010 2010 2018 among contemporary samples, and thus do not exploit the potential diversity present in current infections. Recently, the availability of TPA genome-wide data from current infections have increased through the application of new enrichment methods to eliminate background DNA prior to the next generation sequencing (NGS) [ 24, 25 ]. While whole genome sequencing is not yet time and cost effective for the clinical setting, the availability of these new treponemal genomes provide us with the opportunity to more carefully examine diversity patterns and potential new typing loci in order to design the most up-to-date molecular typing method. In this work, we aim to i) identify candidate typing loci based on the analyses of genomewide data from TPA, ii) select loci providing the highest discriminatory power for inclusion in a new MLST typing scheme, and iii) apply the MLST to a set of 120 clinical samples from Switzerland and France collected between 2011 and 2015. Materials and methods Design of the typing scheme We identified candidate loci for TPA typing in a dataset comprising 12 whole or draft genome sequences obtained from rabbit propagated reference TPA strains and 18 draft genomes sequenced directly from clinical specimens. These samples were collected in countries of three continents (Argentina, Austria, Czech Republic, Mexico, The Netherlands, USA, and Switzerland) in the period between 1912 and 2013 [6±10, 24] (S1 Table). The dataset contained sequences from both genetically distinct groups of TPA strains, SS14 clade (n = 21) and Nichols clade (n = 9), respectively. Since the tpr genes were shown to be among the most variable chromosomal loci [ 26 ] and the stability of tprE, tprG, and tprJ genes was recently questioned [ 20 ], we excluded all tpr genes from the analyses. Other paralogous regions, intergenic regions, genes with repetitions and sites with ambiguous data were excluded as well. We selected candidate loci for TPA typing if; i) they contained at least four variable sites with a single nucleotide variant (SNV) density >0.001 (S2 Table), which is 10 times and 50 times higher than the average SNV density inside the Nichols-clade and SS14-clade, respectively, ii) they displayed at least two different haplotypes within each clade, and if iii) the SNVs were not too widespread in order to use as short as possible loci for increasing the amplification efficiency. We used MEGA6 [ 27 ] to visualize the sequence alignment and to identify parsimony-informative sites. Clinical samples, DNA isolation and qPCR We performed the genetic analysis in 120 samples collected from 120 patients diagnosed with early syphilis or syphilis of unknown stage between September 2011 and May 2015 in outpatient clinics in five cities in Switzerland and France (Zurich, n = 70; Geneva, n = 10; Lausanne, n = 2; Lyon, n = 21; and Paris, n = 17). The samples included 119 genital, anal, oral or rectal lesion swabs and 1 tissue sample. DNA isolation and routine real-time PCR for diagnosis have been performed as described previously [ 28, 29 ] (detailed in the Supplementary Methods). We had clinical data about the patient age, gender, site of infection, HIV status and syphilis stage for 93/120 patients. Nested PCR and DNA sequencing Five candidate loci (TP0136, TP0462, TP0548, TP0705 and TP0865) and the 23S rRNA genes were amplified by a nested-PCR approach with a touchdown protocol in the first PCR reactions (primer sequences and conditions detailed in S1 Supplementary Methods and S3 Table). 3 / 16 PCR products were purified using a QIAquick PCR Purification Kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions and sequenced on an automated capillary DNA sequencing system (GATC-Biotech AG, Constance, Germany). Sequence analyses were performed using the Lasergene software (DNASTAR v. 7.1.0.; DNASTAR, Madison, WI). Sequences representing allelic variants in different typing loci were deposited in GenBank under the following accession numbers: MG894082ÐMG894119 (S4 Table). For the 23S rRNA genes, positions corresponding to 2058 and 2059 in the Escherichia coli 23S rRNA gene (accession no. V00331) were checked for A!G mutations indicative of macrolide resistance [15±17, 30±31]. Comparison of typing schemes We analyzed the number of 60 bp-long repetitions in the arp gene and amplified the tprE, G, and J genes as described previously [ 11, 13 ] with minor modifications for the tpr genes, including an increase of the extension temperature from 68ÊC to 72ÊC and the use of GoTaq Flexi DNA polymerase (Promega, Mannheim, Germany). The information revealed from the RFLP of the tpr genes, determination of repeats in the arp gene and sequencing of TP0136, TP0548 and 23S rDNA loci allowed us to compare the typing efficiency of the currently used typing schemes with the MLST scheme proposed in this study. Typing stability In order to test the typing stability of the new typing scheme, we amplified all loci in eleven passages of the TPA reference strain DAL-1, which was continuously propagated in rabbits during 142 days (representing 114 TPA generations). Moreover, we tested the concordance in five samples collected in the Czech Republic during 2004±2010 and pertaining to five epidemiologically related patients. Phylogenetic analyses We generated maximum likelihood (ML) phylogenetic trees with MEGA 6 [ 27 ] using the Tamura Nei model and 1000 pseudorandom bootstrap replicates. Median-joining (MJ) networks were generated with Network version 4 [ 32 ]. We created genome-wide data MJ networks for the 30 samples published in [ 24 ] using 768 variable positions (including putative recombinant loci and excluding positions with missing data). The MJ networks for the candidate typing loci included 87 variable positions (5 candidate loci) and 74 variable positions (3 candidate loci). Deletions and insertions were counted as single events. The software Sequence Matrix 1.8 was used for the sequence concatenations [ 33 ]. Statistical methods We had a fixed number of samples (n = 120) allowing us to detect medium effect size associations between some genetic and clinical characteristics. Genetic characteristics of clinical isolates (macrolide resistance, full allelic profiles, allelic variants per gene and clade) were examined to test for associations with clinical characteristics such as patient's origin, MSM status, age, gender, HIV infection, serological results and syphilis stage (as detailed in the S1 Supplementary Methods). Statistical analyses were carried out in RStudio (version 1.1.383) using Fischer's exact tests for categorical variables and Kruskal Wallis test for continuous variables. Mean nucleotide diversity (π) was calculated with MEGA 6 [ 27 ]. 4 / 16 Ethics statement This study was approved by the ethics committees of Zurich (2016±01518), Lausanne (380/ 11), Geneva (11±151) and Lyon (L11-150). All patients provided informed written consent. Results Candidate typing loci and their discriminatory power Based on the selection criteria, we identified five genes as candidates for TPA typing: TP0136, TP0462, TP0548, TP0705 and TP0865. TP0136 had the largest SNV density across all genes examined, followed by TP0462, TP0548, TP0865 and TP0705 (S2 Table). Four of these loci code for outer membrane proteins and one for penicillin-binding protein. All but one (TP0705) have been previously identified as putative recombinant genes [ 24, 34, 35 ]. To assess the resolution power of the five candidate loci, we compared the MJ network built from their concatenated sequences (incorporating 87 variable positions) to the genome-wide data network (incorporating 768 variable positions) for 30 draft and complete genome sequences (S1 Table). The concatenated sequences of the five loci were able to distinguish 8/26 (30.8%) haplotypes identified using genome-wide data (data not shown). By only using three of the five candidate loci (TP0136, TP0548 and TP0705), the same resolution power as with all five loci was achieved. We further analyzed the five candidate loci in a combined dataset of all (n = 63) previously published sequences [ 24, 25, 36 ] and of the sequences of the 120 qPCR-positive clinical samples from this study. In this combined sequence dataset, 16, 13, 10, 6 and 5 different allele variants were found in the TP0548, TP0136, TP0705, TP0462 and TP0865 loci, respectively (Figs 1 and 2). For the SS14-clade, the highest variability was found in TP0705, followed by TP0548, TP0136, TP0462 and TP0865 loci, while for the Nichols-clade, the highest variability was found in TP0548 followed by TP0136, TP0462, TP0865 and TP0705 loci. A much higher nucleotide mean diversity was found within the Nichols-clade (π = 0.0121) when compared to SS14-clade strains (π = 0.0002). The TP0136 and TP0548 loci clearly distinguished between the SS14 and Nichols genetic groups with 41 and 30 parsimony-informative sites, respectively. In contrast, the TP0462 and TP0865 loci revealed different haplotypes for the two genetic groups but no parsimony-informative sites distinguishing them. The TP0705 locus shared identical alleles among strains from SS14- and Nichols-clades. A new multilocus sequence typing (MLST) scheme for molecular typing To keep the number of typing loci minimal, we propose MLST of the three loci TP0136, TP0548 and TP0705 as a new typing scheme of TPA. The MLST scheme is based on different allelic profiles (haplotypes) described by a three-letter code, where the first number corresponds to the TP0136 allele, the second to the TP0548 allele and the third to the TP0705 allele (e.g. 1.3.1). MLST and prevalence of macrolide resistance in clinical samples from France and Switzerland We analyzed 120 clinical samples with the new MLST. Clinical characteristics of the patients are depicted in Table 2. Ninety-seven (80.8%) samples were PCR-positive for at least one of the three typing loci, and 54 (45%) of them could be fully typed. We identified eight allele variants for TP0136, 11 for TP0548 and ten for TP0705 (Fig 1B, shown in red). In total, we found 23 different allelic profiles (18 within the SS14-clade, four within the Nichols-clade and one 5 / 16 Fig 1. Allelic variants of all available strain sequences for TP0136, TP0462, TP0548, TP0705 and TP0865. Alleles found in 120 clinical samples examined in this study (samples from France and Switzerland) and in 63 samples published in [ 24, 25, 36 ] (total n = 183) are shown. A. Structural arrangement of candidate loci on the genome. The TP0705 and TP0865 genes lie on the complementary DNA chain. B. Allelic variants. Sequence alignments show only positions containing nucleotide variants. Numbers correspond to nucleotide positions in TPASS_0136, TPASS_0462, TPASS_0548, TPASS_0705, TPASS_0865 genes (TPA SS14; CP004011.1) or to TPANIC_0136, TPANIC_0462, TPANIC_0548, TPANIC_0705, TPANIC_0865 genes (TPA Nichols; CP004010.2). Deletions are shown with dashes or ªdelº. Allele variants that were found among clinical samples in this study (clinical samples from France and Switzerland) are shown in red. For TP0136 and TP0548, the sequence variants were translated to both current molecular typing systems (SBMT and ECDCT) and new sequence variants, which were not identified in previous typing studies, are shown in red. Sequence variant TP0136_11 represent TPA strain Dallas (DAL-1), which contains 30 SNVs and a 58 bp long deletion in TP0136 locus. undetermined) (Table 3). The allelic profiles 1.1.1 and 1.3.1 (pertaining to the SS14-clade) were the most frequent profiles and together accounted for 67% of fully typed samples with 1.3.1 being the only profile detected every year (2011±2015) in all five cities involved in this study. Only 9/97 samples (9.3%) belonged to the Nichols-clade. Apart from the candidate loci, we analyzed the 23S rRNA gene operons to check for macrolide susceptibility/resistance. The operons were successfully amplified in 79 samples. In 63 (79.7%) of them, we either found the A2058G (60 samples, 95.2%) or the A2059G (3 samples, 4.8%) mutation indicating resistance. No strain harbored both mutations. Macrolide resistance was present in both the SS14-clade and Nichols-clade (81% and 60%, respectively) and was relatively high in all tested locations including Lausanne (100%), Paris (85.7%), Lyon (83.3%), Zurich (80%) and Geneva (62.5%). The allelic variant 1.3.1, which is part of the SS14-clade, was associated with macrolide resistance (P = 0.0006). We further search for associations between allelic profiles and clinical characteristics of patients. Besides the association between clade and city (all Nichols-clade strains were found in Zurich), no others were found. 6 / 16 Fig 2. Median joining networks of all available strain sequences for TP0136, TP0462, TP0548, TP0705 and TP0865. Alleles found in 120 clinical samples examined in this study (samples from France and Switzerland) and in 63 samples published in [ 24, 25, 36 ] (total n = 183) are shown. Median joining networks per locus show the different allelic variants and the number of mutational differences among them. Number of mutations, when above one, is given close to branches. Red branch indicates the mutation that did not result in amino acid replacements. Inferred allelic variants (median vectors) are shown as black connecting dots. If contiguous, indels were considered as a single event only. Yellow circles represent the SS14-clade and brown circles represents the Nichols-clade. The number found among the 183 samples is shown inside the circles. Allelic variant numbering as proposed for the MLST scheme are shown in the grey fields. PLOS ONE | https://doi.org/10.1371/journal.pone.0200773 7 / 16 aInformation on clinical characteristics was unavailable for 27 clinical samples. bThere is no information about sexual orientation for 18 samples. cThere is no information about HIV status for 17 samples. dThere is no information about material for 11 samples. Comparison of MLST with other typing schemes While MLST revealed 23 different allelic profiles among 97 typeable clinical samples (Table 3), SBMT revealed 19 different genotypes among 99 typeable samples, of which 12 have not been observed in previous studies (S5 Table). With the ECDCT scheme, 89 samples could be typed revealing 15 different subtypes (S6 Table). When comparing only fully typed samples, we identified 15 allelic profiles, 11 genotypes and eight subtypes with MLST, SBMT and ECDCT, respectively. The ECDCT subtype 14d/g was further divided into five allelic profiles with MLST. Based on arp repeat sizes and tpr restriction fragment length polymorphisms, the MLST allelic profiles 1.1.1 and 1.3.1 were divided into two and four different ECDCT subtypes, respectively. The three typing schemes MLST, ECDCT and SBMT share an 83 bp overlapping fragment of the TP0548, enabling comparative analyses of samples typed across the globe with these schemes. As illustrated in Fig 3, a phylogenetic tree revealed that most of the Swiss and French samples grouped together with seven out of 19 previously identified subtypes. Interestingly, two Swiss samples and two French samples produced two new subtypes within the SS14 clade, respectively, not observed in other studies. These two new subtypes are defined as MLST TP0548 alleles 5 and 11 (and designated as Y and Z in the ECDCT nomenclature). Typing stability and epidemiological concordance In order to determine the typing stability of MLST, we analyzed eleven passages of the TPA reference strain DAL-1 that was continuously propagated in rabbits during 142 days (representing 114 TPA generations). All passages revealed the same allelic profile. In addition, samples of five epidemiologically related patients were typed with MLST; we did not find any discrepancies in allelic profiles. 8 / 16 aAllelic profiles are based on a three-number code: the first number corresponds to the allelic variant in TP0136 locus, the second corresponds to the allelic variant in the TP0548 locus and the third corresponds to the allelic variant in the TP0705 locus. We used X to denote that the allelic variant was not determined. bThe SNVs that determine the allelic variant are shown in the Fig 1. c SÐno mutation in 23S rDNA (sensitive), R8 ±A2058G mutation in 23S rDNA (resistance), R9 ±A2059G mutation in 23S rDNA (resistance). When sensitive and resistant cases were found for one profile, the frequency is given in parenthesis. d NA, not available. e For partially typed samples the no. of samples is not shown. Discussion Syphilis is an important public health problem and an increasing incidence has been noted in recent years. Efficient and affordable typing methods that can be performed in any laboratory around the world with access to sequencing facilities are warranted. These methods would provide information on the frequency of strain types and changes throughout time, on whether particular strain types are associated with specific patient groups, and on patterns among antibiotic resistant strains. With the introduction of TPA enrichment and WGS directly from clinical samples in 2016 [ 24, 25 ] the number of available TPA whole genomes increased, opening an opportunity to identify new genes suitable for typing. We analyzed 30 whole or draft genome TPA sequences to propose a new, efficient and easy-to-perform molecular typing scheme. Taking into account that TPA is an uncultured pathogen and that clinical samples contain low quantities of pathogen DNA, our aim was to keep the number of typing loci minimal. At the same time, considering that TPA is classified 9 / 16 Fig 3. Phylogeny for sequences obtained in typing studies. Maximum likelihood tree produced in MEGA 6 for the 83 bp overlapping fragment of the TP0548 gene typed with MLST, ECDCT and SBMT. The letters A-S stand for the classification in ECDCT and SBMT typing studies. Rabbit passed strains BAL9 was included as representative for type ªHº, which was not found among clinical samples yet. For subtypeºLº the number of identified samples is not known [ 37 ]. SubtypeºJº does not correspond to TPA but was found in T. pallidum subsp. endemicum and T. pallidum subsp. 10 / 16 pertenue [ 35, 38 ]. As observed, TP0548 does not distinguish TPA from TEN or TPE, reflecting potential recombination events between these three. All samples were obtained directly from patients except the BAL9 (passed through rabbits) and isolates from Tanzania (isolated from flies). The numbers in the brackets represent the number of samples identified in the given country. The numbers in red represent the Swiss and French samples from this study. as a monomorphic bacterium and each SNV may be informative, we targeted loci with the highest SNV density to ensure sufficient variability to capture overall population diversity. Our results reveal three variable genes suitable for typing (TP0136, TP0548 and TP0705) which, together with the 23S rRNA genes, are here proposed as a new molecular typing method for syphilis. The proposed TPA MLST shares two typing loci (TP0136, and TP0548) with a recently published MLST for a close relative of TPAÐTreponema pallidum subsp. pertenue [ 38 ], as well as with the previously published TPA SBMT [ 12, 21 ]. Thus, MLST can be considered as enhanced SBMT with additional analysis of TP0705. This latter gene codes for a penicillin-binding protein. TP0136 and TP0548 code for outer membrane proteins, with TP0136 products binding to host fibronectin in the extracellular matrix and playing an important role in local replication [34]. While both TP0136 and TP0548 have been identified as putative recombinant loci, the clustering patterns within TPA recapitulate the distinction between the SS14 and Nichols clade observed with genome-wide data, and are thus congruent with vertical relationships among TPA strains. This congruence may result from recombination not occurring within TPA, but between TPA and Treponema pallidum subsp. endemicum (TEN) or Treponema pallidum subsp. pertenue (TPE) (for TP0548 see Fig 3) [ 24, 35 ]. Nonetheless, the proposed MLST is based on allelic profiles in order to account for potential horizontal gene exchange. As with other MLST schemes, the allelic profiles can easily be converted to the corresponding nucleotide sequence data. The new MLST reveals high resolution power (more than 30% of the resolution achieved with genome-wide data), allowing the distinction of TPA from other treponemal subspecies (S1 Fig), the distinction between the two TPA clades SS14 and Nichols (TP0136, TP0548), and the differentiation of strains within each of these clades (TP0548 for the Nichols clade, and TP0705 for the SS14 clade). When designing a new MLST it is important to find a balance between increasing discriminatory power and limiting the number of samples from which one can get complete typing data. The amplification efficiency of typing schemes is dependent on several factors, e.g. type of material, treatment history of patients, time between sampling and DNA isolation, DNA extraction method, length of the amplification product, and amplification protocol. In order to optimize the resolution power, we opted for relatively long loci that resulted in an overall moderate rate of 45% of fully typed samples. This number is comparable to the amplification efficiency of recently published TP typing studies (ranging from 15.3% to 97% with a median of 46.9% [19, 23, 39±43]. During the course of our study, we further optimized the PCR protocol (see Supplementary Methods) to increase the amplification efficiency (data not shown). Unfortunately, we did not have enough material left in many of the samples, thus, we did not have the opportunity to apply this new enhanced protocol to the whole set of samples. However, when the enhanced protocol was applied to a subset of 30 samples, the amplification efficiency reached 77%. For this reason, we are expecting higher amplification efficiency in the future. The CDCT system has revealed different arp and tpr subtypes in parallel samples of whole blood and genital smears taken from the same patient [ 20 ] indicating inconsistent typing stability. Furthermore, CDCT subtypes, which are based on sequence length differences, do not provide phylogenetic information that could be used to track syphilis transmission dynamics. 11 / 16 The application of the different typing schemes to the Swiss and French dataset showed that MLST yields the highest discriminatory power, with 23 different allelic profiles among 97 typeable samples. Eighteen belonged to the SS14-clade, four to the Nichols-clade and one could not be assigned because in this case only TP0705 data were available. Only nine samples (9.3%) were part of the Nichols-clade, all of which were detected in Zurich. The relative frequency of Nichols and SS14 strains is consistent with that found in other studies [reviewed in 19, 24]. However, a much higher prevalence of samples belonging to the Nichols-clade was reported in Peru (21.4%) [ 39 ], Taiwan (20.2%) [ 40 ], Madagascar (100%) [ 13 ] and Argentina (26.8%) [ 19 ]. We found an alarmingly high percentage of macrolide resistance strains (80% of samples), consistent with other studies performed in Australia [ 41 ], Cuba [ 42 ], China [ 43 ], USA [ 14 ] and Europe including the Czech Republic [ 22 ], Ireland [ 44 ], England [ 45 ], Belgium [ 23 ], Portugal [ 25 ], the Netherlands and Austria [ 24 ]. Accordingly, azithromycin should no longer be recommended as second line therapy for patients with penicillin allergy or bleeding disorders [ 46 ]. Of the two dominant allelic profiles found in this study, 1.3.1 was associated with macrolide resistance. Only a limited number of associations between genetic variants and patient characteristics have been identified so far. A study by Marra and colleagues published in 2010 [ 13 ] found that 50% of patients in Seattle with ECDC subtype 14d/f had neurosyphilis. Other studies have found correlations between genotypes and patient's age, origin, MSM-status and serofast status [ 19, 22, 47 ]. However, we did not observe any associations for the Swiss and French samples other than the location of Nichols clade samples, which were all collected in Zurich. In this study, we present MLST as a new molecular typing scheme for syphilis based on the sequencing of three loci (TP0136, TP0548 and TP0705), and the additional analysis of the 23S rDNA genes to determine macrolide resistance/sensitivity. As demonstrated in this study, the novel MLST scheme yields higher discriminatory power compared to previous typing schemes. We encourage the scientific community and public health authorities to use this MLST method, which can be easily and quickly performed in standard laboratories We expect this tool to open new opportunities in epidemiology allowing longitudinal studies of TPA allelic profiles in different locations, better tracking of syphilis infections, and finding associations of particular strains with specific patient groups, all in favor of a better understanding of the epidemiology of syphilis. Supporting information S1 Supplementary Methods. (DOCX) S1 Table. Draft and complete genome sequences (n = 30) used for identification of candidate loci. (DOCX) S2 Table. Genes with the highest SNV density among samples listed in S1 Table (n = 30). Intergenic regions, paralogous genes (including tpr genes), genes with repetitions and sites with ambiguous data were not included in this analysis. Only genes containing at least four variable sites resulting in the SNV density > 0.001 are shown. Candidate loci are shown in bold. (DOCX) S3 Table. Primers used for the nested-PCR amplification for the candidate loci and the 23S rRNA gene. (DOCX) 12 / 16 S4 Table. GenBank accession numbers corresponding to the particular allelic variants. (DOCX) S5 Table. Nineteen different genotypes found among 99 typeable Swiss and French clinical samples by Sequencing-based molecular typing (SBMT). (DOCX) S6 Table. Fifteen different subtypes found among 89 typeable Swiss and French clinical samples by ECDCT. (DOCX) S1 Fig. Phylogeny of MLST allelic profiles obtained from complete genomes representing different Treponema pallidum subspecies. Maximum likelihood tree produced in MEGA 6 for concatenated sequences of typing loci (TP_0136, TP_0548, TP_0705) in available complete genomes of reference strains representing different Treponema pallidum subspecies (TPA± Treponema pallidum subsp. pallidum; TPEÐTreponema pallidum subsp. pertenue; TENÐTreponema pallidum subsp. endemicum). (TIF) Acknowledgments We thank Nada Juricevic (University Hospital Zurich) for technical assistance as well as Jacques Schrenzel, Gisela Getaz-Jimenez and Manuela Tangomo (University of Geneva) and HeÂlène Salord, FreÂdeÂric Laurent and Fatima Oria (Hospices civils de Lyon) for the performance of DNA isolation and routine real-time PCR on 60 samples included in current study and recruitment of patients. Author Contributions Conceptualization: Linda GrillovaÂ, Fernando GonzaÂlez-Candelas, Homayoun C. Bagheri, David SÏmajs, Natasha Arora, Philipp P. Bosshard. Formal analysis: Linda GrillovaÂ, Kay Nieselt, Marta Pla-DÂõaz, Fernando GonzaÂlez-Candelas. Funding acquisition: David SÏmajs, Natasha Arora. Investigation: Linda GrillovaÂ, Tanika Bawa, Lenka MikalovaÂ. Methodology: Linda GrillovaÂ. Project administration: Michael KruÈtzen. Resources: Angèle Gayet-Ageron, Michal Strouhal, Patrice Sednaoui, Tristan Ferry, Matthias Cavassini, Stephan Lautenschlager, Fabrizio Dutly. Supervision: Linda GrillovaÂ, Philipp P. Bosshard. Writing ± original draft: Linda GrillovaÂ. Writing ± review & editing: Linda GrillovaÂ, David SÏmajs, Natasha Arora, Philipp P. Bosshard. 13 / 16 14 / 16 15 / 16 multicenter study. J Clin Microbiol. 2012; 50(7): 2299±2304. https://doi.org/10.1128/JCM.00341-12 PMID: 22518868 1. World Health Organization. WHO guidelines for the treatment of Treponema pallidum (syphilis ). 2016 ; http://www.who.int/reproductivehealth/publications/rtis/syphilis-treatment-guidelines/en/. 2. Cameron CE , Lukehart SA . Current status of syphilis vaccine development: need, challenges, prospects . Vaccine . 2014 ; 32 : 1602 ± 1609 . https://doi.org/10.1016/j.vaccine. 2013 . 09 .053 PMID: 24135571 3. Grange PA , Mikalova L , Gaudin C , Strouhal M , Janier M , Benhaddou N , et al. Treponema pallidum 11qj Subtype May Correspond to a Treponema pallidum Subsp . Endemicum Strain. Sex Transm Dis . 2016 ; 43 : 517 ± 518 . https://doi.org/10.1097/OLQ.0000000000000474 PMID: 27419817 4. Noda AA , Grillova L , Lienhard R , Blanco O , RodrÂõguez I , SÏ majs D. Bejel in Cuba: molecular identification of Treponema pallidum subsp. endemicum in patients diagnosed with venereal syphilis . Clin Microbiol Infect. Forthcoming 2018 . 5. Fraser CM , Norris SJ , Weinstock GM , White O , Sutton GG , Dodson R , et al. Complete genome sequence of Treponema pallidum, the syphilis spirochete . Science . 1998 ; 281 ( 5375 ): 375 ± 388 . PMID: 9665876 6. Giacani L , Jeffrey BM , Molini BJ , Le HT , Lukehart SA , Centurion-Lara A , et al. Complete genome sequence and annotation of the Treponema pallidum subsp. pallidum Chicago strain . J Bacteriol . 2010 ; 192 ( 10 ): 2645 ± 2646 . https://doi.org/10.1128/JB.00159-10 PMID: 20348263 7. PětrosÏova H , ZobanÂõkova M , Čejkova D , Mikalova L , PospÂõsÏilova P , Strouhal M , et al. Whole genome sequence of Treponema pallidum ssp. pallidum, strain Mexico A, suggests recombination between yaws and syphilis strains . PLoS Negl Trop Dis . 2012 ; 6 ( 9 ):e1832. https://doi.org/10.1371/journal.pntd. 0001832 PMID: 23029591 8. ZobanÂõkova M , Mikolka P , Čejkova D , PospÂõsÏilova P , Chen L , Strouhal M , et al. Complete genome sequence of Treponema pallidum strain DAL-1 . Stand Genomic Sci. 2012 ; 7 ( 1 ): 12 ± 21 . https://doi.org/ 10.4056/sigs.2615838 PMID: 23449808 9. PětrosÏova H , PospÂõsÏilova P , Strouhal M , Čejkova D , ZobanÂõkova M , Mikalova L , et al. Resequencing of Treponema pallidum ssp. pallidum strains Nichols and SS14: correction of sequencing errors resulted in increased separation of syphilis treponeme subclusters . PLoS ONE . 2013 ; 8 ( 9 ):e74319. https://doi. org/10.1371/journal.pone. 0074319 PMID: 24058545 10. Giacani L , Iverson-Cabral SL , King JCK , Molini BJ , Lukehart SA , Centurion-Lara A . Complete Genome Sequence of the Treponema pallidum subsp . pallidum Sea81-4 Strain. Genome Announc . 2014 ; 2 ( 2 ). 11. Pillay A , Liu H , Chen CY , Holloway B , Sturm AW , Steiner B , et al. Molecular subtyping of Treponema pallidum subspecies pallidum . Sex Transm Dis . 1998 ; 25 : 408 ± 414 . PMID: 9773432 12. Flasarova M , SÏmajs D , Matějkova P , Woznicova V , HeroldovaÂ-DvořaÂkova M , Votava M. Molecular detection and subtyping of Treponema pallidum subsp. pallidum in clinical specimens . Epidemiol Mikrobiol Imunol . 2006 ; 55 : 105 ± 111 . PMID: 16970074 13. Marra CM , Sahi SK , Tantalo LC , Godornes C , Reid T , Behets F , et al. Enhanced Molecular Typing of Treponema pallidum: Geographical Distribution of Strain Types and Association with Neurosyphilis . J Infect Dis . 2010 ; 202 : 1380 ± 1388 . https://doi.org/10.1086/656533 PMID: 20868271 14. Katz KA , Pillay A , Ahrens K , Kohn RP , Hermanstyne K , Bernstein KT , et al. Molecular epidemiology of syphilisÐSan Francisco , 2004 ± 2007 . Sex Transm Dis. 2010 ; 37 : 660 ± 663 . https://doi.org/10.1097/ OLQ.0b013e3181e1a77a PMID: 20601928 15. Stamm LV , Bergen HL . A point mutation associated with bacterial macrolide resistance is present in both 23S rRNA genes of an erythromycin-resistant Treponema pallidum clinical isolate . Antimicrob Agents Chemother . 2000 ; 44 ( 3 ): 806 ± 807 . PMID: 10755994 16. Matějkova P , Flasarova M , ZaÂkoucka H , Borek M , Kremenova S , Arenberger P , et al. Macrolide treatment failure in a case of secondary syphilis: a novel A2059G mutation in the 23S rRNA gene of Treponema pallidum subsp . pallidum. J Med Microbiol . 2009 ; 58 ( 6 ): 832 ± 836 . 17. Molini BJ , Tantalo LC , Sahi SK , Rodriguez VI , Brandt SL , Fernandez MC , et al. Macrolide Resistance in Treponema pallidum Correlates With 23S rDNA Mutations in Recently Isolated Clinical Strains . Sex Transm Dis . 2016 ; 43 : 579 ± 583 . https://doi.org/10.1097/OLQ.0000000000000486 PMID: 27513385 18. Peng RR , Wang AL , Li J , Tucker JD , Yin YP , Chen XS . Molecular typing of Treponema pallidum: a systematic review and meta-analysis . PLoS Negl Trop Dis . 2011 ; 5 ( 11 ):e1273. https://doi.org/10.1371/ journal.pntd. 0001273 PMID: 22087340 19. Gallo Vaulet L , Grillova L , Mikalova L , Casco R , Fermepin MR , Pando MA , et al. Molecular typing of Treponema pallidum isolates from Buenos Aires, Argentina: frequent Nichols-like isolates and low levels of macrolide resistance . PLOS ONE . 2017 ; 12 ( 2 ): e0172905. https://doi.org/10.1371/journal.pone. 0172905 PMID: 28235102 20. Mikalova L , PospÂõsÏilova P , Woznicova V , Kuklova I , ZaÂkoucka H , SÏmajs D. Comparison of CDC and sequence-based molecular typing of syphilis treponemes: tpr and arp loci are variable in multiple samples from the same patient . BMC Microbiol . 2013 ; 13 : 178 . https://doi.org/10.1186/ 1471 -2180-13-178 PMID: 23898829 21. Flasarova M , PospÂõsÏilova P , Mikalova L , ValisÏova Z , Dastychova E , Strnadel R , et al. Sequencing based molecular typing of Treponema pallidum strains in the Czech Republic: all identified genotypes are related to the sequence of the SS14 strain . Acta Derm Venereol . 2012 ; 92 : 669 ± 674 . https://doi. org/10.2340/ 00015555 -1335 PMID: 22434073 22. Grillova L , PětrosÏova H , Mikalova L , Strnadel R , Dastychova E , Kuklova I , et al. Molecular typing of Treponema pallidum in the Czech Republic during 2011 to 2013: increased prevalence of identified genotypes and of isolates with macrolide resistance . J Clin Microbiol 2014 ; 52 : 3693 ± 3700 . https://doi. org/10.1128/JCM.01292-14 PMID: 25100820 23. Mikalova L , Grillova L , Osbak K , Strouhal M , Kenyon C , Crucitti T , et al. Molecular typing of syphiliscausing strains among HIV-positive patients in Antwerp, Belgium . Sex Transm Dis . 2017 ; 44 ( 6 ): 376 ± 379 . https://doi.org/10.1097/OLQ.0000000000000600 PMID: 28499290 24. Arora N , Schuenemann VJ , JaÈger G , Peltzer A , Seitz A , Herbig A , et al. Origin of modern syphilis and emergence of a pandemic Treponema pallidum cluster . Nat Microbiol . 2016 ; 2 : 16245 . https://doi.org/ 10.1038/nmicrobiol. 2016 .245 PMID: 27918528 25. Pinto M , Borges V , Antelo M , Pinheiro M , Nunes A , Azevedo J , et al. Genome-scale analysis of the noncultivable Treponema pallidum reveals extensive within-patient genetic variation . Nat Microbiol . 2016 ; 2 : 16190 . https://doi.org/10.1038/nmicrobiol. 2016 .190 PMID: 27748767 26. Centurion-Lara A , Giacani L , Godornes C , Molini BJ , Brinck Reid T , Lukehart SA . Fine analysis of genetic diversity of the tpr gene family among treponemal species, subspecies and strains . PLoS Negl Trop Dis . 2013 ; 7: e2222 . https://doi.org/10.1371/journal.pntd. 0002222 PMID: 23696912 27. Tamura K , Stecher G , Peterson D , Filipski A , Kumar S . MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Mol Biol Evol . 2013 ; 30 : 2725 ± 2729 . https://doi.org/10.1093/molbev/mst197 PMID: 24132122 28. Glatz M , Juricevic N , Altwegg M , Bruisten S , Komericki P , Lautenschlager S , et al. A multicenter prospective trial to asses a new real-time polymerase chain reaction for detection of Treponema pallidum, herpes simplex-1/2 and Haemophilus ducreyi in genital, anal and oropharyngeal ulcers . Clin Microbiol Infect . 2014 ; 20 : 1020 ± 1027 . 29. Gayet-Ageron A , Ninet B , Toutous-Trellu L , Lautenschlager S , Furrer H , Piguet V , et al. Assessment of the real-time PCR test to diagnose syphilis from diverse biological samples . Sex Transm Infect . 2009 ; 85 ( 4 ): 264 ± 269 . https://doi.org/10.1136/sti. 2008 .034314 PMID: 19155240 30. Woznicova V , SÏ majs D , Wechsler D , Matĕjkova P , Flasarova M. Detection of Treponema pallidum subsp. pallidum from skin lesions, serum, and cerebrospinal fluid in an infant with congenital syphilis after clindamycin treatment of the mother during pregnancy . J Clin Microbiol . 2007 ; 45 ( 2 ): 659 ± 661 . https://doi.org/10.1128/JCM.02209-06 PMID: 17151205 31. Lukehart SA , Godornes C , Molini BJ , Sonnett P , Hopkins S , Mulcahy F , et al. Macrolide resistance in Treponema pallidum in the United States and Ireland . N Engl J Med 2004 ; 351 ( 2 ): 154 ± 158 . https://doi. org/10.1056/NEJMoa040216 PMID: 15247355 32. Bandelt HJ , Forster P , RoÈhl A . Median-joining networks for inferring intraspecific phylogenies . Mol Biol Evol . 1999 ; 16 : 37 ± 48 . https://doi.org/10.1093/oxfordjournals.molbev. a026036 PMID: 10331250 33. Vaidya G , Lohman DJ , Meier R. SequenceMatrix: concatenation software for the fast assembly of multigene datasets with character set and codon information . Cladistics 2011 ; 27 ( 2 ): 171 ± 180 . 34. Ke W , Molini BJ , Lukehart SA , Giacani L. Treponema pallidum subsp. pallidum TP0136 protein is heterogeneous among isolates and binds cellular and plasma fibronectin via its NH2-terminal end . PLoS Negl Trop Dis . 2015 ; 9 ( 3 ): e0003662. https://doi.org/10.1371/journal.pntd. 0003662 PMID: 25793702 35. Mikalova L , Strouhal M , Oppelt J , Grange PA , Janier M , Benhaddou N , et al. Human Treponema pallidum 11q/j isolate belongs to subsp. endemicum but contains two loci with a sequence in TP0548 and TP0488 similar to subsp. pertenue and subsp. pallidum, respectively . PLoS Negl Trop Dis . 2017 ; 11 ( 3 ): e0005434. https://doi.org/10.1371/journal.pntd. 0005434 PMID: 28263990 36. Sun J , Meng Z , Wu K , Liu B , Zhang S , Liu Y , et al., Tracing the origin of Treponema pallidum in China using next-generation sequencing . Oncotarget . 2016 ; 7 ( 28 ): 42904 ± 42918 . https://doi.org/10.18632/ oncotarget.10154 PMID: 27344187 37. Tian H , Li Z , Li Z , Hou J , Zheng R , Li F , et al. Molecular typing of Treponema pallidum: identification of a new sequence of TP0548 gene in Shandong, China . Sex. Transm. Dis . 2014 ; 41 ( 9 ): 551 . https://doi. org/10.1097/OLQ.0000000000000155 PMID: 25118968 38. Godornes C , Giacani L , Barry AE , Mitja O , Lukehart SA . Development of a Multilocus Sequence Typing (MLST) scheme for Treponema pallidum subsp. pertenue: Application to yaws in Lihir Island, Papua New Guinea . PLoS Negl Trop Dis . 2017 ; 11 ( 12 ): e0006113. https://doi.org/10.1371/journal.pntd. 0006113 PMID: 29281641 39. Flores JA , Vargas SK , Leon SR , Perez DG , Ramos LB , Chow J , et al. Treponema pallidum pallidum Genotypes and Macrolide Resistance Status in Syphilitic Lesions among Patients at 2 Sexually Transmitted Infection Clinics in Lima, Peru . Sex Transm Dis . 2016 ; 43 ( 7 ): 465 ± 466 . https://doi.org/10.1097/ OLQ.0000000000000465 PMID: 27322050 40. Wu H , Chang SY , Lee NY , Huang WC , Wu BR , Yang CJ , et al. Evaluation of macrolide resistance and enhanced molecular typing of Treponema pallidum in patients with syphilis in Taiwan: a prospective 41. Read P , Tagg KA , Jeoffreys N , Guy RJ , Gilbert GL , Donovan B . Treponema pallidum Strain Types and Association with Macrolide Resistance in Sydney, Australia: New TP0548 Gene Types Identified . J Clin Microbiol . 2016 ; 54 ( 8 ): 2172 ± 2174 . https://doi.org/10.1128/JCM.00959-16 PMID: 27194693 42. Noda AA , Matos N , Blanco O , RodrÂõguez I , Stamm LV . First Report of the 23S rRNA Gene A2058G Point Mutation Associated With Macrolide Resistance in Treponema pallidum From Syphilis Patients in Cuba . Sex Transm Dis . 2016 ; 43 ( 5 ): 332 ± 334 . https://doi.org/10.1097/OLQ.0000000000000440 PMID: 27100771 43. Xiao Y , Liu S , Liu Z , Xie Y , Jiang C , Xu M , et al. Molecular Subtyping and Surveillance of Resistance Genes In Treponema pallidum DNA From Patients With Secondary and Latent Syphilis in Hunan, China . Sex Transm Dis . 2016 ; 43 ( 5 ): 310 ± 316 . https://doi.org/10.1097/OLQ.0000000000000445 PMID: 27100768 44. Muldoon EG , Walsh A , Crowley B , Mulcahy F . Treponema pallidum azithromycin resistance in Dublin, Ireland . Sex Transm Dis . 2012 ; 39 ( 10 ): 784 ± 786 . https://doi.org/10.1097/OLQ.0b013e318269995f PMID: 23001265 45. Tipple C , Taylor GP . Syphilis testing, typing, and treatment follow-up: a new era for an old disease . Curr Opin Infect Dis . 2015 ; 28 ( 1 ): 53 ± 60 . https://doi.org/10.1097/QCO.0000000000000124 PMID: 25485649 46. Janier M , Hegyi V , Dupin N , Unemo M , Tiplica GS , Potočnik M , et al. 2014 European guideline on the management of syphilis . J Eur Acad Dermatol Venereol . 2014 ; 28 ( 12 ): 1581 ± 1593 . https://doi.org/10. 1111/jdv.12734 PMID: 25348878 47. Zhang RL , Wang QQ , Zhang JP , Yang LJ . Molecular typing of Treponema pallidum and associated factors of serofast status in early syphilis patients: Identified novel genotype and cytokine marker . PLoS One 2017 ; 12 ( 4 ): e0175477. https://doi.org/10.1371/journal.pone. 0175477 PMID: 28410389


This is a preview of a remote PDF: http://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0200773&type=printable

Linda Grillová, Tanika Bawa, Lenka Mikalová, Angèle Gayet-Ageron, Kay Nieselt, Michal Strouhal, Patrice Sednaoui, Tristan Ferry, Matthias Cavassini, Stephan Lautenschlager, Fabrizio Dutly, Marta Pla-Díaz, Michael Krützen, Fernando González-Candelas, Homayoun C. Bagheri, David Šmajs, Natasha Arora, Philipp P. Bosshard. Molecular characterization of Treponema pallidum subsp. pallidum in Switzerland and France with a new multilocus sequence typing scheme, PLOS ONE, 2018, DOI: 10.1371/journal.pone.0200773