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% of broad 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 broad 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. Author summary The incidence of venereal syphilis has greatly increased in the last years, however endemic syphilis (bejel) which have been considered as a disease restricted to dry arid areas such as the Sahel and the Middle East, remain as a neglected disease. In Cuba, which is a tropical country, several bejel cases were unexpectedly detected few years ago in “syphilis” patients with no records of travel abroad or sex with foreign partners. In this study, we explored the whole genome sequences from four of the Cuban Treponema pallidum subsp. endemicum (TEN) isolates and the substantial genetic diversity detected among them suggests a persistent infection of TEN within the human population rather than a single outbreak of a TEN isolate introduced from an area where it is typically endemic. This finding has significant implications on this neglected and also possibly tropical disease in terms of geographical/temporal distribution, and highlights the importance of keeping in mind neglected diseases in apparently non-endemic areas.

[1]  Kelika A. Konda,et al.  Treponema pallidum genome sequencing from six continents reveals variability in vaccine candidate genes and dominance of Nichols clade strains in Madagascar , 2021, PLoS neglected tropical diseases.

[2]  D. Šmajs,et al.  Genetic recombination in Treponema pallidum: Implications for diagnosis, epidemiology, and vaccine development. , 2021, Sexually Transmitted Diseases.

[3]  M. Strouhal,et al.  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 , 2020, PloS one.

[4]  T. Kawahata,et al.  Bejel, a Nonvenereal Treponematosis, among Men Who Have Sex with Men, Japan , 2019, Emerging infectious diseases.

[5]  A. Mechaly,et al.  Directly Sequenced Genomes of Contemporary Strains of Syphilis Reveal Recombination-Driven Diversity in Genes Encoding Predicted Surface-Exposed Antigens , 2019, Front. Microbiol..

[6]  O. Blanco,et al.  Multilocus Sequence Typing of Treponema pallidum subsp. pallidum in Cuba From 2012 to 2017 , 2018, The Journal of infectious diseases.

[7]  O. Blanco,et al.  Bejel in Cuba: molecular identification of Treponema pallidum subsp. endemicum in patients diagnosed with venereal syphilis. , 2018, Clinical microbiology and infection : the official publication of the European Society of Clinical Microbiology and Infectious Diseases.

[8]  S. Knauf,et al.  Complete genome sequences of two strains of Treponema pallidum subsp. pertenue from Indonesia: Modular structure of several treponemal genes , 2018, PLoS neglected tropical diseases.

[9]  G. Russo,et al.  Sequencing of Treponema pallidum subsp. pallidum from isolate UZ1974 using Anti-Treponemal Antibodies Enrichment: First complete whole genome sequence obtained directly from human clinical material , 2018, PloS one.

[10]  A. Gayet-Ageron,et al.  Molecular characterization of Treponema pallidum subsp. pallidum in Switzerland and France with a new multilocus sequence typing scheme , 2018, PloS one.

[11]  S. Knauf,et al.  Genetics of human and animal uncultivable treponemal pathogens. , 2018, Infection, genetics and evolution : journal of molecular epidemiology and evolutionary genetics in infectious diseases.

[12]  S. Bruisten,et al.  Complete genome sequences of two strains of Treponema pallidum subsp. pertenue from Ghana, Africa: Identical genome sequences in samples isolated more than 7 years apart , 2017, PLoS neglected tropical diseases.

[13]  N. Dupin,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 , 2017, PLoS neglected tropical diseases.

[14]  J. Krause,et al.  Origin of modern syphilis and emergence of a pandemic Treponema pallidum cluster , 2016, Nature Microbiology.

[15]  Sudhir Kumar,et al.  MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. , 2016, Molecular biology and evolution.

[16]  D. Raoult,et al.  Incidental Syphilis Diagnosed by Real-Time PCR Screening of Urine Samples , 2015, Journal of Clinical Microbiology.

[17]  L. Fulton,et al.  Whole Genome Sequence of the Treponema pallidum subsp. endemicum Strain Bosnia A: The Genome Is Related to Yaws Treponemes but Contains Few Loci Similar to Syphilis Treponemes , 2014, PLoS neglected tropical diseases.

[18]  L. Pham,et al.  Selective Microbial Genomic DNA Isolation Using Restriction Endonucleases , 2014, PloS one.

[19]  Heng Li,et al.  Toward better understanding of artifacts in variant calling from high-coverage samples , 2014, Bioinform..

[20]  L. Giacani,et al.  The Endemic Treponematoses , 2014, Clinical Microbiology Reviews.

[21]  Yunlong Liu,et al.  NGSUtils: a software suite for analyzing and manipulating next-generation sequencing datasets , 2013, Bioinform..

[22]  R. Tsang,et al.  Local Transmission of Imported Endemic Syphilis, Canada, 2011 , 2012, Emerging infectious diseases.

[23]  G. Weinstock,et al.  Genetic diversity in Treponema pallidum: implications for pathogenesis, evolution and molecular diagnostics of syphilis and yaws. , 2012, Infection, genetics and evolution : journal of molecular epidemiology and evolutionary genetics in infectious diseases.

[24]  L. Fulton,et al.  Whole Genome Sequences of Three Treponema pallidum ssp. pertenue Strains: Yaws and Syphilis Treponemes Differ in Less than 0.2% of the Genome Sequence , 2012, PLoS neglected tropical diseases.

[25]  Marcel Martin Cutadapt removes adapter sequences from high-throughput sequencing reads , 2011 .

[26]  M. Golden,et al.  Enhanced molecular typing of treponema pallidum: geographical distribution of strain types and association with neurosyphilis. , 2010, The Journal of infectious diseases.

[27]  M. DePristo,et al.  The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. , 2010, Genome research.

[28]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[29]  M. Nei,et al.  Prospects for inferring very large phylogenies by using the neighbor-joining method. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[30]  B. Molini,et al.  Gene conversion: a mechanism for generation of heterogeneity in the tprK gene of Treponema pallidum during infection , 2004, Molecular microbiology.

[31]  B. Spratt Exploring the concept of clonality in bacteria. , 2004, Methods in molecular biology.

[32]  G. Weinstock,et al.  BAC library of T. pallidum DNA in E. coli. , 2002, Genome research.

[33]  G. Weinstock,et al.  From microbial genome sequence to applications. , 2000, Research in microbiology.

[34]  M. Nei,et al.  Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. , 1993, Molecular biology and evolution.

[35]  J. Felsenstein CONFIDENCE LIMITS ON PHYLOGENIES: AN APPROACH USING THE BOOTSTRAP , 1985, Evolution; international journal of organic evolution.

[36]  E. Grin Endemic syphilis in Bosnia; clinical and epidemiological observations on a successful mass-treatment campaign. , 1952, Bulletin of the World Health Organization.