Diagnostics for Yaws Eradication: Insights From Direct Next-Generation Sequencing of Cutaneous Strains of Treponema pallidum.

Background: Yaws-like chronic ulcers can be caused by Treponema pallidum subspecies pertenue, Haemophilus ducreyi, or other, still-undefined bacteria. To permit accurate evaluation of yaws elimination efforts, programmatic use of molecular diagnostics is required. The accuracy and sensitivity of current tools remain unclear because our understanding of T. pallidum diversity is limited by the low number of sequenced genomes. Methods: We tested samples from patients with suspected yaws collected in the Solomon Islands and Ghana. All samples were from patients whose lesions had previously tested negative using the Centers for Disease Control and Prevention (CDC) diagnostic assay in widespread use. However, some of these patients had positive serological assays for yaws on blood. We used direct whole-genome sequencing to identify T. pallidum subsp pertenue strains missed by the current assay. Results: From 45 Solomon Islands and 27 Ghanaian samples, 11 were positive for T. pallidum DNA using the species-wide quantitative polymerase chain reaction (PCR) assay, from which we obtained 6 previously undetected T. pallidum subsp pertenue whole-genome sequences. These show that Solomon Islands sequences represent distinct T. pallidum subsp pertenue clades. These isolates were invisible to the CDC diagnostic PCR assay, due to sequence variation in the primer binding site. Conclusions: Our data double the number of published T. pallidum subsp pertenue genomes. We show that Solomon Islands strains are undetectable by the PCR used in many studies and by health ministries. This assay is therefore not adequate for the eradication program. Next-generation genome sequence data are essential for these efforts.

Direct Whole-Genome Sequencing of Cutaneous Strains of Haemophilus ducreyi.

Haemophilus ducreyi, which causes chancroid, has emerged as a cause of pediatric skin disease. Isolation of H. ducreyi in low-income settings is challenging, limiting phylogenetic investigation. Next-generation sequencing demonstrates that cutaneous strains arise from class I and II H. ducreyi clades and that class II may represent a distinct subspecies.

Functional analysis of Salmonella Typhi adaptation to survival in water.

Contaminated water is a major risk factor associated with the transmission of Salmonella enterica serovar Typhi (S. Typhi), the aetiological agent of human typhoid. However, little is known about how this pathogen adapts to living in the aqueous environment. We used transcriptome analysis (RNA-seq) and transposon mutagenesis (TraDIS) to characterize these adaptive changes and identify multiple genes that contribute to survival. Over half of the genes in the S. Typhi genome altered expression level within the first 24 h following transfer from broth culture to water, although relatively few did so in the first 30 min. Genes linked to central metabolism, stress associated with arrested proton motive force and respiratory chain factors changed expression levels. Additionally, motility and chemotaxis genes increased expression, consistent with a scavenging lifestyle. The viaB-associated gene tviC encoding a glcNAc epimerase that is required for Vi polysaccharide biosynthesis was, along with several other genes, shown to contribute to survival in water. Thus, we define regulatory adaptation operating in S. Typhi that facilitates survival in water.

Correction: Elucidation of the RamA Regulon in Klebsiella pneumoniae Reveals a Role in LPS Regulation.

[This corrects the article DOI: 10.1371/journal.ppat.1004627.].

Patterns of genome evolution that have accompanied host adaptation in Salmonella.

Many bacterial pathogens are specialized, infecting one or few hosts, and this is often associated with more acute disease presentation. Specific genomes show markers of this specialization, which often reflect a balance between gene acquisition and functional gene loss. Within Salmonella enterica subspecies enterica, a single lineage exists that includes human and animal pathogens adapted to cause infection in different hosts, including S. enterica serovar Enteritidis (multiple hosts), S. Gallinarum (birds), and S. Dublin (cattle). This provides an excellent evolutionary context in which differences between these pathogen genomes can be related to host range. Genome sequences were obtained from ∼ 60 isolates selected to represent the known diversity of this lineage. Examination and comparison of the clades within the phylogeny of this lineage revealed signs of host restriction as well as evolutionary events that mark a path to host generalism. We have identified the nature and order of events for both evolutionary trajectories. The impact of functional gene loss was predicted based upon position within metabolic pathways and confirmed with phenotyping assays. The structure of S. Enteritidis is more complex than previously known, as a second clade of S. Enteritidis was revealed that is distinct from those commonly seen to cause disease in humans or animals, and that is more closely related to S. Gallinarum. Isolates from this second clade were tested in a chick model of infection and exhibited a reduced colonization phenotype, which we postulate represents an intermediate stage in pathogen-host adaptation.

Mycoplasma genitalium: whole genome sequence analysis, recombination and population structure.

BACKGROUND: Although Mycoplasma genitalium is a common sexually transmitted pathogen causing clinically distinct diseases both in male and females, few genomes have been sequenced up to now, due mainly to its fastidious nature and slow growth. Hence, we lack a robust phylogenetic framework to provide insights into the population structure of the species. Currently our understanding of the nature and diversity of M. genitalium relies on molecular tests targeting specific genes or regions of the genome and knowledge is limited by a general under-testing internationally. This is set against a background of drug resistance whereby M. genitalium has developed resistance to mainly all therapeutic antimicrobials. RESULTS: We sequenced 28 genomes of Mycoplasma genitalium from temporally (1980-2010) and geographically (Europe, Japan, Australia) diverse sources. All the strain showed essentially the same genomic content without any accessory regions found. However, we identified extensive recombination across their genomes with a total of 25 regions showing heightened levels of SNP density. These regions include the MgPar loci, associated with host interactions, as well as other genes that could also be involved in this role. Using these data, we generated a robust phylogeny which shows that there are two main clades with differing degrees of genomic variability. SNPs found in region V of 23S rRNA and parC were consistent with azithromycin/erythromycin and fluoroquinolone resistances, respectively, and with their phenotypic MIC data. CONCLUSIONS: The sequence data here generated is essential for designing rational approaches to type and track Mycoplasma genitalium as antibiotic resistance increases. It represents a first approach to its population genetics to better appreciate the role of this organism as a sexually transmitted pathogen.

Elucidation of the RamA regulon in Klebsiella pneumoniae reveals a role in LPS regulation.

Klebsiella pneumoniae is a significant human pathogen, in part due to high rates of multidrug resistance. RamA is an intrinsic regulator in K. pneumoniae established to be important for the bacterial response to antimicrobial challenge; however, little is known about its possible wider regulatory role in this organism during infection. In this work, we demonstrate that RamA is a global transcriptional regulator that significantly perturbs the transcriptional landscape of K. pneumoniae, resulting in altered microbe-drug or microbe-host response. This is largely due to the direct regulation of 68 genes associated with a myriad of cellular functions. Importantly, RamA directly binds and activates the lpxC, lpxL-2 and lpxO genes associated with lipid A biosynthesis, thus resulting in modifications within the lipid A moiety of the lipopolysaccharide. RamA-mediated alterations decrease susceptibility to colistin E, polymyxin B and human cationic antimicrobial peptide LL-37. Increased RamA levels reduce K. pneumoniae adhesion and uptake into macrophages, which is supported by in vivo infection studies, that demonstrate increased systemic dissemination of ramA overexpressing K. pneumoniae. These data establish that RamA-mediated regulation directly perturbs microbial surface properties, including lipid A biosynthesis, which facilitate evasion from the innate host response. This highlights RamA as a global regulator that confers pathoadaptive phenotypes with implications for our understanding of the pathogenesis of Enterobacter, Salmonella and Citrobacter spp. that express orthologous RamA proteins.

Parallel independent evolution of pathogenicity within the genus Yersinia.

The genus Yersinia has been used as a model system to study pathogen evolution. Using whole-genome sequencing of all Yersinia species, we delineate the gene complement of the whole genus and define patterns of virulence evolution. Multiple distinct ecological specializations appear to have split pathogenic strains from environmental, nonpathogenic lineages. This split demonstrates that contrary to hypotheses that all pathogenic Yersinia species share a recent common pathogenic ancestor, they have evolved independently but followed parallel evolutionary paths in acquiring the same virulence determinants as well as becoming progressively more limited metabolically. Shared virulence determinants are limited to the virulence plasmid pYV and the attachment invasion locus ail. These acquisitions, together with genomic variations in metabolic pathways, have resulted in the parallel emergence of related pathogens displaying an increasingly specialized lifestyle with a spectrum of virulence potential, an emerging theme in the evolution of other important human pathogens.

The Type III Secretion System Effector SeoC of Salmonella enterica subsp. salamae and S. enterica subsp. arizonae ADP-Ribosylates Src and Inhibits Opsonophagocytosis.

Salmonella species utilize type III secretion systems (T3SSs) to translocate effectors into the cytosol of mammalian host cells, subverting cell signaling and facilitating the onset of gastroenteritis. In this study, we compared a draft genome assembly of Salmonella enterica subsp. salamae strain 3588/07 against the genomes of S. enterica subsp. enterica serovar Typhimurium strain LT2 and Salmonella bongori strain 12419. S. enterica subsp. salamae encodes the Salmonella pathogenicity island 1 (SPI-1), SPI-2, and the locus of enterocyte effacement (LEE) T3SSs. Though several key S Typhimurium effector genes are missing (e.g., avrA, sopB, and sseL), S. enterica subsp. salamae invades HeLa cells and contains homologues of S. bongori sboK and sboC, which we named seoC SboC and SeoC are homologues of EspJ from enteropathogenic and enterohemorrhagic Escherichia coli (EPEC and EHEC, respectively), which inhibit Src kinase-dependent phagocytosis by ADP-ribosylation. By screening 73 clinical and environmental Salmonella isolates, we identified EspJ homologues in S. bongori, S. enterica subsp. salamae, and Salmonella enterica subsp. arizonae The ß-lactamase TEM-1 reporter system showed that SeoC is translocated by the SPI-1 T3SS. All the Salmonella SeoC/SboC homologues ADP-ribosylate Src E310 in vitro Ectopic expression of SeoC/SboC inhibited phagocytosis of IgG-opsonized beads into Cos-7 cells stably expressing green fluorescent protein (GFP)-FcγRIIa. Concurrently, S. enterica subsp. salamae infection of J774.A1 macrophages inhibited phagocytosis of beads, in a seoC-dependent manner. These results show that S. bongori, S. enterica subsp. salamae, and S. enterica subsp. arizonae share features of the infection strategy of extracellular pathogens EPEC and EHEC and shed light on the complexities of the T3SS effector repertoires of Enterobacteriaceae.

Travel- and Community-Based Transmission of Multidrug-Resistant Shigella sonnei Lineage among International Orthodox Jewish Communities.

Shigellae are sensitive indicator species for studying trends in the international transmission of antimicrobial-resistant Enterobacteriaceae. Orthodox Jewish communities (OJCs) are a known risk group for shigellosis; Shigella sonnei is cyclically epidemic in OJCs in Israel, and sporadic outbreaks occur in OJCs elsewhere. We generated whole-genome sequences for 437 isolates of S. sonnei from OJCs and non-OJCs collected over 22 years in Europe (the United Kingdom, France, and Belgium), the United States, Canada, and Israel and analyzed these within a known global genomic context. Through phylogenetic and genomic analysis, we showed that strains from outbreaks in OJCs outside of Israel are distinct from strains in the general population and relate to a single multidrug-resistant sublineage of S. sonnei that prevails in Israel. Further Bayesian phylogenetic analysis showed that this strain emerged approximately 30 years ago, demonstrating the speed at which antimicrobial drug-resistant pathogens can spread widely through geographically dispersed, but internationally connected, communities.

Roary: rapid large-scale prokaryote pan genome analysis.

UNLABELLED: A typical prokaryote population sequencing study can now consist of hundreds or thousands of isolates. Interrogating these datasets can provide detailed insights into the genetic structure of prokaryotic genomes. We introduce Roary, a tool that rapidly builds large-scale pan genomes, identifying the core and accessory genes. Roary makes construction of the pan genome of thousands of prokaryote samples possible on a standard desktop without compromising on the accuracy of results. Using a single CPU Roary can produce a pan genome consisting of 1000 isolates in 4.5 hours using 13 GB of RAM, with further speedups possible using multiple processors. AVAILABILITY AND IMPLEMENTATION: Roary is implemented in Perl and is freely available under an open source GPLv3 license from http://sanger-pathogens.github.io/Roary CONTACT: roary@sanger.ac.uk SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Dominant role of nucleotide substitution in the diversification of serotype 3 pneumococci over decades and during a single infection.

Streptococcus pneumoniae of serotype 3 possess a mucoid capsule and cause disease associated with high mortality rates relative to other pneumococci. Phylogenetic analysis of a complete reference genome and 81 draft sequences from clonal complex 180, the predominant serotype 3 clone in much of the world, found most sampled isolates belonged to a clade affected by few diversifying recombinations. However, other isolates indicate significant genetic variation has accumulated over the clonal complex's entire history. Two closely related genomes, one from the blood and another from the cerebrospinal fluid, were obtained from a patient with meningitis. The pair differed in their behaviour in a mouse model of disease and in their susceptibility to antimicrobials, with at least some of these changes attributable to a mutation that up-regulated the patAB efflux pump. This indicates clinically important phenotypic variation can accumulate rapidly through small alterations to the genotype.

RNA-seq analysis of the influence of anaerobiosis and FNR on Shigella flexneri.

BACKGROUND: Shigella flexneri is an important human pathogen that has to adapt to the anaerobic environment in the gastrointestinal tract to cause dysentery. To define the influence of anaerobiosis on the virulence of Shigella, we performed deep RNA sequencing to identify transcriptomic differences that are induced by anaerobiosis and modulated by the anaerobic Fumarate and Nitrate Reduction regulator, FNR. RESULTS: We found that 528 chromosomal genes were differentially expressed in response to anaerobic conditions; of these, 228 genes were also influenced by FNR. Genes that were up-regulated in anaerobic conditions are involved in carbon transport and metabolism (e.g. ptsG, manX, murQ, cysP, cra), DNA topology and regulation (e.g. ygiP, stpA, hns), host interactions (e.g. yciD, nmpC, slyB, gapA, shf, msbB) and survival within the gastrointestinal tract (e.g. shiA, ospI, adiY, cysP). Interestingly, there was a marked effect of available oxygen on genes involved in Type III secretion system (T3SS), which is required for host cell invasion and pathogenesis. These genes, located on the large Shigella virulence plasmid, were down regulated in anaerobiosis in an FNR-dependent manner. We also confirmed anaerobic induction of csrB and csrC small RNAs in an FNR-independent manner. CONCLUSIONS: Anaerobiosis promotes survival and adaption strategies of Shigella, while modulating virulence plasmid genes involved in T3SS-mediated host cell invasion. The influence of FNR on this process is more extensive than previously appreciated, although aside from the virulence plasmid, this transcriptional regulator does not govern expression of genes on other horizontally acquired sequences on the chromosome such as pathogenicity islands.

Using 16s rRNA sequencing to characterize the microbiome of tropical cutaneous ulcer disease: insights into the microbial landscape and implications for diagnosis and treatment.

Cutaneous ulcers are common in yaws-endemic areas. Although often attributed to 'Treponema pallidum subsp. pertenue' and Haemophilus ducreyi, quantitative PCR has highlighted a significant proportion of these ulcers are negative for both pathogens and are considered idiopathic. This is a retrospective analysis utilising existing 16S rRNA sequencing data from two independent yaws studies that took place in Ghana and the Solomon Islands. We characterized bacterial diversity in 38 samples to identify potential causative agents for idiopathic cutaneous ulcers. We identified a diverse bacterial profile, including Arcanobacterium haemolyticum, Campylobacter concisus, Corynebacterium diphtheriae, Staphylococcus spp. and Streptococcus pyogenes, consistent with findings from previous cutaneous ulcer microbiome studies. No single bacterial species was universally present across all samples. The most prevalent bacterium, Campylobacter ureolyticus, appeared in 42% of samples, suggesting a multifactorial aetiology for cutaneous ulcers in yaws-endemic areas. This study emphasizes the need for a nuanced understanding of potential causative agents. The findings prompt further exploration into the intricate microbial interactions contributing to idiopathic yaw-like ulcers, guiding future research toward comprehensive diagnostic and therapeutic strategies.

Salmonella bongori provides insights into the evolution of the Salmonellae.

The genus Salmonella contains two species, S. bongori and S. enterica. Compared to the well-studied S. enterica there is a marked lack of information regarding the genetic makeup and diversity of S. bongori. S. bongori has been found predominantly associated with cold-blooded animals, but it can infect humans. To define the phylogeny of this species, and compare it to S. enterica, we have sequenced 28 isolates representing most of the known diversity of S. bongori. This cross-species analysis allowed us to confidently differentiate ancestral functions from those acquired following speciation, which include both metabolic and virulence-associated capacities. We show that, although S. bongori inherited a basic set of Salmonella common virulence functions, it has subsequently elaborated on this in a different direction to S. enterica. It is an established feature of S. enterica evolution that the acquisition of the type III secretion systems (T3SS-1 and T3SS-2) has been followed by the sequential acquisition of genes encoding secreted targets, termed effectors proteins. We show that this is also true of S. bongori, which has acquired an array of novel effector proteins (sboA-L). All but two of these effectors have no significant S. enterica homologues and instead are highly similar to those found in enteropathogenic Escherichia coli (EPEC). Remarkably, SboH is found to be a chimeric effector protein, encoded by a fusion of the T3SS-1 effector gene sopA and a gene highly similar to the EPEC effector nleH from enteropathogenic E. coli. We demonstrate that representatives of these new effectors are translocated and that SboH, similarly to NleH, blocks intrinsic apoptotic pathways while being targeted to the mitochondria by the SopA part of the fusion. This work suggests that S. bongori has inherited the ancestral Salmonella virulence gene set, but has adapted by incorporating virulence determinants that resemble those employed by EPEC.

Citrobacter rodentium is an unstable pathogen showing evidence of significant genomic flux.

Citrobacter rodentium is a natural mouse pathogen that causes attaching and effacing (A/E) lesions. It shares a common virulence strategy with the clinically significant human A/E pathogens enteropathogenic E. coli (EPEC) and enterohaemorrhagic E. coli (EHEC) and is widely used to model this route of pathogenesis. We previously reported the complete genome sequence of C. rodentium ICC168, where we found that the genome displayed many characteristics of a newly evolved pathogen. In this study, through PFGE, sequencing of isolates showing variation, whole genome transcriptome analysis and examination of the mobile genetic elements, we found that, consistent with our previous hypothesis, the genome of C. rodentium is unstable as a result of repeat-mediated, large-scale genome recombination and because of active transposition of mobile genetic elements such as the prophages. We sequenced an additional C. rodentium strain, EX-33, to reveal that the reference strain ICC168 is representative of the species and that most of the inactivating mutations were common to both isolates and likely to have occurred early on in the evolution of this pathogen. We draw parallels with the evolution of other bacterial pathogens and conclude that C. rodentium is a recently evolved pathogen that may have emerged alongside the development of inbred mice as a model for human disease.

Salmonella Typhi sense host neuroendocrine stress hormones and release the toxin haemolysin E.

Salmonella enterica serovar Typhi (S. typhi) causes typhoid fever. We show that exposure of S. typhi to neuroendocrine stress hormones results in haemolysis, which is associated with the release of haemolysin E in membrane vesicles. This effect is attributed to increased expression of the small RNA micA and RNA chaperone Hfq, with concomitant downregulation of outer membrane protein A. Deletion of micA or the two-component signal-transduction system, CpxAR, abolishes the phenotype. The hormone response is inhibited by the ß-blocker propranolol. We provide mechanistic insights into the basis of neuroendocrine hormone-mediated haemolysis by S. typhi, increasing our understanding of inter-kingdom signalling.

Structure, diversity, and mobility of the Salmonella pathogenicity island 7 family of integrative and conjugative elements within Enterobacteriaceae.

Integrative and conjugative elements (ICEs) are self-mobile genetic elements found in the genomes of some bacteria. These elements may confer a fitness advantage upon their host bacteria through the cargo genes that they carry. Salmonella pathogenicity island 7 (SPI-7), found within some pathogenic strains of Salmonella enterica, possesses features indicative of an ICE and carries genes implicated in virulence. We aimed to identify and fully analyze ICEs related to SPI-7 within the genus Salmonella and other Enterobacteriaceae. We report the sequence of two novel SPI-7-like elements, found within strains of Salmonella bongori, which share 97% nucleotide identity over conserved regions with SPI-7 and with each other. Although SPI-7 within Salmonella enterica serovar Typhi appears to be fixed within the chromosome, we present evidence that these novel elements are capable of excision and self-mobility. Phylogenetic analyses show that these Salmonella mobile elements share an ancestor which existed approximately 3.6 to 15.8 million years ago. Additionally, we identified more distantly related ICEs, with distinct cargo regions, within other strains of Salmonella as well as within Citrobacter, Erwinia, Escherichia, Photorhabdus, and Yersinia species. In total, we report on a collection of 17 SPI-7 related ICEs within enterobacterial species, of which six are novel. Using comparative and mutational studies, we have defined a core of 27 genes essential for conjugation. We present a growing family of SPI-7-related ICEs whose mobility, abundance, and cargo variability indicate that these elements may have had a large impact on the evolution of the Enterobacteriaceae.

The TCA cycle is not required for selection or survival of multidrug-resistant Salmonella.

OBJECTIVES: The initial aim of this study was to use a systems biology approach to analyse a ciprofloxacin-selected multidrug-resistant (MDR) Salmonella enterica serotype Typhimurium, L664. METHODS: The whole genome sequence and transcriptome of L664 were analysed. Site-directed mutagenesis to recreate each mutation was carried out, followed by phenotypic characterization and mutation frequency analysis. As a mutation in the TCA cycle was detected we tested the controversial hypothesis regarding the bacterial response to bactericidal antibiotics, put forward by Kohanski et al. (Cell 2007; 130: 797-810 and Mol Cell 2010; 37: 311-20), that exposure of bacteria to agents such as ciprofloxacin produces reactive oxygen species (ROS), which transiently increase the mutation rate giving rise to MDR bacteria. RESULTS: L664 contained a mutation in ramR that conferred MDR. A mutation in tctA affected the TCA cycle and conferred the inability to grow on minimal agar. The virulence of L664 was not attenuated. Ciprofloxacin exposure produced ROS in L664 and SL1344 (tctA::aph), but it was reduced and occurred later. There were no significant differences in the rates of killing or mutations per generation to antibiotic resistance between the strains. CONCLUSIONS: Whilst we confirm production of ROS in response to ciprofloxacin, we have no data to support the hypothesis that this leads to selection of MDR strains. Our results indicate that the mutations in tctA and glgA were random as they did not pre-exist in the parental strain, and that the mutation in tctA did not provide a survival advantage or disadvantage in the presence of antibiotic.

A strand-specific RNA-Seq analysis of the transcriptome of the typhoid bacillus Salmonella typhi.

High-density, strand-specific cDNA sequencing (ssRNA-seq) was used to analyze the transcriptome of Salmonella enterica serovar Typhi (S. Typhi). By mapping sequence data to the entire S. Typhi genome, we analyzed the transcriptome in a strand-specific manner and further defined transcribed regions encoded within prophages, pseudogenes, previously un-annotated, and 3'- or 5'-untranslated regions (UTR). An additional 40 novel candidate non-coding RNAs were identified beyond those previously annotated. Proteomic analysis was combined with transcriptome data to confirm and refine the annotation of a number of hpothetical genes. ssRNA-seq was also combined with microarray and proteome analysis to further define the S. Typhi OmpR regulon and identify novel OmpR regulated transcripts. Thus, ssRNA-seq provides a novel and powerful approach to the characterization of the bacterial transcriptome.