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.

Genetic Markers in S. Paratyphi C Reveal Primary Adaptation to Pigs.

Salmonella enterica with the identical antigenic formula 6,7:c:1,5 can be differentiated biochemically and by disease syndrome. One grouping, Salmonella Paratyphi C, is currently considered a typhoidal serovar, responsible for enteric fever in humans. The human-restricted typhoidal serovars (S. Typhi and Paratyphi A, B and C) typically display high levels of genome degradation and are cited as an example of convergent evolution for host adaptation in humans. However, S. Paratyphi C presents a different clinical picture to S. Typhi/Paratyphi A, in a patient group with predisposition, raising the possibility that its natural history is different, and that infection is invasive salmonellosis rather than enteric fever. Using whole genome sequencing and metabolic pathway analysis, we compared the genomes of 17 S. Paratyphi C strains to other members of the 6,7:c:1,5 group and to two typhoidal serovars: S. Typhi and Paratyphi A. The genome degradation observed in S. Paratyphi C was much lower than S. Typhi/Paratyphi A, but similar to the other 6,7:c:1,5 strains. Genomic and metabolic comparisons revealed little to no overlap between S. Paratyphi C and the other typhoidal serovars, arguing against convergent evolution and instead providing evidence of a primary adaptation to pigs in accordance with the 6,7:c:1.5 strains.

Whole genome sequence of Vibrio cholerae directly from dried spotted filter paper.

BACKGROUND: Global estimates for cholera annually approximate 4 million cases worldwide with 95,000 deaths. Recent outbreaks, including Haiti and Yemen, are reminders that cholera is still a global health concern. Cholera outbreaks can rapidly induce high death tolls by overwhelming the capacity of health facilities, especially in remote areas or areas of civil unrest. Recent studies demonstrated that stool specimens preserved on filter paper facilitate molecular analysis of Vibrio cholerae in resource limited settings. Specimens preserved in a rapid, low-cost, safe and sustainable manner for sequencing provides previously unavailable data about circulating cholera strains. This may ultimately contribute new information to shape public policy response on cholera control and elimination. METHODOLOGY/PRINCIPAL FINDINGS: Whole genome sequencing (WGS) recovered close to a complete sequence of the V. cholerae O1 genome with satisfactory genome coverage from stool specimens enriched in alkaline peptone water (APW) and V. cholerae culture isolates, both spotted on filter paper. The minimum concentration of V. cholerae DNA sufficient to produce quality genomic information was 0.02 ng/µL. The genomic data confirmed the presence or absence of genes of epidemiological interest, including cholera toxin and pilus loci. WGS identified a variety of diarrheal pathogens from APW-enriched specimen spotted filter paper, highlighting the potential for this technique to explore the gut microbiome, potentially identifying co-infections, which may impact the severity of disease. WGS demonstrated that these specimens fit within the current global cholera phylogenetic tree, identifying the strains as the 7th pandemic El Tor. CONCLUSIONS: WGS results allowed for mapping of short reads from APW-enriched specimen and culture isolate spotted filter papers. This provided valuable molecular epidemiological sequence information on V. cholerae strains from remote, low-resource settings. These results identified the presence of co-infecting pathogens while providing rare insight into the specific V. cholerae strains causing outbreaks in cholera-endemic areas.

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.

Morphological, genomic and transcriptomic responses of Klebsiella pneumoniae to the last-line antibiotic colistin.

Colistin remains one of the few antibiotics effective against multi-drug resistant (MDR) hospital pathogens, such as Klebsiella pneumoniae. Yet resistance to this last-line drug is rapidly increasing. Characterized mechanisms of colR in K. pneumoniae are largely due to chromosomal mutations in two-component regulators, although a plasmid-mediated colR mechanism has recently been uncovered. However, the effects of intrinsic colistin resistance are yet to be characterized on a whole-genome level. Here, we used a genomics-based approach to understand the mechanisms of adaptive colR acquisition in K. pneumoniae. In controlled directed-evolution experiments we observed two distinct paths to colistin resistance acquisition. Whole genome sequencing identified mutations in two colistin resistance genes: in the known colR regulator phoQ which became fixed in the population and resulted in a single amino acid change, and unstable minority variants in the recently described two-component sensor crrB. Through RNAseq and microscopy, we reveal the broad range of effects that colistin exposure has on the cell. This study is the first to use genomics to identify a population of minority variants with mutations in a colR gene in K. pneumoniae.

Dynamics of cholera epidemics from Benin to Mauritania.

BACKGROUND: The countries of West Africa are largely portrayed as cholera endemic, although the dynamics of outbreaks in this region of Africa remain largely unclear. METHODOLOGY/PRINCIPAL FINDINGS: To understand the dynamics of cholera in a major portion of West Africa, we analyzed cholera epidemics from 2009 to 2015 from Benin to Mauritania. We conducted a series of field visits as well as multilocus variable tandem repeat analysis and whole-genome sequencing analysis of V. cholerae isolates throughout the study region. During this period, Ghana accounted for 52% of the reported cases in the entire study region (coastal countries from Benin to Mauritania). From 2009 to 2015, we found that one major wave of cholera outbreaks spread from Accra in 2011 northwestward to Sierra Leone and Guinea in 2012. Molecular epidemiology analysis confirmed that the 2011 Ghanaian isolates were related to those that seeded the 2012 epidemics in Guinea and Sierra Leone. Interestingly, we found that many countries deemed "cholera endemic" actually suffered very few outbreaks, with multi-year lulls. CONCLUSIONS/SIGNIFICANCE: This study provides the first cohesive vision of the dynamics of cholera epidemics in a major portion of West Africa. This epidemiological overview shows that from 2009 to 2015, at least 54% of reported cases concerned populations living in the three urban areas of Accra, Freetown, and Conakry. These findings may serve as a guide to better target cholera prevention and control efforts in the identified cholera hotspots in West Africa.

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.

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.

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.

Transcriptome and proteome analysis of Salmonella enterica serovar Typhimurium systemic infection of wild type and immune-deficient mice.

Salmonella enterica are a threat to public health. Current vaccines are not fully effective. The ability to grow in infected tissues within phagocytes is required for S. enterica virulence in systemic disease. As the infection progresses the bacteria are exposed to a complex host immune response. Consequently, in order to continue growing in the tissues, S. enterica requires the coordinated regulation of fitness genes. Bacterial gene regulation has so far been investigated largely using exposure to artificial environmental conditions or to in vitro cultured cells, and little information is available on how S. enterica adapts in vivo to sustain cell division and survival. We have studied the transcriptome, proteome and metabolic flux of Salmonella, and the transcriptome of the host during infection of wild type C57BL/6 and immune-deficient gp91-/-phox mice. Our analyses advance the understanding of how S. enterica and the host behaves during infection to a more sophisticated level than has previously been reported.

The secondary resistome of multidrug-resistant Klebsiella pneumoniae.

Klebsiella pneumoniae causes severe lung and bloodstream infections that are difficult to treat due to multidrug resistance. We hypothesized that antimicrobial resistance can be reversed by targeting chromosomal non-essential genes that are not responsible for acquired resistance but essential for resistant bacteria under therapeutic concentrations of antimicrobials. Conditional essentiality of individual genes to antimicrobial resistance was evaluated in an epidemic multidrug-resistant clone of K. pneumoniae (ST258). We constructed a high-density transposon mutant library of >430,000 unique Tn5 insertions and measured mutant depletion upon exposure to three clinically relevant antimicrobials (colistin, imipenem or ciprofloxacin) by Transposon Directed Insertion-site Sequencing (TraDIS). Using this high-throughput approach, we defined three sets of chromosomal non-essential genes essential for growth during exposure to colistin (n = 35), imipenem (n = 1) or ciprofloxacin (n = 1) in addition to known resistance determinants, collectively termed the "secondary resistome". As proof of principle, we demonstrated that inactivation of a non-essential gene not previously found linked to colistin resistance (dedA) restored colistin susceptibility by reducing the minimum inhibitory concentration from 8 to 0.5 µg/ml, 4-fold below the susceptibility breakpoint (S ≤ 2 µg/ml). This finding suggests that the secondary resistome is a potential target for developing antimicrobial "helper" drugs that restore the efficacy of existing antimicrobials.

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.

What's in a Name? Species-Wide Whole-Genome Sequencing Resolves Invasive and Noninvasive Lineages of Salmonella enterica Serotype Paratyphi B.

UNLABELLED: For 100 years, it has been obvious that Salmonella enterica strains sharing the serotype with the formula 1,4,[5],12:b:1,2-now known as Paratyphi B-can cause diseases ranging from serious systemic infections to self-limiting gastroenteritis. Despite considerable predicted diversity between strains carrying the common Paratyphi B serotype, there remain few methods that subdivide the group into groups that are congruent with their disease phenotypes. Paratyphi B therefore represents one of the canonical examples in Salmonella where serotyping combined with classical microbiological tests fails to provide clinically informative information. Here, we use genomics to provide the first high-resolution view of this serotype, placing it into a wider genomic context of the Salmonella enterica species. These analyses reveal why it has been impossible to subdivide this serotype based upon phenotypic and limited molecular approaches. By examining the genomic data in detail, we are able to identify common features that correlate with strains of clinical importance. The results presented here provide new diagnostic targets, as well as posing important new questions about the basis for the invasive disease phenotype observed in a subset of strains. IMPORTANCE: Salmonella enterica strains carrying the serotype Paratyphi B have long been known to possess Jekyll and Hyde characteristics; some cause gastroenteritis, while others cause serious invasive disease. Understanding what makes up the population of strains carrying this serotype, as well as the source of their invasive disease, is a 100-year-old puzzle that we address here using genomics. Our analysis provides the first high-resolution view of this serotype, placing strains carrying serotype Paratyphi B into the wider genomic context of the Salmonella enterica species. This work reveals a history of disease dating back to the middle ages, caused by a group of distinct lineages with various abilities to cause invasive disease. By quantifying the key genomic differences between the invasive and noninvasive populations, we are able to identify key virulence-related targets that can form the basis of simple, rapid, point-of-care tests.

Distinct Salmonella Enteritidis lineages associated with enterocolitis in high-income settings and invasive disease in low-income settings.

An epidemiological paradox surrounds Salmonella enterica serovar Enteritidis. In high-income settings, it has been responsible for an epidemic of poultry-associated, self-limiting enterocolitis, whereas in sub-Saharan Africa it is a major cause of invasive nontyphoidal Salmonella disease, associated with high case fatality. By whole-genome sequence analysis of 675 isolates of S. Enteritidis from 45 countries, we show the existence of a global epidemic clade and two new clades of S. Enteritidis that are geographically restricted to distinct regions of Africa. The African isolates display genomic degradation, a novel prophage repertoire, and an expanded multidrug resistance plasmid. S. Enteritidis is a further example of a Salmonella serotype that displays niche plasticity, with distinct clades that enable it to become a prominent cause of gastroenteritis in association with the industrial production of eggs and of multidrug-resistant, bloodstream-invasive infection in Africa.

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.

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.].

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.

Whole genome analysis to detect potential vaccine-induced changes on Shigella sonnei genome.

Shigellosis or bacillary dysentery is endemic worldwide and is a significant cause of death in children less than five years of age in developing countries. There are no licensed Shigella vaccines and glycoconjugates are among the leading candidate vaccines against shigellosis today. We used whole genome sequence analysis (WGA) to find out whether immunization, with an investigational Shigella sonnei glycoconjugate, could induce selective pressure leading to changes in the genome of S. sonnei. An outbreak of culture-proven S. sonnei shigellosis which occurred immediately after vaccination in one of the cohorts of volunteers participating in a phase III trial of the vaccine in Israel created a unique condition in which the epidemic agent "co-existed" with the developing immune responses induced by the vaccine and natural infection among vaccinees who developed S. sonnei shigellosis. By comparing the whole genomes of S. sonnei isolated from vaccinees and from volunteers in the control group, we show at a very high sensitivity that a potent S. sonnei glycoconjugate that conferred 74% protective efficacy against the homologous disease did not induce changes in the genome of S. sonnei and in particular on the O-antigen gene cluster.