High frequency of Nichols-like strains and increased levels of macrolide resistance in Treponema pallidum isolates from Buenos Aires, Argentina

doi:10.21203/rs.3.rs-1514085/v1

Globally, 94% of Treponema pallidum subsp. pallidum (TPA) clinical isolates belong to the SS14-like group and 6% to the Nichols-like group, with a prevalence of macrolide resistance of 90%. Our goal was to determine whether local TPA strain distribution and macrolide resistance frequency have changed significantly since our last report, which revealed that Argentina had a high frequency of Nichols-like strains (27%) and low levels of macrolide resistance (14%). Swab samples from patients with suspected syphilis were collected during 2015-2019 and loci TP0136, TP0548, TP0705 were sequenced in order to perform multilocus sequence typing. Strains were classified as Nichols-like or SS14-like. The presence of macrolide resistance-associated mutations was determined by examination of the 23S rDNA gene sequence. Of 46 typable samples, 37.0% were classified as Nichols-like and 63.0% as SS14-like. Macrolide resistance prevalence was 45.7%. Seven allelic profiles were found, five were SS14-like and two were Nichols-like. The frequency of Nichols-like strains increased between studies (26.8% vs. 37.0%, p=0.36). A dramatic increase was found in the frequency of macrolide resistant isolates between studies (14.3% vs. 45.7%, p=0.005). Our results are in agreement with international tendencies and underscore the need to pursue further TPA molecular typing studies in South America.

Treponema pallidum subsp. pallidum (TPA) is the bacterial pathogen that causes syphilis, a chronic venereal disease in humans¹. Despite the availability of several diagnostic tools as well as the effective treatment with penicillin, syphilis prevalence and incidence remain high around the world. The World Health Organization estimated that approximately 22.3 million people globally between the ages of 15 and 49 had syphilis in 2020, with 7.1 million new syphilis cases per year globally, and 1.1 million new cases in the Region of the Americas².

In Argentina, syphilis incidence has been steadily on the rise since 2010, with an increase from 21.2 new cases per 100,000 habitants in 2015 to 56.1 per 100,000 habitants in 2019, and 1.55 cases of congenital syphilis per 1,000 live births in 2019³. Several epidemiological studies have identified groups that are disproportionally vulnerable to syphilis, such as men who have sex with men, female sex workers and male to female transgender people, with syphilis prevalences between 17–50%^4–7.

The molecular characterization of TPA has been of increasing interest in recent years as a useful tool for improving understanding of syphilis network transmission and epidemiology dynamics. Whole genome sequencing (WGS) of clinical isolates from around the world reveals that the global syphilis population comprises just two deeply branching lineages, Nichols and SS14, each with multiple sublineages. The dominant TPA sublineages circulating today are believed to be the result of a rapid population expansion in the 2000s, following a population bottleneck in the 1990s^8,9. The first typing scheme was designed by the Center for Disease Control and Prevention (CDC) in 1998¹⁰ and subsequently improved^11,12. An alternative sequencing-based molecular typing (SBMT) scheme was introduced in 2006¹³ and recently enhanced¹⁴. The latest multilocus sequence typing (MLST) scheme analyzes the TP0136, TP0548 and TP0705 loci and provides a simple alternative to WGS, retaining around 30% of WGS’s resolution power. Additionally, MLST allows for the distinction between the two TPA clades, SS14 and Nichols (based on TP0136 and TP0548 sequences), as well as the differentiation of strains within each clade (based on TP0136, TP0548 and TP0705 sequences). MLST also has the advantage of distinguishing TPA from other treponemes such as Treponema pallidum subsp. pertenue and Treponema pallidum subsp. endemicum¹⁴. MLST of TPA isolates has been expanding since its proposal: the TPA PubMLST database¹⁵ currently contains 1,572 isolate records from 40 different countries. However, over 70% of records come from Europe. Both SBMT and MLST can be paired with the analysis of 23S rDNA to determine the presence of mutations associated with macrolide resistance.

During 2017 we conducted the first characterization of circulating TPA strains in Argentina using SBMT. We detected the presence of nine distinct allelic profiles, five of which had not been described before. Out of these nine genotypes, 77% were related to SS14 and 23% to Nichols strains¹⁶. These results suggested that the syphilis epidemic in Argentina differed considerably from the global distribution pattern, according to which Nichols-like strains are much rarer¹⁷.

Even though penicillin is the antibiotic of choice for the treatment of syphilis, there are other antibiotics which can be prescribed in certain situations (e.g. penicillin allergy). Macrolides in particular are a desirable treatment option because of their oral administration, which makes them more readily acceptable by patients. Two point mutations in the 23S rRNA locus (A2058G and A2059G) have been linked to macrolide resistance in TPA strains¹⁸. Global prevalence of macrolide resistance mutations is 90%, with wide regional variations: Cuba, the US and several European countries report high prevalence rates (61%-100%), whereas Canada, Russia, South Africa, Madagascar, Taiwan and Peru report little to no macrolide resistance (0%-12%). Our first study seemed to indicate that the situation in Argentina was closer to the latter, with 14% prevalence of macrolide resistance mutations¹⁶.

Several studies have shown that both the predominant TPA strains and the frequency of mutations associated with macrolide resistance in a particular population can change over time¹⁹. With the aim of studying the dynamics of TPA isolates in Argentina, the purpose of the present study was to characterize the syphilis epidemic in Buenos Aires by analyzing TPA clinical isolates collected during 2015–2019 in terms of their strain and genotype distribution as well as the presence of antibiotic resistance-associated mutations, in order to update and expand on our previous report using the newer MLST methodology.

Characteristics of patients and clinical samples. A total of 99 swabs were collected from patients with suspected syphilis lesions (Fig. 1). In 28 cases, active syphilis was ruled out due to negative serology (non reactive Venereal Disease Research Laboratory –VDRL– test). DNA isolation was unsuccessful in 14 cases as evidenced by negative beta-actin PCR. Treponemal DNA (polA and/or tmpC) amplification was unsuccessful in 11 cases.

Clinical characteristics of patients with syphilis diagnosis (based on serology and dark-field microscopy –DFM– results) and successful DNA isolation are presented in Table 1. Briefly, 83% of patients were men with a median age of 25 years and 46% of those with known HIV status were HIV-positive. Most patients were Argentine nationals (74%) and unmarried (88%), with approximately half of them residing in the city of Buenos Aires (52%) or the Province of Buenos Aires (48%); 52% were students or unemployed. Lesions were mainly located in the genitals (53%), some were oral (28%) and only one was an anal lesion. Most patients (83%) had positive VDRL serology, 65% with high titres (> 1:32). Although data on syphilis stage could not be recovered for almost half of the patients, among those with a determined clinical stage, 45% were primary syphilis cases and 55% were secondary syphilis cases. No significant differences were found between PCR-positive and PCR-negative patients in terms of sociodemographic or clinical characteristics.

Table 1

Characteristics of 57 syphilis-positive patients with successful DNA isolation enrolled in this study during 2015–2019 in Buenos Aires, Argentina.
Clinical characteristics of patients	All patients (n = 57)^a	PCR-positive patients (n = 46)^a	PCR-negative patients (n = 11)^a
Gender, n (%)	M 47 (82.5); W 10 (17.5)	M 37 (80.4); W 9 (19.6)	M 10 (90.9); W 1 (9.1)
Median age (men/women)	25 (16–69) / 27 (18–35)	25 (16–59) / 26 (18–35)	29 (16–69) / 27
HIV positive^b, n (%)	10/22 (45.5)	9/20 (45.0)	1/2 (50.0)
DFM^c (%)	10 P (17.5); 38 N (66.7); 9 n.d. (15.8)	9 P (19.6); 31 N (67.4); 6 n.d. (13.0)	1 P (9.1); 7 N (63.6); 3 n.d. (27.3)
Serology
FTA/TPHA^d (%)	20 P (35.1); 37 N/n.d. (64.9)	17 P (37.0); 29 N/n.d. (63.0)	3 P (27.3); 8 N/n.d. (72.7)
VDRL (%)	47 P (82.5); 9 N (15.8); 1 n.d. (1.8)	37 P (80.4); 8 N (17.4); 1 n.d. (2.2)	10 P (90.9); 1 N (9.1); 0 n.d. (0)
≥ 1:32	37 (64.9)	29 (63.0)	9 (81.8)
Syphilis stage (%)
Primary	13 (22.8)	12 (26.1)	1 (9.1)
Secondary	16 (28.1)	14 (30.4)	2 (18.2)
Undetermined	28 (49.1)	20 (43.5)	8 (72.7)
Lesion type	Samples (n = 57)	Samples (n = 46)	Samples (n = 11)
Genital	30 (52.6%)	25 (54.3%)	5 (45.5%)
Oral	16 (28.1%)	12 (26.1%)	4 (36.4%)
Anal	1 (1.8%)	1 (2.2%)	0 (0)
Oral-genital	1 (1.8%)	1 (2.2%)	0 (0)
Undetermined	9 (15.8%)	7 (15.2%)	2 (18.2%)

^a M, men; W, women; P, positive; N, negative; n.d., not determined.

^b HIV status was known only in 22 patients (38.6%).

^c DFM, dark-field microscopy.

^d FTA, Fluorescent Treponemal Antibody Absorption; TPHA, Treponema pallidum Particle Hemagglutination Assay.

Multilocus sequence typing (MLST). All 46 PCR-positive samples were typeable using MLST. Samples were considered typeable when at least one sequence of a typing locus (either TP0136, TP0548 or TP0705) was obtained. By comparing the TP0136 and/or TP0548 sequences from our samples with reference sequences we determined that 63.0% (29/46) of isolates belonged to the SS14-like group and 37.0% (17/46) to the Nichols-like group. In all cases where both TP0136 and TP0548 were successfully sequenced (n = 37), clade assignment was consistent between loci.

Samples from 34 patients were completely typed (including the TP0136, TP0548 and TP0705 gene sequences), 29 also including the 23S rRNA gene sequence. Altogether, seven different TPA genotypes were found among the fully typed: five were related to the sequence from the TPA strain SS14 (1.1.1., 1.3.1., 1.1.8., 1.32.1. and 1.48.3.) and two were related to the sequence from the TPA strain Nichols (9.7.3. and 9.20.3.). All identified genotypes and their frequencies are shown in Table 2. Briefly, the most frequent allelic profile was 9.7.3. (35.3% among fully-typed samples) followed by 1.1.1. and 1.3.1. (each 26.5% among fully-typed samples), and finally by 1.1.8., 1.32.1, 1.48.3 and 9.20.3. with only one case each. Among the partially typed samples, only one SS14-like strain showed a genotype that was not consistent with any of the fully typed profiles (X.1.3.) All alleles found in this study have been described previously; however, we newly describe the allelic profile 1.48.3. A phylogenetic analysis of allelic profiles found among fully typed samples is shown in Fig. 2.

Table 2

MLST allelic profiles for 34 fully typed and 12 partially typed lesion swabs samples from patients attending a sexually transmitted infections clinic in Buenos Aires, 2015–2019.
Sequence type^a	Allelic profile^b	Typing	23S rDNA^c (no. of samples)	Genetic group	No. of samples (%)^d
26	9.7.3.	Complete	S(7)/R8(2)/X(3)	Nichols	12 (26.1)
2	1.1.1.	Complete	S(6)/R8(2)/X(1)	SS14	9 (19.6)
1	1.3.1.	Complete	S(1)/R8(8)	SS14	9 (19.6)
31	9.20.3.	Complete	R9 (1)	Nichols	1 (2.2)
3	1.1.8.	Complete	X(1)	SS14	1 (2.2)
50	1.32.1.	Complete	R8(1)	SS14	1 (2.2)
117^e	1.48.3.^e	Complete	S(1)	SS14	1 (2.2)
-	1.X.1.	Partial	R8(1)/X(1)	SS14	2 (4.3)
-	1.X.X.	Partial	X(2)	SS14	2 (4.3)
-	9.7.X.	Partial	S(2)	Nichols	2 (4.3)
-	X.1.1.	Partial	X(1)	SS14	1 (2.2)
-	X.1.3.	Partial	S(1)	SS14	1 (2.2)
-	X.3.1.	Partial	R8(1)	SS14	1 (2.2)
-	X.3.X.	Partial	X(1)	SS14	1 (2.2)
-	9.20.X.	Partial	S(1)	Nichols	1 (2.2)
-	9.X.X.	Partial	X(1)	Nichols	1 (2.2)

^a According to the PubMLST database of Treponema pallidum subsp. pallidum¹⁵.

^b Allelic 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.

^c S, no mutation in 23S rDNA (sensitive); R8, A2058G mutation in 23S rDNA (resistance); R9, A2059G mutation in 23S rDNA (resistance); X, not determined. The frequency (number) of each case is given in parenthesis.

^d Percentages calculated over the total number of typed samples (no. of fully typed samples + no. of partially typed samples).

^e Newly identified allelic profile, with alleles known from previous studies.

Prevalence of mutations associated with macrolide resistance. The 23S rRNA gene locus was amplified in 35 samples, and in 16 of them (45.7%), either the A2058G (15 patients, 93.7%) or the A2059G (1 patient, 6.3%) mutation was found. No strains harbored both A2058G and A2059G mutations. Macrolide resistance was detected in both SS14-like (13/22) and Nichols-like (3/13) strains. The frequency of resistance was higher in the SS14 clade, however without statistical significance (59.1% vs. 23.1%, p = 0.08).

Associations of clade and allelic profiles with patient characteristics. We found no statistically significant associations of clade or allelic profile with patient age, gender, residence, nationality, occupation, marital status, sexually transmitted infection (STI) history, location of lesion, VDRL titer, DFM results, syphilis stage or HIV status. Patients carrying Nichols-like strains tended to be slightly older than patients carrying SS14-like strains (median age: 31 vs. 25 years, p = 0.07) and to have higher VDRL titers (median titer: 128 vs. 32, p = 0.09).

At the local level, this study is the second to analyze the genetic characteristics of circulating TPA in Argentina, serving as both a continuation and an update of the previous report by our research group¹⁶. While the first study used the SBMT methodology, this study is the first to apply the enhanced MLST scheme in Argentina. At the regional level, this study constitutes as well as one of the first initiatives to characterize the diversity of TPA in South America, where only a few molecular epidemiology studies of TPA have been conducted^20,21, and even fewer using the MLST methodology²².

In this study, as in most reports worldwide, the SS14 clade was predominant over the Nichols clade. However, the frequency of Nichols-like strains observed (37%) is considerably higher than is usually the case. Clade SS14 is reported to constitute the vast majority (94.1%) of worldwide TPA clinical isolates, though the relative frequencies of these clades show some geographical variation: while in most countries the prevalence of Nichols-like strains is lower than 10%, this has been found to be 21.4% in Peru, 20.2% in Taiwan and 100% in Madagascar²³. In Argentina, our first study¹⁶ showed a frequency of Nichols-like strains of 26.8%, which would make Argentina the country with the highest reported frequency at the time after Madagascar. In the present study we found an even higher prevalence of the Nichols clade which might indicate that this prevalence is currently on the rise. Moreover, our results, together with those of the report from Peru²¹, suggest that South America is a region with high frequency of Nichols-like circulating strains, though admittedly more data is required to support this claim.

Among fully typed samples, we found seven distinct genotypes, five of which were related to the SS14 strain (1.1.1., 1.3.1., 1.1.8., 1.32.1. and 1.48.3.), while the other two were related to the Nichols strain (9.7.3. and 9.20.3.). Genotypes 9.7.3. (35%), 1.1.1. (27%) and 1.3.1. (27%) accounted for most samples. In our previous study, a different SBMT nomenclature was used and also that study did not include the TP0705 locus, making comparisons difficult. However, a rough comparison clearly reveals that the most common genotypes have changed. In our previous study, the SS14-like strains SSS, SSR9 and SSR8 (now named 1.1.X.) constituted the most prevalent (43.9%), followed by the Nichols-like strains U3U6S and U3U6R8 (now named 9.7.X.) representing 17% of the cases. However, in this new study, strain 9.7.X. (specifically 9.7.3.) constitutes the most prevalent (35.3%) followed by strain 1.1.X. (29.4%, including 1.1.1., 26.5%, and 1.1.8., 2.9%) and strain 1.3.X. (specifically 1.3.1., 26.5%). Thus, although SS14-like strains as a whole remained dominant, the most prevalent allelic profile changed from an SS14-related genotype to a Nichols-related genotype, reflecting the increase in the overall frequency of Nichols strains. The difference in genotype distribution between the Argentinean studies may be due to the fact that samples were collected at different times; several studies show that TPA strains in a particular population can vary over time¹⁹.

A comparison of the frequencies of TPA allelic profiles observed in this study with the average frequencies reported in other continents (Table 3) shows that, not only are the relative frequencies of the SS14 and Nichols clades observed in Argentina unusual, but so are the allelic profiles detected. Genotype 9.7.3., which was found in 35% of our samples, is only present in less than 5% of samples in other continents. Meanwhile genotype 1.3.1., which accounts for roughly half of clinical samples in Europe and 60% of samples in North America, was found in approximately a fourth of Argentinian samples. Genotype 1.1.1., though not rare in other continents, was detected in our samples with a considerably higher frequency than usually reported. These differences underscore the need for further TPA molecular typing studies in South America, since it cannot be assumed that clade and allelic profile frequencies among local clinical samples will be similar to those in other regions.

Table 3

Comparison of reported TPA allelic profile frequencies between this study and studies conducted in other continents.
	1.3.1.	1.1.1.	9.7.3.	1.1.8.	9.20.3.	1.32.1.	1.48.3.
Argentina (this study)	26.5%	26.5%	35.3%	2.9%	2.9%	2.9%	2.9%
Europe^a	48.9%	15.1%	4.4%	7.7%	0.15%	0.15%	0%
North America^a	60%	9.6%	2.2%	0.74%	0%	0%	0%
Asia^a	3.9%	3.9%	0.8%	75.8%	0%	0%	0%
Total^a	44.2%	12.8%	3.5%	15.9%	0.11%	0.11%	0%

^a Data according to the PubMLST database¹⁵. Only the frequencies of fully typed isolates were used for calculations.

The observed prevalence of macrolide resistance-associated mutations in the present study was nearly 50%, the A2058G mutation accounting for all but one of these cases. The presence of this mutation was detected in SS14-like as well as Nichols-like isolates, with a tendency towards higher mutation frequency among SS14-like strains (59% vs. 23%), which was also the case in our first study¹⁶. The frequency of macrolide resistance has tripled between our first study and this study (14.3% vs 45.7%, p < 0.005) and it places Argentina closer to other countries with high prevalence of macrolide resistance such as Australia²⁴, Cuba²⁵, Czech Republic^19,22, China^27–30, Ireland^31,32, the UK³³ and the US^11,34−36. As shown in Table 4, even though the presence of macrolide resistance mutations increased with time in both clades, this increase showed statistical significance only in the SS14 clade. The observed increase in the frequency of macrolide-resistant TPA strains in our study population coincides with the global trend³⁷.

Table 4

Comparison of clade and macrolide resistance frequencies between studies.
	Gallo Vaulet ML, 2017	This study	p-value
Clade SS14 Nichols	73.2% (30/41) 26.8% (11/41)	63.0% (29/46) 37.0% (17/46)	0.36
Macrolide resistance mutations All SS14 Nichols	14.3% (6/42) 16.7% (5/30) 9.1% (1/11)	45.7% (16/35) 59.1% (13/22) 23.1% (3/13)	0.005 0.003 0.60

This rise might reflect a change in the consumption of macrolide antibiotics in Argentina. However, data on antibiotic use in our country is limited: a study on the evolution of the purchase of antibiotic drugs from Argentine pharmacies between 2015 and 2017 found no substantial change in the purchase of macrolides³⁸. Another explanation could be an increased veterinary use of macrolides, as is the case with other pathogens³⁹. More publicly available information on medical and veterinary uses of macrolides in Argentina would be necessary to shed light on the observed increase in macrolide resistance-associated mutations in circulating TPA strains.

A total of 99 swabs were collected from patients with suspected syphilis lesions, 28 of which were ruled out as syphilis-negative. Among the remaining 71 samples, DNA isolation was unsuccessful in 14 cases, while treponemal DNA amplification was unsuccessful in 11 cases. These failures in DNA isolation and amplification may be due to low amounts of material collected during swabbing or to the presence of PCR inhibitors.

In this study, no statistically significant associations were found between clade or allelic profile and patient characteristics. However, patients carrying Nichols-like strains tended to be slightly older than patients carrying SS14-like strains. A significant association between the Nichols clade and older patients has been reported in clinical samples from France⁴⁰. The Nichols clade also tended to be associated with higher VDRL titers, but this association has not, to our knowledge, been found in other studies.

Molecular typing studies of human pathogenic bacteria can help elucidate the evolutionary history and origin of old diseases such as syphilis¹⁷, while also contributing to the understanding of transmission mechanisms, epidemic diversity and dynamics, as well as to the improvement of clinical practice and public policies of control and prevention⁴¹. Particularly in the case of TPA, a bacterium which so far has resisted any attempt at long-term in vitro culture⁴², molecular typing techniques prove especially useful. Global TPA typing data has facilitated the characterization of syphilis outbreaks, the discovery of associations between subtypes and neurosyphilis, the surveillance of macrolide resistance, the distinction between re-infection and reactivation, and the comprehension of the geographical, temporal and population distribution of TPA. However, in spite of all its benefits, only a limited number of studies from a few countries have focused on molecular typing of TPA since the emergence of the first typing assays, especially in South America⁴¹. Future molecular epidemiology studies of syphilis should aim at providing useful information for the development of prevention and control programs, focusing for instance on the identification of at-risk populations and sources of infection, as well as the verification of links between certain subtypes and pathology.

Clinical samples. Clinical samples were collected between May 2015 and November 2019 from patients with suspected syphilis at an STI clinic in the city of Buenos Aires (Programa de Enfermedades de Transmisión Sexual –PETS-, Hospital de Clínicas “José de San Martín”, University of Buenos Aires, Argentina). Although it is possible to detect treponemal DNA in non-lesion samples, lesions tend to have a higher treponemal load and thus usually yield higher DNA amplification efficiencies, therefore swab samples were taken only from lesions⁴³. Swab samples from genital, anal, oral or rectal lesions were taken without cleaning the lesion. The swab samples remained at room temperature for no more than a week before being stored at -20°C until DNA isolation.

Patients suspected of having syphilis based on clinical findings were diagnosed with laboratory tests that included the examination of lesion material by DFM and serological tests. Serology included VDRL, FTA and TPHA. Patients were considered to have syphilis when one of the tests (DFM and/or serology) was positive. Patients gave their written, informed consent and received recommendations of treatment, care, referral to the laboratory for serological diagnosis, and a new visit for follow-up. Clinical data included patients’ age, gender, STI history, location of lesion, results of serology and DFM, syphilis stage and HIV status. Socio-demographic data included residence, nationality, occupation and marital status.

DNA isolation. DNA isolation from swabs was performed using a commercial kit (QIAamp DNA Blood Mini Kit, Qiagen, Germany). First, swabs were submerged in 600 µl of phosphate buffered saline. The eluate was then processed according to the manufacturer’s recommendations. The efficacy of the extraction process was verified by PCR amplification of the human beta-actin gene. The PCR mixture (final volume of 25 µl) contained 16.3 µl of water, 0.5 µl of a 10 mM deoxynucleotide triphosphate (dNTP) mixture, 2.5 µl of 10x buffer, 3 µl of 25 mM MgCl₂, 0.75 µl of each primer (10 µM), 0.25 µl of Taq DNA polymerase (5,000 U/ml), and 1 µl of isolated total DNA. PCR amplification was performed under the following conditions: 95˚C (2 min); 94˚C (30 s), 60˚C (30 s), 72˚C (30 s) for 35 cycles; and 72˚C (7 min). Primers used can be found in Supplementary Table S1 online.

Detection of treponemal DNA. In order to screen for TPA-positive samples, amplification of loci TP0105 (polA) and TP0319 (tmpC) was performed by nested PCR, as described previously²⁶. The PCR mixture (final volume of 25 µl) contained in the first step 15.3 µl of water, 0.1 µl of a 10 mM dNTP mixture, 2.5 µl of 10x buffer, 1.5 µl of 25 mM MgCl₂, 0.25 µl of each primer (10 µM), 0.1 µl of Taq DNA polymerase (5,000 U/ml), and 5 µl of isolated DNA. PCR amplification was performed under the following conditions: 94˚C (1 min); 94˚C (30 s), 58˚C (30 s), 72˚C (1 min) for 30 cycles; and 72˚C (10 min). The mixture for the second step was the same except that it contained 19.3 µl of water and 1 µl of PCR product from the first step. PCR amplification was performed under the following conditions: 94˚C (1 min); 94˚C (30 s), 58˚C (30 s), 72˚C (1 min) for 40 cycles; and 72˚C (10 min). Samples were considered PCR-positive for treponemal DNA when PCR amplification was positive for at least one treponemal locus. Primers used can be found in Supplementary Table S1 online.

PCR amplification for MLST. MLST was performed to type all TPA-positive samples as described previously⁴⁴. Four loci, including TP0136, TP0548, TP0705 and 23S rRNA genes, were amplified using the nested PCR protocol. The PCR mixture (final volume of 25 µl) contained in both steps 20.5 µl of water, 0.5 µl of a 10 mM dNTP mixture, 2.5 µl of 10x buffer, 0.25 µl of each primer (100 µM), 0.05 µl of Taq DNA polymerase (5,000 U/ml), and 1 µl of isolated DNA. In the first step, PCR amplification was performed under the following conditions: 95˚C (1 min); 94˚C (30 s), 50˚C (30 s), 72˚C (1 min, 45 s) for 40 cycles; and 72˚C (7 min). In the second step, PCR amplification was performed under the following conditions: 95˚C (1 min); 94˚C (30 s), 55˚C (30 s), 72˚C (1 min, 45 s) for 40 cycles; and 72˚C (7 min). The second amplification step was repeated an additional three times, in order to quadruple the volume of PCR product per sample prior to DNA purification. PCR products were purified with polyethylene glycol (PEG), using an in-house protocol (see Supplementary Information). Primers used can be found in Supplementary Table S1 online.

DNA sequencing and sequence analysis. Sequencing reactions were performed on the purified PCR products using a commercial kit (BigDye™ Terminator v3.1 Cycle Sequencing Kit, Applied Biosystems, USA). The products were sequenced on an automated capillary DNA sequencing system. Sequence analyses were performed using Sequencher software (Sequencher® Version 5.4.6 DNA Sequence Analysis Software, Gene Codes Corporation, Ann Arbor, MI, USA). TP0136 and TP0548 sequences were aligned to reference sequences TPA Nichols (CP004010.2) and TPA SS14 (CP004011.1) in order to assign samples to the corresponding clade (clade assignment was possible where at least one of these loci was successfully sequenced). 23S rDNA sequences were evaluated at positions 2058 and 2059 in the 23S rDNA gene of Escherichia coli(accession no. V00331), where A→G mutations are associated with macrolide resistance. Alleles encoding resistance were marked A2058G or A2059G depending on the site of substitution.

Phylogenetic analyses. Clade assignment was confirmed through phylogenetic analysis: the software MEGA X⁴⁵ was used to construct multiple sequence alignments for each gene separately via the ClustalW algorithm⁴⁶, as well as for a concatenation of TP0136, TP0548 and TP0705 sequences, which were then used for the construction of phylogenetic trees, using the maximum likelihood method⁴⁷ with bootstrapping⁴⁸ (Tamura-Nei model). MLST of T. pallidum is based on the different allelic profiles (haplotypes) described by a 3-letter code, where the first number corresponds to the TP0136, the second to the TP0548 allele, and the third to the TP0705 allele (e.g. 1.3.1.).

Statistical analyses. Comparisons were performed using Fisher's exact test, for categorical variables, and the Mann–Whitney U test, for numerical variables, via the SPSS software (IBM Corp. Released 2011. IBM SPSS Statistics for Windows, Version 20.0. Armonk, NY: IBM Corp).

Ethics statement. This study was approved by the Ethics Committee of the Hospital de Clínicas “José de San Martín” (University of Buenos Aires). All patients gave their written, informed consent. The study was conducted according to the principles expressed in the Declaration of Helsinki⁴⁹, as well as national regulations regarding research.

Acknowledgments

We would like to thank all participants who made this study possible. This work was supported through a grant provided by the Universidad de Buenos Aires (UBACYT 20720130200010BA) awarded to Pando MA. The funding source had no role in the design or carrying out of the study, data collection, analysis, or interpretation, the writing of the manuscript, or the decision to submit it for publication.

Author Contributions

D.S. and M.A.P. conceived the study. A.M. enrolled participants. N.M. performed experiments. N.M., E.V., D.S., R.D.R. and M.A.P. analysed the data. N.M. and M.A.P. drafted the report. All authors read the manuscript, provided input into the report and approved the final version.

Data availability statement

All alleles found in this study were identical to previously reported sequences, which can be found in the PubMLST database¹⁵.

Competing Interests

The authors declare no competing interests.

Radolf, J. D. et al. Treponema pallidum, the syphilis spirochete: making a living as a stealth pathogen. Nat. Rev. Microbiol. 14, 744–759 (2016).
World Health Organization. Global progress report on HIV, viral hepatitis and sexually transmitted infections, 2021: accountability for the global health sector strategies 2016–2021: actions for impact: web annex 2: data methods. https://apps.who.int/iris/bitstream/handle/10665/342813/9789240030992-eng.pdf?sequence=1&isAllowed=y (2021).
Dirección de Sida y ETS & Ministerio de Salud de la Nación. Boletín N° 38: Respuesta al VIH y las ITS en la Argentina. https://bancos.salud.gob.ar/recurso/boletin-ndeg-38-respuesta-al-vih-y-las-its-en-la-argentina (2021).
Pando, M. A. et al. Prevalence of HIV and other sexually transmitted infections among female commercial sex workers in Argentina. Am. J. Trop. Med. Hyg. 74, 233–238 (2006).
Pando, M. A. et al. HIV and Other Sexually Transmitted Infections among Men Who Have Sex with Men Recruited by RDS in Buenos Aires, Argentina: High HIV and HPV Infection. PLoS One 7, e39834 (2012).
Pando, M. A. et al. Violence as a Barrier for HIV Prevention among Female Sex Workers in Argentina. PLoS One 8, e54147 (2013).
dos Ramos Farías, M. S. et al. First report on sexually transmitted infections among trans (male to female transvestites, transsexuals, or transgender) and male sex workers in Argentina: high HIV, HPV, HBV, and syphilis prevalence. Int. J. Infect. Dis. 15, e635–e640 (2011).
Beale, M. A. et al. Global phylogeny of Treponema pallidum lineages reveals recent expansion and spread of contemporary syphilis. Nat. Microbiol. 6, 1549–1560 (2021).
Pla-Díaz, M. et al. Evolutionary Processes in the Emergence and Recent Spread of the Syphilis Agent, Treponema pallidum. Mol. Biol. Evol. 39, msab318 (2022).
Pillay, A. et al. Molecular Subtyping of Treponema pallidum Subspecies pallidum. Sex. Transm. Dis. 25, 408–414 (1998).
Katz, K. A. et al. Molecular Epidemiology of Syphilis—San Francisco, 2004-2007. Sex. Transm. Dis. 37, 660–663 (2010).
Marra, C. M. et al. Enhanced Molecular Typing of Treponema pallidum: Geographical Distribution of Strain Types and Association with Neurosyphilis. J. Infect. Dis. 202, 1380–1388 (2010).
Flasarová, M. et al. Molekulární detekce a typizace Treponema pallidum subsp. pallidum v klinickém materiálu. Epidemiol. Mikrobiol. Imunol. 55, 105–111 (2006).
Grillová, L. et al. Molecular characterization of Treponema pallidum subsp. pallidum in Switzerland and France with a new multilocus sequence typing scheme. PLoS One 13, e0200773 (2018).
Grillová, L., Jolley, K., Šmajs, D. & Picardeau, M. A public database for the new MLST scheme for Treponema pallidum subsp. pallidum: surveillance and epidemiology of the causative agent of syphilis. PeerJ 6, e6182 (2019).
Gallo Vaulet, L. et al. Molecular typing of Treponema pallidum isolates from Buenos Aires, Argentina: Frequent Nichols-like isolates and low levels of macrolide resistance. PLoS One 12, e0172905 (2017).
Arora, N. et al. Origin of modern syphilis and emergence of a pandemic Treponema pallidum cluster. Nat. Microbiol. 2, 16245 (2017).
Stamm, L. V. Syphilis: antibiotic treatment and resistance. Epidemiol. Infect. 143, 1567–1574 (2015).
Grillová, L. 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. 52, 3693–3700 (2014).
Cruz, A. R. et al. Secondary Syphilis in Cali, Colombia: New Concepts in Disease Pathogenesis. PLoS Negl. Trop. Dis. 4, e690 (2010).
Flores, J. A. 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. 43, 465–466 (2016).
Lieberman, N. A. P. et al. Treponema pallidum genome sequencing from six continents reveals variability in vaccine candidate genes and dominance of Nichols clade strains in Madagascar. PLoS Negl. Trop. Dis. 15, e0010063 (2021).
Šmajs, D., Strouhal, M. & Knauf, S. Genetics of human and animal uncultivable treponemal pathogens. Infect. Genet. Evol. 61, 92–107 (2018).
Read, P. et al. Treponema pallidum Strain Types and Association with Macrolide Resistance in Sydney, Australia: New TP0548 Gene Types Identified. J. Clin. Microbiol. 54, 2172–2174 (2016).
Noda, A. A., Matos, N., Blanco, O., Rodríguez, I. & Stamm, L. V. 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. 43, 332–334 (2016).
Flasarová, M. 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. Venerol. 92, 669–674 (2012).
Xiao, Y. 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. 43, 310–316 (2016).
Martin, I. E., Gu, W., Yang, Y. & Tsang, R. S. W. Macrolide Resistance and Molecular Types of Treponema pallidum Causing Primary Syphilis in Shanghai, China. Clin. Infect. Dis. 49, 515–521 (2009).
Li, Z. et al. Two Mutations Associated with Macrolide Resistance in Treponema pallidum in Shandong, China. J. Clin. Microbiol. 51, 4270–4271 (2013).
Chen, X.-S. et al. High prevalence of azithromycin resistance to Treponema pallidum in geographically different areas in China. Clin. Microbiol. Infect. 19, 975–979 (2013).
Lukehart, S. A. et al. Macrolide Resistance in Treponema pallidum in the United States and Ireland. N. Engl. J. Med. 351, 154–158 (2004).
Muldoon, E. G., Walsh, A., Crowley, B. & Mulcahy, F. Treponema pallidum Azithromycin Resistance in Dublin, Ireland. Sex. Transm. Dis. 39, 784–786 (2012).
Tipple, C., McClure, M. O. & Taylor, G. P. High prevalence of macrolide resistant Treponema pallidum strains in a London centre. Sex. Transm. Infect. 87, 486–488 (2011).
Grimes, M. et al. Two Mutations Associated with Macrolide Resistance in Treponema pallidum: Increasing Prevalence and Correlation with Molecular Strain Type in Seattle, Washington. Sex. Transm. Dis. 39, 954–958 (2012).
Marra, C. M. et al. Antibiotic Selection May Contribute to Increases in Macrolide‐Resistant Treponema pallidum. J. Infect. Dis. 194, 1771–1773 (2006).
Su, J. R. Prevalence of the 23S rRNA A2058G Point Mutation and Molecular Subtypes in Treponema pallidum in the United States, 2007 to 2009. Sex. Transm. Dis. 39, 794–798 (2012).
Šmajs, D., Grillová, L. & Paštěková, L. Macrolide Resistance in the Syphilis Spirochete, Treponema pallidum ssp. pallidum: Can We Also Expect Macrolide-Resistant Yaws Strains? Am. J. Trop. Med. Hyg. 93, 678–683 (2015).
Maito, M. A. Use of commercialized antibiotics in Argentine Republic pharmacies (2015-2017). Revista Ciencia Reguladora 3, 30-34 (2018).
Kenyon, C. Positive association between the use of macrolides in food-producing animals and pneumococcal macrolide resistance: a global ecological analysis. Int. J. Infect. Dis. 116, 344–347 (2022).
Pospíšilová, P. et al. Multi-locus sequence typing of Treponema pallidum subsp. pallidum present in clinical samples from France: Infecting treponemes are genetically diverse and belong to 18 allelic profiles. PLoS One 13, e0201068 (2018).
Peng, R.-R. et al. Molecular Typing of Treponema pallidum: A Systematic Review and Meta-Analysis. PLoS Negl. Trop. Dis. 5, e1273 (2011).
Edmondson, D. G., Hu, B. & Norris, S. J. Long-Term In Vitro Culture of the Syphilis Spirochete Treponema pallidum subsp. pallidum. mBio 9, e01153-18 (2018).
Towns, J. M. et al. Treponema pallidum detection in lesion and non-lesion sites in men who have sex with men with early syphilis: a prospective, cross-sectional study. Lancet Infect. Dis. 21, 1324–1331 (2021).
Vrbová, E. et al. MLST typing of Treponema pallidum subsp. pallidum in the Czech Republic during 2004-2017: Clinical isolates belonged to 25 allelic profiles and harbored 8 novel allelic variants. PLoS One 14, e0217611 (2019).
Kumar, S., Stecher, G., Li, M., Knyaz, C. & Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol. Biol. Evol. 35, 1547–1549 (2018).
Thompson, J. D., Higgins, D. G. & Gibson, T. J. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22, 4673–4680 (1994).
Tamura, K., Nei, M. & Kumar, S. Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc. Natl. Acad. Sci. U.S.A. 101, 11030–11035 (2004).
Felsenstein, J. Confidence limits on phylogenies: An Approach using the bootstrap. Evolution 39, 783–791 (1985).
World Medical Association. World Medical Association Declaration of Helsinki: Ethical Principles for Medical Research Involving Human Subjects. JAMA 310, 2191–2194 (2013).

No competing interests reported.

SupplementaryInformation.pdf

High frequency of Nichols-like strains and increased levels of macrolide resistance in Treponema pallidum isolates from Buenos Aires, Argentina

Status:

Version 1

Abstract

Figures

Introduction

Results

Discussion

Methods

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1