Species-wide whole genome sequencing reveals historical global spread and recent local persistence in Shigella flexneri

Protein-specific immune response elicited by the Shigella sonnei 1790GAHB GMMA-based candidate vaccine in adults with varying exposure to Shigella

Shigella is a leading cause of diarrheal morbidity and mortality in young children from low- and middle-income countries. Here, we aimed to verify the ability of the generalized modules for membrane antigens (GMMA)-based Shigella sonnei candidate vaccine 1790GAHB to elicit an anti-protein antibody response. Serum samples from previous clinical trials in adults (a dose-escalation study and its extension in France, a vaccine efficacy study after human challenge in the United States, and a study in Kenya) were investigated using pan-proteome microarrays consisting of 3,150 full-length or fragmented Shigella proteins. Pre-/post-vaccination comparisons identified subsets of proteins that were highly immunoreactive and largely overlapped across all trials; the T3SS lipochaperone family protein (expressed on GMMA) was the most reactive in all studies. Responses to several microarray antigens correlated well with S. sonnei LPS serum IgG antibody levels. Overall, we confirmed the ability of GMMA to elicit an anti-protein IgG/IgA response; however, no association with protection against shigellosis was identified. In the challenge study, IgG response to seven antigens (IpaC, IpaB, IpaA, IpaD, IpaH, IpgC, and MxiD; not expressed on GMMA) was associated with a decreased risk of shigellosis. These antigens were observed to also have high IgG responses at baseline in individuals naturally exposed to Shigella and could constitute targets for future vaccine development.IMPORTANCEShigella remains a major cause of diarrheal disease, especially in children aged under 5 years from low-to-middle-income countries. No vaccine against shigellosis is yet widely available despite the high public health need. An ideal vaccine would provide protection against the most prevalent species, Shigella flexneri and Shigella sonnei; therefore, it could be relevant to identify common antigens. We developed a microarray containing 3,150 full-length or fragmented proteins selected across Shigella species. Sera collected in four clinical trials conducted in three countries of varying endemicity to evaluate a S. sonnei GMMA-based candidate vaccine were tested against these proteins. We identified several Shigella proteins (IpaC, IpaB, IpaA, IpaD, IpaH, IpgC, MxiD) that induced robust antibody response following experimental challenge or natural infection. These proteins correlated with a reduced risk of shigellosis after the S. sonnei challenge. We found no apparent role for anti-GMMA proteins' IgG or IgA response in protection against shigellosis.

Complete genome sequence of Shigella phage Moonfish isolated from Mid-Michigan

A Shigella podophage, Moonfish, was isolated from a longhorn cattle ranch in Barry County, MI, USA. The 69,166-bp genome was annotated and places Moonfish in a similar category to N4-like phages. Host range studies show that Moonfish has a narrow host profile relative to similar phages.

Shigella sonnei: epidemiology, evolution, pathogenesis, resistance and host interactions

Shigella sonnei is a major cause of diarrhoea globally and is increasing in prevalence relative to other Shigella because of multiple demographic and environmental influences. This single-serotype species has traditionally received less attention in comparison to Shigella flexneri and Shigella dysenteriae, which were more common in low-income countries and more tractable in the laboratory. In recent years, we have learned that Shigella are highly complex and highly susceptible to environmental change, as exemplified by epidemiological trends and increasing relevance of S. sonnei. Ultimately, methods, tools and data generated from decades of detailed research into S. flexneri have been used to gain new insights into the epidemiology, microbiology and pathogenesis of S. sonnei. In parallel, widespread adoption of genomic surveillance has yielded insights into antimicrobial resistance, evolution and organism transmission. In this Review, we provide an overview of current knowledge of S. sonnei, highlighting recent insights into this globally disseminated antimicrobial-resistant pathogen and assessing how novel data may impact future vaccine development and implementation.

Bacterial lysis or survival after infection with phage Sf14 depends on combined nutrient and temperature conditions

Bacteriophage Sf14 is a Moogle-like myovirus that infects Shigella flexneri. S. flexneri is a human pathogen that replicates intracellularly in the intestine, but it persists in a low metabolic state in environmental fresh water sources. Though closely related to FelixO1, Moogleviruses were more recently discovered within the last 10 years; thus, mechanistic knowledge of their infection cycles is still being gathered. This work investigated the combined effects of temperature and nutrient concentration on both host growth and phage replication. In combination, a total of 16 different conditions were analyzed. Results indicate that nutrient-rich media facilitate shorter infection cycles and support phage production at all temperatures. As nutrient content decreased, temperature significantly affected both host cell replication and phage production. Results indicate phage genomes are entering the cells and genes are actively expressed; however, there is a significant delay in expression, which could allow bacterial populations to outpace phage growth.

Molecular epidemiology of a multidrug-resistant Shigella sonnei outbreak in Tunisia (2022-2023) using whole-genome sequencing

Purpose. The prevalence of multidrug-resistant (MDR) Shigella sonnei is increasing globally, raising concerns for public health. In 2022, an outbreak of MDR S. sonnei was observed in Tunisia. We aimed to evaluate the genetic profile of S. sonnei isolates during the outbreak, including their clonal relationship, antimicrobial determinants and connection to international strains.Methods. In this study, we sequenced the whole genome of 24 S. sonnei strains collected from South Tunisia between July 2022 and November 2023. Bioinformatic analysis was conducted to confirm species identification, assign sequence types, determine core genome sequence types, analyse phylogenetic relationships and identify antimicrobial resistance determinants. Phylodynamic and phylogeographic analyses were performed to trace the spatiotemporal spread of the outbreak genotype.Results. Our investigation revealed that 23 out of 24 isolates were grouped into the HC10-20662 genotype within the 3.6.3 subclade. All isolates carried the blaCTX-M-15 gene associated with extended-spectrum beta-lactamase production, as well as the dfrA1 and qnrS1 genes, along with the D87G mutation in gyrA. Additionally, the sul2, tet(A) and mph(A) resistance genes were present in most isolates (96%, 96 and 83, respectively). Phylogeographic analysis suggested that the outbreak genotype likely spread in Europe before being introduced into Tunisia.Conclusion. To the best of our knowledge, this is the first MDR S. sonnei outbreak in the country. The HC10-20662 genotype appears to be responsible for a multi-country outbreak, affecting both Tunisia and Europe. Continued genomic surveillance efforts, both nationally and internationally, are essential for monitoring the dynamic evolution and global spread of MDR S. sonnei.

Understanding and addressing β-lactam resistance mechanisms in gram-negative bacteria in Lebanon: A scoping review

Background

A growing threat to public health is the worldwide problem of antimicrobial resistance (AMR), in which gram-negative organisms are playing a significant role. Antibiotic abuse and misuse, together with inadequate monitoring and control protocols, have contributed to the emergence of resistant strains. This global scenario prepares us to look more closely at the situation in Lebanon. The aim of this review is to investigate in detail the resistance mechanisms and related genes that are displayed by gram-negative organisms in Lebanon.

Methods

A comprehensive analysis was carried out to pinpoint and gather information regarding gram-negative bacteria displaying resistance to antibiotics. To contribute to a complete understanding of the current state of antibiotic resistance in gram-negative strains, it was intended to collect and evaluate data on these organisms' resistance patterns in a comprehensive manner.

Results

Several studies have emphasized the prevalence of carbapenem-resistant Enterobacteriaceae (CRE) in Lebanon, specifically noting Escherichia coli and Klebsiella pneumoniae as the most frequent culprits, with OXA-48 and NDM-1 being the primary carbapenemases discovered. Furthermore, the TEM β-lactamase families are the primary source of extended-spectrum β-lactamases (ESBLs) in Shigella and Salmonella. Additionally, resistant strains of Acinetobacter baumannii and Pseudomonas aeruginosa have been linked to nosocomial infections in the country.

Conclusion

There is a considerable frequency of antibiotic overuse and misuse in Lebanon, based to the limited data available on antibiotic consumption. In conclusion, antibiotic stewardship initiatives and additional research beyond the confines of single-center studies in Lebanon are needed.

Molecular epidemiology of extended-spectrum beta-lactamases and carbapenemases-producing Shigella in Africa: a systematic review and meta-analysis

BACKGROUND

The treatment of Shigella infections has become a major challenge due to the emergence of multidrug-resistant Shigella. There is however insufficient knowledge regarding the molecular epidemiology of Shigella strains producing beta-lactamases in Africa. This systematic review investigated the scientific literature on the molecular epidemiology of extended-spectrum beta-lactamase (ESBL) and carbapenemases producing Shigella in Africa.

METHODS

Papers published in English and French from African countries on the molecular epidemiology of ESBL and carbapenemase producing Shigella from January 1999 to July 5, 2024 were reviewed. An extensive literature search was conducted through electronic databases including PubMed, Scopus, Web of Sciences, African Journals Online (AJOL) and Google scholar using specific keywords. The meta-analysis and forest plots of Shigella species, ESBL and carbapenemases genes were done using the comprehensive Meta-Analysis software. All data were analyzed using a binary random-effects model by the DerSimonian-Laird method at a 95% confidence interval.

RESULTS

Out of the 583 research articles, only 18 (3.1%) articles representing eleven countries were included in the meta-analysis. The overall pooled prevalence of ESBL and Carbapenem Resistant (CR)-producing Shigella was estimated as 41.2% (95% CI: 22.8-62.4; I2 = 93.7%, p < 0.05). The leading ESBL and CR-producing Shigella species reported in this review was Shigella flexneri, 34.5% (95% CI: 16.6-58.2; I2 = 94.9%, p < 0.05). ESBL and CR-producing Shigella sonnei was the least reported with estimated prevalence of 6.7% (95% CI: 3.4-13.2; I2 = 80.7%, p < 0.05). In this review, blaTEM, blaOXA-1, and blaCTX-M were the most prevalent genes in Africa with prevalence of 25.9% (95% CI: 13.9-43.2; I2 = 90.9%, p < 0.05), 25.7% (95% CI: 14.9-43.0; I2 = 93.7%, p < 0.05), and 10.8% (95% CI: 4.5-23.4; I2 = 85.6%, p < 0.05) respectively. The prevalence of Carbapenemases genes on the other hand was low, reported as 0.8% (95% CI: 0.2-10.3; I2 = 51.4%; p < 0.05) for blaNDM, 1.1% for blaKPC, and 0.5% for blaIMP.

CONCLUSION

This study highlighted Shigella flexneri as the most prevalent ESBL and CR-producing Shigella species found in various African countries. The findings from this review indicate blaTEM, blaOXA-1 and blaCTX-M as the most prevalent genes in Africa while carbapenemases are least prevalent. The findings from the study suggest that ESBL and CR-producing Shigella pose a significant public health threat in Africa. Effective antimicrobial resistance (AMR) surveillance strategies are needed to improve the management and treatment of Shigella infections in Africa.

Characterization of Shigella flexneri serotype 6 strains from geographically diverse low- and middle-income countries

Shigella flexneri serotype 6 (Sf6) is one of the most common serotypes recovered from surveillance studies of moderate to severe diarrhea. Despite the clinical significance of Sf6, this serotype is understudied. In this work, we have performed both serotype-specific genomic and phenotypic comparisons of Sf6 isolates to one another and non-S. flexneri serotypes. Comparative genomic analyses identified significant nucleotide homology between Sf6 strains (n = 325), despite a broad range of collection timeframes and geographic locations. We identified Sf6 specific factors, including a potential novel Shigella virulence factor (type II secretion system). Additionally, we identified established Shigella virulence genes (ospG) and metabolic genes (rutABCDEFGR) that were absent in Sf6 strains while present in the majority of 728 non-Sf6 strains. Complete sequencing of 11 clinical Sf6 strains, demonstrated that the Sf6 virulence plasmid (pINV) is ~38 kb smaller than the average non-Sf6 pINV (~228 kb). Comparisons of S. flexneri species level antibiotic susceptibility highlighted that clinical Sf6 isolates from Africa in the Global Enteric Multicenter Study (GEMS) and Vaccine Impact on Diarrhea in Africa (VIDA) study demonstrated geographic, serotype-specific susceptibility pattern. Phenotypic analyses of Sf6 identified reduced intracellular invasion and cytokine induction from HT-29 cells, as well as reduced Ipa protein effector secretion, compared with S. flexneri serotype 2a strain 2457T. Together these data highlight conserved and unique serotype-specific genotypic and phenotypic features for Sf6. This level of conservation has not been noted for other S. flexneri serotypes and is promising for vaccine and diagnostic assays to provide global Sf6-specific coverage.IMPORTANCEShigellosis is an ongoing global public health crisis with >270 million annual episodes among all age groups; however, the greatest disease burden is among children in low- and middle-income countries (LMIC). The lack of a licensed Shigella vaccine and the observed rise in antimicrobial-resistant Shigella spp. highlights the urgency for effective preventative and interventional strategies. The inclusion of S. flexneri serotype 6 (Sf6) is a necessary component of a multivalent vaccine strategies based on its clinical and epidemiological importance. Given the genomic diversity of Sf6 compared with other S. flexneri serotypes and Sf6 unique O-antigen core structure, serotype-specific characterization of Sf6 is a critical step to inform Shigella-directed vaccine and alternative therapeutic designs. Herein, we identified conserved genomic content among a large collection of temporally and geographically diverse Sf6 clinical isolates and characterized genotypic and phenotypic properties that separate Sf6 from non-Sf6 S. flexneri serotypes.

Characterization and application of LysSGF2 and HolSGF2 as potential biocontrol agents against planktonic and biofilm cells of common pathogenic bacteria

Antimicrobial resistance represents a global health emergency, necessitating the introduction of novel antimicrobial agents. In the present study, lysozyme and holin from Shigella flexneri 1.1868 phage SGF2, named LysSGF2 and HolSGF2, respectively, were cloned, expressed, and characterized. LysSGF2 and HolSGF2 showed lytic activities against S. flexneri 1.1868 cells at 4-55 °C and pH 3.1-10.3. LysSGF2 exhibited antimicrobial activity against five gram-negative and two gram-positive bacteria. HolSGF2 showed antimicrobial activity against four gram-negative and one gram-positive species. The antibacterial activities of LysSGF2 and HolSGF2 were determined in liquid beverages, including bottled water and milk. The relative lytic activity of LysSGF2 combined with HolSGF2 against the tested bacteria was approximately 46-77 % in water. Furthermore, the combination markedly decreased the viable counts of tested bacteria by approximately 3-5 log CFU/mL. LysSGF2 and HolSGF2 could efficiently remove biofilms on polystyrene, glass, and stainless-steel. The efficacy of the LysSGF2 and HolSGF2 combination against the tested bacteria on polystyrene was 58-71 %. Combination treatment effectively killed biofilm cells formed on stainless-steel and glass by 1-4 log CFU/mL. ese results indicate that LysSGF2 and HolSGF2 can successfully control both the planktonic and biofilm cells of common pathogenic bacteria, suggesting that the combined or single use of LysSGF2 and HolSGF2 may be of great value in food processing.

Whole genome sequencing and antimicrobial resistance among clinical isolates of Shigella sonnei in Addis Ababa, Ethiopia

BACKGROUND

Shigellosis is an acute gastroenteritis infection and one of Ethiopia's most common causes of morbidity and mortality, especially in children under five. Antimicrobial resistance (AMR) has spread quickly among Shigella species due to inappropriate antibiotic use, inadequacies of diagnostic facilities, and unhygienic conditions. This study aimed to characterize Shigella sonnei (S. sonnei) using whole genome sequence (WGS) analysis in Addis Ababa, Ethiopia.

METHODS

The raw reads were quality-filtered and trimmed, and a minimum length of 50bp was retained and taxonomically classified using MiniKraken version 1. The whole genome data were aligned with Antibiotic Resistance Gene (ARG) sequences of the Comprehensive Antibiotic Resistance Database (CARD) by Resistance Gene Identifier (RGI). Plasmids were analyzed using the PlasmidFinder tool version 2.1. Additionally, AMR and virulence genes were screened at the Centre for Genomic Epidemiology (CGE) web-based server.

RESULTS

All isolates in our investigation contained genes encoding blaEC-8 and blaZEG-1. Here, 60.7% of the isolates were phenotypically sensitive to cefoxitin among the blaEC-8 genes detected in the genotyping analysis, whereas all isolates were completely resistant to amoxicillin and erythromycin phenotypically. The study also identified genes that conferred resistance to trimethoprim (dfrA). Plasmid Col156 and Col (BS512) types were found in all isolates, while IncFII and Col (MG828) plasmids were only identified in one isolate.

CONCLUSION

This study found that many resistant genes were present, confirming the high variety in S. sonnei strains and hence a divergence in phylogenetic relationships. Thus, combining WGS methods for AMR prediction and strain identification into active surveillance may be beneficial for monitoring the spread of AMR in S. sonnei and detecting the potential emergence of novel variations.

Genome and antibiotic resistance characteristics of Shigella clinical isolates in Fujian Province, Southeast China, 2005-2019

Shigellosis is a serious public health issue in many developing countries. The emergence of multidrug-resistant (MDR) Shigella isolates has deepened the treatment difficulty and health burden of shigellosis. China is the largest developing country in the world, but so far, the genome of MDR Shigella isolates has not been well characterized. In this study, 60 clinical isolates of Shigella spp. in Fujian Province, southeast China, from 2005 to 2019 were characterized for drug resistance phenotype, whole-genome sequencing and bioinformatics analysis. The results showed that the MDR rate of Shigella isolates was 100%, among which the resistance rates of cefotaxime, ciprofloxacin and azithromycin were 36.67, 21.67 and 10.00 %, respectively. The positive rate of extended-spectrum beta-lactamase (ESBL)-producing strains was 23.33%. The resistance profiles of Shigella flexneri and Shigella sonnei to some antimicrobials differed. The MDR isolates carried multiple antimicrobial resistance genes, among which blaCTX-M-14 and blaCTX-M-15 mediated ESBL resistance. 'ISEcp1 -blaCTX-M -IS903' (type I) and 'ISEcp1 -blaCTX-M' (type II) were the most common genetic environments around the blaCTX-M genes, and plasmids containing these structures included IncFII, IncI1, IncI2 and IncN. The double gene mutation pattern of gyrA and parC resulted in a significant decrease in the sensitivity of S. flexneri to ciprofloxacin. The overall resistance phenotype and genotype concordance rate was 88.50%, and the sensitivity and specificity of the genotype antimicrobial susceptibility test (AST) were 93.35 and 82.53 %, respectively. However, inconsistency occurred between phenotypic and genotype profiles for a few antibiotics. Phylogenomic investigation with core genome multi-locus sequence typing (cgMLST) and SNPs were used to characterize the endemic transmission of these infections in Fujian and their evolutionary origin within the global context. For S. flexneri, Fujian isolates were all limited to PG3 and could be divided into three phylogenetic clusters. The ciprofloxacin-resistant strains were mainly F2a and FXv and assigned to the three clusters with different quinolone resistance-determining region mutation patterns. For S. sonnei, most Fujian strains belonged to Lineage III with genotype 3.7.6, except three isolates of Lineage I with genotype 1.3. The strains carrying the blaCTX-M genes were dispersed, indicating different origins of gene acquisition. Most of the circulating isolates in Fujian Province were not related to major international outbreak lineages and were only endemic to the country. In conclusion, multi-drug resistance of Shigella isolates in Fujian Province was serious, and genome-based laboratory surveillance will be crucial to the clinical treatment and public health measures for shigellosis.

Foodborne Bacterial Pathogens in Animal, Food, and Environmental Samples Collected From the Physical Exposure of Children With Diarrhea in Ethiopia: A One Health Approach

Foodborne pathogens (FBPs) are transmitted principally through the consumption of contaminated food or drinking water and pose a remarkable public health risk, particularly in low-income countries. A cross-sectional study was conducted between November 2021 and January 2023 to estimate the prevalence, co-occurrence, and monthly patterns of FBPs in the physical exposures of children with diarrhea in Harar town and Kersa district, Ethiopia. Animal, food, and environmental samples were collected from direct or indirect contact sites of children with diarrhea. The isolation and identification of FBPs, including nontyphoidal Salmonella (NTS), diarrheagenic E. coli (DEC), and Shigella, was performed using selective and deferential culture media and a series of biochemical tests. Among the 438 analyzed samples, the overall prevalence of these pathogens was 18.3%, with 3.9% co-occurrence and 14.4% single pathogen occurrence rates. The highest prevalence was observed in wastewater (40.9%; AOR = 3.3; 95% CI: 1.1-10.1). The pathogen detection rate in food was 17.9% (AOR = 1.2; 95% CI: 0.4-3.6), with no significant difference between animal-sourced and other food categories. The occurrence rates of NTS, DEC, and Shigella in the meat samples were 13.9%, 5.4%, and 6.5%, respectively. Interestingly, DEC and Shigella were detected in cooked food. Moreover, Shigella was detected in drinking water (5%) and other water sources (10%). A significantly higher prevalence of FBPs was detected in poultry than in cattle and camel feces. This study revealed fluctuations in the monthly occurrence patterns of FBPs, with a peak of 37.1% during the dry season. In conclusion, the study revealed a high prevalence of FBPs, with no significant differences between rural and urban areas or food and water sources, highlighting the need for food safety measures in both settings. Further studies with larger sample sizes and advanced diagnostics are recommended to determine the relative contribution of each source.

Trends in shigellosis notifications in England, January 2016 to March 2023

We reviewed all diagnoses of Shigella species notified to the UK Health Security Agency from January 2016 to March 2023. An overall increase in notifications of shigellosis was seen between 2016 (n = 415/quarter) and 2023 (n = 1 029/quarter). However, notifications dramatically declined between March 2020 and September 2021 during the COVID-19 pandemic (n = 208/quarter) highlighting the impact of travel and social distancing restrictions on transmission. S. sonnei diagnoses were more affected by lockdown restrictions than S. flexneri, most likely due to a combination of species-specific characteristics and host attributes. Azithromycin resistance continued to be associated with epidemics of sexually transmissible S. flexneri (adult males = 45.6% vs. adult females = 8.7%) and S. sonnei (adult males = 59.5% vs. adult females = 14.6%). We detected resistance to ciprofloxacin in S. sonnei from adult male cases not reporting travel at a higher frequency (79.4%) than in travel-associated cases (61.7%). Extensively drug-resistant Shigella species associated with sexual transmission among men almost exclusively had ESBL encoded by blaCTX-M-27, whereas those associated with returning travellers had blaCTX-M-15. Given the increasing incidence of infections and AMR, we recommend that enhanced surveillance is used to better understand the impact of travel and sexual transmission on the acquisition and spread of MDR and XDR Shigella species.

Synthetic extremophiles via species-specific formulations improve microbial therapeutics

Microorganisms typically used to produce food and pharmaceuticals are now being explored as medicines and agricultural supplements. However, maintaining high viability from manufacturing until use remains an important challenge, requiring sophisticated cold chains and packaging. Here we report synthetic extremophiles of industrially relevant gram-negative bacteria (Escherichia coli Nissle 1917, Ensifer meliloti), gram-positive bacteria (Lactobacillus plantarum) and yeast (Saccharomyces boulardii). We develop a high-throughput pipeline to define species-specific materials that enable survival through drying, elevated temperatures, organic solvents and ionizing radiation. Using this pipeline, we enhance the stability of E. coli Nissle 1917 by more than four orders of magnitude over commercial formulations and demonstrate its capacity to remain viable while undergoing tableting and pharmaceutical processing. We further show, in live animals and plants, that synthetic extremophiles remain functional against enteric pathogens and as nitrogen-fixing plant supplements even after exposure to elevated temperatures. This synthetic, material-based stabilization enhances our capacity to apply microorganisms in extreme environments on Earth and potentially during exploratory space travel.

Lipidomic approach to identify Escherichia coli and Shigella spp. by matrix-assisted laser desorption/ionization mass spectrometry

PURPOSE

Escherichia coli (E. coli) and Shigella species, being highly similar, present a challenge for differentiation using classical methods such as phenotyping, 16S rRNA sequencing, or protein profiling using matrix-assisted laser desorption/ionization mass spectrometry (MALDI MS). The paper proposes a method for identifying E. coli, S. flexneri, S. sonnei, and S. boydii by augmenting the Bruker Biotyper database with reference spectra of lipid profiles obtained using MALDI MS in the positive mode.

MATERIALS/METHODS

Lipid extracts were made from cultured E. coli, S. flexneri, S. sonnei and S. boydii using the Bligh & Dyer protocol. MALDI MS spectra in positive ion mode were performed for the extracts. Reference spectra were created from 30 spectra for each bacterium and added to the Bruker Biotyper database.

RESULTS

Identification of bacteria based on lipid profiles in the Biotyper database gave correct results with scores above 2.49. Statistical analysis of the results by Partial Least Squares-Discriminant Analysis (PLS-DA) showed that it is possible to correctly differentiate the microorganisms studied using the lipidomic approach. A panel of six m/z values was proposed for which the value of the area under the ROC curve is 1, thus enabling the identification of E. coli and S. flexneri with 100 ​% accuracy.

CONCLUSIONS

Identification of bacteria from lipid fingerprints obtained by the MALDI MS technique is possible and may become a useful tool in the future, especially for microorganisms that are difficult to distinguish by other methods.

The Tomato Variety Affects the Survival of Shigella flexneri 2a in Fruit Pericarp

The presence of enteric pathogens in produce can serve as a significant means of transmitting infections to consumers. Notably, tomatoes, as a type of produce, have been implicated in outbreaks caused by various human pathogens, such as Salmonella enterica and pathogenic Escherichia coli. However, the survival characteristics of Shigella spp. in tomatoes have not been thoroughly investigated. In this study, we assess the survival of S. flexneri 2a in two distinct varieties of post-harvested tomatoes. S. flexneri 2a was used to inoculate both regular-sized Vine tomatoes and cherry-type Mini Plum tomatoes. Our findings reveal no significant difference in Shigella survival in the pericarp of both varieties on day 2 post-inoculation. However, a significant disparity emerges on day 6, where all recovered Shigella colonies exclusively belong to the Mini Plum variety, with none associated with the Vine type. When Shigella was inoculated into the locular cavity (deep inoculation), no significant difference between varieties was observed. Additionally, we investigate the potential role of the SRL pathogenicity island (SRL PAI) in the survival and fitness of S. flexneri 2a in post-harvested tomatoes. Our results indicate that while the SRL PAI is not linked to the survival of the strains in tomato, it does impact their fitness. These findings underscore the variability in Shigella strains' survival capabilities depending on the tomato variety, highlighting the importance of understanding Shigella ecology beyond the human host and identifying molecular determinants influencing bacterial survival to mitigate the risk of future outbreaks. The significance of this data on Shigella persistence in fresh vegetables should not be underestimated, as even a small number of Shigella cells can pose a threat to the health of individuals.

Microbiological Methods Used in the Enterics for Global Health Shigella Surveillance Study

Background

Shigella is a major cause of diarrhea in young children worldwide. Multiple vaccines targeting Shigella are in development, and phase 3 clinical trials are imminent to determine efficacy against shigellosis.

Methods

The Enterics for Global Health (EFGH) Shigella surveillance study is designed to determine the incidence of medically attended shigellosis in 6- to 35-month-old children in 7 resource-limited settings. Here, we describe the microbiological methods used to isolate and identify Shigella. We developed a standardized laboratory protocol for isolation and identification of Shigella by culture. This protocol was implemented across all 7 sites, ensuring consistency and comparability of results. Secondary objectives of the study are to determine the antibiotic resistance profiles of Shigella, compare isolation of Shigella from rectal swabs versus whole stool, and compare isolation of Shigella following transport of rectal swabs in Cary-Blair versus a modified buffered glycerol saline transport medium.

Conclusions

Data generated from EFGH using culture methods described herein can potentially be used for microbiological endpoints in future phase 3 clinical trials to evaluate vaccines against shigellosis and for other clinical and public health studies focused on these organisms.

Genomic, transcriptomic, and phenotypic differences among archetype Shigella flexneri strains of serotypes 2a, 3a, and 6

IMPORTANCE

Given the genomic diversity between S. flexneri serotypes and the paucity of data to support serotype-specific phenotypic differences, we applied in silico and in vitro functional analyses of archetype strains of 2457T (Sf2a), J17B (Sf3a), and CH060 (Sf6). These archetype strains represent the three leading S. flexneri serotypes recommended for inclusion in multivalent vaccines. Characterizing the genomic and phenotypic variation among these clinically prevalent serotypes is an important step toward understanding serotype-specific host-pathogen interactions to optimize the efficacy of multivalent vaccines and therapeutics. This study underpins the importance for further large-scale serotype-targeted analyses.

The L-DBF vaccine cross protects mice against different Shigella serotypes after prior exposure to the pathogen

IMPORTANCE

Shigellosis is endemic to low- and middle-income regions of the world where children are especially vulnerable. In many cases, there are pre-existing antibodies in the local population and the effect of prior exposure should be considered in the development and testing of vaccines against Shigella infection. Our study shows that L-DBF-induced immune responses are not adversely affected by prior exposure to this pathogen. Moreover, somewhat different cytokine profiles were observed in the lungs of vaccinated mice not having been exposed to Shigella, suggesting that the immune responses elicited by Shigella infection and L-DBF vaccination follow different pathways.

Genomic analysis of Shigella isolates from Lebanon reveals marked genetic diversity and antimicrobial resistance

In this study, we characterized 54 clinical isolates of Shigella collected in North Lebanon between 2009 and 2017 through phenotypic and genomic analyses. The most prevalent serogroup was S. sonnei, accounting for 46.3 % (25/54) of the isolates, followed by S. flexneri (27.8 %, 15/54), S. boydii (18.5 %, 10/54) and S. dysenteriae (7.4 %, 4/54). Only three isolates were pan-susceptible, and 87 % (47/54) of the isolates had multidrug resistance phenotypes. Notably, 27.8 % (15/54) of the isolates were resistant to third-generation cephalosporins (3GCs) and 77.8 % (42/54) were resistant to nalidixic acid. 3GC resistance was mediated by the extended-spectrum beta-lactamase genes bla CTX-M-15 and bla CTX-M-3, which were present on various plasmids. Quinolone resistance was conferred by single point mutations in the gyrA DNA gyrase gene, leading to GyrA S83L, GyrA D87Y or GyrA S83A amino acid substitutions. This is the first study, to our knowledge, to provide genomic insights into the serotypes of Shigella circulating in Lebanon and the various antimicrobial resistance determinants carried by these strains.

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.

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.

Characterization of Shigella flexneri in northern Vietnam in 2012-2016

Introduction

Shigellosis remains a considerable public health concern in developing countries. Shigella flexneri and Shigella sonnei are prevalent worldwide and S. sonnei has been replacing S. flexneri .

Gap Statement

S. flexneri still causes outbreaks of shigellosis in northern Vietnam but limited information is available on its genetic characteristics.

Aim

This study aimed to characterize the genetic characteristics of S. flexneri strains from northern Vietnam.

Methodology

This study used 17 isolates from eight incidents, collected in northern Vietnam between 2012 and 2016. The samples were subjected to whole genome sequencing, molecular serotyping, cluster analysis and identification of antimicrobial resistance genes. Additionally, phylogenetic analysis was performed including isolates from previous studies.

Results

Clusters were identified according to spatiotemporal backgrounds. The results suggested that two incidents in Yen Bai province in 2015 and 2016 were derived from a very recent common ancestor. All isolates belonged to phylogroup (PG) 3, which was divided into two sub-lineages. Thirteen of 17 isolates, including those from the Yen Bai incidents, belonged to sub-lineage Sub-1 and were serotyped as 1a. The remaining four isolates belonged to sub-lineage Sub-2 and were the globally predominant serotype 2a. The Sub-1 S. flexneri isolates possessed the gtrI gene, which encodes the glycosyl transferase that determines serotype 1a, with bacteriophage elements in the vicinity.

Conclusion

This study revealed two PG3 sub-lineages of S. flexneri in northern Vietnam, of which Sub-1 might be specific to the region.

Genomic reconstruction and directed interventions in a multidrug-resistant Shigellosis outbreak in Seattle, WA, USA: a genomic surveillance study

BACKGROUND

Shigella spp have been associated with community-wide outbreaks in urban settings. We analysed a sustained shigellosis outbreak in Seattle, WA, USA, to understand its origins and mechanisms of antimicrobial resistance, define ongoing transmission patterns, and optimise strategies for treatment and infection control.

METHODS

We did a retrospective study of all Shigella isolates identified from stool samples at the clinical laboratories at Harborview Medical Center and University of Washington Medical Center (Seattle, WA, USA) from May 1, 2017, to Feb 28, 2022. We characterised isolates by species identification, phenotypic susceptibility testing, and whole-genome sequencing. Demographic characteristics and clinical outcomes of the patients were retrospectively examined.

FINDINGS

171 cases of shigellosis were included. 78 (46%) patients were men who have sex with men (MSM), and 88 (52%) were people experiencing homelessness (PEH). Although 84 (51%) isolates were multidrug resistant, 100 (70%) of 143 patients with data on antimicrobial therapy received appropriate empirical therapy. Phylogenomic analysis identified sequential outbreaks of multiple distinct lineages of Shigella flexneri and Shigella sonnei. Discrete clonal lineages (ten in S flexneri and nine in S sonnei) and resistance traits were responsible for infection in different at-risk populations (ie, MSM, PEH), enabling development of effective guidelines for empirical treatment. The most prevalent lineage in Seattle was probably introduced to Washington State via international travel, with subsequent domestic transmission between at-risk groups.

INTERPRETATION

An outbreak in Seattle was driven by parallel emergence of multidrug-resistant strains involving international transmission networks and domestic transmission between at-risk populations. Genomic analysis elucidated not only outbreak origin, but directed optimal approaches to testing, treatment, and public health response. Rapid diagnostics combined with detailed knowledge of local epidemiology can enable high rates of appropriate empirical therapy even in multidrug-resistant infection.

FUNDING

None.

Characterization of Shigella flexneri Serotype 6 Strains Isolated from Bangladesh and Identification of a New Phylogenetic Cluster

A significant cause of shigellosis in Bangladesh and other developing countries is Shigella flexneri serotype 6. This serotype has been subtyped, on the basis of the absence or presence of a group-specific antigen, E1037, into S. flexneri 6a and 6b, respectively. Here, we provided rationales for the subclassification, using several phenotypic and molecular tools. A set of S. flexneri 6a and 6b strains isolated between 1997 and 2015 were characterized by analyzing their biochemical properties, plasmid profiles, virulence markers, pulsed-field gel electrophoresis (PFGE) results, and ribotype. Additionally, the genomic relatedness of these subserotypes was investigated with global isolates of serotype 6 using publicly available genomes. Both subserotypes of S. flexneri 6 agglutinated with monoclonal antiserum against S. flexneri (MASF) B and type VI-specific antiserum (MASF VI) and were PCR positive for O-antigen flippase-specific genes and virulence markers (ipaH, ial, sen, and sigA). Unlike S. flexneri 6a strains, S. flexneri 6b strains seroagglutinated with anti-E1037 antibodies, MASF IV-I. Notably, these two antigenically distinct subserotypes were clonally diverse, showing two distinct PFGE patterns following the digestion of chromosomal DNA with either XbaI or IceuI. In addition, hybridization of a 16S rRNA gene probe with HindIII-digested genomic DNA yielded two distinguishing ribotypes. Genomic comparison of S. flexneri subserotype 6a and 6b strains from Bangladesh indicated that, although these strains were in genomic synteny, the majority of them formed a unique phylogroup (PG-4) that was missing for the global isolates. This study supports the subserotyping and emphasizes the need for global monitoring of the S. flexneri subserotypes 6a and 6b. IMPORTANCE Shigella flexneri serotype 6 is one of the predominant serotypes among shigellosis cases in Bangladesh. Characterization of a novel subserotype of S. flexneri 6 (VI:E1037), agglutinated with type 6-specific antibody and anti-E1037, indicates a unique evolutionary ancestry. PFGE genotyping supports the finding that these two antigenically distinct subserotypes are clonally diverse. A phylogenetic study based on single-nucleotide polymorphism (SNP) data revealed that these two subserotypes were in genomic synteny, although their genomes were reduced. Interestingly, a majority of the S. flexneri 6 strains isolated from Bangladesh form a novel phylogenetic cluster. Therefore, this report underpins the global monitoring and tracking of the novel subserotype.

Influence of Physicochemical Factors on Adsorption of Ten Shigella&nbsp;flexneri Phages

Bacterial viruses known as bacteriophages have been demonstrated to be effective in killing foodborne pathogens such as Shigella flexneri. Adsorption is the first step in the phage-host interaction. In the present work, 10 Shigella phages were used to characterize the adsorption process on Shigella flexneri ATCC12022 in several physicochemical conditions related to food and in a food matrix. One-step growth curves were drawn for all the Shigella-phages evaluated. Furthermore, the adsorption rate for each of the 10 phages was determined. In addition, the influence of temperature, Na⁺, Mg²⁺, pH, sucrose and glycerol on phage adsorption was investigated. Two phages (Shi22 and Shi30) showed higher burst sizes values (67 and 64 PFU cell^-1, respectively) and burst times of 25 min to 30 min, while the other eight phages exhibited burst sizes ranging from 14 to 17 PFU cell^-1 with slower burst times. Furthermore, most phages achieved a high adsorption rate, and the adsorption constants (k) ranged from ~10^-9 to 10^-10 mL min^-1. Regarding the influence of temperature, cations and pH, a high or moderate percentage of adsorption was observed for most of the phages evaluated. The adsorption decreased at increasing concentrations of Na⁺, sucrose and glycerol, although at different levels, since adsorption was more affected by sucrose than by glycerol and Na⁺ for most phages. The adsorption obtained in Triptein soy broth (TSB) for most of the phages/strain systems evaluated was moderate or high, as well as those observed in a food matrix. Thus, our phages could potentially be used to improve food safety under a wide range of environmental conditions against foodborne pathogens.

Evolutionary and genomic insights into the long-term colonization of in animals

The enteroinvasive bacterium Shigella flexneri is known as a highly host-adapted human pathogen. There had been no known other reservoirs reported until recently. Here 34 isolates obtained from animals (yaks, dairy cows and beef cattle) from 2016-2017 and 268 human S. flexneri isolates from China were sequenced to determine the relationships between animal and human isolates and infer the evolutionary history of animal-associated S. flexneri. The 18 animal isolates (15 yak and 3 beef cattle isolates) in PG1 were separated into 4 lineages, and the 16 animal isolates (1 yak, 5 beef cattle and 10 dairy cow isolates) in PG3 were clustered in 8 lineages. The most recent human isolates from China belonged to PG3 whereas Chinese isolates from the 1950s-1960s belonged to PG1. PG1 S. flexneri may has been transmitted to the yaks during PG1 circulation in the human population in China and has remained in the yak population since, while PG3 S. flexneri in animals were likely recent transmissions from the human population. Increased stability of the large virulence plasmid and acquisition of abundant antimicrobial resistance determinants may have enabled PG3 to expand globally and replaced PG1 in China. Our study confirms that animals may act as a reservoir for S. flexneri. Genomic analysis revealed the evolutionary history of multiple S. flexneri lineages in animals and humans in China. However, further studies are required to determine the public health threat of S. flexneri from animals.

Host Range Expansion of Shigella Phage Sf6 Evolves through Point Mutations in the Tailspike

The first critical step in a virus’s infection cycle is attachment to its host. This interaction is precise enough to ensure the virus will be able to productively infect the cell, but some flexibility can be beneficial to enable coevolution and host range switching or expansion. Bacteriophage Sf6 utilizes a two-step process to recognize and attach to its host Shigella flexneri. Sf6 first recognizes the lipopolysaccharide (LPS) of S. flexneri and then binds outer membrane protein (Omp) A or OmpC. This phage infects serotype Y strains but can also form small, turbid plaques on serotype 2a₂; turbid plaques appear translucent rather than transparent, indicating greater survival of bacteria. Reduced plating efficiency further suggested inefficient infection. To examine the interactions between Sf6 and this alternate host, phages were experimentally evolved using mixed populations of S. flexneri serotypes Y and 2a₂. The recovered mutants could infect serotype 2a₂ with greater efficiency than the ancestral Sf6, forming clear plaques on both serotypes. All mutations mapped to two distinct regions of the receptor-binding tailspike protein: (i) adjacent to the LPS binding site near the N terminus; and (ii) at the distal, C-terminal tip of the protein. Although we anticipated interactions between the Sf6 tailspike and 2a₂ O-antigen to be weak, LPS of this serotype appears to inhibit infection through strong binding of particles, effectively removing them from the environment. The mutations of the evolved strains reduce the inhibitory effect by either reducing electrostatic interactions with the O-antigen or increasing reliance on the Omp secondary receptors.

IMPORTANCE Viruses depend on host cells to propagate themselves. In mixed populations and communities of host cells, finding these susceptible host cells may have to be balanced with avoiding nonhost cells. Alternatively, being able to infect new cell types can increase the fitness of the virus. Many bacterial viruses use a two-step process to identify their hosts, binding first to an LPS receptor and then to a host protein. For Shigella virus Sf6, the tailspike protein was previously known to bind the LPS receptor. Genetic data from this work imply the tailspike also binds to the protein receptor. By experimentally evolving Sf6, we also show that point mutations in this protein can dramatically affect the binding of one or both receptors. This may provide Sf6 flexibility in identifying host cells and the ability to rapidly alter its host range under selective pressure.

Genomic and proteomic characterization of two strains of Shigella flexneri 2 isolated from infants' stool samples in Argentina

Background

Shigella specie is a globally important intestinal pathogen disseminated all over the world. In this study we analyzed the genome and the proteomic component of two Shigella flexneri 2a clinical isolates, collected from pediatric patients with gastroenteritis of the Northwest region of Argentina (NWA) in two periods of time, with four years of difference. Our goal was to determine putative changes at molecular levels occurred during these four years, that could explain the presence of this Shigella`s serovar as the prevalent pathogen in the population under study.

Results

As previously reported, our findings support the idea of Shigella has a conserved “core” genome, since comparative studies of CI133 and CI172 genomes performed against 80 genomes obtained from the NCBI database, showed that there is a large number of genes shared among all of them. However, we observed that CI133 and CI172 harbors a small number of strain-specific genes, several of them present in mobile genetic elements, supporting the hypothesis that these isolates were established in the population by horizontal acquisition of genes. These differences were also observed at proteomic level, where it was possible to detect the presence of certain secreted proteins in a culture medium that simulates the host environment.

Conclusion

Great similarities were observed between the CI133 and CI172 strains, confirming the high percentage of genes constituting the “core” genome of S. flexneri 2. However, numerous strain specific genes were also determined. The presence of the here identified molecular elements into other strain of our culture collation, is currently used to develop characteristic markers of local pathogens. In addition, the most outstanding result of this study was the first description of a S. flexneri 2 producing Colicin E, as one of the characteristics that allows S. flexneri 2 to persist in the microbial community. These findings could also contribute to clarify the mechanism and the evolution strategy used by this pathogen to specifically colonize, survive, and cause infection within the NWA population.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08711-5.

Exploiting genomics to mitigate the public health impact of antimicrobial resistance

Antimicrobial resistance (AMR) is a major global public health threat, which has been largely driven by the excessive use of antimicrobials. Control measures are urgently needed to slow the trajectory of AMR but are hampered by an incomplete understanding of the interplay between pathogens, AMR encoding genes, and mobile genetic elements at a microbial level. These factors, combined with the human, animal, and environmental interactions that underlie AMR dissemination at a population level, make for a highly complex landscape. Whole-genome sequencing (WGS) and, more recently, metagenomic analyses have greatly enhanced our understanding of these processes, and these approaches are informing mitigation strategies for how we better understand and control AMR. This review explores how WGS techniques have advanced global, national, and local AMR surveillance, and how this improved understanding is being applied to inform solutions, such as novel diagnostic methods that allow antimicrobial use to be optimised and vaccination strategies for better controlling AMR. We highlight some future opportunities for AMR control informed by genomic sequencing, along with the remaining challenges that must be overcome to fully realise the potential of WGS approaches for international AMR control.

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.

Shigella Outer Membrane Vesicles as Promising Targets for Vaccination

The clinical symptoms of shigellosis, a gastrointestinal infection caused by Shigella spp. range from watery diarrhea to fulminant dysentery. Endemic infections, particularly among children in developing countries, represent the majority of clinical cases. The situation is aggravated due to the high mortality rate of shigellosis, the rapid dissemination of multi-resistant Shigella strains and the induction of only serotype-specific immunity. Thus, infection prevention due to vaccination, encompassing as many of the circulating serotypes as possible, has become a topic of interest. However, vaccines have turned out to be ineffective so far. Outer membrane vesicles (OMVs) are promising novel targets for vaccination. OMVs are constitutively secreted by Gram-negative bacteria including Shigella during growth. They are composed of soluble luminal portions and an insoluble membrane and can contain toxins, bioactive periplasmic and cytoplasmic (lipo-) proteins, (phospho-) lipids, nucleic acids and/or lipopolysaccharides. Thus, OMVs play an important role in bacterial cell-cell communication, growth, survival and pathogenesis. Furthermore, they modulate the secretion and transport of biomolecules, the stress response, antibiotic resistance and immune responses of the host. Thus, OMVs serve as novel secretion machinery. Here, we discuss the current literature and highlight the properties of OMVs as potent vaccine candidates because of their immunomodulatory, antigenic and adjuvant properties.

Linkage of Whole Genome Sequencing, Epidemiological, and Clinical Data to Understand the Genetic Diversity and Clinical Outcomes of Shigella flexneri among Men Who Have Sex with Men in England

The public health value of whole genome sequencing (WGS) for Shigella spp. in England has been limited by a lack of information on sexual identity and behavior. We combined WGS data with other data sources to better understand Shigella flexneri transmission in men who have sex with men (MSM). WGS data for all S. flexneri isolates referred to the national reference laboratory were linked to i) clinical and behavioral data collected in seven of 21 health regions in England using a standardized exposure questionnaire and, ii) national HIV surveillance data. We included 926 S. flexneri isolates, of which 43.0% (n = 398) fell phylogenetically within two domestically circulating clades associated with genotypic markers of azithromycin resistance. Approximately one third of isolates in these clades were from people living with HIV, primarily acquired through sex between men. 182 (19.7%) isolates had linked questionnaire data; 88% (84/95) of MSM isolates fell phylogenetically within the domestically circulating clades, while 92% (72/78) of isolates from other cases fell within lineages linked with travel to high-risk regions. There was no evidence of sustained transmission between networks of MSM and the wider community. MSM were more likely to be admitted to hospital and receive antimicrobials. Our study emphasizes the importance of sex between men as a major route of transmission for S. flexneri. Combined WGS, epidemiological and clinical data provide unique insights that can inform contact tracing, clinical management and the delivery of targeted prevention activities. Future studies should investigate why MSM experience more severe clinical outcomes. IMPORTANCE Within the last 2 decades there have been an increasing number of Shigella spp. outbreaks among men who have sex with men (MSM) worldwide. In 2015, Public Health England (PHE) introduced routine whole genome sequencing (WGS) for the national surveillance of Shigella spp. However, the lack of information on sexual identity and behavior has hindered interpretation. Our study illustrates the power of linking WGS data with epidemiological, behavioral, and clinical data. We provide unique population-level insights into different transmission networks that can inform the delivery of appropriate public health interventions and patient management. Furthermore, we describe and compare clinical characteristics and outcomes of S. flexneri infection in MSM and other exposure groups. We found that MSM were more likely to be admitted to hospital and receive antimicrobials, indicating that their infections were potentially more severe. The exact reasons for this are unclear and require further exploration.

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.

Cluster-specific gene markers enhance Shigella and enteroinvasive Escherichia coli in silico serotyping

Shigella and enteroinvasive Escherichia coli (EIEC) cause human bacillary dysentery with similar invasion mechanisms and share similar physiological, biochemical and genetic characteristics. Differentiation of Shigella from EIEC is important for clinical diagnostic and epidemiological investigations. However, phylogenetically, Shigella and EIEC strains are composed of multiple clusters and are different forms of E. coli, making it difficult to find genetic markers to discriminate between Shigella and EIEC. In this study, we identified 10 Shigella clusters, seven EIEC clusters and 53 sporadic types of EIEC by examining over 17000 publicly available Shigella and EIEC genomes. We compared Shigella and EIEC accessory genomes to identify cluster-specific gene markers for the 17 clusters and 53 sporadic types. The cluster-specific gene markers showed 99.64% accuracy and more than 97.02% specificity. In addition, we developed a freely available in silico serotyping pipeline named Shigella EIEC Cluster Enhanced Serotype Finder (ShigEiFinder) by incorporating the cluster-specific gene markers and established Shigella and EIEC serotype-specific O antigen genes and modification genes into typing. ShigEiFinder can process either paired-end Illumina sequencing reads or assembled genomes and almost perfectly differentiated Shigella from EIEC with 99.70 and 99.74% cluster assignment accuracy for the assembled genomes and read mapping respectively. ShigEiFinder was able to serotype over 59 Shigella serotypes and 22 EIEC serotypes and provided a high specificity of 99.40% for assembled genomes and 99.38% for read mapping for serotyping. The cluster-specific gene markers and our new serotyping tool, ShigEiFinder (installable package: https://github.com/LanLab/ShigEiFinder, online tool: https://mgtdb.unsw.edu.au/ShigEiFinder/), will be useful for epidemiological and diagnostic investigations.

Historical, current, and emerging tools for identification and serotyping of Shigella

The Shigella genus includes serious foodborne disease etiologic agents, with 4 species and 54 serotypes. Identification at species and serotype levels is a crucial task in microbiological laboratories. Nevertheless, the genetic similarity between Shigella spp. and Escherichia coli challenges the correct identification and serotyping of Shigella spp., with subsequent negative repercussions on surveillance, epidemiological investigations, and selection of appropriate treatments. For this purpose, multiple techniques have been developed historically ranging from phenotype-based methods and single or multilocus molecular techniques to whole-genome sequencing (WGS). To facilitate the selection of the most relevant method, we herein provide a global overview of historical and emerging identification and serotyping techniques with a particular focus on the WGS-based approaches. This review highlights the excellent discriminatory power of WGS to more accurately elucidate the epidemiology of Shigella spp., disclose novel promising genomic targets for surveillance methods, and validate previous well-established methods.

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.

Large metabolic rewiring from small genomic changes between strains of Shigella flexneri

The instability of Shigella genomes has been described, but how this instability causes phenotypic differences within the Shigella flexneri species is largely unknown and likely variable. We describe herein the genome of S. flexneri strain PE577, originally a clinical isolate, which exhibits several phenotypic differences compared to the model strain 2457T. Like many previously described strains of S. flexneri, PE577 lacks discernible, functional CRISPR and restriction-modification systems. Its phenotypic differences when compared to 2457T include lower transformation efficiency, higher oxygen sensitivity, altered carbon metabolism, and greater susceptibility to a wide variety of lytic bacteriophage isolates. Since relatively few Shigella phages have been isolated on 2457T or the previously characterized strain M90T, developing a more universal model strain for isolating and studying Shigella phages is critical to understanding both phages and phage-host interactions. In addition to phage biology, the genome sequence of PE577 was used to generate and test hypotheses of how pseudogenes in this strain-whether interrupted by degraded prophages, transposases, frameshifts, or point mutations-have led to metabolic rewiring compared to the model strain 2457T. Results indicate that PE577 can utilise the less-efficient pyruvate oxidase/acetyl-CoA synthetase (PoxB/Acs) pathway to produce acetyl-CoA, while strain 2457T cannot due to a nonsense mutation in acs, rendering it a pseudogene in this strain. Both strains also utilize pyruvate-formate lyase to oxidize formate but cannot survive with this pathway alone, possibly because a component of the formate-hydrogen lyase (fdhF) is a pseudogene in both strains.Importance Shigella causes millions of dysentery cases worldwide, primarily affecting children under five years old. Despite active research in developing vaccines and new antibiotics, relatively little is known about the variation of physiology or metabolism across multiple isolates. In this work, we investigate two strains of S. flexneri that share 98.9% genetic identity but exhibit drastic differences in metabolism, ultimately affecting the growth of the two strains. Results suggest additional strains within the S. flexneri species utilize different metabolic pathways to process pyruvate. Metabolic differences between these closely-related isolates suggest an even wider variety of differences in growth across S. flexneri and Shigella in general. Exploring this variation further may assist the development or application of vaccines and therapeutics to combat Shigella infections.

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.

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.

Diarrhea in an infant due to Shigella flexneri 1 carrying multiple cephalosporinase-encoding genes

Background

Infections caused by multidrug-resistant shigellae resistant to broad-spectrum cephalosporins are becoming more prevalent in the Middle East. We report a case of severe diarrhea due to a multiresistant Shigella flexneri 1 strain carrying four different ß-lactamase genes.

Case presentation

A one-year-old Syrian infant presented with severe acute diarrhea, vomiting and dehydration. She did not respond to empirical treatment with amoxicillin-clavulanic acid followed by cefotaxime. Later, stool culture revealed S. flexneri 1 resistant to both these drugs. The patient was successfully treated with meropenem to which S. flexneri 1 was susceptible. The isolate was resistant to eight classes of antibiotics, and the whole genome sequence (WGS) identified four ß-lactamase genes (bla_CTX-M-15, bla_EC-8, bla_OXA-1, and bla_TEM-1) along with genes mediating resistance to seven other antibiotic classes. The WGS also identified several virulence genes including senA that encodes ShET-2 which induces watery diarrhea. Phylogenetically, the isolate was closely related to isolates from South Asia.

Conclusions

This report highlights the emergence of extremely resistant Shigella that has acquired multiple resistance genes to cephalosporins rendering these drugs ineffective.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13099-021-00413-9.

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.

Biodiversity of New Lytic Bacteriophages Infecting Shigella spp. in Freshwater Environment

Bacteriophages, viruses that infect and replicate within prokaryotic cells are the most abundant life forms in the environment, yet the vast majority of them have not been properly reported or even discovered. Almost all reported bacteriophages infecting the Enterobacteriaceae family, with Escherichia coli being the major subject of studies, have been isolated from wastewater, sewage, and effluent resources. In the present study, we focused on the distribution and biodiversity of Shigella phages in an aquatic ecosystem. While no Shigella bacteria was recovered from the Yangtze River, three lytic phages were isolated from this ecosystem and were subjected to biological, morphological, and genomic characteristics. Comparative genomics and phylogenetic analyses demonstrated that vB _SflM_004 isolate belongs to Myoviridae family, Felixounavirus genus of Ounavirinae subfamily, vB_SdyM_006 was classified under the same family, however, it is suggested to be in a new genus under Tevenvirinae subfamily with some other related bacteriophages. vB_SsoS_008 phage belongs to the Siphoviridae family, Tunavirus genus, Tunavirinae subfamily. The phages did not harbor any genes involved in the lysogenic cycles and showed a high temperature and pH stability. The biodiversity of the isolated phages highly suggests that continued isolation on non-model members of Enterobacteriaceae family is necessary to fully understand bacteriophage diversity in aquatic environments.

Outbreak of multi-drug-resistant (MDR) Shigella flexneri in northern Australia due to an endemic regional clone acquiring an IncFII plasmid

Epidemiological surveillance of Shigella spp. in Australia is conducted to inform public health response. Multi-drug resistance has recently emerged as a contributing factor to sustained local transmission of Shigella spp. All data were collected as part of routine public health surveillance, and strains were whole-genome sequenced for further molecular characterisation. 108 patients with an endemic regional Shigella flexneri strain were identified between 2016 and 2019. The S. flexneri phylogroup 3 strain endemic to northern Australia acquired a multi-drug resistance conferring blaDHA plasmid, which has an IncFII plasmid backbone with virulence and resistance elements typically found in IncR plasmids. This is the first report of multi-drug resistance in Shigella sp. in Australia that is not associated with men who have sex with men. This strain caused an outbreak of multi-drug-resistant S. flexneri in northern Australia that disproportionality affects Aboriginal and Torres Strait Islander children. Community controlled public health action is recommended.

Evaluation of Molecular Serotyping Assays for Shigella flexneri Directly on Stool Samples

Shigella flexneri is prevalent worldwide and is the most common Shigella species in many countries. At least 19 S. flexneri serotypes exist, and serotype information is important for epidemiologic and vaccine development purposes.

Shigella flexneri is prevalent worldwide and is the most common Shigella species in many countries. At least 19 S. flexneri serotypes exist, and serotype information is important for epidemiologic and vaccine development purposes. We evaluated the performance of real-time PCR assays for O-antigen modification genes to identify the major serotypes on isolates and direct stool samples. The assays were formulated into two multiplex panels: one panel included gtrII, gtrV, gtrX, oac, and wzx6 to identify S. flexneri serotypes 2a, 2b, 3a, 5a, 5b, 6, and X, and the other panel included ipaH, gtrI, gtrIc, and gtrIV to confirm Shigella detection and further identify S. flexneri serotypes 1a, 1b, 1d, 3b, 4a, 4b, 7a, and 7b. We first evaluated 283 Shigella isolates, and PCR serotyping demonstrated 97.0% (95% confidence interval, 93.0% to 99.0%) sensitivity and 99.9% (99.9% to 100%) specificity compared to conventional serotyping. The assays then were utilized on direct stool specimens. A quantitative detection algorithm was developed with a validation set of 174 Shigella culture-positive stool samples and further tested with a derivation set of 164 samples. The PCR serotyping on stool achieved 93% (89% to 96%) sensitivity and 99% (99% to 100%) specificity compared to serotyping. Most discrepancies were genotypic-phenotypic discordance, not genotypic failure. These real-time PCR assays provide an efficient and novel tool for S. flexneri serotype identification.

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

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.