Characterization of the core and accessory genomes of Pseudomonas aeruginosa using bioinformatic tools Spine and AGEnt

Background

Pseudomonas aeruginosa is an important opportunistic pathogen responsible for many infections in hospitalized and immunocompromised patients. Previous reports estimated that approximately 10% of its 6.6 Mbp genome varies from strain to strain and is therefore referred to as “accessory genome”. Elements within the accessory genome of P. aeruginosa have been associated with differences in virulence and antibiotic resistance. As whole genome sequencing of bacterial strains becomes more widespread and cost-effective, methods to quickly and reliably identify accessory genomic elements in newly sequenced P. aeruginosa genomes will be needed.

Results

We developed a bioinformatic method for identifying the accessory genome of P. aeruginosa. First, the core genome was determined based on sequence conserved among the completed genomes of twelve reference strains using Spine, a software program developed for this purpose. The core genome was 5.84 Mbp in size and contained 5,316 coding sequences. We then developed an in silico genome subtraction program named AGEnt to filter out core genomic sequences from P. aeruginosa whole genomes to identify accessory genomic sequences of these reference strains. This analysis determined that the accessory genome of P. aeruginosa ranged from 6.9-18.0% of the total genome, was enriched for genes associated with mobile elements, and was comprised of a majority of genes with unknown or unclear function. Using these genomes, we showed that AGEnt performed well compared to other publically available programs designed to detect accessory genomic elements. We then demonstrated the utility of the AGEnt program by applying it to the draft genomes of two previously unsequenced P. aeruginosa strains, PA99 and PA103.

Conclusions

The P. aeruginosa genome is rich in accessory genetic material. The AGEnt program accurately identified the accessory genomes of newly sequenced P. aeruginosa strains, even when draft genomes were used. As P. aeruginosa genomes become available at an increasingly rapid pace, this program will be useful in cataloging the expanding accessory genome of this bacterium and in discerning correlations between phenotype and accessory genome makeup. The combination of Spine and AGEnt should be useful in defining the accessory genomes of other bacterial species as well.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-737) contains supplementary material, which is available to authorized users.     

Core and accessory genome architecture in a group of Pseudomonas aeruginosa Mu-like phages

Background

Bacteriophages that infect the opportunistic pathogen Pseudomonas aeruginosa have been classified into several groups. One of them, which includes temperate phage particles with icosahedral heads and long flexible tails, bears genomes whose architecture and replication mechanism, but not their nucleotide sequences, are like those of coliphage Mu. By comparing the genomic sequences of this group of P. aeruginosa phages one could draw conclusions about their ontogeny and evolution.

Results

Two newly isolated Mu-like phages of P. aeruginosa are described and their genomes sequenced and compared with those available in the public data banks. The genome sequences of the two phages are similar to each other and to those of a group of P. aeruginosa transposable phages. Comparing twelve of these genomes revealed a common genomic architecture in the group. Each phage genome had numerous genes with homologues in all the other genomes and a set of variable genes specific for each genome. The first group, which comprised most of the genes with assigned functions, was named “core genome”, and the second group, containing mostly short ORFs without assigned functions was called “accessory genome”. Like in other phage groups, variable genes are confined to specific regions in the genome.

Conclusion

Based on the known and inferred functions for some of the variable genes of the phages analyzed here, they appear to confer selective advantages for the phage survival under particular host conditions. We speculate that phages have developed a mechanism for horizontally acquiring genes to incorporate them at specific loci in the genome that help phage adaptation to the selective pressures imposed by the host.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1146) contains supplementary material, which is available to authorized users.     

Complete genome sequence of Pseudomonas aeruginosa lytic bacteriophage PA1O which resembles temperate bacteriophage D3112

A novel Pseudomonas aeruginosa lytic bacteriophage (phage), PA1Ø, was isolated, and its genome was sequenced completely. This phage is able to lyse not only P. aeruginosa but also Staphylococcus aureus. Genome analysis of PA1Ø showed that it is similar to a P. aeruginosa temperate phage, D3112, with the exception of the absence of a c repressor-encoding gene, which is known to play a critical role in the maintenance of the lysogenic state of D3112 in P. aeruginosa.     

The Electrophysiological Underpinnings of Processing Gender Stereotypes in Language

Despite the widely documented influence of gender stereotypes on social behaviour, little is known about the electrophysiological substrates engaged in the processing of such information when conveyed by language. Using event-related brain potentials (ERPs), we examined the brain response to third-person pronouns (lei “she” and lui “he”) that were implicitly primed by definitional (passeggera _FEM “passenger”, pensionato _MASC “pensioner”), or stereotypical antecedents (insegnante “teacher”, conducente “driver”). An N400-like effect on the pronoun emerged when it was preceded by a definitionally incongruent prime (passeggera _FEM – lui; pensionato _MASC – lei), and a stereotypically incongruent prime for masculine pronouns only (insegnante – lui). In addition, a P300-like effect was found when the pronoun was preceded by definitionally incongruent primes. However, this effect was observed for female, but not male participants. Overall, these results provide further evidence for on-line effects of stereotypical gender in language comprehension. Importantly, our results also suggest a gender stereotype asymmetry in that male and female stereotypes affected the processing of pronouns differently.     

Highly selective end-tagged antimicrobial peptides derived from PRELP

Background

Antimicrobial peptides (AMPs) are receiving increasing attention due to resistance development against conventional antibiotics. Pseudomonas aeruginosa and Staphylococcus aureus are two major pathogens involved in an array of infections such as ocular infections, cystic fibrosis, wound and post-surgery infections, and sepsis. The goal of the study was to design novel AMPs against these pathogens.

Methodology and Principal Findings

Antibacterial activity was determined by radial diffusion, viable count, and minimal inhibitory concentration assays, while toxicity was evaluated by hemolysis and effects on human epithelial cells. Liposome and fluorescence studies provided mechanistic information. Protease sensitivity was evaluated after subjection to human leukocyte elastase, staphylococcal aureolysin and V8 proteinase, as well as P. aeruginosa elastase. Highly active peptides were evaluated in ex vivo skin infection models. C-terminal end-tagging by W and F amino acid residues increased antimicrobial potency of the peptide sequences GRRPRPRPRP and RRPRPRPRP, derived from proline arginine-rich and leucine-rich repeat protein (PRELP). The optimized peptides were antimicrobial against a range of Gram-positive S. aureus and Gram-negative P. aeruginosa clinical isolates, also in the presence of human plasma and blood. Simultaneously, they showed low toxicity against mammalian cells. Particularly W-tagged peptides displayed stability against P. aeruginosa elastase, and S. aureus V8 proteinase and aureolysin, and the peptide RRPRPRPRPWWWW-NH₂ was effective against various “superbugs” including vancomycin-resistant enterococci, multi-drug resistant P. aeruginosa, and methicillin-resistant S. aureus, as well as demonstrated efficiency in an ex vivo skin wound model of S. aureus and P. aeruginosa infection.

Conclusions/Significance

Hydrophobic C-terminal end-tagging of the cationic sequence RRPRPRPRP generates highly selective AMPs with potent activity against multiresistant bacteria and efficiency in ex vivo wound infection models. A precise “tuning” of toxicity and proteolytic stability may be achieved by changing tag-length and adding W- or F-amino acid tags.     

Chemical analysis of cellular and extracellular carbohydrates of a biofilm-forming strain Pseudomonas aeruginosa PA14

Background

Pseudomonas aeruginosa is a Gram-negative bacterium and an opportunistic pathogen, which causes persisting life-threatening infections in cystic fibrosis (CF) patients. Biofilm mode of growth facilitates its survival in a variety of environments. Most P. aeruginosa isolates, including the non-mucoid laboratory strain PA14, are able to form a thick pellicle, which results in a surface-associated biofilm at the air-liquid (A–L) interface in standing liquid cultures. Exopolysaccharides (EPS) are considered as key components in the formation of this biofilm pellicle. In the non-mucoid P. aeruginosa strain PA14, the “scaffolding” polysaccharides of the biofilm matrix, and the molecules responsible for the structural integrity of rigid A–L biofilm have not been identified. Moreover, the role of LPS in this process is unclear, and the chemical structure of the LPS O-antigen of PA14 has not yet been elucidated.

Principal Findings

In the present work we carried out a systematic analysis of cellular and extracellular (EC) carbohydrates of P. aeruginosa PA14. We also elucidated the chemical structure of the LPS O-antigen by chemical methods and 2-D NMR spectroscopy. Our results showed that it is composed of linear trisaccharide repeating units, identical to those described for P. aeruginosa Lanýi type O:2a,c (Lanýi-Bergman O-serogroup 10a, 10c; IATS serotype 19) and having the following structure: -4)-α-L-GalNAcA-(1–3)-α-D-QuiNAc-(1–3)- α-L-Rha-(1-. Furthermore, an EC O-antigen polysaccharide (EC O-PS) and the glycerol-phosphorylated cyclic β-(1,3)-glucans were identified in the culture supernatant of PA14, grown statically in minimal medium. Finally, the extracellular matrix of the thick biofilm formed at the A-L interface contained, in addition to eDNA, important quantities (at least ∼20% of dry weight) of LPS-like material.

Conclusions

We characterized the chemical structure of the LPS O-antigen and showed that the O-antigen polysaccharide is an abundant extracellular carbohydrate of PA14. We present evidence that LPS-like material is found as a component of a biofilm matrix of P. aeruginosa.     

The Susceptibility of Pseudomonas aeruginosa Strains from Cystic Fibrosis Patients to Bacteriophages

Phage therapy may become a complement to antibiotics in the treatment of chronic Pseudomonas aeruginosa infection. To design efficient therapeutic cocktails, the genetic diversity of the species and the spectrum of susceptibility to bacteriophages must be investigated. Bacterial strains showing high levels of phage resistance need to be identified in order to decipher the underlying mechanisms. Here we have selected genetically diverse P. aeruginosa strains from cystic fibrosis patients and tested their susceptibility to a large collection of phages. Based on plaque morphology and restriction profiles, six different phages were purified from “pyophage”, a commercial cocktail directed against five different bacterial species, including P. aeruginosa. Characterization of these phages by electron microscopy and sequencing of genome fragments showed that they belong to 4 different genera. Among 47 P. aeruginosa strains, 13 were not lysed by any of the isolated phages individually or by pyophage. We isolated two new phages that could lyse some of these strains, and their genomes were sequenced. The presence/absence of a CRISPR-Cas system (Clustered Regularly Interspaced Short Palindromic Repeats and Crisper associated genes) was investigated to evaluate the role of the system in phage resistance. Altogether, the results show that some P. aeruginosa strains cannot support the growth of any of the tested phages belonging to 5 different genera, and suggest that the CRISPR-Cas system is not a major defence mechanism against these lytic phages.     

Non Mycobacterial Virulence Genes in the Genome of the Emerging Pathogen Mycobacterium abscessus

Mycobacterium abscessus is an emerging rapidly growing mycobacterium (RGM) causing a pseudotuberculous lung disease to which patients with cystic fibrosis (CF) are particularly susceptible. We report here its complete genome sequence. The genome of M. abscessus (CIP 104536T) consists of a 5,067,172-bp circular chromosome including 4920 predicted coding sequences (CDS), an 81-kb full-length prophage and 5 IS elements, and a 23-kb mercury resistance plasmid almost identical to pMM23 from Mycobacterium marinum. The chromosome encodes many virulence proteins and virulence protein families absent or present in only small numbers in the model RGM species Mycobacterium smegmatis. Many of these proteins are encoded by genes belonging to a “mycobacterial” gene pool (e.g. PE and PPE proteins, MCE and YrbE proteins, lipoprotein LpqH precursors). However, many others (e.g. phospholipase C, MgtC, MsrA, ABC Fe(3+) transporter) appear to have been horizontally acquired from distantly related environmental bacteria with a high G+C content, mostly actinobacteria (e.g. Rhodococcus sp., Streptomyces sp.) and pseudomonads. We also identified several metabolic regions acquired from actinobacteria and pseudomonads (relating to phenazine biosynthesis, homogentisate catabolism, phenylacetic acid degradation, DNA degradation) not present in the M. smegmatis genome. Many of the “non mycobacterial” factors detected in M. abscessus are also present in two of the pathogens most frequently isolated from CF patients, Pseudomonas aeruginosa and Burkholderia cepacia. This study elucidates the genetic basis of the unique pathogenicity of M. abscessus among RGM, and raises the question of similar mechanisms of pathogenicity shared by unrelated organisms in CF patients.     

Genes but not genomes reveal bacterial domestication of Lactococcus lactis

Background

The population structure and diversity of Lactococcus lactis subsp. lactis, a major industrial bacterium involved in milk fermentation, was determined at both gene and genome level. Seventy-six lactococcal isolates of various origins were studied by different genotyping methods and thirty-six strains displaying unique macrorestriction fingerprints were analyzed by a new multilocus sequence typing (MLST) scheme. This gene-based analysis was compared to genomic characteristics determined by pulsed-field gel electrophoresis (PFGE).

Methodology/Principal Findings

The MLST analysis revealed that L. lactis subsp. lactis is essentially clonal with infrequent intra- and intergenic recombination; also, despite its taxonomical classification as a subspecies, it displays a genetic diversity as substantial as that within several other bacterial species. Genome-based analysis revealed a genome size variability of 20%, a value typical of bacteria inhabiting different ecological niches, and that suggests a large pan-genome for this subspecies. However, the genomic characteristics (macrorestriction pattern, genome or chromosome size, plasmid content) did not correlate to the MLST-based phylogeny, with strains from the same sequence type (ST) differing by up to 230 kb in genome size.

Conclusion/Significance

The gene-based phylogeny was not fully consistent with the traditional classification into dairy and non-dairy strains but supported a new classification based on ecological separation between “environmental” strains, the main contributors to the genetic diversity within the subspecies, and “domesticated” strains, subject to recent genetic bottlenecks. Comparison between gene- and genome-based analyses revealed little relationship between core and dispensable genome phylogenies, indicating that clonal diversification and phenotypic variability of the “domesticated” strains essentially arose through substantial genomic flux within the dispensable genome.     

Trehalose Biosynthesis Promotes Pseudomonas aeruginosa Pathogenicity in Plants

Pseudomonas aeruginosa strain PA14 is a multi-host pathogen that infects plants, nematodes, insects, and vertebrates. Many PA14 factors are required for virulence in more than one of these hosts. Noting that plants have a fundamentally different cellular architecture from animals, we sought to identify PA14 factors that are specifically required for plant pathogenesis. We show that synthesis by PA14 of the disaccharide trehalose is required for pathogenesis in Arabidopsis, but not in nematodes, insects, or mice. In-frame deletion of two closely-linked predicted trehalose biosynthetic operons, treYZ and treS, decreased growth in Arabidopsis leaves about 50 fold. Exogenously co-inoculated trehalose, ammonium, or nitrate, but not glucose, sulfate, or phosphate suppressed the phenotype of the double ΔtreYZΔtreS mutant. Exogenous trehalose or ammonium nitrate does not suppress the growth defect of the double ΔtreYZΔtreS mutant by suppressing the plant defense response. Trehalose also does not function intracellularly in P. aeruginosa to ameliorate a variety of stresses, but most likely functions extracellularly, because wild-type PA14 rescued the in vivo growth defect of the ΔtreYZΔtreS in trans. Surprisingly, the growth defect of the double ΔtreYZΔtreS double mutant was suppressed by various Arabidopsis cell wall mutants that affect xyloglucan synthesis, including an xxt1xxt2 double mutant that completely lacks xyloglucan, even though xyloglucan mutants are not more susceptible to pathogens and respond like wild-type plants to immune elicitors. An explanation of our data is that trehalose functions to promote the acquisition of nitrogen-containing nutrients in a process that involves the xyloglucan component of the plant cell wall, thereby allowing P. aeruginosa to replicate in the intercellular spaces in a leaf. This work shows how P. aeruginosa, a multi-host opportunistic pathogen, has repurposed a highly conserved “house-keeping” anabolic pathway (trehalose biosynthesis) as a potent virulence factor that allows it to replicate in the intercellular environment of a leaf.     

The Periplasmic Protein TolB as a Potential Drug Target in Pseudomonas aeruginosa

The Gram-negative bacterium Pseudomonas aeruginosa is one of the most dreaded pathogens in the hospital setting, and represents a prototype of multi-drug resistant “superbug” for which effective therapeutic options are very limited. The identification and characterization of new cellular functions that are essential for P. aeruginosa viability and/or virulence could drive the development of anti-Pseudomonas compounds with novel mechanisms of action. In this study we investigated whether TolB, the periplasmic component of the Tol-Pal trans-envelope protein complex of Gram-negative bacteria, represents a potential drug target in P. aeruginosa. By combining conditional mutagenesis with the analysis of specific pathogenicity-related phenotypes, we demonstrated that TolB is essential for P. aeruginosa growth, both in laboratory and clinical strains, and that TolB-depleted P. aeruginosa cells are strongly defective in cell-envelope integrity, resistance to human serum and several antibiotics, as well as in the ability to cause infection and persist in an insect model of P. aeruginosa infection. The essentiality of TolB for P. aeruginosa growth, resistance and pathogenicity highlights the potential of TolB as a novel molecular target for anti-P. aeruginosa drug discovery.     

The leukemia-specific fusion gene ETV6/RUNX1 perturbs distinct key biological functions primarily by gene repression

Background

ETV6/RUNX1 (E/R) (also known as TEL/AML1) is the most frequent gene fusion in childhood acute lymphoblastic leukemia (ALL) and also most likely the crucial factor for disease initiation; its role in leukemia propagation and maintenance, however, remains largely elusive. To address this issue we performed a shRNA-mediated knock-down (KD) of the E/R fusion gene and investigated the ensuing consequences on genome-wide gene expression patterns and deducible regulatory functions in two E/R-positive leukemic cell lines.

Findings

Microarray analyses identified 777 genes whose expression was substantially altered. Although approximately equal proportions were either up- (KD-UP) or down-regulated (KD-DOWN), the effects on biological processes and pathways differed considerably. The E/R KD-UP set was significantly enriched for genes included in the “cell activation”, “immune response”, “apoptosis”, “signal transduction” and “development and differentiation” categories, whereas in the E/R KD-DOWN set only the “PI3K/AKT/mTOR signaling” and “hematopoietic stem cells” categories became evident. Comparable expression signatures obtained from primary E/R-positive ALL samples underline the relevance of these pathways and molecular functions. We also validated six differentially expressed genes representing the categories “stem cell properties”, “B-cell differentiation”, “immune response”, “cell adhesion” and “DNA damage” with RT-qPCR.

Conclusion

Our analyses provide the first preliminary evidence that the continuous expression of the E/R fusion gene interferes with key regulatory functions that shape the biology of this leukemia subtype. E/R may thus indeed constitute the essential driving force for the propagation and maintenance of the leukemic process irrespective of potential consequences of associated secondary changes. Finally, these findings may also provide a valuable source of potentially attractive therapeutic targets.     

Genomic signatures of strain selection and enhancement in Bacillus atrophaeus var. globigii, a historical biowarfare simulant

Background

Despite the decades-long use of Bacillus atrophaeus var. globigii (BG) as a simulant for biological warfare (BW) agents, knowledge of its genome composition is limited. Furthermore, the ability to differentiate signatures of deliberate adaptation and selection from natural variation is lacking for most bacterial agents. We characterized a lineage of BGwith a long history of use as a simulant for BW operations, focusing on classical bacteriological markers, metabolic profiling and whole-genome shotgun sequencing (WGS).

Results

Archival strains and two “present day” type strains were compared to simulant strains on different laboratory media. Several of the samples produced multiple colony morphotypes that differed from that of an archival isolate. To trace the microevolutionary history of these isolates, we obtained WGS data for several archival and present-day strains and morphotypes. Bacillus-wide phylogenetic analysis identified B. subtilis as the nearest neighbor to B. atrophaeus. The genome of B. atrophaeus is, on average, 86% identical to B. subtilis on the nucleotide level. WGS of variants revealed that several strains were mixed but highly related populations and uncovered a progressive accumulation of mutations among the “military” isolates. Metabolic profiling and microscopic examination of bacterial cultures revealed enhanced growth of “military” isolates on lactate-containing media, and showed that the “military” strains exhibited a hypersporulating phenotype.

Conclusions

Our analysis revealed the genomic and phenotypic signatures of strain adaptation and deliberate selection for traits that were desirable in a simulant organism. Together, these results demonstrate the power of whole-genome and modern systems-level approaches to characterize microbial lineages to develop and validate forensic markers for strain discrimination and reveal signatures of deliberate adaptation.     

Isolation and Characterization of a “phiKMV-Like” Bacteriophage and Its Therapeutic Effect on Mink Hemorrhagic Pneumonia

The objective of this study was to investigate the potential of using phages as a therapy against hemorrhagic pneumonia in mink both in vitro and in vivo. Five Pseudomonas aeruginosa (P. aeruginosa) strains were isolated from lungs of mink with suspected hemorrhagic pneumonia and their identity was confirmed by morphological observation and 16S rDNA sequence analysis. Compared to P. aeruginosa strains isolated from mink with hemorrhagic pneumonia in 2002, these isolates were more resistant to antibiotics selected. A lytic phage vB_PaeP_PPA-ABTNL (PPA-ABTNL) of the Podoviridae family was isolated from hospital sewage using a P. aeruginosa isolate as host, showing broad host range against P. aeruginosa. A one-step growth curve analysis of PPA-ABTNL revealed eclipse and latent periods of 20 and 35 min, respectively, with a burst size of about 110 PFU per infected cell. Phage PPA-ABTNL significantly reduced the growth of P. aeruginosa isolates in vitro. The genome of PPA-ABTNL was 43,227 bp (62.4% G+C) containing 54 open reading frames and lacked regions encoding known virulence factors, integration-related proteins and antibiotic resistance determinants. Genome architecture analysis showed that PPA-ABTNL belonged to the “phiKMV-like Viruses” group. A repeated dose inhalational toxicity study using PPA-ABTNL crude preparation was conducted in mice and no significantly abnormal histological changes, morbidity or mortality were observed. There was no indication of any potential risk associated with using PPA-ABTNL as a therapeutic agent. The results of a curative treatment experiment demonstrated that atomization by ultrasonic treatment could efficiently deliver phage to the lungs of mink and a dose of 10 multiplicity of infection was optimal for treating mink hemorrhagic pneumonia. Our work demonstrated the potential for phage to fight P. aeruginosa involved in mink lung infections when administered by means of ultrasonic nebulization.     

Genome Analysis of Environmental and Clinical P. aeruginosa Isolates from Sequence Type-1146

The genomes of Pseudomonas aeruginosa isolates of the new sequence type ST-1146, three environmental (P37, P47 and P49) and one clinical (SD9) isolates, with differences in their antibiotic susceptibility profiles have been sequenced and analysed. The genomes were mapped against P. aeruginosa PAO1-UW and UCBPP-PA14. The allelic profiles showed that the highest number of differences were in “Related to phage, transposon or plasmid” and “Secreted factors” categories. The clinical isolate showed a number of exclusive alleles greater than that for the environmental isolates. The phage Pf1 region in isolate SD9 accumulated the highest number of nucleotide substitutions. The ORF analysis of the four genomes assembled de novo indicated that the number of isolate-specific genes was higher in isolate SD9 (132 genes) than in isolates P37 (24 genes), P47 (16 genes) and P49 (21 genes). CRISPR elements were found in all isolates and SD9 showed differences in the spacer region. Genes related to bacteriophages F116 and H66 were found only in isolate SD9. Genome comparisons indicated that the isolates of ST-1146 are close related, and most genes implicated in pathogenicity are highly conserved, suggesting a genetic potential for infectivity in the environmental isolates similar to the clinical one. Phage-related genes are responsible of the main differences among the genomes of ST-1146 isolates. The role of bacteriophages has to be considered in the adaptation processes of isolates to the host and in microevolution studies.     

Molecular typing of Treponema pallidum: a systematic review and meta-analysis

Background

Syphilis is resurgent in many regions of the world. Molecular typing is a robust tool for investigating strain diversity and epidemiology. This study aimed to review original research on molecular typing of Treponema pallidum (T. pallidum) with three objectives: (1) to determine specimen types most suitable for molecular typing; (2) to determine T. pallidum subtype distribution across geographic areas; and (3) to summarize available information on subtypes associated with neurosyphilis and macrolide resistance.

Methodology/Principal Findings

Two researchers independently searched five databases from 1998 through 2010, assessed for eligibility and study quality, and extracted data. Search terms included “Treponema pallidum,” or “syphilis,” combined with the subject headings “molecular,” “subtyping,” “typing,” “genotype,” and “epidemiology.” Sixteen eligible studies were included. Publication bias was not statistically significant by the Begg rank correlation test. Medians, inter-quartile ranges, and 95% confidence intervals were determined for DNA extraction and full typing efficiency. A random-effects model was used to perform subgroup analyses to reduce obvious between-study heterogeneity. Primary and secondary lesions and ear lobe blood specimens had an average higher yield of T. pallidum DNA (83.0% vs. 28.2%, χ² = 247.6, p<0.001) and an average higher efficiency of full molecular typing (80.9% vs. 43.1%, χ² = 102.3, p<0.001) compared to plasma, whole blood, and cerebrospinal fluid. A pooled analysis of subtype distribution based on country location showed that 14d was the most common subtype, and subtype distribution varied across geographic areas. Subtype data associated with macrolide resistance and neurosyphilis were limited.

Conclusions/Significance

Primary lesion was a better specimen for obtaining T. pallidum DNA than blood. There was wide geographic variation in T. pallidum subtypes. More research is needed on the relationship between clinical presentation and subtype, and further validation of ear lobe blood for obtaining T. pallidum DNA would be useful for future molecular studies of syphilis.