The genomic epidemiology of shigellosis in South Africa.

Shigellosis, a leading cause of diarrhoeal mortality and morbidity globally, predominantly affects children under five years of age living in low- and middle-income countries. While whole genome sequence analysis (WGSA) has been effectively used to further our understanding of shigellosis epidemiology, antimicrobial resistance, and transmission, it has been under-utilised in sub-Saharan Africa. In this study, we applied WGSA to large sub-sample of surveillance isolates from South Africa, collected from 2011 to 2015, focussing on Shigella flexneri 2a and Shigella sonnei. We find each serotype is epidemiologically distinct. The four identified S. flexneri 2a clusters having distinct geographical distributions, and antimicrobial resistance (AMR) and virulence profiles, while the four sub-Clades of S. sonnei varied in virulence plasmid retention. Our results support serotype specific lifestyles as a driver for epidemiological differences, show AMR is not required for epidemiological success in S. flexneri, and that the HIV epidemic may have promoted Shigella population expansion.

The evolution and international spread of extensively drug resistant Shigella sonnei.

Shigella sonnei causes shigellosis, a severe gastrointestinal illness that is sexually transmissible among men who have sex with men (MSM). Multidrug resistance in S. sonnei is common including against World Health Organisation recommended treatment options, azithromycin, and ciprofloxacin. Recently, an MSM-associated outbreak of extended-spectrum ß-lactamase producing, extensively drug resistant S. sonnei was reported in the United Kingdom. Here, we aimed to identify the genetic basis, evolutionary history, and international dissemination of the outbreak strain. Our genomic epidemiological analyses of 3,304 isolates from the United Kingdom, Australia, Belgium, France, and the United States of America revealed an internationally connected outbreak with a most recent common ancestor in 2018 carrying a low-fitness cost resistance plasmid, previously observed in travel associated sublineages of S. flexneri. Our results highlight the persistent threat of horizontally transmitted antimicrobial resistance and the value of continuing to work towards early and open international sharing of genomic surveillance data.

The genomic epidemiology of Escherichia albertii infecting humans and birds in Great Britain.

Escherichia albertii is a recently identified gastrointestinal bacterial pathogen of humans and animals which is typically misidentified as pathotypes of diarrhoeagenic Escherichia coli or Shigella species and is generally only detected during genomic surveillance of other Enterobacteriaceae. The incidence of E. albertii is likely underestimated, and its epidemiology and clinical relevance are poorly characterised. Here, we whole genome sequenced E. albertii isolates from humans (n = 83) and birds (n = 79) isolated in Great Britain between 2000 and 2021 and analysed these alongside a broader public dataset (n = 475) to address these gaps. We found human and avian isolates typically (90%; 148/164) belonged to host-associated monophyletic groups with distinct virulence and antimicrobial resistance profiles. Overlaid patient epidemiological data suggested that human infection was likely related to travel and possibly foodborne transmission. The Shiga toxin encoding stx2f gene was associated with clinical disease (OR = 10.27, 95% CI = 2.98-35.45 p = 0.0002) in finches. Our results suggest that improved future surveillance will further elucidate disease ecology and public and animal health risks associated with E. albertii.

Comparative genomics of Staphylococcus capitis reveals species determinants.

Staphylococcus capitis is primarily described as a human skin commensal but is now emergent as an opportunistic pathogen isolated from the bloodstream and prosthetic joint infections, and neonatal intensive care unit (NICU)-associated sepsis. We used comparative genomic analyses of S. capitis to provide new insights into commensal scalp isolates from varying skin states (healthy, dandruff lesional, and non-lesional), and to expand our current knowledge of the species populations (scalp isolates, n = 59; other skin isolates, n = 7; publicly available isolates, n = 120). A highly recombinogenic population structure was revealed, with genomes including the presence of a range of previously described staphylococcal virulence factors, cell wall-associated proteins, and two-component systems. Genomic differences between the two described S. capitis subspecies were explored, which revealed the determinants associated exclusively with each subspecies. The subspecies ureolyticus was distinguished from subspecies capitis based on the differences in antimicrobial resistance genes, ß-lactam resistance genes, and ß-class phenol soluble modulins and gene clusters linked to biofilm formation and survival on skin. This study will aid further research into the classification of S. capitis and virulence-linked phylogroups to monitor the spread and evolution of S. capitis.

A tale of two plasmids: contributions of plasmid associated phenotypes to epidemiological success among Shigella.

Dissemination of antimicrobial resistance (AMR) genes by horizontal gene transfer (HGT) mediated through plasmids is a major global concern. Genomic epidemiology studies have shown varying success of different AMR plasmids during outbreaks, but the underlying reasons for these differences are unclear. Here, we investigated two Shigella plasmids (pKSR100 and pAPR100) that circulated in the same transmission network but had starkly contrasting epidemiological outcomes to identify plasmid features that may have contributed to the differences. We used plasmid comparative genomics to reveal divergence between the two plasmids in genes encoding AMR, SOS response alleviation and conjugation. Experimental analyses revealed that these genomic differences corresponded with reduced conjugation efficiencies for the epidemiologically successful pKSR100, but more extensive AMR, reduced fitness costs, and a reduced SOS response in the presence of antimicrobials, compared with the less successful pAPR100. The discrepant phenotypes between the two plasmids are consistent with the hypothesis that plasmid-associated phenotypes contribute to determining the epidemiological outcome of AMR HGT and suggest that phenotypes relevant in responding to antimicrobial pressure and fitness impact may be more important than those around conjugation in this setting. Plasmid phenotypes could thus be valuable tools in conjunction with genomic epidemiology for predicting AMR dissemination.

Pathogenomic analyses of Shigella isolates inform factors limiting shigellosis prevention and control across LMICs.

Shigella spp. are the leading bacterial cause of severe childhood diarrhoea in low- and middle-income countries (LMICs), are increasingly antimicrobial resistant and have no widely available licenced vaccine. We performed genomic analyses of 1,246 systematically collected shigellae sampled from seven countries in sub-Saharan Africa and South Asia as part of the Global Enteric Multicenter Study (GEMS) between 2007 and 2011, to inform control and identify factors that could limit the effectiveness of current approaches. Through contemporaneous comparison among major subgroups, we found that S. sonnei contributes ≥6-fold more disease than other Shigella species relative to its genomic diversity, and highlight existing diversity and adaptative capacity among S. flexneri that may generate vaccine escape variants in <6 months. Furthermore, we show convergent evolution of resistance against ciprofloxacin, the current WHO-recommended antimicrobial for the treatment of shigellosis, among Shigella isolates. This demonstrates the urgent need to integrate existing genomic diversity into vaccine and treatment plans for Shigella, providing a framework for the focused application of comparative genomics to guide vaccine development, and the optimization of control and prevention strategies for other pathogens relevant to public health policy considerations.

Temporal, Spatial, and Genomic Analyses of Enterobacteriaceae Clinical Antimicrobial Resistance in Companion Animals Reveals Phenotypes and Genotypes of One Health Concern.

BACKGROUND: Antimicrobial resistance (AMR) is a globally important one health threat. The impact of resistant infections on companion animals, and the potential public health implications of such infections, has not been widely explored, largely due to an absence of structured population-level data. OBJECTIVES: We aimed to efficiently capture and repurpose antimicrobial susceptibility test (AST) results data from several veterinary diagnostic laboratories (VDLs) across the United Kingdom to facilitate national companion animal clinical AMR surveillance. We also sought to harness and genotypically characterize isolates of potential AMR importance from these laboratories. METHODS: We summarized AST results for 29,330 canine and 8,279 feline Enterobacteriaceae isolates originating from companion animal clinical practice, performed between April 2016 and July 2018 from four VDLs, with submissions from 2,237 United Kingdom veterinary practice sites. RESULTS: Escherichia coli (E. coli) was the most commonly isolated Enterobacteriaceae in dogs (69.4% of AST results, 95% confidence interval, CI, 68.7-70.0) and cats (90.5%, CI 89.8-91.3). Multi-drug resistance was reported in 14.1% (CI 13.5-14.8) of canine and 12.0% (CI 11.1-12.9) of feline E. coli isolates. Referral practices were associated with increased E. coli 3rd generation ≤ cephalosporin resistance odds (dogs: odds ratio 2.0, CI 1.2-3.4). We selected 95 E. coli isolates for whole genome analyses, of which seven belonged to sequence type 131, also carrying the plasmid-associated extended spectrum ß-lactamase gene bla CTX-M- 15. The plasmid-mediated colistin resistance gene mcr-9 was also identified for the first time in companion animals. CONCLUSIONS: Linking clinical AMR data with genotypic characterization represents an efficient means of identifying important resistance trends in companion animals on a national scale.

Whole genome sequence analysis of Shigella from Malawi identifies fluoroquinolone resistance.

Increasing antimicrobial resistance and limited alternative treatments have led to fluoroquinolone-resistant Shigella strain inclusion on the WHO global priority pathogens list. In this study we characterized multiple Shigella isolates from Malawi with whole genome sequence analysis, identifying the acquirable fluoroquinolone resistance determinant qnrS1.

Global population structure and genotyping framework for genomic surveillance of the major dysentery pathogen, Shigella sonnei.

Shigella sonnei is the most common agent of shigellosis in high-income countries, and causes a significant disease burden in low- and middle-income countries. Antimicrobial resistance is increasingly common in all settings. Whole genome sequencing (WGS) is increasingly utilised for S. sonnei outbreak investigation and surveillance, but comparison of data between studies and labs is challenging. Here, we present a genomic framework and genotyping scheme for S. sonnei to efficiently identify genotype and resistance determinants from WGS data. The scheme is implemented in the software package Mykrobe and tested on thousands of genomes. Applying this approach to analyse >4,000 S. sonnei isolates sequenced in public health labs in three countries identified several common genotypes associated with increased rates of ciprofloxacin resistance and azithromycin resistance, confirming intercontinental spread of highly-resistant S. sonnei clones and demonstrating the genomic framework can facilitate monitoring the spread of resistant clones, including those that have recently emerged, at local and global scales.

Accessory Genome Dynamics and Structural Variation of Shigella from Persistent Infections.

Shigellosis is a diarrheal disease caused mainly by Shigella flexneri and Shigella sonnei Infection is thought to be largely self-limiting, with short- to medium-term and serotype-specific immunity provided following clearance. However, cases of men who have sex with men (MSM)-associated shigellosis have been reported where Shigella of the same serotype were serially sampled from individuals between 1 and 1,862 days apart, possibly due to persistent carriage or reinfection with the same serotype. Here, we investigate the accessory genome dynamics of MSM-associated S. flexneri and S. sonnei isolates serially sampled from individual patients at various days apart to shed light on the adaptation of these important pathogens during infection. We find that pairs likely associated with persistent infection/carriage and with a smaller single nucleotide polymorphism (SNP) distance, demonstrated significantly less variation in accessory genome content than pairs likely associated with reinfection, and with a greater SNP distance. We observed antimicrobial resistance acquisition during Shigella carriage, including the gain of an extended-spectrum beta-lactamase gene during carriage. Finally, we explored large chromosomal structural variations and rearrangements in seven (five chronic and two reinfection associated) pairs of S. flexneri 3a isolates from an MSM-associated epidemic sublineage, which revealed variations at several common regions across isolate pairs, mediated by insertion sequence elements and comprising a distinct predicted functional profile. This study provides insight on the variation of accessory genome dynamics and large structural genomic changes in Shigella during persistent infection/carriage. In addition, we have also created a complete reference genome and biobanked isolate of the globally important pathogen, S. flexneri 3a.IMPORTANCEShigella spp. are Gram-negative bacteria that are the etiological agent of shigellosis, the second most common cause of diarrheal illness among children under the age of five in low-income countries. In high-income countries, shigellosis is also a sexually transmissible disease among men who have sex with men. Within the latter setting, we have captured prolonged and/or recurrent infection with shigellae of the same serotype, challenging the belief that Shigella infection is short lived and providing an early opportunity to study the evolution of the pathogen over the course of infection. Using this recently emerged transmission scenario, we comprehensively characterize the genomic changes that occur over the course of individual infection with Shigella and uncover a distinct functional profile of variable genomic regions, findings that have relevance for other Enterobacteriaceae.

Stepwise evolution of Salmonella Typhimurium ST313 causing bloodstream infection in Africa.

Bloodstream infections caused by nontyphoidal Salmonella are a major public health concern in Africa, causing ~49,600 deaths every year. The most common Salmonella enterica pathovariant associated with invasive nontyphoidal Salmonella disease is Salmonella Typhimurium sequence type (ST)313. It has been proposed that antimicrobial resistance and genome degradation has contributed to the success of ST313 lineages in Africa, but the evolutionary trajectory of such changes was unclear. Here, to define the evolutionary dynamics of ST313, we sub-sampled from two comprehensive collections of Salmonella isolates from African patients with bloodstream infections, spanning 1966 to 2018. The resulting 680 genome sequences led to the discovery of a pan-susceptible ST313 lineage (ST313 L3), which emerged in Malawi in 2016 and is closely related to ST313 variants that cause gastrointestinal disease in the United Kingdom and Brazil. Genomic analysis revealed degradation events in important virulence genes in ST313 L3, which had not occurred in other ST313 lineages. Despite arising only recently in the clinic, ST313 L3 is a phylogenetic intermediate between ST313 L1 and L2, with a characteristic accessory genome. Our in-depth genotypic and phenotypic characterization identifies the crucial loss-of-function genetic events that occurred during the stepwise evolution of invasive S. Typhimurium across Africa.

Metagenomic sequencing of clinical samples reveals a single widespread clone of Lawsonia intracellularis responsible for porcine proliferative enteropathy.

Lawsonia intracellularis is a Gram-negative obligate intracellular bacterium that is the aetiological agent of proliferative enteropathy (PE), a common intestinal disease of major economic importance in pigs and other animal species. To date, progress in understanding the biology of L. intracellularis for improved disease control has been hampered by the inability to culture the organism in vitro. In particular, our understanding of the genomic diversity and population structure of clinical L. intercellularis is very limited. Here, we utilized a metagenomic shotgun approach to directly sequence and assemble 21 L. intracellularis genomes from faecal and ileum samples of infected pigs and horses across three continents. Phylogenetic analysis revealed a genetically monomorphic clonal lineage responsible for infections in pigs, with distinct subtypes associated with infections in horses. The genome was highly conserved, with 94 % of genes shared by all isolates and a very small accessory genome made up of only 84 genes across all sequenced strains. In part, the accessory genome was represented by regions with a high density of SNPs, indicative of recombination events importing novel gene alleles. In summary, our analysis provides the first view of the population structure for L. intracellularis, revealing a single major lineage associated with disease of pigs. The limited diversity and broad geographical distribution suggest the recent emergence and clonal expansion of an important livestock pathogen.

Lawsonia intracellularis exploits ß-catenin/Wnt and Notch signalling pathways during infection of intestinal crypt to alter cell homeostasis and promote cell proliferation.

Lawsonia intracellularis is an obligate intracellular bacterial pathogen that causes proliferative enteropathy (PE) in pigs. L. intracellularis infection causes extensive intestinal crypt cell proliferation and inhibits secretory and absorptive cell differentiation. However, the affected host upstream cellular pathways leading to PE are still unknown. ß-catenin/Wnt signalling is essential in maintaining intestinal stem cell (ISC) proliferation and self-renewal capacity, while Notch signalling governs differentiation of secretory and absorptive lineage specification. Therefore, in this report we used immunofluorescence (IF) and quantitative reverse transcriptase PCR (RTqPCR) to examine ß-catenin/Wnt and Notch-1 signalling levels in uninfected and L. intracellularis infected pig ileums at 3, 7, 14, 21 and 28 days post challenge (dpc). We found that while the significant increase in Ki67+ nuclei in crypts at the peak of L. intracellularis infection suggested enhanced cell proliferation, the expression of c-MYC and ASCL2, promoters of cell growth and ISC proliferation respectively, was down-regulated. Peak infection also coincided with enhanced cytosolic and membrane-associated ß-catenin staining and induction of AXIN2 and SOX9 transcripts, both encoding negative regulators of ß-catenin/Wnt signalling and suggesting a potential alteration to ß-catenin/Wnt signalling levels, with differential regulation of the expression of its target genes. We found that induction of HES1 and OLFM4 and the down-regulation of ATOH1 transcript levels was consistent with the increased Notch-1 signalling in crypts at the peak of infection. Interestingly, the significant down-regulation of ATOH1 transcript levels coincided with the depletion of MUC2 expression at 14 dpc, consistent with the role of ATOH1 in promoting goblet cell maturation. The lack of significant change to LGR5 transcript levels at the peak of infection suggested that the crypt hyperplasia was not due to the expansion of ISC population. Overall, simultaneous induction of Notch-1 signalling and the attenuation of ß-catenin/Wnt pathway appear to be associated with the inhibition of goblet cell maturation and enhanced crypt cell proliferation at the peak of L. intracellularis infection. Moreover, the apparent differential regulation of apoptosis between crypt and lumen cells together with the strong induction of Notch-1 signalling and the enhanced SOX9 expression along crypts 14 dpc suggest an expansion of actively dividing transit amplifying and/or absorptive progenitor cells and provide a potential basis for understanding the development and maintenance of PE.

Lawsonia intracellularis infection of intestinal crypt cells is associated with specific depletion of secreted MUC2 in goblet cells.

The expression patterns of secreted (MUC2 and MUC5AC) and membrane-tethered (MUC1, MUC4, MUC12 and MUC13) mucins were monitored in healthy pigs and pigs challenged orally with Lawsonia intracellularis. These results showed that the regulation of mucin gene expression is distinctive along the GI tract of the healthy pig, and may reflect an association between the function of the mucin subtypes and different physiological demands at various sites. We identified a specific depletion of secreted MUC2 from goblet cells in infected pigs that correlated with the increased level of intracellular bacteria in crypt cells. We concluded that L. intracellularis may influence MUC2 production, thereby altering the mucus barrier and enabling cellular invasion.