Complete genome sequences of two strains of Treponema pallidum subsp. pertenue from Indonesia: Modular structure of several treponemal genes

Redefining the treponemal history through pre-Columbian genomes from Brazil

The origins of treponemal diseases have long remained unknown, especially considering the sudden onset of the first syphilis epidemic in the late 15th century in Europe and its hypothesized arrival from the Americas with Columbus' expeditions1,2. Recently, ancient DNA evidence has revealed various treponemal infections circulating in early modern Europe and colonial-era Mexico3-6. However, there has been to our knowledge no genomic evidence of treponematosis recovered from either the Americas or the Old World that can be reliably dated to the time before the first trans-Atlantic contacts. Here, we present treponemal genomes from nearly 2,000-year-old human remains from Brazil. We reconstruct four ancient genomes of a prehistoric treponemal pathogen, most closely related to the bejel-causing agent Treponema pallidum endemicum. Contradicting the modern day geographical niche of bejel in the arid regions of the world, the results call into question the previous palaeopathological characterization of treponeme subspecies and showcase their adaptive potential. A high-coverage genome is used to improve molecular clock date estimations, placing the divergence of modern T. pallidum subspecies firmly in pre-Columbian times. Overall, our study demonstrates the opportunities within archaeogenetics to uncover key events in pathogen evolution and emergence, paving the way to new hypotheses on the origin and spread of treponematoses.

Low genetic diversity of Treponema pallidum ssp. pertenue (TPE) isolated from patients' ulcers in Namatanai District of Papua New Guinea: Local human population is infected by three TPE genotypes

Yaws is an endemic disease caused by Treponema pallidum subsp. pertenue (TPE) that primarily affects children in rural regions of the tropics. The endemic character of yaws infections and the expected exclusive reservoir of TPE in humans opened a new opportunity to start a yaws eradication campaign. We have developed a multi-locus sequence typing (MLST) scheme for TPE isolates combining the previously published (TP0548, TP0488) and new (TP0858) chromosomal loci, and we compared this typing scheme to the two previously published MLST schemes. We applied this scheme to TPE-containing clinical isolates obtained during a mass drug administration study performed in the Namatanai District of Papua New Guinea between June 2018 and December 2019. Of 1081 samples collected, 302 (28.5%) tested positive for TPE DNA, from which 255 (84.4%) were fully typed. The TPE PCR-positivity in swab samples was higher in younger patients, patients with single ulcers, first ulcer episodes, and with ulcer duration less than six months. Non-treponemal serological test positivity correlated better with PCR positivity compared to treponema-specific serological tests. The MLST revealed a low level of genetic diversity among infecting TPE isolates, represented by just three distinct genotypes (JE11, SE22, and TE13). Two previously used typing schemes revealed similar typing resolutions. Two new alleles (one in TP0858 and one in TP0136) were shown to arise by intragenomic recombination/deletion events. Compared to samples genotyped as JE11, the minor genotypes (TE13 and SE22) were more frequently detected in samples from patients with two or more ulcers and patients with higher values of specific TP serological tests. Moreover, the A2058G mutation in the 23S rRNA genes of three JE11 isolates was found, resulting in azithromycin resistance.

The genomes of the yaws bacterium, Treponema pallidum subsp. pertenue, of nonhuman primate and human origin are not genomically distinct

BACKGROUND

Treponema pallidum subsp. pertenue (TPE) is the causative agent of human yaws. Yaws is currently reported in 13 endemic countries in Africa, southern Asia, and the Pacific region. During the mid-20th century, a first yaws eradication effort resulted in a global 95% drop in yaws prevalence. The lack of continued surveillance has led to the resurgence of yaws. The disease was believed to have no animal reservoirs, which supported the development of a currently ongoing second yaws eradication campaign. Concomitantly, genetic evidence started to show that TPE strains naturally infect nonhuman primates (NHPs) in sub-Saharan Africa. In our current study we tested hypothesis that NHP- and human-infecting TPE strains differ in the previously unknown parts of the genomes.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we determined complete (finished) genomes of ten TPE isolates that originated from NHPs and compared them to TPE whole-genome sequences from human yaws patients. We performed an in-depth analysis of TPE genomes to determine if any consistent genomic differences are present between TPE genomes of human and NHP origin. We were able to resolve previously undetermined TPE chromosomal regions (sequencing gaps) that prevented us from making a conclusion regarding the sequence identity of TPE genomes from NHPs and humans. The comparison among finished genome sequences revealed no consistent differences between human and NHP TPE genomes.

CONCLUSION/SIGNIFICANCE

Our data show that NHPs are infected with strains that are not only similar to the strains infecting humans but are genomically indistinguishable from them. Although interspecies transmission in NHPs is a rare event and evidence for current spillover events is missing, the existence of the yaws bacterium in NHPs is demonstrated. While the low risk of spillover supports the current yaws treatment campaign, it is of importance to continue yaws surveillance in areas where NHPs are naturally infected with TPE even if yaws is successfully eliminated in humans.

Whole genome sequences of Treponema pallidum subsp. endemicum isolated from Cuban patients: The non-clonal character of isolates suggests a persistent human infection rather than a single outbreak

Bejel (endemic syphilis) is a neglected non-venereal disease caused by Treponema pallidum subsp. endemicum (TEN). Although it is mostly present in hot, dry climates, a few cases have been found outside of these areas. The aim of this work was the sequencing and analysis of TEN isolates obtained from “syphilis patients” in Cuba, which is not considered an endemic area for bejel. Genomes were obtained by pool segment genome sequencing or direct sequencing methods, and the bioinformatics analysis was performed according to an established pipeline. We obtained four genomes with 100%, 81.7%, 52.6%, and 21.1% breadth of coverage, respectively. The sequenced genomes revealed a non-clonal character, with nucleotide variability ranging between 0.2–10.3 nucleotide substitutions per 100 kbp among the TEN isolates. Nucleotide changes affected 27 genes, and the analysis of the completely sequenced genome also showed a recombination event between tprC and tprI, in TP0488 as well as in the intergenic region between TP0127–TP0129. Despite limitations in the quality of samples affecting breadth of sequencing coverage, the determined non-clonal character of the isolates suggests a persistent infection in the Cuban population rather than a single outbreak caused by imported case.

Development of multiplex HRM-based loop-mediated isothermal amplification method for specific and sensitive detection of Treponema pallidum

Syphilis is a sexually transmitted disease caused by the spirochaete bacterium Treponema pallidum. This study has developed a multiplex High-Resolution Melt-curve Loop-mediated isothermal amplification (multiplex HRM-LAMP) assay targeting the marker genes polA and tprL to detect T. pallidum. The multiplex HRM-LAMP assay conditions were optimized at 65 °C for 45 min. Real-time melt-curve analysis of multiplex HRM-LAMP shows two melt-curve peaks corresponding to polA and tprL with a Tm value of 80 ± 0.5 °C and 87 ± 0.5 °C, respectively. The detection limit of multiplex HRM-LAMP was found to be 6.4 × 10^-4 ng/μL (3.79 copies/μL) of T. pallidum. The specificity was evaluated using seven different bacterial species, and the developed method was 100% specific in detecting T. pallidum. A total of 64 blood samples of T. pallidum suspected cases were used to validate the assay method. The clinical validation showed that the assay was 96.43% sensitive and 100% specific in detecting syphilis. Thus, the developed method was more rapid and sensitive than other available methods and provides a multigene-based diagnostic approach to detect T. pallidum.

Evolutionary processes in the emergence and recent spread of the syphilis agent, Treponema pallidum

The incidence of syphilis has risen worldwide in the last decade in spite of being an easily treated infection. The causative agent of this sexually transmitted disease is the bacterium Treponema pallidum subspecies pallidum (TPA), very closely related to subsp. pertenue (TPE) and endemicum (TEN), responsible for the human treponematoses yaws and bejel, respectively. Although much focus has been placed on the question of the spatial and temporary origins of TPA, the processes driving the evolution and epidemiological spread of TPA since its divergence from TPE and TEN are not well understood. Here, we investigate the effects of recombination and selection as forces of genetic diversity and differentiation acting during the evolution of T. pallidum subspecies. Using a custom-tailored procedure, named phylogenetic incongruence method, with 75 complete genome sequences, we found strong evidence for recombination among the T. pallidum subspecies, involving 12 genes and 21 events. In most cases, only one recombination event per gene was detected and all but one event corresponded to intersubspecies transfers, from TPE/TEN to TPA. We found a clear signal of natural selection acting on the recombinant genes, which is more intense in their recombinant regions. The phylogenetic location of the recombination events detected and the functional role of the genes with signals of positive selection suggest that these evolutionary processes had a key role in the evolution and recent expansion of the syphilis bacteria and significant implications for the selection of vaccine candidates and the design of a broadly protective syphilis vaccine.

Characterisation of the novel clinical isolate X-4 containing a new tp0548 sequence-type

OBJECTIVES

A novel tp0548 sequence-type was identified in one clinical isolate (X-4) from a patient diagnosed with primary syphilis in Xiamen, China. To precisely define and characterise a new clinical isolate, we performed further genome-scale molecular analysis.

METHODS

The pooled segment genome sequencing method followed by Illumina sequencing was performed.

RESULTS

This novel sequence-type contained a unique nucleotide substitution 'T' at position 167 and belonged to the SS14-like clade of TPA strains, as determined by phylogenetic analysis. Multi-locus sequence analysis of nine chromosomal loci demonstrated that the X-4 isolate was clustered within a monophyletic group of TPA strains. Whole-genome phylogenetic analysis subsequently corroborated the TPA strain classification of the X-4 isolate and revealed that the isolate was closely related to the SS14 strain, with 42 single-nucleotide variations and 12 insertions/deletions. In addition, high intrastrain heterogeneity in the length of the poly G/C tracts was found in the TPAChi_0347 locus, which might indicate that this gene of the X-4 isolate is likely involved in phase variation events. The length heterogeneity of the poly A/T tracts was lower than the genetic variability of the poly G/C tracts, and all the observed intrastrain variations fell within coding regions.

CONCLUSION

The novel tp0548 sequence-type was determined to belong to a new TPA isolate, X-4. The identification of variable length in homopolymetic tracts (G/C and A/T) could provide a snapshot of the genes that potentially involved in genotype-phenotype variations. These findings provide an unequivocal characterisation for better understanding the molecular variation of this emerging isolate.

Yaws re-emergence and bacterial drug resistance selection after mass administration of azithromycin: a genomic epidemiology investigation

BACKGROUND

In a longitudinal study assessing the WHO strategy for yaws eradication using mass azithromycin treatment, we observed resurgence of yaws cases with dominance of a single JG8 sequence type and emergence of azithromycin-resistant Treponema pallidum subspecies pertenue (T p pertenue). Here, we analyse genomic changes in the bacterial population using samples collected during the study.

METHODS

We did whole bacterial genome sequencing directly on DNA extracted from 37 skin lesion swabs collected from patients on Lihir Island, Papua New Guinea, between April 1, 2013, and Nov 1, 2016. We produced phylogenies and correlated these with spatiotemporal information to investigate the source of new cases and the emergence of five macrolide-resistant cases. We used deep amplicon sequencing of surveillance samples to assess the presence of minority macrolide-resistant populations.

FINDINGS

We recovered 20 whole T p pertenue genomes, and phylogenetic analysis showed that the re-emerging JG8 sequence type was composed of three bacterial sublineages characterised by distinct spatiotemporal patterns. Of five patients with resistant T p pertenue, all epidemiologically linked, we recovered genomes from three and found no variants. Deep sequencing showed that before treatment, the index patient had fixed macrolide-sensitive T p pertenue, whereas the post-treatment sample had a fixed resistant genotype, as did three of four contact cases.

INTERPRETATION

In this study, re-emergence of yaws cases was polyphyletic, indicating multiple epidemiological sources. However, given the genomic and epidemiological linkage of resistant cases and the rarity of resistance alleles in the general population, azithromycin resistance is likely to have evolved only once in this study, followed by onward dissemination.

FUNDING

Wellcome and Provincial Deputation of Barcelona.

A retrospective study on nested PCR detection of syphilis treponemes in clinical samples: PCR detection contributes to the diagnosis of syphilis in patients with seronegative and serodiscrepant results

Syphilis, caused by Treponema pallidum ssp. pallidum (TPA), is a persisting global health problem. Although syphilis diagnostics relies mainly on serology, serological tests have some limitations, and it is recommended that the final diagnosis be supported by additional tests. The purpose of this study was to analyze the relationship between serology and PCR in syphilis diagnostics. From the year 2004 to May 2019, a total of 941 samples were taken from 833 patients suspected of having syphilis, in Czech Republic. In all these samples, both nested PCR detection of TPA and serology testing were performed. Of the 941 samples, 126 were seronegative, 651 were seropositive, and 164 were serodiscrepant. Among seronegative samples (n = 126), 11 were PCR-positive (8.7%). Among seropositive samples (n = 651; i.e., samples positive for both non-treponemal and treponemal serology tests), 368 samples were PCR-positive (56.5%). The remaining 164 serodiscrepant samples included RPR negative and treponemal serological test-positive samples (n = 154) and a set of 10 RPR-positive samples negative in treponemal serological tests. While the first group revealed 73 PCR-positive samples (47.4%), the second revealed 5 PCR positive samples (50.0%). PCR detection rates were highest in primary syphilis, with lower rates in the secondary and undetermined syphilis stages. As shown here, the nested PCR can improve diagnostics of syphilis, especially in seronegative patients and in patients with discrepant serology.

Ancient Bacterial Genomes Reveal a High Diversity of Treponema pallidum Strains in Early Modern Europe

Syphilis is a globally re-emerging disease, which has marked European history with a devastating epidemic at the end of the 15^th century. Together with non-venereal treponemal diseases, like bejel and yaws, which are found today in subtropical and tropical regions, it currently poses a substantial health threat worldwide. The origins and spread of treponemal diseases remain unresolved, including syphilis' potential introduction into Europe from the Americas. Here, we present the first genetic data from archaeological human remains reflecting a high diversity of Treponema pallidum in early modern Europe. Our study demonstrates that a variety of strains related to both venereal syphilis and yaws-causing T. pallidum subspecies were already present in Northern Europe in the early modern period. We also discovered a previously unknown T. pallidum lineage recovered as a sister group to yaws- and bejel-causing lineages. These findings imply a more complex pattern of geographical distribution and etiology of early treponemal epidemics than previously understood.

Whole genome sequence of the Treponema pallidum subsp. endemicum strain Iraq B: A subpopulation of bejel treponemes contains full-length tprF and tprG genes similar to those present in T. p. subsp. pertenue strains

Treponema pallidum subsp. endemicum (TEN) is the causative agent of endemic syphilis (bejel). Until now, only a single TEN strain, Bosnia A, has been completely sequenced. The only other laboratory TEN strain available, Iraq B, was isolated in Iraq in 1951 by researchers from the US Centers for Disease Control and Prevention. In this study, the complete genome of the Iraq B strain was amplified as overlapping PCR products and sequenced using the pooled segment genome sequencing method and Illumina sequencing. Total average genome sequencing coverage reached 3469×, with a total genome size of 1,137,653 bp. Compared to the genome sequence of Bosnia A, a set of 37 single nucleotide differences, 4 indels, 2 differences in the number of tandem repetitions, and 18 differences in the length of homopolymeric regions were found in the Iraq B genome. Moreover, the tprF and tprG genes that were previously found deleted in the genome of the TEN Bosnia A strain (spanning 2.3 kb in length) were present in a subpopulation of TEN Iraq B and Bosnia A microbes, and their sequence was highly similar to those found in T. p. subsp. pertenue strains, which cause the disease yaws. The genome sequence of TEN Iraq B revealed close genetic relatedness between both available bejel-causing laboratory strains (i.e., Iraq B and Bosnia A) and also genetic variability within the bejel treponemes comparable to that found within yaws- or syphilis-causing strains. In addition, genetic relatedness to TPE strains was demonstrated by the sequence of the tprF and tprG genes found in subpopulations of both TEN Iraq B and Bosnia A. The loss of the tprF and tprG genes in most TEN microbes suggest that TEN genomes have been evolving via the loss of genomic regions, a phenomenon previously found among the treponemes causing both syphilis and rabbit syphilis.

Directly Sequenced Genomes of Contemporary Strains of Syphilis Reveal Recombination-Driven Diversity in Genes Encoding Predicted Surface-Exposed Antigens

Syphilis, caused by Treponema pallidum subsp. pallidum (TPA), remains an important public health problem with an increasing worldwide prevalence. Despite recent advances in in vitro cultivation, genetic variability of this pathogen during infection is poorly understood. Here, we present contemporary and geographically diverse complete treponemal genome sequences isolated directly from patients using a methyl-directed enrichment prior to sequencing. This approach reveals that approximately 50% of the genetic diversity found in TPA is driven by inter- and/or intra-strain recombination events, particularly in strains belonging to one of the defined genetic groups of syphilis treponemes: Nichols-like strains. Recombinant loci were found to encode putative outer-membrane proteins and the recombination variability was almost exclusively found in regions predicted to be at the host-pathogen interface. Genetic recombination has been considered to be a rare event in treponemes, yet our study unexpectedly showed that it occurs at a significant level and may have important impacts in the biology of this pathogen, especially as these events occur primarily in the outer membrane proteins. This study reveals the existence of strains with different repertoires of surface-exposed antigens circulating in the current human population, which should be taken into account during syphilis vaccine development.

Identification of positively selected genes in human pathogenic treponemes: Syphilis-, yaws-, and bejel-causing strains differ in sets of genes showing adaptive evolution

BACKGROUND

Pathogenic treponemes related to Treponema pallidum are both human (causing syphilis, yaws, bejel) and animal pathogens (infections of primates, venereal spirochetosis in rabbits). A set of 11 treponemal genome sequences including those of five Treponema pallidum ssp. pallidum (TPA) strains (Nichols, DAL-1, Mexico A, SS14, Chicago), four T. p. ssp. pertenue (TPE) strains (CDC-2, Gauthier, Samoa D, Fribourg-Blanc), one T. p. ssp. endemicum (TEN) strain (Bosnia A) and one strain (Cuniculi A) of Treponema paraluisleporidarum ecovar Cuniculus (TPeC) were tested for the presence of positively selected genes.

METHODOLOGY/PRINCIPAL FINDINGS

A total of 1068 orthologous genes annotated in all 11 genomes were tested for the presence of positively selected genes using both site and branch-site models with CODEML (PAML package). Subsequent analyses with sequences obtained from 62 treponemal draft genomes were used for the identification of positively selected amino acid positions. Synthetic biotinylated peptides were designed to cover positively selected protein regions and these peptides were tested for reactivity with the patient's syphilis sera. Altogether, 22 positively selected genes were identified in the TP genomes and TPA sets of positively selected genes differed from TPE genes. While genetic variability among TPA strains was predominantly present in a number of genetic loci, genetic variability within TPE and TEN strains was distributed more equally along the chromosome. Several syphilitic sera were shown to react with some peptides derived from the protein sequences evolving under positive selection.

CONCLUSIONS/SIGNIFICANCE

The syphilis-, yaws-, and bejel-causing strains differed relative to sets of positively selected genes. Most of the positively selected chromosomal loci were identified among the TPA treponemes. The local accumulation of genetic variability suggests that the diversification of TPA strains took place predominantly in a limited number of genomic regions compared to the more dispersed genetic diversity differentiating TPE and TEN strains. The identification of positively selected sites in tpr genes and genes encoding outer membrane proteins suggests their role during infection of human and animal hosts. The driving force for adaptive evolution at these loci thus appears to be the host immune response as supported by observed reactivity of syphilitic sera with some peptides derived from protein sequences showing adaptive evolution.

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

A recently introduced Multilocus Sequence Typing scheme for Treponema pallidum subsp. pallidum was applied to clinical samples collected from 2004 to 2017 from the two largest cities (Prague and Brno) in the Czech Republic. Altogether, a total of 675 samples were tested in this study and 281 of them were found PCR-positive for treponemal DNA and typeable. Most of the typed samples (n = 281) were swabs from primary or secondary syphilis lesions (n = 231), and only a minority were whole blood or tissue samples (n = 50). Swab samples from patients with rapid plasma regain (RPR) values of 1–1024 were more frequently PCR-positive (84.6%) compared to samples from patients with non-reactive RPR test (46.5%; p-value = 0.0001). Out of 281 typeable samples, 136 were fully-typed at all TP0136, TP0548, and TP0705 loci. Among the fully and partially typed samples, 25 different allelic profiles were identified. Altogether, eight novel allelic variants were found among fully (n = 5) and partially (n = 3) typed samples. The distribution of TPA allelic profiles identified in the Czech Republic from 2004 to 2017 revealed a dynamic character with allelic profiles disappearing and emerging over time. While the number of samples with the A2058G mutation was seen to increase (86.7% in 2016/2017), the number of samples harboring the A2059G mutation was found to have decreased over time (3.3% in 2016/2017). In addition, we found several allelic profile associations with macrolide resistance or susceptibility, the gender of patients, as well as patient residence.