'Add, stir and reduce': Yersinia spp. as model bacteria for pathogen evolution

In the study of molecular microbiology and bacterial genetics, pathogenic species of the Yersinia genus have been pillars for research aimed at understanding how bacteria evolve into mammalian pathogens. The advent of large-scale population genomic studies has hugely accelerated progress in this field, and the pathogenic Yersinia species have re-emerged as model organisms to help shape our understanding of the evolutionary processes involved in pathogenesis. In this review, we highlight how microbial genomic studies of the yersiniae have revealed distinct features marking the evolutionary path towards pathogenesis, which are changing our understanding of pathogen evolution. As these features are also found in the genomes of other members of the Enterobacteriaceae, they provide a blueprint for the evolution of enteropathogenic bacteria

Directional gene flow and ecological separation in Yersinia enterocolitica

Yersinia enterocolitica is a common cause of food-borne gastroenteritis worldwide. Recent work defining the phylogeny of the genus Yersinia subdivided Y. enterocolitica into six distinct phylogroups. Here, we provide detailed analyses of the evolutionary processes leading to the emergence of these phylogroups. The dominant phylogroups isolated from human infections, PG3–5, show very little diversity at the sequence level, but do present marked patterns of gain and loss of functions, including those involved in pathogenicity and metabolism, including the acquisition of phylogroup-specific O-antigen loci. We tracked gene flow across the species in the core and accessory genome, and show that the non-pathogenic PG1 strains act as a reservoir for diversity, frequently acting as donors in recombination events. Analysis of the core and accessory genome also suggested that the different Y. enterocolitica phylogroups may be ecologically separated, in contrast to the long-held belief of common shared ecological niches across the Y. enterocolitica species

Use of whole-genome sequencing to identify clusters of Shigella flexneri associated with sexual transmission in men who have sex with men in England: a validation study using linked behavioural data

Since the 1970s, shigellosis has been reported as a sexually transmissible infection, and in recent years, genomic data have revealed the breadth of Shigella spp. transmission among global networks of men who have sex with men (MSM). In 2015, Public Health England (PHE) introduced routine whole-genome sequencing (WGS) of Shigella spp. to identify transmission clusters. However, limited behavioural information for the cases hampers interpretation. We investigated whether WGS can distinguish between clusters representing sexual transmission in MSM and clusters representing community (non-sexual) transmission to inform infection control. WGS data for Shigella flexneri from August 2015 to July 2017 were aggregated into single linkage clusters based on SNP typing using a range of SNP distances (the standard for Shigella surveillance at PHE is 10 SNPs). Clusters were classified as 'adult male', 'household', 'travel-associated' or 'community' using routine demographic data submitted alongside laboratory cultures. From August 2015 to March 2017, PHE contacted those with shigellosis as part of routine public-health follow-up and collected exposure data on a structured questionnaire, which for the first time included questions about sexual identity and behaviour. The questionnaire data were used to determine whether clusters classified as 'adult male' represented likely sexual transmission between men, thereby validating the use of the SNP clustering tool for informing appropriate public-health responses. Overall, 1006 S. flexneri cases were reported, of which 563 clustered with at least one other case (10-SNP threshold). Linked questionnaire data were available for 106 clustered cases, of which 84.0 % belonged to an 'adult male' cluster. At the 10-SNP threshold, 95.1 % [95 % confidence interval (CI) 88.0-98.1%] of MSM belonged to an 'adult male' cluster, while 73.2 % (95 % CI 49.1-87.5%) of non-MSM belonged to a 'community' or 'travel-associated' cluster. At the 25-SNP threshold, all MSM (95 % CI 96.0-100%) belonged to an 'adult male' cluster and 77.8 % (95 % CI 59.2-89.4%) of non-MSM belonged to a 'community' or 'travel-associated' cluster. Within one phylogenetic clade of S. flexneri, 9 clusters were identified (7 'adult male'; 2 'community') using a 10-SNP threshold, while a single 'adult male' cluster was identified using a 25-SNP threshold. Genotypic markers of azithromycin resistance were detected in 84.5 % (294/348) of 'adult male' cases and 20.9 % (9/43) of cases in other clusters (10-SNP threshold), the latter of which contained gay-identifying men who reported recent same-sex sexual contact. Our study suggests that SNP clustering can be used to identify Shigella clusters representing likely sexual transmission in MSM to inform infection control. Defining clusters requires a flexible approach in terms of genetic relatedness to ensure a clear understanding of underlying transmission networks

Genomic and phenotypic analyses of recent Acinetobacter baumannii isolates from tertiary care hospitals in Thailand

Antibiotic resistant strains of Acinetobacter baumannii are responsible for a large and increasing burden of nosocomial infections in Thailand and other countries of Southeast Asia. New approaches to their control and treatment are urgently needed and we are actively seeking biological agents that remove the polysaccharide capsules that protect these pathogens from the host’s immune system. To examine phylogenetic relationships, distribution of capsule chemotypes, acquired antibiotic resistance determinants, susceptibility to complement and other traits associated with systemic infection, we sequenced 191 recent isolates from three tertiary referral hospitals in Thailand and used phenotypic assays to characterise key aspects of infectivity. Several distinct lineages were circulating in three hospitals and the majority belonged to global clonal group 2 (GC2). Very high levels of resistance to carbapenems and other front-line antibiotics were found, as were a number of widespread plasmid replicons. A high diversity of capsule genotypes were encountered with only three (KL6, KL10 and KL47) above 10% frequency. Almost 90% of GC2 isolates belonged to the most common capsule genotypes and were fully resistant to the bactericidal action of human serum complement; we attribute this trait to the presence of a substantial protective capsule and for this to represent a key determinant of virulence for systemic infection. We conclude that current Thai nosocomial isolates represent potential targets for therapeutic strategies designed to remove the polysaccharide capsule from extensively drug-resistant A. baumanii during the course of systemic infection

The impact of xerostomia on oral-health-related quality of life among younger adults

BACKGROUND: Recent research has suggested that chronic dry mouth affects the day-to-day lives of older people living in institutions. The condition has usually been considered to be a feature of old age, but recent work by our team produced the somewhat surprising finding that 10% of people in their early thirties are affected. This raises the issue of whether dry mouth is a trivial condition or a more substantial threat to quality of life among younger people. The objective of this study was to examine the association between xerostomia and oral-health-related quality of life among young adults while controlling for clinical oral health status and other potential confounding factors. METHODS: Cross-sectional analysis of data from a longstanding prospective observational study of a Dunedin (New Zealand) birth cohort: clinical dental examinations and questionnaires were used at age 32. The main measures were xerostomia (the subjective feeling of dry mouth, measured with a single question) and oral-health-related quality of life (OHRQoL) measured using the short-form Oral Health Impact Profile (OHIP-14). RESULTS: Of the 923 participants (48.9% female), one in ten were categorised as 'xerostomic', with no apparent gender difference. There was a strong association between xerostomia and OHRQoL (across all OHIP-14 domains) which persisted after multivariate analysis to control for clinical characteristics, gender, smoking status and personality characteristics (negative emotionality and positive emotionality). CONCLUSION: Xerostomia is not a trivial condition; it appears to have marked and consistent effects on sufferers' day-to-day lives

Spatial and genomic data to characterize endemic typhoid transmission

BACKGROUND: Diverse environmental exposures and risk factors have been implicated in the transmission of Salmonella Typhi, however, the dominant transmission pathways through the environment to susceptible humans remain unknown. Here, we utilize spatial, bacterial genomic, and hydrological data to refine our view of Typhoid transmission in an endemic setting. METHODS: 546 patients presenting to Queen Elizabeth Central Hospital in Blantyre, Malawi with blood culture-confirmed typhoid fever between April 2015 and January 2017 were recruited to a cohort study. The households of a subset of these patients were geolocated, and 256 S. Typhi isolates were whole genome sequenced. Pairwise single nucleotide variant (SNV) distances were incorporated into a geostatistical modeling framework using multidimensional scaling. RESULTS: Typhoid fever was not evenly distributed across Blantyre, with estimated minimum incidence ranging across the city from less than 15 to over 100 cases/100,000/year. Pairwise SNV distance and physical household distances were significantly correlated (p=0.001). We evaluated the ability of river catchment to explain the spatial patterns of genomics observed, finding that it significantly improved the fit of the model (p=0.003). We also found spatial correlation at a smaller spatial scale, of households living <192 meters apart. CONCLUSIONS: These findings reinforce the emerging view that hydrological systems play a key role in the transmission of typhoid fever. By combining genomic and spatial data, we show how multi-faceted data can be used to identify high incidence areas, understand the connections between them, and inform targeted environmental surveillance, all of which will be critical to shape local and regional typhoid control strategies

Population-based analysis of ocular Chlamydia trachomatis in trachoma-endemic West African communities identifies genomic markers of disease severity.

BACKGROUND: Chlamydia trachomatis (Ct) is the most common infectious cause of blindness and bacterial sexually transmitted infection worldwide. Ct strain-specific differences in clinical trachoma suggest that genetic polymorphisms in Ct may contribute to the observed variability in severity of clinical disease. METHODS: Using Ct whole genome sequences obtained directly from conjunctival swabs, we studied Ct genomic diversity and associations between Ct genetic polymorphisms with ocular localization and disease severity in a treatment-naïve trachoma-endemic population in Guinea-Bissau, West Africa. RESULTS: All Ct sequences fall within the T2 ocular clade phylogenetically. This is consistent with the presence of the characteristic deletion in trpA resulting in a truncated non-functional protein and the ocular tyrosine repeat regions present in tarP associated with ocular tissue localization. We have identified 21 Ct non-synonymous single nucleotide polymorphisms (SNPs) associated with ocular localization, including SNPs within pmpD (odds ratio, OR = 4.07, p* = 0.001) and tarP (OR = 0.34, p* = 0.009). Eight synonymous SNPs associated with disease severity were found in yjfH (rlmB) (OR = 0.13, p* = 0.037), CTA0273 (OR = 0.12, p* = 0.027), trmD (OR = 0.12, p* = 0.032), CTA0744 (OR = 0.12, p* = 0.041), glgA (OR = 0.10, p* = 0.026), alaS (OR = 0.10, p* = 0.032), pmpE (OR = 0.08, p* = 0.001) and the intergenic region CTA0744-CTA0745 (OR = 0.13, p* = 0.043). CONCLUSIONS: This study demonstrates the extent of genomic diversity within a naturally circulating population of ocular Ct and is the first to describe novel genomic associations with disease severity. These findings direct investigation of host-pathogen interactions that may be important in ocular Ct pathogenesis and disease transmission

MLVA Subtyping of Genovar E Chlamydia trachomatis Individualizes the Swedish Variant and Anorectal Isolates from Men who Have Sex with Men

This study describes a new multilocus variable number tandem-repeat (VNTR) analysis (MLVA) typing system for the discrimination of Chlamydia trachomatis genovar D to K isolates or specimens. We focused our MLVA scheme on genovar E which predominates in most populations worldwide. This system does not require culture and therefore can be performed directly on DNA extracted from positive clinical specimens. Our method was based on GeneScan analysis of five VNTR loci labelled with fluorescent dyes by multiplex PCR and capillary electrophoresis. This MLVA, called MLVA-5, was applied to a collection of 220 genovar E and 94 non-E genovar C. trachomatis isolates and specimens obtained from 251 patients and resulted in 38 MLVA-5 types. The genetic stability of the MLVA-5 scheme was assessed for results obtained both in vitro by serial passage culturing and in vivo using concomitant and sequential isolates and specimens. All anorectal genovar E isolates from men who have sex with men exhibited the same MLVA-5 type, suggesting clonal spread. In the same way, we confirmed the clonal origin of the Swedish new variant of C. trachomatis. The MLVA-5 assay was compared to three other molecular typing methods, ompA gene sequencing, multilocus sequence typing (MLST) and a previous MLVA method called MLVA-3, on 43 genovar E isolates. The discriminatory index was 0.913 for MLVA-5, 0.860 for MLST and 0.622 for MLVA-3. Among all of these genotyping methods, MLVA-5 displayed the highest discriminatory power and does not require a time-consuming sequencing step. The results indicate that MLVA-5 enables high-resolution molecular epidemiological characterisation of C. trachomatis genovars D to K infections directly from specimens

European Chlamydia abortus livestock isolate genomes reveal unusual stability and limited diversity, reflected in geographical signatures.

BACKGROUND: Chlamydia abortus (formerly Chlamydophila abortus) is an economically important livestock pathogen, causing ovine enzootic abortion (OEA), and can also cause zoonotic infections in humans affecting pregnancy outcome. Large-scale genomic studies on other chlamydial species are giving insights into the biology of these organisms but have not yet been performed on C. abortus. Our aim was to investigate a broad collection of European isolates of C. abortus, using next generation sequencing methods, looking at diversity, geographic distribution and genome dynamics. RESULTS: Whole genome sequencing was performed on our collection of 57 C. abortus isolates originating primarily from the UK, Germany, France and Greece, but also from Tunisia, Namibia and the USA. Phylogenetic analysis of a total of 64 genomes shows a deep structural division within the C. abortus species with a major clade displaying limited diversity, in addition to a branch carrying two more distantly related Greek isolates, LLG and POS. Within the major clade, seven further phylogenetic groups can be identified, demonstrating geographical associations. The number of variable nucleotide positions across the sampled isolates is significantly lower than those published for C. trachomatis and C. psittaci. No recombination was identified within C. abortus, and no plasmid was found. Analysis of pseudogenes showed lineage specific loss of some functions, notably with several Pmp and TMH/Inc proteins predicted to be inactivated in many of the isolates studied. CONCLUSIONS: The diversity within C. abortus appears to be much lower compared to other species within the genus. There are strong geographical signatures within the phylogeny, indicating clonal expansion within areas of limited livestock transport. No recombination has been identified within this species, showing that different species of Chlamydia may demonstrate different evolutionary dynamics, and that the genome of C. abortus is highly stable

Shotgun sequencing of Yersinia enterocolitica strain W22703 (biotype 2, serotype O:9): genomic evidence for oscillation between invertebrates and mammals

<p>Abstract</p> <p>Background</p> <p><it>Yersinia enterocolitica </it>strains responsible for mild gastroenteritis in humans are very diverse with respect to their metabolic and virulence properties. Strain W22703 (biotype 2, serotype O:9) was recently identified to possess nematocidal and insecticidal activity. To better understand the relationship between pathogenicity towards insects and humans, we compared the W22703 genome with that of the highly pathogenic strain 8081 (biotype1B; serotype O:8), the only <it>Y. enterocolitica </it>strain sequenced so far.</p> <p>Results</p> <p>We used whole-genome shotgun data to assemble, annotate and analyse the sequence of strain W22703. Numerous factors assumed to contribute to enteric survival and pathogenesis, among them osmoregulated periplasmic glucan, hydrogenases, cobalamin-dependent pathways, iron uptake systems and the <it>Yersinia </it>genome island 1 (YGI-1) involved in tight adherence were identified to be common to the 8081 and W22703 genomes. However, sets of ~550 genes revealed to be specific for each of them in comparison to the other strain. The plasticity zone (PZ) of 142 kb in the W22703 genome carries an ancient flagellar cluster Flg-2 of ~40 kb, but it lacks the pathogenicity island YAPI_Ye, the secretion system <it>ysa </it>and <it>yts1</it>, and other virulence determinants of the 8081 PZ. Its composition underlines the prominent variability of this genome region and demonstrates its contribution to the higher pathogenicity of biotype 1B strains with respect to W22703. A novel type three secretion system of mosaic structure was found in the genome of W22703 that is absent in the sequenced strains of the human pathogenic <it>Yersinia </it>species, but conserved in the genomes of the apathogenic species. We identified several regions of differences in W22703 that mainly code for transporters, regulators, metabolic pathways, and defence factors.</p> <p>Conclusion</p> <p>The W22703 sequence analysis revealed a genome composition distinct from other pathogenic <it>Yersinia enterocolitica </it>strains, thus contributing novel data to the <it>Y. enterocolitica </it>pan-genome. This study also sheds further light on the strategies of this pathogen to cope with its environments.</p