Shigella sonnei genome sequencing and phylogenetic analysis indicate recent global dissemination from Europe

Investigation of the vaccine potential of an in silico designed FepA peptide vaccine against Shigella flexneri in mice model

Background

Shigellosis is one of the significant causes of diarrhea in Bangladesh. It is a global health problem; approximately 1.3 million people die yearly from Shigellosis. The current treatment method, using different antibiotics against Shigellosis is ineffective. Moreover, it becomes a worrying situation due to the emergence of antibiotic-resistant pathogenic microbes responsible for these diarrheal diseases.

Methodology

Previous immunoinformatics study predicted a potential peptide from the Ferric enterobactin protein (FepA) of Shigella spp. In this study, we have chemically synthesized the FepA peptide. As a highly immunogenic, FepA peptide conjugated with KLH has been tested in mice model with complete and incomplete adjuvants as a vaccine candidate.

Results

Immunological analysis showed that all vaccinated mice were immunologically boosted, which was statistically significant (P-value 0.0325) compared to control mice. Immunological analysis for bacterial neutralization test result was also statistically significant (P-value 0.0468), where each ELISA plate was coated with 1 × 107S. flexneri cells. The Challenge test with 1 × 1012S. flexneri cells to each vaccinated and controlled mice showed that 37.5 % of control (non-vaccinated) mice died within seven days after the challenge was given while 100 % of vaccinated mice remained strong and stout. The analyses of the post-challenge weight loss of the mice were also significant (P-value 0.0367) as the weight loss percentage in control mice was much higher than in the vaccinated mice. The pathological and phenotypic appearances of vaccinated mice were also clearly differentiable compared with control mice. Thus all these immunological analysis and pathological appearances directly supported our FepA peptide as a potential immune booster.

Conclusion

This study provides evidence that the FepA peptide is a highly immunogenic vaccine candidate against S. flexneri. Therefore, these findings inspire future trials for the evaluation of the suitability of this vaccine candidate against Shigellosis.

Adaptation of a fluoroquinolone-sensitive Shigella sonnei to norfloxacin exposure

Shigella causes shigellosis that requires antibiotic treatment in severe cases. Sublethal antibiotic concentrations can promote resistance, but their effect on antibiotic-sensitive bacteria before resistance development is unclear. This study investigated the effects of sublethal norfloxacin (NOR) challenges on a NOR-sensitive strain, Shigella sonnei UKMCC1015. Firstly, the whole genome of S. sonnei UKMCC1015 was assembled, and 45 antimicrobial resistance (AMR) genes were identified. Interestingly, transcriptomic analysis showed that low NOR levels do not change either the expression of the AMR genes or NOR targets such as gyrA. Instead, multiple ribosomal protein genes were downregulated, which could be attributed to decreased ribosomal protein promoter activity, modulated by elevated guanosine pentaphosphate and tetraphosphate (ppGpp) levels. This alarmone is involved in the bacterial stringent response during environmental stress, and it is mainly produced from the ppGpp synthetase (relA). Additionally, we observed that a relA overexpression (prolonged period of elevated ppGpp levels) may negatively affect the NOR tolerance of the bacteria. In conclusion, this study revealed that a NOR-sensitive strain responds differently to sublethal NOR than commonly reported in resistant strains.

Shigella sonnei utilises colicins during inter-bacterial competition

The mammalian colon is one of the most densely populated habitats currently recognised, with 1011-1013 commensal bacteria per gram of colonic contents. Enteric pathogens must compete with the resident intestinal microbiota to cause infection. Among these enteric pathogens are Shigella species which cause approximately 125 million infections annually, of which over 90 % are caused by Shigella flexneri and Shigella sonnei. Shigella sonnei was previously reported to use a Type VI Secretion System (T6SS) to outcompete E. coli and S. flexneri in in vitro and in vivo experiments. S. sonnei strains have also been reported to harbour colicinogenic plasmids, which are an alternative anti-bacterial mechanism that could provide a competitive advantage against the intestinal microbiota. We sought to determine the contribution of both T6SS and colicins to the anti-bacterial killing activity of S. sonnei. We reveal that whilst the T6SS operon is present in S. sonnei, there is evidence of functional degradation of the system through SNPs, indels and IS within key components of the system. We created strains with synthetically inducible T6SS operons but were still unable to demonstrate anti-bacterial activity of the T6SS. We demonstrate that the anti-bacterial activity observed in our in vitro assays was due to colicin activity. We show that S. sonnei no longer displayed anti-bacterial activity against bacteria that were resistant to colicins, and removal of the colicin plasmid from S. sonnei abrogated anti-bacterial activity of S. sonnei. We propose that the anti-bacterial activity demonstrated by colicins may be sufficient for niche competition by S. sonnei within the gastrointestinal environment.

Pathogenicity and virulence of Shigella sonnei: A highly drug-resistant pathogen of increasing prevalence

Shigella spp. are the causative agent of shigellosis (or bacillary dysentery), a diarrhoeal disease characterized for the bacterial invasion of gut epithelial cells. Among the 4 species included in the genus, Shigella flexneri is principally responsible for the disease in the developing world while Shigella sonnei is the main causative agent in high-income countries. Remarkably, as more countries improve their socioeconomic conditions, we observe an increase in the relative prevalence of S. sonnei. To date, the reasons behind this change in aetiology depending on economic growth are not understood. S. flexneri has been widely used as a model to study the pathogenesis of the genus, but as more research data are collected, important discrepancies with S. sonnei have come to light. In comparison to S. flexneri, S. sonnei can be differentiated in numerous aspects; it presents a characteristic O-antigen identical to that of one serogroup of the environmental bacterium Plesiomonas shigelloides, a group 4 capsule, antibacterial mechanisms to outcompete and displace gut commensal bacteria, and a poorer adaptation to an intracellular lifestyle. In addition, the World Health Organization (WHO) have recognized the significant threat posed by antibiotic-resistant strains of S. sonnei, demanding new approaches. This review gathers knowledge on what is known about S. sonnei within the context of other Shigella spp. and aims to open the door for future research on understanding the increasing spread of this pathogen.

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.

Escherichia coli killing by epidemiologically successful sublineages of Shigella sonnei is mediated by colicins

BACKGROUND

Shigella sp. are enteric pathogens which causes >125 million cases of shigellosis annually. S. sonnei accounts for about a quarter of those cases and is increasingly prevalent in industrialising nations. Being an enteric pathogen, S. sonnei benefits from outcompeting gut commensals such as Escherichia coli to establish itself and cause disease. There are numerous mechanisms that bacterial pathogens use to outcompete its rivals including molecules called colicins. A Type 6 Secretion System (T6SS) was recently described as contributing to E. coli killing in S. sonnei.

METHODS

We used Bulk Phenotyping of Epidemiological Replicates (BPER) which combined bacterial Genome Wide Association Studies (bGWAS) and high throughput phenotyping on a collection of S. sonnei surveillance isolates to identify the genetic features associated with E. coli killing and explore their relationship with epidemiological behaviour. We further explored the presence of colicins and T6SS components in the isolates using genomics, laboratory experimentation, and proteomics.

FINDINGS

Our bGWAS analysis returned known and novel colicin and colicin related genes as significantly associated with E. coli killing. In silico analyses identified key colicin clusters responsible for the killing phenotype associated with epidemiologically successful sub-lineages. The killing phenotype was not associated with the presence of a T6SS. Laboratory analyses confirmed the presence of the key colicin clusters and that killing was contact-independent.

INTERPRETATION

Colicins are responsible for E. coli killing by S. sonnei, not a T6SS. This phenotype contributes to shaping the observed epidemiology of S. sonnei and may contribute to its increasing prevalence globally. BPER is an epidemiologically relevant approach to phenotypic testing that enables the rapid identification of genetic drivers of phenotypic changes, and assessment of their relevance to epidemiology in natural settings.

FUNDING

Biotechnology and Biological Sciences Research Council, Biotechnology and Biological Sciences Research Council Doctoral Training Partnership studentship, Wellcome Trust, Medical Research Council (UK), French National Research Agency.

Shigella Serotypes Associated With Carriage in Humans Establish Persistent Infection in Zebrafish

Shigella represents a paraphyletic group of enteroinvasive Escherichia coli. More than 40 Shigella serotypes have been reported. However, most cases within the men who have sex with men (MSM) community are attributed to 3 serotypes: Shigella sonnei unique serotype and Shigella flexneri 2a and 3a serotypes. Using the zebrafish model, we demonstrate that Shigella can establish persistent infection in vivo. Bacteria are not cleared by the immune system and become antibiotic tolerant. Establishment of persistent infection depends on the O-antigen, a key constituent of the bacterial surface and a serotype determinant. Representative isolates associated with MSM transmission persist in zebrafish, while representative isolates of a serotype not associated with MSM transmission do not. Isolates of a Shigella serotype establishing persistent infections elicited significantly less macrophage death in vivo than isolates of a serotype unable to persist. We conclude that zebrafish are a valuable platform to illuminate factors underlying establishment of Shigella persistent infection in humans.

A Tale about Shigella: Evolution, Plasmid, and Virulence

Shigella spp. cause hundreds of millions of intestinal infections each year. They target the mucosa of the human colon and are an important model of intracellular bacterial pathogenesis. Shigella is a pathovar of Escherichia coli that is characterized by the presence of a large invasion plasmid, pINV, which encodes the characteristic type III secretion system and icsA used for cytosol invasion and cell-to-cell spread, respectively. First, we review recent advances in the genetic aspects of Shigella, shedding light on its evolutionary history within the E. coli lineage and its relationship to the acquisition of pINV. We then discuss recent insights into the processes that allow for the maintenance of pINV. Finally, we describe the role of the transcription activators VirF, VirB, and MxiE in the major virulence gene regulatory cascades that control the expression of the type III secretion system and icsA. This provides an opportunity to examine the interplay between these pINV-encoded transcriptional activators and numerous chromosome-encoded factors that modulate their activity. Finally, we discuss novel chromosomal genes icaR, icaT, and yccE that are regulated by MxiE. This review emphasizes the notion that Shigella and E. coli have walked the fine line between commensalism and pathogenesis for much of their history.

Concordance between Genotypic and Phenotypic Drug-Resistant Profiles of Shigella Isolates from Taiyuan City, Shanxi Province, China, 2005 to 2016

Antimicrobial resistance in Shigella spp. is a global public health concern. In this study, the AMR phenotypic profiles of 10 kinds of antibiotics were compared with the genotypic profiles using genomic analysis of 218 Shigella isolates from Taiyuan City, Shanxi Province, China, 2005 to 2016. Core genome Multilocus Sequence Typing (cgMLST) based on the EnteroBase Escherichia/Shigella scheme was used to obtain the genetic relatedness of Shigella isolates. Multiple-drug resistance was observed in 96.79% Shigella spp., and the resistance to antimicrobial agents varied between S. flexneri and S. sonnei. The genotypic results correlated well with the phenotypic profiles with concordance rates of 96.42% and 94.50% in S. flexneri and S. sonnei isolates, respectively, from Taiyuan City, Shanxi Province. The sensitivity and specificity of the genotypic antimicrobial susceptibility testing (AST) were 97.56% and 95.34% for S. flexneri, and 95.65% and 93.31% for S. sonnei isolates, respectively. A discrepancy of genotypic and phenotypic AST results existed in some cephalosporin- and azithromycin-resistant Shigella isolates; there were no clear resistance patterns to predict ciprofloxacin resistance. There were major discrepancies between genotypic and phenotypic AST in the genotypically resistant but phenotypically susceptible isolates. The drug-resistance patterns and essential drug-resistance genes to predict the phenotypic drug-resistant profiles were the discrepancies between S. flexneri and S. sonnei isolates. Phylogenetic analysis showed that isolates of the same cluster but with different antibiotic-resistance gene patterns occurred because of the loss or gain of antibiotic-resistance genes located in the plasmids and multidrug-resistance islands. IMPORTANCE Antimicrobial resistance in Shigella spp. has become a global public health concern. In this study, we identified the antimicrobial susceptibility testing (AST) characteristics based on genomic sequences of 218 Shigella isolates and analyzed the correlation between genotypic and phenotypic antibiotic resistance profiles of Shigella spp., especially for fluoroquinolone, macrolides, and third-generation cephalosporins. Our results show that the genotypic results correlated with the phenotypic profiles with concordance rates of 96.42% and 94.50% in S. flexneri and S. sonnei isolates, respectively. The drug-resistance patterns and essential drug-resistance genes to predict the phenotypic drug-resistant profiles of S. flexneri and S. sonnei isolates in Taiyuan city were distinct. The discrepancy between genotypic and phenotypic AST was considerable in the genotypically resistant but phenotypically susceptible isolates. The information on drug resistance and resistance genes in this study can offer more details on the prevalence of drug resistance of Shigella spp.

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.

Shigellosis in Southeast Asia: A systematic review and meta-analysis

BACKGROUND

Southeast Asia is attractive for tourism. Unfortunately, travelers to this region are at risk of becoming infected with Shigella. We conducted a meta-analysis to provide updates on Shigella prevalence in Southeast Asia, along with their serogroups and serotypes.

METHODS

We conducted a systematic search using PubMed, EMBASE, and Web of Science for peer-reviewed studies from 2000 to November 2022. We selected studies that detected Shigella in stools by culture or polymerase chain reaction (PCR). Two reviewers extracted the data using a standardized form and performed quality assessments using the Joanna Briggs Institute checklist. The random effects model was used to estimate the pooled prevalence of Shigella.

RESULTS

During our search, we identified 4376 studies. 29 studies (from six Southeast Asian countries) were included in the systematic review, 21 each in the meta-analysis of the prevalence of Shigella (Sample size: 109545) and the prevalence of Shigella serogroups. The pooled prevalence of Shigella was 4% (95% CI: 4-5%) among diarrhea cases. Shigella sonnei was the most abundant serogroup in Thailand (74%) and Vietnam (57%), whereas Shigella flexneri was dominant in Indonesia (72%) and Cambodia (71%). Shigella dysenteriae and Shigella boydii were uncommon (pooled prevalence of 1% each). The pooled prevalence of Shigella was 5% (95% CI: 4-6%) in children aged <5 years. The pooled prevalence showed a decreasing trend comparing data collected between 2000-2013 (5%; 95% CI: 4-6%) and between 2014-2022 (3%; 95% CI: 2-4%). Shigella prevalence was 6% in studies that included participants with mixed pathogens versus 3% in those without. Shigella flexneri serotype 2a was the most frequently isolated (33%), followed by 3a (21%), 1b (10%), 2b (3%), and 6 (3%).

CONCLUSIONS

This study provides compelling evidence for the development of effective Shigella vaccines for residents of endemic regions and travellers to these areas.

Rapid emergence of extensively drug-resistant Shigella sonnei in France

Shigella sonnei, the main cause of bacillary dysentery in high-income countries, has become increasingly resistant to antibiotics. We monitored the antimicrobial susceptibility of 7121 S. sonnei isolates collected in France between 2005 and 2021. We detected a dramatic increase in the proportion of isolates simultaneously resistant to ciprofloxacin (CIP), third-generation cephalosporins (3GCs) and azithromycin (AZM) from 2015. Our genomic analysis of 164 such extensively drug-resistant (XDR) isolates identified 13 different clusters within CIP-resistant sublineage 3.6.1, which was selected in South Asia ∼15 years ago. AZM resistance was subsequently acquired, principally through IncFII (pKSR100-like) plasmids. The last step in the development of the XDR phenotype involved various extended-spectrum beta-lactamase genes (bla_CTX-M-3, bla_CTX-M-15, bla_CTX-M-27, bla_CTX-M-55, and bla_CTX-M-134) carried by different plasmids (IncFII, IncI1, IncB/O/K/Z) or even integrated into the chromosome, and encoding resistance to 3GCs. This rapid emergence of XDR S. sonnei, including an international epidemic strain, is alarming, and good laboratory-based surveillance of shigellosis will be crucial for informed decision-making and appropriate public health action.

Using a combination of short- and long-read sequencing to investigate the diversity in plasmid- and chromosomally encoded extended-spectrum beta-lactamases (ESBLs) in clinical Shigella and Salmonella isolates in Belgium

For antimicrobial resistance (AMR) surveillance, it is important not only to detect AMR genes, but also to determine their plasmidic or chromosomal location, as this will impact their spread differently. Whole-genome sequencing (WGS) is increasingly used for AMR surveillance. However, determining the genetic context of AMR genes using only short-read sequencing is complicated. The combination with long-read sequencing offers a potential solution, as it allows hybrid assemblies. Nevertheless, its use in surveillance has so far been limited. This study aimed to demonstrate its added value for AMR surveillance based on a case study of extended-spectrum beta-lactamases (ESBLs). ESBL genes have been reported to occur also on plasmids. To gain insight into the diversity and genetic context of ESBL genes detected in clinical isolates received by the Belgian National Reference Center between 2013 and 2018, 100 ESBL-producing Shigella and 31 ESBL-producing Salmonella were sequenced with MiSeq and a representative selection of 20 Shigella and six Salmonella isolates additionally with MinION technology, allowing hybrid assembly. The bla _CTX-M-15 gene was found to be responsible for a rapid rise in the ESBL Shigella phenotype from 2017. This gene was mostly detected on multi-resistance-carrying IncFII plasmids. Based on clustering, these plasmids were determined to be distinct from the circulating plasmids before 2017. They were spread to different Shigella species and within Shigella sonnei between multiple genotypes. Another similar IncFII plasmid was detected after 2017 containing bla _CTX-M-27 for which only clonal expansion occurred. Matches of up to 99 % to plasmids of various bacterial hosts from all over the world were found, but global alignments indicated that direct or recent ESBL-plasmid transfers did not occur. It is most likely that travellers introduced these in Belgium and subsequently spread them domestically. However, a clear link to a specific country could not be made. Moreover, integration of bla _CTX-M in the chromosome of two Shigella isolates was determined for the first time, and shown to be related to ISEcp1. In contrast, in Salmonella, ESBL genes were only found on plasmids, of which bla _CTX-M-55 and IncHI2 were the most prevalent, respectively. No matching ESBL plasmids or cassettes were detected between clinical Shigella and Salmonella isolates. The hybrid assembly data allowed us to check the accuracy of plasmid prediction tools. MOB-suite showed the highest accuracy. However, these tools cannot replace the accuracy of long-read and hybrid assemblies. This study illustrates the added value of hybrid assemblies for AMR surveillance and shows that a strategy where even just representative isolates of a collection used for hybrid assemblies could improve international AMR surveillance as it allows plasmid tracking.

Characterization of MxiE- and H-NS-Dependent Expression of ipaH7.8, ospC1, yccE, and yfdF in Shigella flexneri

Shigella flexneri uses a type 3 secretion system (T3SS) apparatus to inject virulence effector proteins into the host cell cytosol. Upon host cell contact, MxiE, an S. flexneri AraC-like transcriptional regulator, is required for the expression of a subset of T3SS effector genes encoded on the large virulence plasmid. Here, we defined the MxiE regulon using RNA-seq. We identified virulence plasmid- and chromosome-encoded genes that are activated in response to type 3 secretion in a MxiE-dependent manner. Bioinformatic analysis revealed that similar to previously known MxiE-dependent genes, chromosome-encoded genes yccE and yfdF contain a regulatory element known as the MxiE box, which is required for their MxiE-dependent expression. The significant AT enrichment of MxiE-dependent genes suggested the involvement of H-NS. Using a dominant negative H-NS system, we demonstrate that H-NS silences the expression of MxiE-dependent genes located on the virulence plasmid (ipaH7.8 and ospC1) and the chromosome (yccE and yfdF). Furthermore, we show that MxiE is no longer required for the expression of ipaH7.8, ospC1, yccE, and yfdF when H-NS silencing is relieved. Finally, we show that the H-NS anti-silencer VirB is not required for ipaH7.8 and yccE expression upon MxiE/IpgC overexpression. Based on these genetic studies, we propose a model of MxiE-dependent gene regulation in which MxiE counteracts H-NS-mediated silencing. IMPORTANCE The expression of horizontally acquired genes, including virulence genes, is subject to complex regulation involving xenogeneic silencing proteins, and counter-silencing mechanisms. The pathogenic properties of Shigella flexneri mainly rely on the acquisition of the type 3 secretion system (T3SS) and cognate effector proteins, whose expression is repressed by the xenogeneic silencing protein H-NS. Based on previous studies, releasing H-NS-mediated silencing mainly relies on two mechanisms involving (i) a temperature shift leading to the release of H-NS at the virF promoter, and (ii) the virulence factor VirB, which dislodges H-NS upon binding to specific motifs upstream of virulence genes, including those encoding the T3SS. In this study, we provide genetic evidence supporting the notion that, in addition to VirB, the AraC family member MxiE also contributes to releasing H-NS-mediated silencing in S. flexneri.

A Shigella sonnei clone with extensive drug resistance associated with waterborne outbreaks in China

Antimicrobial resistance of Shigella sonnei has become a global concern. Here, we report a phylogenetic group of S. sonnei with extensive drug resistance, including a combination of multidrug resistance, coresistance to ceftriaxone and azithromycin (cef^Razi^R), reduced susceptibility to fluoroquinolones, and even colistin resistance (col^R). This distinct clone caused six waterborne shigellosis outbreaks in China from 2015 to 2020. We collect 155 outbreak isolates and 152 sporadic isolates. The cef^Razi^R isolates, including outbreak strains, are mainly distributed in a distinct clade located in global Lineage III. The outbreak strains form a recently derived monophyletic group that may have emerged circa 2010. The cef^Razi^R and col^R phenotypes are attributed to the acquisition of different plasmids, particularly the IncB/O/K/Z plasmid coharboring the bla_CTX-M-14, mphA, aac(3)-IId, dfrA17, aadA5, and sul1 genes and the IncI2 plasmid with an mcr-1 gene. Genetic analyses identify 92 accessory genes and 60 single-nucleotide polymorphisms associated with the cef^Razi^R phenotype. Surveillance of this clone is required to determine its dissemination and threat to global public health.

Antimicrobial resistance in shigellosis: A surveillance study among urban and rural children over 20 years in Bangladesh

Antimicrobial resistance against shigellosis is increasingly alarming. However, evidence-based knowledge gaps regarding the changing trends of shigellosis in Bangladesh exist due to the scarcity of longitudinal data on antimicrobial resistance. Our study evaluated the last 20 years antimicrobial resistance patterns against shigellosis among under-5 children in the urban and rural sites of Bangladesh. Data were extracted from the Diarrheal Disease Surveillance System (DDSS) of Dhaka Hospital (urban site) and Matlab Hospital (rural site) of the International Centre for Diarrheal Disease Research, Bangladesh (icddr,b) between January 2001 and December 2020. We studied culture-confirmed shigellosis cases from urban Dhaka Hospital (n = 883) and rural Matlab Hospital (n = 1263). Since 2001, a declining percentage of shigellosis in children observed in urban and rural sites. Moreover, higher isolation rates of Shigella were found in the rural site [1263/15684 (8.1%)] compared to the urban site [883/26804 (3.3%)] in the last 20 years. In both areas, S. flexneri was the predominant species. The upward trend of S. sonnei in both the study sites was statistically significant after adjusting for age and sex. WHO-recommended 1st line antibiotic ciprofloxacin resistance gradually reached more than 70% in both the urban and rural site by 2020. In multiple logistic regression after adjusting for age and sex, ciprofloxacin, azithromycin, mecillinam, ceftriaxone, and multidrug resistance (resistance to any two of these four drugs) among under-5 children were found to be increasing significantly (p<0.01) in the last 20 years in both sites. The study results underscore the importance of therapeutic interventions for shigellosis by appropriate drugs based on their current antibiogram for under-5 children. These observations may help policymakers in formulating better case management strategies for shigellosis.

The international and intercontinental spread and expansion of antimicrobial-resistant Salmonella Typhi: a genomic epidemiology study

Background

The emergence of increasingly antimicrobial-resistant Salmonella enterica serovar Typhi (S Typhi) threatens to undermine effective treatment and control. Understanding where antimicrobial resistance in S Typhi is emerging and spreading is crucial towards formulating effective control strategies.

Methods

In this genomic epidemiology study, we sequenced the genomes of 3489 S Typhi strains isolated from prospective enteric fever surveillance studies in Nepal, Bangladesh, Pakistan, and India (between 2014 and 2019), and combined these with a global collection of 4169 S Typhi genome sequences isolated between 1905 and 2018 to investigate the temporal and geographical patterns of emergence and spread of antimicrobial-resistant S Typhi. We performed non-parametric phylodynamic analyses to characterise changes in the effective population size of fluoroquinolone-resistant, extensively drug-resistant (XDR), and azithromycin-resistant S Typhi over time. We inferred timed phylogenies for the major S Typhi sublineages and used ancestral state reconstruction methods to estimate the frequency and timing of international and intercontinental transfers.

Findings

Our analysis revealed a declining trend of multidrug resistant typhoid in south Asia, except for Pakistan, where XDR S Typhi emerged in 2016 and rapidly replaced less-resistant strains. Mutations in the quinolone-resistance determining region (QRDR) of S Typhi have independently arisen and propagated on at least 94 occasions, nearly all occurring in south Asia. Strains with multiple QRDR mutations, including triple mutants with high-level fluoroquinolone resistance, have been increasing in frequency and displacing strains with fewer mutations. Strains containing acrB mutations, conferring azithromycin resistance, emerged in Bangladesh around 2013 and effective population size of these strains has been steadily increasing. We found evidence of frequent international (n=138) and intercontinental transfers (n=59) of antimicrobial-resistant S Typhi, followed by local expansion and replacement of drug-susceptible clades.

Interpretation

Independent acquisition of plasmids and homoplastic mutations conferring antimicrobial resistance have occurred repeatedly in multiple lineages of S Typhi, predominantly arising in south Asia before spreading to other regions.

Funding

Bill & Melinda Gates Foundation.

Pyroptosis in host defence against bacterial infection

Pyroptosis, a regulated form of pro-inflammatory cell death, is characterised by cell lysis and by the release of cytokines, damage- and pathogen-associated molecular patterns. It plays an important role during bacterial infection, where it can promote an inflammatory response and eliminate the replicative niche of intracellular pathogens. Recent work, using a variety of bacterial pathogens, has illuminated the versatility of pyroptosis, revealing unexpected and important concepts underlying host defence. In this Review, we overview the molecular mechanisms underlying pyroptosis and discuss their role in host defence, from the single cell to the whole organism. We focus on recent studies using three cellular microbiology paradigms - Mycobacterium tuberculosis, Salmonella Typhimurium and Shigella flexneri - that have transformed the field of pyroptosis. We compare insights discovered in tissue culture, zebrafish and mouse models, highlighting the advantages and disadvantages of using these complementary infection models to investigate pyroptosis and for modelling human infection. Moving forward, we propose that in-depth knowledge of pyroptosis obtained from complementary infection models can better inform future studies using higher vertebrates, including humans, and help develop innovative host-directed therapies to combat bacterial infection.

Virulence plasmid pINV as a genetic signature for Shigella flexneri phylogeny

Shigella flexneri is a major health burden in low- and middle-income countries, where it is a leading cause of mortality associated with diarrhoea in children, and shows an increasing incidence among travellers and men having sex with men. Like all Shigella spp., S. flexneri has evolved from commensal Escherichia coli following the acquisition of a large plasmid pINV, which contains genes essential for virulence. Current sequence typing schemes of Shigella are based on combinations of chromosomal genetic loci, since pINV-encoded virulence genes are often lost during growth in the laboratory, making these elements inappropriate for sequence typing. By performing comparative analysis of pINVs from S. flexneri strains isolated from different geographical regions and belonging to different serotypes, we found that in contrast to plasmid-encoded virulence genes, plasmid maintenance genes are highly stable pINV-encoded elements. For the first time, to our knowledge, we have developed a S. flexneri plasmid multilocus sequence typing (pMLST) method based on different combinations of alleles of the vapBC and yacAB toxin–antitoxin (TA) systems, and the parAB partitioning system. This enables typing of S. flexneri pINV plasmids into distinct ‘virulence sequence types’ (vSTs). Furthermore, the phylogenies of vST alleles and bacterial host core genomes suggests an intimate co-evolution of pINV with the chromosome of its bacterial host, consistent with previous findings. This work demonstrates the potential of plasmid maintenance loci as genetic characteristics to study as well as to trace the molecular phylogenesis of S. flexneri pINV and the phylogenetic relationship of this plasmid with its bacterial host.

Phylogenetic analysis and antibiotic resistance of Shigella sonnei isolates

Shigellosis is one of the most important gastric infections caused by different species of Shigella and has been regarded as a serious threat to public health. Lineage/sublineage profile of S. sonnei is strongly associated with the antibiotic resistance and population structure of this pathogen. In this study, we determined the phylogeny and antibiotic resistance profiles of S. sonnei strains, isolated from 1246 stool and 580 food samples, using multiplex PCR-HRMA genotyping and Kirby-Bauer disk diffusion methods, respectively. A total of 64 S. sonnei strains were isolated (13 food and 51 clinical isolates). Multiplex PCR-HMR assay was able to differentiate the lineages II and III, and sublineages IIIb and IIIc strains successfully considering the definite melting curves and temperatures. Lineage I and sublineage IIIa strain were not isolated in this study. We also demonstrated that most of the S. sonnei strains isolated from both food and clinical samples clustered within the lineage III and sublineage IIIc. Resistance against trimethoprim-sulfamethoxazole, tetracycline, chloramphenicol and streptomycin antibiotics were the most prevalent phenotypes among the S. sonnei lineage III and sublineage IIIc strains.

A fatal outbreak of neonatal sepsis caused by mcr-10 carrying Enterobacter kobei in a tertiary care hospital in Nepal

BACKGROUND

Enterobacter kobei is an emerging cause of outbreak of nosocomial infections in neonatal intensive care units (NICUs). Between July and September of 2016, an NICU in a tertiary care hospital of Nepal observed an abrupt increase in the number of neonatal sepsis cases caused by Enterobacter spp. infecting 11 of 23 admitted neonates, 5 of whom died of an exacerbated sepsis.

AIM

Main aims of this study were to confirm the suspected outbreak, identify environmental source of infection, and characterize genetic determinants of antimicrobial resistance (AMR) and virulence of the pathogen.

METHODS

We performed whole genome sequencing of all Enterobacter spp. isolated from blood cultures of septic neonates admitted to NICU between May 2016 and December 2017. Also, an environmental sampling was intensified from fortnightly to weekly during the outbreak.

FINDINGS

The genomic analysis revealed that 10 of 11 non-duplicated E. kobei isolated from neonatal blood cultures between July and September 2016 were clonal, confirming the outbreak. The isolates carried AMR genes including bla_AmpC and mcr-10 conferring reduced susceptibility to carbapenem and colistin respectively. The environmental sampling however failed to isolate any Enterobacter spp. Reinforcement of aseptic protocols in invasive procedures, hand hygiene, environmental decontamination, fumigation, and secluded care of culture positive cases successfully terminated the outbreak.

CONCLUSION

Our study underscored the need to implement stringent infection control measures to prevent infection outbreaks. Further, for the first time, we report the emergence of carbapenem and colistin non-susceptible E. kobei carrying mcr-10 gene as an important cause of nosocomial neonatal sepsis in an NICU.

Rates of evolutionary change of resident Escherichia coli O157:H7 differ within the same ecological niche

Background

Shiga toxin-producing Escherichia coli (STEC) O157:H7 is a pathogen known to reside in cattle feedlots. This retrospective study examined 181 STEC O157:H7 strains collected over 23 years from a closed-system feedlot. All strains were subjected to short-read sequencing, with a subset of 36 also subjected to long-read sequencing.

Results

Over 96% of the strains fell into four phylogenetically distinct clades. Clade membership was associated with multiple factors including stx composition and the alleles of a well-characterized polymorphism (tir 255 T > A). Small plasmids (2.7 to 40 kb) were found to be primarily clade specific. Within each clade, chromosomal rearrangements were observed along with a core phageome and clade specific phages. Across both core and mobile elements of the genome, multiple SNP alleles were in complete linkage disequilibrium across all strains within specific clades. Clade evolutionary rates varied between 0.9 and 2.8 SNP/genome/year with two tir A allele clades having the lowest evolutionary rates. Investigation into possible causes of the differing rates was not conclusive but revealed a synonymous based mutation in the DNA polymerase III of the fastest evolving clade. Phylogenetic trees generated through our bioinformatic pipeline versus the NCBI’s pathogen detection project were similar, with the two tir A allele clades matching individual NCBI SNP clusters, and the two tir T allele clades assigned to multiple closely-related SNP clusters.

Conclusions

In one ecological niche, a diverse STEC O157:H7 population exhibited different rates of evolution that associated with SNP alleles in linkage disequilibrium in the core genome and mobile elements, including tir 255 T > A.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08497-6.

Long-term spatiotemporal variation of antimicrobial resistance genes within the Serratia marcescens population and transmission of S. marcescens revealed by public whole-genome datasets

The development of antimicrobial resistance (AMR) is accelerated by the selective pressure exerted by the widespread use of antimicrobial drugs, posing an increasing danger to public health. However, long-term spatiotemporal variation in AMR genes in microorganisms, particularly in bacterial pathogens in response to antibiotic consumption, is not fully understood. Here, we used the NCBI RefSeq database to collect 478 whole-genome sequences for Serratia marcescens ranging from 1961 up to 2019, to document global long-term AMR trends in S. marcescens populations. In total, 100 AMR gene subtypes (16 AMR gene types) were detected in the genomes of S. marcescens populations. We identified 3 core resistance genes in S. marcescens genomes, and a high diversity of AMR genes was observed in S. marcescens genomes after corresponding antibiotics were discovered and introduced into clinical practice, suggesting the adaptation of S. marcescens populations to challenges with therapeutic antibiotics. Our findings indicate spatiotemporal variation of AMR genes in S. marcescens populations in relation to antibiotic consumption and suggest the potential transmission of S. marcescens isolates harboring AMR genes among countries and between the environment and the clinic, representing a public health threat that necessitates international solidarity to overcome.

Proteome based mapping and molecular docking revealed DnaA as a potential drug target against Shigella sonnei

Shigella sonnei is one of the major causes of diarrhea and remained a critical microbe responsible for higher morbidity and mortality rates resulting from dysentery every year across the world. Antibiotic therapy of Shigella diseases plays a critical role in decreasing the prevalence as well as the fatality rate of this infection. However, the management of these diseases remains challenging, owing to the overall increase in resistance against many antimicrobials. The situation necessitates the rapid development of effective and feasible S. sonnei treatments. In the present study, the subtractive genomics approach was utilized to find the potential drug targets for S. sonnei strain Ss046. Various tools of bioinformatics were implemented to remove the human-specific homologous and pathogen-specific paralogous sequences from the bacterial proteome. Then, metabolic pathway and subcellular location analysis were performed of essential bacterial proteins to describe their role in various cellular processes. Only one essential protein i-e Chromosomal replication initiator protein DnaA was found in the proteome of the pathogen that could be used as a potent target for designing new drugs. 3D structure prediction of DnaA protein was carried out using Phyre 2. Molecular docking of 5000 phytochemicals was performed against DnaA to identify four top-ranked phytochemicals (Riccionidin A, Dothistromin, Fustin, and Morin) based on scoring functions and interaction with the active site. This study suggests that these phytochemicals could be used as antibacterial drugs to treat S. sonnei infections in the future. To confirm their efficacy and evaluate their drug potency, further in vitro analyses are required.

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.

Population structure analysis and laboratory monitoring of Shigella by core-genome multilocus sequence typing

The laboratory surveillance of bacillary dysentery is based on a standardised Shigella typing scheme that classifies Shigella strains into four serogroups and more than 50 serotypes on the basis of biochemical tests and lipopolysaccharide O-antigen serotyping. Real-time genomic surveillance of Shigella infections has been implemented in several countries, but without the use of a standardised typing scheme. Here, we study over 4000 reference strains and clinical isolates of Shigella, covering all serotypes, with both the current serotyping scheme and the standardised EnteroBase core-genome multilocus sequence typing scheme (cgMLST). The Shigella genomes are grouped into eight phylogenetically distinct clusters, within the E. coli species. The cgMLST hierarchical clustering (HC) analysis at different levels of resolution (HC2000 to HC400) recognises the natural population structure of Shigella. By contrast, the serotyping scheme is affected by horizontal gene transfer, leading to a conflation of genetically unrelated Shigella strains and a separation of genetically related strains. The use of this cgMLST scheme will facilitate the transition from traditional phenotypic typing to routine whole-genome sequencing for the laboratory surveillance of Shigella infections.

Lab-based surveillance of Shigella has traditionally been based on serotyping but increasing availability of whole genome sequencing could enable higher resolution typing. Here, the authors apply a core genome multilocus sequence typing scheme to Shigella sequence data and describe its population structure.

Maintenance of the Shigella sonnei virulence plasmid is dependent on its repertoire and amino acid sequence of toxin:antitoxin systems

Shigella sonnei is a major cause of bacillary dysentery and an increasing concern due to the spread of multidrug resistance. S. sonnei harbors pINV, an ∼210 kb plasmid that encodes a type III secretion system (T3SS), which is essential for virulence. During growth in the laboratory, avirulence arises spontaneously in S. sonnei at high frequency, hampering studies on and vaccine development against this important pathogen. Here, we investigated the molecular basis for the emergence of avirulence in S. sonnei and showed that avirulence mainly results from pINV loss, which is consistent with previous findings. Ancestral deletions have led to the loss from S. sonnei pINV of two toxin-antitoxin (TA) systems involved in plasmid maintenance, CcdAB and GmvAT, which are found on pINV in Shigella flexneri. We showed that the introduction of these TA systems into S. sonnei pINV reduced but did not eliminate pINV loss, while the single amino acid polymorphisms found in the S. sonnei VapBC TA system compared with S. flexneri VapBC also contributed to pINV loss. Avirulence also resulted from deletions of T3SS-associated genes in pINV through recombination between insertion sequences (ISs) on the plasmid. These events differed from those observed in S. flexneri due to the different distribution and repertoire of ISs. Our findings demonstrated that TA systems and ISs influenced plasmid dynamics and loss in S. sonnei and could be exploited for the design and evaluation of vaccines.

IMPORTANCEShigella sonnei is the major cause of shigellosis in high-income and industrializing countries and is an emerging, multidrug-resistant pathogen. A significant challenge when studying this bacterium is that it spontaneously becomes avirulent during growth in the laboratory through loss of its virulence plasmid (pINV). Here, we deciphered the mechanisms leading to avirulence in S. sonnei and how the limited repertoire and amino acid sequences of plasmid-encoded toxin-antitoxin (TA) systems make the maintenance of pINV in this bacterium less efficient compared with Shigella flexneri. Our findings highlighted how subtle differences in plasmids in closely related species have marked effects and could be exploited to reduce plasmid loss in S. sonnei. This should facilitate research on this bacterium and vaccine development.

Genomic insights into the circulation of pandemic fluoroquinolone-resistant extra-intestinal pathogenic Escherichia coli ST1193 in Vietnam

Extra-intestinal pathogenic Escherichia coli (ExPEC) ST1193, a globally emergent fluoroquinolone-resistant clone, has become an important cause of bloodstream infections (BSIs) associated with significant morbidity and mortality. Previous studies have reported the emergence of fluoroquinolone-resistant ExPEC ST1193 in Vietnam; however, limited data exist regarding the genetic structure, antimicrobial resistance (AMR) determinants and transmission dynamics of this pandemic clone. Here, we performed genomic and phylogenetic analyses of 46 ST1193 isolates obtained from BSIs and healthy individuals in Ho Chi Minh City, Vietnam, to investigate the pathogen population structure, molecular mechanisms of AMR and potential transmission patterns. We further examined the phylogenetic structure of ST1193 isolates in a global context. We found that the endemic E. coli ST1193 population was heterogeneous and highly dynamic, largely driven by multiple strain importations. Several well-supported phylogenetic clusters (C1-C6) were identified and associated with distinct bla _CTX-M variants, including bla _CTXM-27 (C1-C3, C5), bla _CTXM-55 (C4) and bla _CTXM-15 (C6). Most ST1193 isolates were multidrug-resistant and carried an extensive array of AMR genes. ST1193 isolates also exhibited the ability to acquire further resistance while circulating in Vietnam. There were phylogenetic links between ST1193 isolates from BSIs and healthy individuals, suggesting these organisms may both establish long-term colonization in the human intestinal tract and induce infections. Our study uncovers factors shaping the population structure and transmission dynamics of multidrug-resistant ST1193 in Vietnam, and highlights the urgent need for local One Health genomic surveillance to capture new emerging ExPEC clones and to better understand the origins and transmission patterns of these pathogens.

The evolutionary history of Shigella flexneri serotype 6 in Asia

Shigella flexneri serotype 6 is an understudied cause of diarrhoeal diseases in developing countries, and has been proposed as one of the major targets for vaccine development against shigellosis. Despite being named as S. flexneri, Shigella flexneri serotype 6 is phylogenetically distinct from other S. flexneri serotypes and more closely related to S. boydii. This unique phylogenetic relationship and its low sampling frequency have hampered genomic research on this pathogen. Herein, by utilizing whole genome sequencing (WGS) and analyses of Shigella flexneri serotype 6 collected from epidemiological studies (1987–2013) in four Asian countries, we revealed its population structure and evolutionary history in the region. Phylogenetic analyses supported the delineation of Asian Shigella flexneri serotype 6 into two phylogenetic groups (PG-1 and −2). Notably, temporal phylogenetic approaches showed that extant Asian S. flexneri serotype 6 could be traced back to an inferred common ancestor arising in the 18^th century. The dominant lineage PG-1 likely emerged in the 1970s, which coincided with the times to most recent common ancestors (tMRCAs) inferred from other major Southeast Asian S. flexneri serotypes. Similar to other S. flexneri serotypes in the same period in Asia, genomic analyses showed that resistance to first-generation antimicrobials was widespread, while resistance to more recent first-line antimicrobials was rare. These data also showed a number of gene inactivation and gene loss events, particularly on genes related to metabolism and synthesis of cellular appendages, emphasizing the continuing role of reductive evolution in the adaptation of the pathogen to an intracellular lifestyle. Together, our findings reveal insights into the genomic evolution of the understudied Shigella flexneri serotype 6, providing a new piece in the puzzle of Shigella epidemiology and evolution.

Genomic epidemiology of persistently circulating MDR Shigella sonnei strains associated with men who have sex with men (MSM) in Belgium (2013-19)

Objectives

Shigella sonnei resistant to first-line antibiotics azithromycin and ciprofloxacin are on the rise globally. The aim of this study was to describe the epidemiology of MDR S. sonnei in Belgium and to identify origins and circulating clusters through WGS.

Methods

We undertook demographic, temporal and geographical analysis of 930 S. sonnei isolates submitted to the Belgian National Reference Centre for Salmonella and Shigella between 2017 and 2019. Phylogenetic analysis of WGS data, genotyping and identification of genetic markers of antimicrobial resistance was performed on 372 Belgian isolates submitted between 2013 and 2019.

Results

S. sonnei was identified in 75% (930/1253) of Belgian Shigella isolates submitted between 2017 and 2019. Overall, 7% (69/930) of isolates were resistant to ciprofloxacin alone, 6% (57/930) showed reduced susceptibility to azithromycin alone, and 24% (223/930) exhibited both. Men were at higher risk of carrying a double resistant S. sonnei strain, compared with women (risk ratio = 8.6, 95% CI = 5.4–13.9). Phylogenetic analysis revealed four independent Belgian clusters of persistently circulating MDR strains, associated with men who have sex with men (MSM) and of the same genotypes as previously described international MSM-related clades. Belgian isolates carried various incompatibility (Inc)-type plasmids, the SpA plasmid and ESBL genes.

Conclusions

In Belgium, S. sonnei isolates from men are much more likely to be resistant to important first-line antibiotics than isolates from women. Multiple co-circulating MDR S. sonnei clusters of different genotypes were identified in the MSM community. Further studies on risk groups are needed for targeted prevention, improved clinical and public health management and antimicrobial stewardship in Belgium.

Acquisition and loss of CTX-M plasmids in Shigella species associated with MSM transmission in the UK

Shigellosis in men who have sex with men (MSM) is caused by multidrug resistant Shigellae, exhibiting resistance to antimicrobials including azithromycin, ciprofloxacin and more recently the third-generation cephalosporins. We sequenced four bla_CTX-M-27-positive MSM Shigella isolates (2018–20) using Oxford Nanopore Technologies; three S. sonnei (identified as two MSM clade 2, one MSM clade 5) and one S. flexneri 3a, to explore AMR context. All S. sonnei isolates harboured Tn7/Int2 chromosomal integrons, whereas S. flexneri 3a contained the Shigella Resistance Locus. All strains harboured IncFII pKSR100-like plasmids (67-83kbp); where present bla_CTX-M-27 was located on these plasmids flanked by IS26 and IS903B, however bla_CTX-M-27 was lost in S. flexneri 3a during storage between Illumina and Nanopore sequencing. IncFII AMR regions were mosaic and likely reorganised by IS26; three of the four plasmids contained azithromycin-resistance genes erm(B) and mph(A) and one harboured the pKSR100 integron. Additionally, all S. sonnei isolates possessed a large IncB/O/K/Z plasmid, two of which carried aph(3’)-Ib/aph(6)-Id/sul2 and tet(A). Monitoring the transmission of mobile genetic elements with co-located AMR determinants is necessary to inform empirical treatment guidance and clinical management of MSM-associated shigellosis.

Association between floods and the risk of dysentery in China: a meta-analysis

The association between floods and the risk of dysentery remain controversial. Therefore, we performed a meta-analysis to clarify this relationship. A literature search was performed in PubMed, Web of science, and Embase for relevant articles published up to November 2019. Random-effects model was used to pool relative risks with 95% confidence intervals. The sensitivity analysis was carried out to evaluate the stability of the results. Publication bias was estimated using Egger's test. Eleven studies from 10 articles evaluated the association between floods and the risk of dysentery in China. The pooled RR (95% CI) of dysentery for the flooded time versus non-flooded period was 1.48 (95% CI: 1.14-1.91). Significant association was found in subgroup analysis stratified by dysentery styles [dysentery: 1.61 (95% CI: 1.34-1.93) and bacillary dysentery: 1.46 (95% CI: 1.06-2.01)]. The pooled RR (95%CI) of sensitivity analysis for dysentery was 1.26 (95% CI: 1.05-1.52). No significant publication bias was found in our meta-analysis. This meta-analysis confirms that floods have significantly increased the risk of dysentery in China. Our findings will provide more evidence to reduce negative health outcomes of floods in China.

Achromobacter spp. genetic adaptation in cystic fibrosis

Achromobacter spp. are emerging pathogens in patients with cystic fibrosis (CF) and Achromobacter spp. caused infections are associated with more severe disease outcomes and high intrinsic antibiotic resistance. While conventional CF pathogens are studied extensively, little is known about the genetic determinants leading to antibiotic resistance and the genetic adaptation in Achromobacter spp. infections. Here, we analysed 101 Achromobacter spp. genomes from 51 patients with CF isolated during the course of up to 20 years of infection to identify within-host adaptation, mutational signatures and genetic variation associated with increased antibiotic resistance. We found that the same regulatory and inorganic ion transport genes were frequently mutated in persisting clone types within and between Achromobacter species, indicating convergent genetic adaptation. Genome-wide association study of six antibiotic resistance phenotypes revealed the enrichment of associated genes involved in inorganic ion transport, transcription gene enrichment in β-lactams, and energy production and translation gene enrichment in the trimethoprim/sulfonamide group. Overall, we provide insights into the pathogenomics of Achromobacter spp. infections in patients with CF airways. Since emerging pathogens are increasingly recognized as an important healthcare issue, our findings on evolution of antibiotic resistance and genetic adaptation can facilitate better understanding of disease progression and how mutational changes have implications for patients with CF.

Circulation of multi-drug-resistant Shigella sonnei and Shigella flexneri among men who have sex with men in Barcelona, Spain, 2015-2019

BACKGROUND

In high-income countries, shigellosis is mainly found in travellers to high-risk regions or in men who have sex with men (MSM). This study investigated the genomic characteristics and the features of antimicrobial resistance of MSM-associated Shigella flexneri and Shigella sonnei circulating in Barcelona, Spain, elucidating their connectivity with contemporaneous Shigella spp. from other countries.

METHODS

Antimicrobial susceptibility, whole-genome sequencing, genomic characterization and phylogenetic analysis were performed in MSM-associated Shigella spp. recovered from 2015 to 2019. Reference genomes of MSM-associated Shigella spp. were included for contextualization and to determine their connection with international outbreaks.

RESULTS

In total, 44 S. flexneri and 26 S. sonnei were identified among MSM. Overall, 80% showed resistance to azithromycin, 65.7% showed resistance to trimethoprim-sulphamethoxazole and 32.8% showed resistance to ciprofloxacin; 27.1% were resistant to all three antimicrobials. mphA and/or ermB, and qnrS and mutations in the quinolone resistance determining regions were found in the azithromycin- and ciprofloxacin-resistant isolates, respectively. Additionally, two isolates carried bla_CTX-M-27. Single-nucleotide-polymorphism-based analysis revealed that the isolates were organized into different lineages, most of which were closely related to dominant MSM-associated lineages described previously in the UK and Australia.

CONCLUSIONS

This study investigated the circulation of lineages of S. flexneri and S. sonnei among MSM in Spain that were mainly resistant to first-/second-line oral treatments, and closely related to dominant MSM-associated lineages described previously in the UK and Australia. These data reinforce the urgent need for the implementation of public health measures focusing on the early detection and prevention of transmission of this emerging pathogen, which is contributing to the antimicrobial resistance crisis in sexually transmitted infections.

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.

Ecological niche adaptation of Salmonella Typhimurium U288 is associated with altered pathogenicity and reduced zoonotic potential

The emergence of new bacterial pathogens is a continuing challenge for agriculture and food safety. Salmonella Typhimurium is a major cause of foodborne illness worldwide, with pigs a major zoonotic reservoir. Two phylogenetically distinct variants, U288 and ST34, emerged in UK pigs around the same time but present different risk to food safety. Here we show using genomic epidemiology that ST34 accounts for over half of all S. Typhimurium infections in people while U288 less than 2%. That the U288 clade evolved in the recent past by acquiring AMR genes, indels in the virulence plasmid pU288-1, and accumulation of loss-of-function polymorphisms in coding sequences. U288 replicates more slowly and is more sensitive to desiccation than ST34 isolates and exhibited distinct pathogenicity in the murine model of colitis and in pigs. U288 infection was more disseminated in the lymph nodes while ST34 were recovered in greater numbers in the intestinal contents. These data are consistent with the evolution of S. Typhimurium U288 adaptation to pigs that may determine their reduced zoonotic potential.

Phylogenomic Investigation of Increasing Fluoroquinolone Resistance among Belgian Cases of Shigellosis between 2013 and 2018 Indicates Both Travel-Related Imports and Domestic Circulation

Shigellosis is an acute enteric infection caused mainly by the species Shigella flexneri and Shigella sonnei. Since surveillance of these pathogens indicated an increase in ciprofloxacin-resistant samples collected in Belgium between 2013 and 2018, a subset of 148 samples was analyzed with whole genome sequencing (WGS) to investigate their dispersion and underlying genomic features associated with ciprofloxacin resistance. A comparison between observed phenotypes and WGS-based resistance prediction to ciprofloxacin revealed perfect correspondence for all samples. Core genome multi-locus sequence typing and single nucleotide polymorphism-typing were used for phylogenomic investigation to characterize the spread of these infections within Belgium, supplemented with data from international reference collections to place the Belgian isolates within their global context. For S. flexneri, substantial diversity was observed with ciprofloxacin-resistant isolates assigned to several phylogenetic groups. Besides travel-related imports, several clusters of highly similar Belgian isolates could not be linked directly to international travel suggesting the presence of domestically circulating strains. For S. sonnei, Belgian isolates were all limited to lineage III, and could often be traced back to travel to countries in Asia and Africa, sometimes followed by domestic circulation. For both species, several clusters of isolates obtained exclusively from male patients were observed. Additionally, we illustrated the limitations of conventional serotyping of S. flexneri, which was impacted by serotype switching. This study contributes to a better understanding of the spread of shigellosis within Belgium and internationally, and highlights the added value of WGS for the surveillance of this pathogen.

Evolutionary histories and antimicrobial resistance in Shigella flexneri and Shigella sonnei in Southeast Asia

Conventional disease surveillance for shigellosis in developing country settings relies on serotyping and low-resolution molecular typing, which fails to contextualise the evolutionary history of the genus. Here, we interrogated a collection of 1,804 Shigella whole genome sequences from organisms isolated in four continental Southeast Asian countries (Thailand, Vietnam, Laos, and Cambodia) over three decades to characterise the evolution of both S. flexneri and S. sonnei. We show that S. sonnei and each major S. flexneri serotype are comprised of genetically diverse populations, the majority of which were likely introduced into Southeast Asia in the 1970s-1990s. Intranational and regional dissemination allowed widespread propagation of both species across the region. Our data indicate that the epidemiology of S. sonnei and the major S. flexneri serotypes were characterised by frequent clonal replacement events, coinciding with changing susceptibility patterns against contemporaneous antimicrobials. We conclude that adaptation to antimicrobial pressure was pivotal to the recent evolutionary trajectory of Shigella in Southeast Asia.

Outbreak of Central American born Shigella sonnei in two youth camps in Belgium in the summer of 2019

In 2019, an outbreak of Shigella sonnei occurred during two youth camps in Belgium. The clustering of isolates from both camps was confirmed by next-generation sequencing, as well as a secondary infection of a technician. The outbreak strain clustered with internationally isolated strains from patients with recent travel history to Central America. This report exemplifies enhanced surveillance and international collaboration between public health institutes by enabling to link local outbreaks to region-specific sublineages circulating abroad.

One is not enough: On the effects of reference genome for the mapping and subsequent analyses of short-reads

Mapping of high-throughput sequencing (HTS) reads to a single arbitrary reference genome is a frequently used approach in microbial genomics. However, the choice of a reference may represent a source of errors that may affect subsequent analyses such as the detection of single nucleotide polymorphisms (SNPs) and phylogenetic inference. In this work, we evaluated the effect of reference choice on short-read sequence data from five clinically and epidemiologically relevant bacteria (Klebsiella pneumoniae, Legionella pneumophila, Neisseria gonorrhoeae, Pseudomonas aeruginosa and Serratia marcescens). Publicly available whole-genome assemblies encompassing the genomic diversity of these species were selected as reference sequences, and read alignment statistics, SNP calling, recombination rates, dN/dS ratios, and phylogenetic trees were evaluated depending on the mapping reference. The choice of different reference genomes proved to have an impact on almost all the parameters considered in the five species. In addition, these biases had potential epidemiological implications such as including/excluding isolates of particular clades and the estimation of genetic distances. These findings suggest that the single reference approach might introduce systematic errors during mapping that affect subsequent analyses, particularly for data sets with isolates from genetically diverse backgrounds. In any case, exploring the effects of different references on the final conclusions is highly recommended.

Mapping consists in the alignment of reads (i.e., DNA fragments) obtained through high-throughput genome sequencing to a previously assembled reference sequence. It is a common practice in genomic studies to use a single reference for mapping, usually the ‘reference genome’ of a species—a high-quality assembly. However, the selection of an optimal reference is hindered by intrinsic intra-species genetic variability, particularly in bacteria. It is known that genetic differences between the reference genome and the read sequences may produce incorrect alignments during mapping. Eventually, these errors could lead to misidentification of variants and biased reconstruction of phylogenetic trees (which reflect ancestry between different bacterial lineages). To our knowledge, this is the first work to systematically examine the effect of different references for mapping on the inference of tree topology as well as the impact on recombination and natural selection inferences. Furthermore, the novelty of this work relies on a procedure that guarantees that we are evaluating only the effect of the reference. This effect has proved to be pervasive in the five bacterial species that we have studied and, in some cases, alterations in phylogenetic trees could lead to incorrect epidemiological inferences. Hence, the use of different reference genomes may be prescriptive to assess the potential biases of mapping.

Models in parasite and pathogen evolution: Genomic analysis reveals predominant clonality and progressive evolution at all evolutionary scales in parasitic protozoa, yeasts and bacteria

The predominant clonal evolution (PCE) model of pathogenic microorganisms postulates that the impact of genetic recombination in those pathogens' natural populations is not enough to erase a persistent phylogenetic signal at all evolutionary scales from microevolution till geological times in the whole ecogeographical range of the species considered. We have tested this model with a set of representative parasitic protozoa, yeasts and bacteria in the light of the most recent genomic data. All surveyed species, including those that were considered as highly recombining, exhibit similar PCE patterns above and under the species level, from macro- to micro-evolutionary scales (Russian doll pattern), suggesting gradual evolution. To our knowledge, it is the first time that such a strong common evolutionary feature among very diverse pathogens has been evidenced. The implications of this model for basic biology and applied research are exposed. These implications include our knowledge on the pathogens' reproductive mode, their population structure, the possibility to type strain and to follow up epidemics (molecular epidemiology) and to revisit pathogens' taxonomy through a flexible use of the phylogenetic species concept (Cracraft, 1983).

Shigella sonnei: virulence and antibiotic resistance

Shigella sonnei is the emerging pathogen globally, as it is the second common infectious species of shigellosis (bloody diarrhoea) in low- and middle-income countries (LMICs) and the leading one in developed world. The multifactorial processes and novel mechanisms have been identified in S. sonnei, that are collectively playing apart a substantial role in increasing its prevalence, while replacing the S. flexneri and other Gram-negative gut pathogens niche occupancy. Recently, studies suggest that due to improvement in sanitation S. sonnei has reduced cross-immunization from Plesiomonas shigelliodes (having same O-antigen as S. sonnei) and also found to outcompete the two major species of Enterobacteriaceae family (Shigella flexneri and Escherichia coli), due to encoding of type VI secretion system (T6SS). This review aimed to highlight S. sonnei as an emerging pathogen in the light of recent research with pondering aspects on its epidemiology, transmission, and pathogenic mechanisms. Additionally, this paper aimed to review S. sonnei disease pattern and related complications, symptoms, and laboratory diagnostic techniques. Furthermore, the available treatment reigns and antibiotic-resistance patterns of S. sonnei are also discussed, as the ciprofloxacin and fluoroquinolone-resistant S. sonnei has already intensified the global spread and burden of antimicrobial resistance. In last, prevention and controlling strategies are briefed to limit and tackle S. sonnei and possible future areas are also explored that needed more research to unravel the hidden mysteries surrounding S. sonnei.

Phylogenetic Methods for Genome-Wide Association Studies in Bacteria

Genome-wide association studies in bacteria have great potential to deliver a better understanding of the genetic basis of many biologically important phenotypes, including antibiotic resistance, pathogenicity, and host adaptation. Such studies need however to account for the specificities of bacterial genomics, especially in terms of population structure, homologous recombination, and genomic plasticity. A powerful way to tackle this challenge is to use a phylogenetic approach, which is based on long-standing methodology for the evolutionary analysis of bacterial genomic data. Here we present both the theoretical and practical aspects involved in the use of phylogenetic methods for bacterial genome-wide association studies.

A Multifactorial Approach for Surveillance of Shigella spp. and Entero-Invasive Escherichia coli Is Important for Detecting (Inter)national Clusters

Shigella spp. and entero-invasive Escherichia coli (EIEC) can cause mild diarrhea to dysentery. In Netherlands, although shigellosis is a notifiable disease, there is no laboratory surveillance for Shigella spp. and EIEC in place. Consequently, the population structure for circulating Shigella spp. and EIEC isolates is not known. This study describes the phenotypic and serological characteristics, the phenotypic and genetic antimicrobial resistance (AMR) profiles, the virulence gene profiles, the classic multi-locus sequence types (MLST) and core genome (cg)MLST types, and the epidemiology of 414 Shigella spp. and EIEC isolates collected during a cross-sectional study in Netherlands in 2016 and 2017. S. sonnei (56%), S. flexneri (25%), and EIEC (15%) were detected predominantly in Netherlands, of which the EIEC isolates were most diverse according to their phenotypical profile, O-types, MLST types, and cgMLST clades. Virulence gene profiling showed that none of the isolates harbored Shiga toxin genes. Most S. flexneri and EIEC isolates possessed nearly all virulence genes examined, while these genes were only detected in approximately half of the S. sonnei isolates, probably due to loss of the large invasion plasmid upon subculturing. Phenotypical resistance correlated well with the resistant genotype, except for the genes involved in resistance to aminoglycosides. A substantial part of the characterized isolates was resistant to antimicrobials advised for treatment, i.e., 73% was phenotypically resistant to co-trimoxazole and 19% to ciprofloxacin. AMR was particularly observed in isolates from male patients who had sex with men (MSM) or from patients that had traveled to Asia. Furthermore, isolates related to international clusters were also circulating in Netherlands. Travel-related isolates formed clusters with isolates from patients without travel history, indicating their emergence into the Dutch population. In conclusion, laboratory surveillance using whole genome sequencing as high-resolution typing technique and for genetic characterization of isolates complements the current epidemiological surveillance, as the latter is not sufficient to detect all (inter)national clusters, emphasizing the importance of multifactorial public health approaches.

Phylogenetic and Evolutionary Analysis Reveals the Recent Dominance of Ciprofloxacin-Resistant Shigella sonnei and Local Persistence of S. flexneri Clones in India

Shigella is the second leading cause of bacterial diarrhea worldwide. This has been categorized as a priority pathogen among enteric bacteria by the Global Antimicrobial Resistance Surveillance System (GLASS) of the World Health Organization (WHO). Recently, S. sonnei seems to be replacing S. flexneri in low- and middle-income countries undergoing economic development. Antimicrobial resistance in S. flexneri and S. sonnei is a growing international concern, specifically with the international dominance of the multidrug-resistant (MDR) lineage. Genomic studies focusing on S. flexneri and S. sonnei in India remain largely unexplored. This study provides information on the introduction and expansion of drug-resistant Shigella strains in India for the first time by comparing the genome sequences of S. flexneri and S. sonnei isolates from India with the publicly available genomes of global strains. The study discusses the key differences between the two dominant species of Shigella at the genomic level to understand the evolutionary trends and genome dynamics of emerging and existing resistance clones. The present work demonstrates evidence for the long-term persistence of all PGs of S. flexneri and the recent dominance of a ciprofloxacin-resistant S. sonnei lineage in India.

Shigella is the second leading cause of bacterial diarrhea worldwide. Recently, Shigella sonnei seems to be replacing Shigella flexneri in low- and middle-income countries undergoing economic development. Despite this, studies focusing on these species at the genomic level remain largely unexplored. Here, we compared the genome sequences of S. flexneri and S. sonnei isolates from India with the publicly available genomes of global strains. Our analysis provides evidence for the long-term persistence of all phylogenetic groups (PGs) of S. flexneri and the recent dominance of the ciprofloxacin-resistant S. sonnei lineage in India. Within S. flexneri PGs, the majority of the study isolates belonged to PG3 within the predominance of serotype 2. For S. sonnei, the current pandemic involves globally distributed multidrug-resistant (MDR) clones that belong to Central Asia lineage III. The presence of such epidemiologically dominant lineages in association with stable antimicrobial resistance (AMR) determinants results in successful survival in the community.

IMPORTANCEShigella is the second leading cause of bacterial diarrhea worldwide. This has been categorized as a priority pathogen among enteric bacteria by the Global Antimicrobial Resistance Surveillance System (GLASS) of the World Health Organization (WHO). Recently, S. sonnei seems to be replacing S. flexneri in low- and middle-income countries undergoing economic development. Antimicrobial resistance in S. flexneri and S. sonnei is a growing international concern, specifically with the international dominance of the multidrug-resistant (MDR) lineage. Genomic studies focusing on S. flexneri and S. sonnei in India remain largely unexplored. This study provides information on the introduction and expansion of drug-resistant Shigella strains in India for the first time by comparing the genome sequences of S. flexneri and S. sonnei isolates from India with the publicly available genomes of global strains. The study discusses the key differences between the two dominant species of Shigella at the genomic level to understand the evolutionary trends and genome dynamics of emerging and existing resistance clones. The present work demonstrates evidence for the long-term persistence of all PGs of S. flexneri and the recent dominance of a ciprofloxacin-resistant S. sonnei lineage in India.

Emergence of NDM-producing Klebsiella pneumoniae and Escherichia coli in Spain: phylogeny, resistome, virulence and plasmids encoding blaNDM-like genes as determined by WGS

OBJECTIVES

NDM carbapenemases have spread worldwide. However, little information exists about the impact of NDM-producing Enterobacteriaceae in Spain. By WGS, we sought to elucidate the population structure of NDM-like-producing Klebsiella pneumoniae and Escherichia coli in Spain and to determine the plasmids harbouring blaNDM-like genes.

METHODS

High-resolution SNP typing, core-genome MLST and plasmid reconstruction (PlasmidID) were performed on 59 NDM-like-producing K. pneumoniae and 8 NDM-like-producing E. coli isolated over an 8 year period in Spain.

RESULTS

Five major epidemic clones of NDM-producing K. pneumoniae caused five important nationwide outbreaks: ST437/NDM-7, ST437/NDM-1, ST147/NDM-1, ST11/NDM-1 and ST101/NDM-1; in contrast, the spread of NDM-producing E. coli was polyclonal. Three blaNDM types were identified: blaNDM-1, 61.2%; blaNDM-7, 32.8%; and blaNDM-5, 6%. Five K. pneumoniae isolates co-produced other carbapenemases (three blaOXA-48 and two blaVIM-1). The average number of acquired resistance genes was higher in K. pneumoniae than in E. coli. The plasmids encoding blaNDM-like genes belonged to IncFII, IncFIB, IncX3, IncR, IncN and IncC types, of which IncF, IncR and IncC were associated with MDR. The genetic surroundings of blaNDM-like genes showed a highly variable region upstream of ISAba125.

CONCLUSIONS

In recent years NDM-producing K. pneumoniae and E. coli have emerged in Spain; the spread of a few high-risk K. pneumoniae clones such as ST437/NDM-7, ST437/NDM-1, ST147/NDM-1, ST11/NDM-1 and ST101/NDM-1 have caused several interregional outbreaks. In contrast, the spread of NDM-producing E. coli has been polyclonal. Plasmid types IncFII, IncFIB, IncX3, IncR, IncN and IncC carried blaNDM, and the same IncX3 plasmid was detected in K. pneumoniae and E. coli.