Defining the phylogenomics of Shigella species: a pathway to diagnostics

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.

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.

Acquisition of a large virulence plasmid (pINV) promoted temperature-dependent virulence and global dispersal of O96:H19 enteroinvasive Escherichia coli

Enteroinvasive Escherichia coli (EIEC) and Shigella are closely related agents of bacillary dysentery. It is widely viewed that EIEC and Shigella species evolved from E. coli via independent acquisitions of a large virulence plasmid (pINV) encoding a type 3 secretion system (T3SS). Sequence Type (ST)99 O96:H19 E. coli is a novel clone of EIEC responsible for recent outbreaks in Europe and South America. Here, we use 92 whole genome sequences to reconstruct a dated phylogeny of ST99 E. coli, revealing distinct phylogenomic clusters of pINV-positive and -negative isolates. To study the impact of pINV acquisition on the virulence of this clone, we developed an EIEC-zebrafish infection model showing that virulence of ST99 EIEC is thermoregulated. Strikingly, zebrafish infection using a T3SS-deficient ST99 EIEC strain and the oldest available pINV-negative isolate reveals a separate, temperature-independent mechanism of virulence, indicating that ST99 non-EIEC strains were virulent before pINV acquisition. Taken together, these results suggest that an already pathogenic E. coli acquired pINV and that virulence of ST99 isolates became thermoregulated once pINV was acquired. IMPORTANCE Enteroinvasive Escherichia coli (EIEC) and Shigella are etiological agents of bacillary dysentery. Sequence Type (ST)99 is a clone of EIEC hypothesized to cause human disease by the recent acquisition of pINV, a large plasmid encoding a type 3 secretion system (T3SS) that confers the ability to invade human cells. Using Bayesian analysis and zebrafish larvae infection, we show that the virulence of ST99 EIEC isolates is highly dependent on temperature, while T3SS-deficient isolates encode a separate temperature-independent mechanism of virulence. These results indicate that ST99 non-EIEC isolates may have been virulent before pINV acquisition and highlight an important role of pINV acquisition in the dispersal of ST99 EIEC in humans, allowing wider dissemination across Europe and South America.

Sequence-based detection and typing procedures for Burkholderia mallei: Assessment and prospects

Although glanders has been eradicated in most of the developed world, the disease still persists in various countries such as Brazil, India, Pakistan, Bangladesh, Nepal, Iran, Bahrain, UAE and Turkey. It is one of the notifiable diseases listed by the World Organization for Animal Health. Occurrence of glanders imposes restriction on equestrian events and restricts equine movement, thus causing economic losses to equine industry. The genetic diversity and global distribution of the causing agent, Burkholderia (B.) mallei, have not been assessed in detail and are complicated by the high clonality of this organism. Among the identification and typing methods, PCR-based methods for distinguishing B. mallei from its close relative B. pseudomallei as well as genotyping using tandem repeat regions (MLVA) are established. The advent and continuous advancement of the sequencing techniques and the reconstruction of closed genomes enable the development of genome guided epidemiological tools. For achieving a higher genomic resolution, genotyping methods based on whole genome sequencing data can be employed, like genome-wide single nucleotide polymorphisms. One of the limitations in obtaining complete genomic sequences for further molecular characterization of B. mallei is its high GC content. In this review, we aim to provide an overview of the widely used detection and typing methods for B. mallei and illustrate gaps that still require development. The genomic features of Burkholderia, their high homology and clonality will be first described from a comparative genomics perspective. Then, the commonly used molecular detection (PCR systems) and typing systems (e.g., multilocus sequence typing, variable number of tandem repeat analysis) will be presented and put in perspective with recently developed genomic methods. Also, the increasing availability of B. mallei genomic sequences and evolution of the sequencing methods offers exciting prospects for further refinement of B. mallei typing, that could overcome the difficulties presently encountered with this particular bacterium.

Genomics and pathotypes of the many faces of Escherichia coli

Escherichia coli is the most researched microbial organism in the world. Its varied impact on human health, consisting of commensalism, gastrointestinal disease, or extraintestinal pathologies, has generated a separation of the species into at least eleven pathotypes (also known as pathovars). These are broadly split into two groups, intestinal pathogenic E. coli (InPEC) and extraintestinal pathogenic E. coli (ExPEC). However, components of E. coli's infinite open accessory genome are horizontally transferred with substantial frequency, creating pathogenic hybrid strains that defy a clear pathotype designation. Here, we take a birds-eye view of the E. coli species, characterizing it from historical, clinical, and genetic perspectives. We examine the wide spectrum of human disease caused by E. coli, the genome content of the bacterium, and its propensity to acquire, exchange, and maintain antibiotic resistance genes and virulence traits. Our portrayal of the species also discusses elements that have shaped its overall population structure and summarizes the current state of vaccine development targeted at the most frequent E. coli pathovars. In our conclusions, we advocate streamlining efforts for clinical reporting of ExPEC, and emphasize the pathogenic potential that exists throughout the entire species.

MALDI-TOF MS Using a Custom-Made Database, Biomarker Assignment, or Mathematical Classifiers Does Not Differentiate Shigella spp. and Escherichia coli

Shigella spp. and E. coli are closely related and cannot be distinguished using matrix-assisted laser desorption-ionization time-of-flight mass spectrometry (MALDI-TOF MS) with commercially available databases. Here, three alternative approaches using MALDI-TOF MS to identify and distinguish Shigella spp., E. coli, and its pathotype EIEC were explored and evaluated using spectra of 456 Shigella spp., 42 E. coli, and 61 EIEC isolates. Identification with a custom-made database resulted in >94% Shigella identified at the genus level and >91% S. sonnei and S. flexneri at the species level, but the distinction of S. dysenteriae, S. boydii, and E. coli was poor. With biomarker assignment, 98% S. sonnei isolates were correctly identified, although specificity was low. Discriminating markers for S. dysenteriae, S. boydii, and E. coli were not assigned at all. Classification models using machine learning correctly identified Shigella in 96% of isolates, but most E. coli isolates were also assigned to Shigella. None of the proposed alternative approaches were suitable for clinical diagnostics for identifying Shigella spp., E. coli, and EIEC, reflecting their relatedness and taxonomical classification. We suggest the use of MALDI-TOF MS for the identification of the Shigella spp./E. coli complex, but other tests should be used for distinction.

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.

Efficient production of immunologically active Shigella invasion plasmid antigens IpaB and IpaH using a cell-free expression system

Abstract

Shigella spp. invade the colonic epithelium and cause bacillary dysentery in humans. Individuals living in areas that lack access to clean water and sanitation are the most affected. Even though infection can be treated with antibiotics, Shigella antimicrobial drug resistance complicates clinical management. Despite decades of effort, there are no licensed vaccines to prevent shigellosis. The highly conserved invasion plasmid antigens (Ipa), which are components of the Shigella type III secretion system, participate in bacterial epithelial cell invasion and have been pursued as vaccine targets. However, expression and purification of these proteins in conventional cell-based systems have been challenging due to solubility issues and extremely low recovery yields. These difficulties have impeded manufacturing and clinical advancement. In this study, we describe a new method to express Ipa proteins using the Xpress^+TM cell-free protein synthesis (CFPS) platform. Both IpaB and the C-terminal domain of IpaH1.4 (IpaH-CTD) were efficiently produced with this technology at yields > 200 mg/L. Furthermore, the expression was linearly scaled in a bioreactor under controlled conditions, and proteins were successfully purified using multimode column chromatography to > 95% purity as determined by SDS-PAGE. Biophysical characterization of the cell-free synthetized IpaB and IpaH-CTD using SEC-MALS analysis showed well-defined oligomeric states of the proteins in solution. Functional analysis revealed similar immunoreactivity as compared to antigens purified from E. coli. These results demonstrate the efficiency of CFPS for Shigella protein production; the practicality and scalability of this method will facilitate production of antigens for Shigella vaccine development and immunological analysis.

Key points

• First report of Shigella IpaB and IpaH produced at high purity and yield using CFPS

• CFPS-IpaB and IpaH perform similarly to E. coli–produced proteins in immunoassays

• CFPS-IpaB and IpaH react with Shigella-specific human antibodies and are immunogenic in mice.

Graphical abstract

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.

Phylum barrier and Escherichia coli intra-species phylogeny drive the acquisition of antibiotic-resistance genes

Escherichia coli is a ubiquitous bacterium that has been widely exposed to antibiotics over the last 70 years. It has adapted by acquiring different antibiotic-resistance genes (ARGs), the census of which we aim to characterize here. To do so, we analysed 70 301 E. coli genomes obtained from the EnteroBase database and detected 1 027 651 ARGs using the AMRFinder, Mustard and ResfinderFG ARG databases. We observed a strong phylogroup and clonal lineage specific distribution of some ARGs, supporting the argument for epistasis between ARGs and the strain genetic background. However, each phylogroup had ARGs conferring a similar antibiotic class resistance pattern, indicating phenotypic adaptive convergence. The G+C content or the type of ARG was not associated with the frequency of the ARG in the database. In addition, we identified ARGs from anaerobic, non-Proteobacteria bacteria in four genomes of E. coli, supporting the hypothesis that the transfer between anaerobic bacteria and E. coli can spontaneously occur but remains exceptional. In conclusion, we showed that phylum barrier and intra-species phylogenetic history are major drivers of the acquisition of a resistome in E. coli.

A Fork Trap in the Chromosomal Termination Area Is Highly Conserved across All Escherichia coli Phylogenetic Groups

Termination of DNA replication, the final stage of genome duplication, is surprisingly complex, and failures to bring DNA synthesis to an accurate conclusion can impact genome stability and cell viability. In Escherichia coli, termination takes place in a specialised termination area opposite the origin. A 'replication fork trap' is formed by unidirectional fork barriers via the binding of Tus protein to genomic ter sites. Such a fork trap system is found in some bacterial species, but it appears not to be a general feature of bacterial chromosomes. The biochemical properties of fork trap systems have been extensively characterised, but little is known about their precise physiological roles. In this study, we compare locations and distributions of ter terminator sites in E. coli genomes across all phylogenetic groups, including Shigella. Our analysis shows that all ter sites are highly conserved in E. coli, with slightly more variability in the Shigella genomes. Our sequence analysis of ter sites and Tus proteins shows that the fork trap is likely to be active in all strains investigated. In addition, our analysis shows that the dif chromosome dimer resolution site is consistently located between the innermost ter sites, even if rearrangements have changed the location of the innermost termination area. Our data further support the idea that the replication fork trap has an important physiological role that provides an evolutionary advantage.

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.

The Evolutionary Conservation of Escherichia coli Drug Efflux Pumps Supports Physiological Functions

Bacteria harness an impressive repertoire of resistance mechanisms to evade the inhibitory action of antibiotics. One such mechanism involves efflux pump-mediated extrusion of drugs from the bacterial cell, which significantly contributes to multidrug resistance. Intriguingly, most drug efflux pumps are chromosomally encoded components of the intrinsic antibiotic resistome. In addition, in terms of xenobiotic detoxification, bacterial efflux systems often exhibit significant levels of functional redundancy. Efflux pumps are also considered to be highly conserved; however, the extent of conservation in many bacterial species has not been reported and the majority of genes that encode efflux pumps appear to be dispensable for growth. These observations, in combination with an increasing body of experimental evidence, imply alternative roles in bacterial physiology. Indeed, the ability of efflux pumps to facilitate antibiotic resistance could be a fortuitous by-product of ancient physiological functions. Using Escherichia coli as a model organism, we here evaluated the evolutionary conservation of drug efflux pumps and we provide phylogenetic analysis of the major efflux families. We show the E. coli drug efflux system has remained relatively stable and the majority (∼80%) of pumps are encoded in the core genome. This analysis further supports the importance of drug efflux pumps in E. coli physiology. In this review, we also provide an update on the roles of drug efflux pumps in the detoxification of endogenously synthesized substrates and pH homeostasis. Overall, gaining insight into drug efflux pump conservation, common evolutionary ancestors, and physiological functions could enable strategies to combat these intrinsic and ancient elements.

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.

The population genetics of pathogenic Escherichia coli

Escherichia coli is a commensal of the vertebrate gut that is increasingly involved in various intestinal and extra-intestinal infections as an opportunistic pathogen. Numerous pathotypes that represent groups of strains with specific pathogenic characteristics have been described based on heterogeneous and complex criteria. The democratization of whole-genome sequencing has led to an accumulation of genomic data that render possible a population phylogenomic approach to the emergence of virulence. Few lineages are responsible for the pathologies compared with the diversity of commensal strains. These lineages emerged multiple times during E. coli evolution, mainly by acquiring virulence genes located on mobile elements, but in a specific chromosomal phylogenetic background. This repeated emergence of stable and cosmopolitan lineages argues for an optimization of strain fitness through epistatic interactions between the virulence determinants and the remaining genome.

Impact of insertion sequences on convergent evolution of Shigella species

Shigella species are specialised lineages of Escherichia coli that have converged to become human-adapted and cause dysentery by invading human gut epithelial cells. Most studies of Shigella evolution have been restricted to comparisons of single representatives of each species; and population genomic studies of individual Shigella species have focused on genomic variation caused by single nucleotide variants and ignored the contribution of insertion sequences (IS) which are highly prevalent in Shigella genomes. Here, we investigate the distribution and evolutionary dynamics of IS within populations of Shigella dysenteriae Sd1, Shigella sonnei and Shigella flexneri. We find that five IS (IS1, IS2, IS4, IS600 and IS911) have undergone expansion in all Shigella species, creating substantial strain-to-strain variation within each population and contributing to convergent patterns of functional gene loss within and between species. We find that IS expansion and genome degradation are most advanced in S. dysenteriae and least advanced in S. sonnei; and using genome-scale models of metabolism we show that Shigella species display convergent loss of core E. coli metabolic capabilities, with S. sonnei and S. flexneri following a similar trajectory of metabolic streamlining to that of S. dysenteriae. This study highlights the importance of IS to the evolution of Shigella and provides a framework for the investigation of IS dynamics and metabolic reduction in other bacterial species.

Staying out or Going in? The Interplay between Type 3 and Type 5 Secretion Systems in Adhesion and Invasion of Enterobacterial Pathogens

Enteric pathogens rely on a variety of toxins, adhesins and other virulence factors to cause infections. Some of the best studied pathogens belong to the Enterobacterales order; these include enteropathogenic and enterohemorrhagic Escherichia coli, Shigella spp., and the enteropathogenic Yersiniae. The pathogenesis of these organisms involves two different secretion systems, a type 3 secretion system (T3SS) and type 5 secretion systems (T5SSs). The T3SS forms a syringe-like structure spanning both bacterial membranes and the host cell plasma membrane that translocates toxic effector proteins into the cytoplasm of the host cell. T5SSs are also known as autotransporters, and they export part of their own polypeptide to the bacterial cell surface where it exerts its function, such as adhesion to host cell receptors. During infection with these enteropathogens, the T3SS and T5SS act in concert to bring about rearrangements of the host cell cytoskeleton, either to invade the cell, confer intracellular motility, evade phagocytosis or produce novel structures to shelter the bacteria. Thus, in these bacteria, not only the T3SS effectors but also T5SS proteins could be considered "cytoskeletoxins" that bring about profound alterations in host cell cytoskeletal dynamics and lead to pathogenic outcomes.

Molecular characterization of Shigella species isolated from diarrheal patients in Tehran, Iran: phylogenetic typing and its association with virulence gene profiles and a novel description of Shigella invasion associated locus

The present study aims to employ a multiplex PCR-based method for phylogenetic typing of Shigella and determine the frequency of several virulence genes among Shigella phylogenetic clades and species. Species identification, phylogenetic typing of 44 previously diagnosed Shigella isolates, and frequency of virulence genes and loci, virA, virB, virF, ipaBCD, ial, sen, and set1A were investigated through performing several PCR assays. Distribution of virulence genes among Shigella phylogenetic clades and species was determined by the statistical analysis. The identities of 40 isolates out of 44 were confirmed as Shigella, and these isolates were classified in four phylogenetic clades, S1 (7.5%), S2 (52.5%), S3 (20%), and S5 (20%) and 4 species, S. sonnei (52.5%), S. flexneri (22.5%), S. dysenteriae (20%), and S. boydii (5%). The prevalence of virA, virB, virF, ipaBCD, ial, sen, and set1A was determined as 67.5%, 72.5%, 72.5%, 65%, 75%, 40%, and 5%, respectively. The presence of sen, uidA, or set1A was found to be statistically correlated with either of Shigella phylogenetic clades or species. A significant statistically association was also determined between set1A and Shigella phylogenetic clades. Furthermore, the nucleotide sequence of invasion-associated locus (ial) was determined and mapped on Shigella genome through in silico analysis. The current study shows the distribution of Shigella isolates and its key virulence genes within the five recently described phylogenetic clades for the first time in the Asia. This is also the first description of ial nucleotide sequence and its exact location on Shigella genome after its initial identification.

DNA sequencing, genomes and genetic markers of microbes on fruits and vegetables

The development of DNA sequencing technology has provided an effective method for studying foodborne and phytopathogenic microorganisms on fruits and vegetables (F & V). DNA sequencing has successfully proceeded through three generations, including the tens of operating platforms. These advances have significantly promoted microbial whole‐genome sequencing (WGS) and DNA polymorphism research. Based on genomic and regional polymorphisms, genetic markers have been widely obtained. These molecular markers are used as targets for PCR or chip analyses to detect microbes at the genetic level. Furthermore, metagenomic analyses conducted by sequencing the hypervariable regions of ribosomal DNA (rDNA) have revealed comprehensive microbial communities in various studies on F & V. This review highlights the basic principles of three generations of DNA sequencing, and summarizes the WGS studies of and available DNA markers for major bacterial foodborne pathogens and phytopathogenic fungi found on F & V. In addition, rDNA sequencing‐based bacterial and fungal metagenomics are summarized under three topics. These findings deepen the understanding of DNA sequencing and its application in studies of foodborne and phytopathogenic microbes and shed light on strategies for the monitoring of F & V microbes and quality control.

Proposal of a Taxonomic Nomenclature for the Bacillus cereus Group Which Reconciles Genomic Definitions of Bacterial Species with Clinical and Industrial Phenotypes

The Bacillus cereus group comprises numerous closely related species, including bioterrorism agent B. anthracis, foodborne pathogen B. cereus, and biopesticide B. thuringiensis Differentiating organisms capable of causing illness or death from those used in industry is essential for risk assessment and outbreak preparedness. However, current species definitions facilitate species-phenotype incongruences, particularly when horizontally acquired genes are responsible for a phenotype. Using all publicly available B. cereus group genomes (n = 2,231), we show that current species definitions lead to overlapping genomospecies clusters, in which 66.2% of genomes belong to multiple genomospecies at a conventional 95 average nucleotide identity (ANI) genomospecies threshold. A genomospecies threshold of ≈92.5 ANI is shown to reflect a natural gap in genome similarity for the B. cereus group, and medoid genomes identified at this threshold are shown to yield resolvable genomospecies clusters with minimal overlap (six of 2,231 genomes assigned to multiple genomospecies; 0.269%). We thus propose a nomenclatural framework for the B. cereus group which accounts for (i) genomospecies using resolvable genomospecies clusters obtained at ≈92.5 ANI, (ii) established lineages of medical importance using a formal collection of subspecies names, and (iii) heterogeneity of clinically and industrially important phenotypes using a formalized and extended collection of biovar terms. We anticipate that the proposed nomenclature will remain interpretable to clinicians, without sacrificing genomic species definitions, which can in turn aid in pathogen surveillance; early detection of emerging, high-risk genotypes; and outbreak preparedness.IMPORTANCE Historical species definitions for many prokaryotes, including pathogens, have relied on phenotypic characteristics that are inconsistent with genome evolution. This scenario forces microbiologists and clinicians to face a tradeoff between taxonomic rigor and clinical interpretability. Using the Bacillus cereus group as a model, a conceptual framework for the taxonomic delineation of prokaryotes which reconciles genomic definitions of species with clinically and industrially relevant phenotypes is presented. The nomenclatural framework outlined here serves as a model for genomics-based bacterial taxonomy that moves beyond arbitrarily set genomospecies thresholds while maintaining congruence with phenotypes and historically important species names.

Incidence, clinical implications and impact on public health of infections with Shigella spp. and entero-invasive Escherichia coli (EIEC): results of a multicenter cross-sectional study in the Netherlands during 2016-2017

Background

Shigella spp. and entero-invasive E. coli (EIEC) use the same invasive mechanism to cause diarrheal diseases. Public health regulations apply only to Shigella spp. infections, but are hampered by the lack of simple methods to distinguish them from EIEC. In the last decades, molecular methods for detecting Shigella spp. and EIEC were implemented in medical microbiological laboratories (MMLs). However, shigellosis cases identified with molecular techniques alone are not notifiable in most countries. Our study investigates the impact of EIEC versus Shigella spp. infections and molecular diagnosed shigellosis versus culture confirmed shigellosis for re-examination of the rationale for the current public health regulations.

Methods

In this multicenter cross-sectional study, fecal samples of patients suspected for gastro-enteritis, referred to 15 MMLs in the Netherlands, were screened by PCR for Shigella spp. or EIEC. Samples were cultured to discriminate between the two pathogens. We compared risk factors, symptoms, severity of disease, secondary infections and socio-economic consequences for (i) culture-confirmed Shigella spp. versus culture-confirmed EIEC cases (ii) culture positive versus PCR positive only shigellosis cases.

Results

In 2016–2017, 777 PCR positive fecal samples with patient data were included, 254 of these were culture-confirmed shigellosis cases and 32 were culture-confirmed EIEC cases. EIEC cases were more likely to report ingestion of contaminated food and were less likely to be men who have sex with men (MSM). Both pathogens were shown to cause serious disease although differences in specific symptoms were observed. Culture-negative but PCR positive cases were more likely report travel or ingestion of contaminated food and were less likely to be MSM than culture-positive cases. Culture-negative cases were more likely to suffer from multiple symptoms. No differences in degree of secondary infections were observed between Shigella spp. and EIEC, and culture-negative and culture-positive cases.

Conclusions

No convincing evidence was found to support the current guidelines that employs different measures based on species or detection method. Therefore, culture and molecular detection methods for Shigella spp. and EIEC should be considered equivalent for case definition and public health regulations regarding shigellosis. Differences were found regarding risks factors, indicating that different prevention strategies may be required.

Comparative genomic analysis and molecular examination of the diversity of enterotoxigenic Escherichia coli isolates from Chile

Enterotoxigenic Escherichia coli (ETEC) is one of the most common diarrheal pathogens in the low- and middle-income regions of the world, however a systematic examination of the genomic content of isolates from Chile has not yet been undertaken. Whole genome sequencing and comparative analysis of a collection of 125 ETEC isolates from three geographic locations in Chile, allowed the interrogation of phylogenomic groups, sequence types and genes specific to isolates from the different geographic locations. A total of 80.8% (101/125) of the ETEC isolates were identified in E. coli phylogroup A, 15.2% (19/125) in phylogroup B, and 4.0% (5/125) in phylogroup E. The over-representation of genomes in phylogroup A was significantly different from other global ETEC genomic studies. The Chilean ETEC isolates could be further subdivided into sub-clades similar to previously defined global ETEC reference lineages that had conserved multi-locus sequence types and toxin profiles. Comparison of the gene content of the Chilean ETEC identified genes that were unique based on geographic location within Chile, phylogenomic classifications or sequence type. Completion of a limited number of genomes provided insight into the ETEC plasmid content, which is conserved in some phylogenomic groups and not conserved in others. These findings suggest that the Chilean ETEC isolates contain unique virulence factor combinations and genomic content compared to global reference ETEC isolates.

Shigella flexneri Adherence Factor Expression in In Vivo-Like Conditions

Bacterial pathogens have evolved to regulate virulence gene expression at critical points in the colonization and infection processes to successfully cause disease. The Shigella species infect the epithelial cells lining the colon to result in millions of cases of diarrhea and a significant global health burden. As antibiotic resistance rates increase, understanding the mechanisms of infection is vital to ensure successful vaccine development. Despite significant gains in our understanding of Shigella infection, it remains unknown how the bacteria initiate contact with the colonic epithelium. Most pathogens harbor multiple adherence factors to facilitate this process, but Shigella was thought to have lost the ability to produce these factors. Interestingly, we have identified conditions that mimic some features of gastrointestinal transit and that enable Shigella to express adherence structural genes. This work highlights aspects of genetic regulation for Shigella adherence factors and may have a significant impact on future vaccine development.

The Shigella species are Gram-negative, facultative intracellular pathogens that invade the colonic epithelium and cause significant diarrheal disease. Despite extensive research on the pathogen, a comprehensive understanding of how Shigella initiates contact with epithelial cells remains unknown. Shigella maintains many of the same Escherichia coli adherence gene operons; however, at least one critical gene component in each operon is currently annotated as a pseudogene in reference genomes. These annotations, coupled with a lack of structures upon microscopic analysis following growth in laboratory media, have led the field to hypothesize that Shigella is unable to produce fimbriae or other traditional adherence factors. Nevertheless, our previous analyses have demonstrated that a combination of bile salts and glucose induces both biofilm formation and adherence to colonic epithelial cells. The goal of this study was to perform transcriptomic and genetic analyses to demonstrate that adherence gene operons in Shigella flexneri strain 2457T are functional, despite the gene annotations. Our results demonstrate that at least three structural genes facilitate S. flexneri 2457T adherence for epithelial cell contact and biofilm formation. Furthermore, our results demonstrate that host factors, namely, glucose and bile salts at their physiological concentrations in the small intestine, offer key environmental stimuli required for adherence factor expression in S. flexneri. This research may have a significant impact on Shigella vaccine development and further highlights the importance of utilizing in vivo-like conditions to study bacterial pathogenesis.

IMPORTANCE Bacterial pathogens have evolved to regulate virulence gene expression at critical points in the colonization and infection processes to successfully cause disease. The Shigella species infect the epithelial cells lining the colon to result in millions of cases of diarrhea and a significant global health burden. As antibiotic resistance rates increase, understanding the mechanisms of infection is vital to ensure successful vaccine development. Despite significant gains in our understanding of Shigella infection, it remains unknown how the bacteria initiate contact with the colonic epithelium. Most pathogens harbor multiple adherence factors to facilitate this process, but Shigella was thought to have lost the ability to produce these factors. Interestingly, we have identified conditions that mimic some features of gastrointestinal transit and that enable Shigella to express adherence structural genes. This work highlights aspects of genetic regulation for Shigella adherence factors and may have a significant impact on future vaccine development.

Improved Species-Level Clinical Identification of Enterobacteriaceae through Broad-Range dnaJ PCR and Sequencing

Enterobacteriaceae represent a diverse and medically important family of bacteria that are difficult to identify to the species level using the standard molecular method of 16S rRNA gene sequencing. Prior work has demonstrated the value of dnaJ gene sequence analysis in resolving different members of the family. However, existing protocols are not optimized for clinical use and exhibit several limitations in practice. Here, we describe an improved assay for dnaJ-based identification of Enterobacteriaceae which boasts increased broad-range specificity across genera, shorter amplicon sizes that are suitable for use with formalin-fixed or direct patient specimens, and enhanced amplification efficiency and assay sensitivity through the incorporation of locked nucleic acid chemistries. Sequence analysis of public databases indicates that the partial dnaJ sequence interrogated by this design retains high discriminatory power among Enterobacteriaceae genera and species, with only particular lineages of Shigella sp. and Escherichia coli proving unresolvable. Limits of detection studies using 8 disparate species indicated that amplification was consistently achievable across organisms and allowed robust dideoxynucleotide chain terminator sequencing from as little as 10 genome equivalents of template, depending on the species interrogated. Retrospective application of the dnaJ assay to patient specimens enabled unambiguous classification of Enterobacteriaceae to the species level in 22 of 27 (81.5%) positive specimens examined, with most remaining cases representing unresolvable calls between closely related Escherichia coli and Shigella species. We expect that this assay will facilitate the accurate molecular identification of species from the Enterobacteriaceae family in a variety of clinical specimens and diagnostic contexts.

First Detection of Carbapenem-Resistant Escherichia fergusonii Strains Harbouring Beta-Lactamase Genes from Clinical Samples

Recently discovered extraintestinal Escherichia fergusonii obtained from non-clinical samples has exhibited the potential for acquiring multiple beta-lactamase genes, just like many extraintestinal Escherichia coli strains. Albeit, they are often omitted or classified as E. coli. This study aimed to, therefore, identify carbapenem-resistant extended-spectrum beta-lactamase (ESBL) producing E. fergusonii isolates from clinical samples, determine their evolutionary relatedness using 16S rRNA sequencing analysis and screen for beta-lactamase genes. A total of 135 septic wound samples were obtained from patients on referral at a General Hospital in Lagos, Nigeria. For the phenotypic identification of isolates from culture-positive samples, morphological, and physiological tests were carried out. Identities of the isolates harbouring beta-lactamase genes were assigned to their genus strains using the 16S rRNA sequencing. The Kirby Bauer disc diffusion technique and double-disc synergy test were used to screen isolates for multidrug resistance and ESBL production. Carbapenem-resistant ESBL producing isolates were screened for beta-lactamase genes in a polymerase chain reaction. Three E. fergusonii isolates (CR11, CR35 and CR49) were obtained during this study. E. fergusonii strains were motile, non-lactose and non-sorbitol fermenting but positive for cellobiose and adonitol fermentation. The I6S rRNA assigned the phenotypically identified isolates to E. fergusonii species. All three isolates were multidrug-resistant, carbapenem-resistant and ESBL producers. Isolates CR11 and CR35 harboured cefotaximase (CTX-M) and temoniera (TEM) beta-lactamase genes while CR49 harboured sulfhydryl variable (SHV) beta-lactamase gene. We herein report the detection of multiple beta-lactamase genes in carbapenem-resistant ESBL producing E. fergusonii from clinical samples.

A tale of two bacterial enteropathogens and one multivalent vaccine

Shigella and enterotoxigenic Escherichia coli (ETEC) are among the top four enteric pathogens that cause diarrheal illness in young children in developing countries and are major etiologic agents of travellers' diarrhoea. A single vaccine that could target both of these pathogens would have significant public health impact. In this review, we highlight the many pivotal contributions of Phillippe Sansonetti to the identification of molecular mechanisms of pathogenesis of Shigella that paved the way for the development of rationally designed, novel vaccines candidates. The CVD developed a series of live attenuated Shigella vaccine strains based on the most prevalent serotypes associated with disease. Shigella vaccine strains were engineered to express critical ETEC antigens to form a broadly protective Shigella-ETEC multivalent vaccine.

Precise, direct, and rapid detection of Shigella Spa gene by a novel unmodified AuNPs-based optical genosensing system

Early detection of infectious bacteria is a necessity for combating infectious diseases. Due to low infectious dose of Shigella, rapid and sensitive detection is needed. Compared to the presented genes, Spa gene can be introduced as a novel sequence for all species of Shigella detection. Herein, the possibility of Spa genes for detection of four species of Shigella was investigated for the first time by AuNPs-based optical genosensing system. In this method, AuNP-DNA probes were hybridized with Spa gene sequence. When the complementary target is present, it prevents the aggregation of the complex under acid environment and the solution remains red whereas in the absence of the specific sequence, it turns to purple. Therefore, visual detection is possible with bare eye. The comparison of this Optical DNA biosensor and PCR-based method showed that the proposed method is simple, cost-effective, rapid operation, with high or comparable detection limit of (LOD and LOQ: 8.14 and 26.6 ng mLl^-1, respectively), without need of any expensive techniques, and equipments compared to the conventional methods. In conclusion, the described method may develop into a platform that could be utilized for detection of various bacterial species with high accuracy and prompt screening of samples.

Current strategy for local- to global-level molecular epidemiological characterisation of global antimicrobial resistance surveillance system pathogens

The prime goal of molecular epidemiology is to identify the origin and evolution of pathogens, which can potentially influence the public health worldwide. Traditional methods provide limited information which is not sufficient for outbreak investigation and studying transmission dynamics. The recent advancement of next-generation sequencing had a major impact on molecular epidemiological studies. Currently, whole-genome sequencing (WGS) has become the gold standard typing method, especially for clinically significant pathogens. Here, we aimed to describe the application of appropriate molecular typing methods for global antimicrobial resistance surveillance system pathogens based on the level of discrimination and epidemiological settings. This shows that sequence-based methods such as multi-locus sequence typing (MLST) are widely used due to cost-effectiveness and database accessibility. However, WGS is the only method of choice for studying Escherichia coli and Shigella spp. WGS is shown to have higher discrimination than other methods in typing Klebsiella pneumoniae, Acinetobacter baumannii and Salmonella spp. due to its changing accessory genome content. For Gram positives such as Streptococcus pneumoniae, WGS would be preferable to understand the evolution of the strains. Similarly, for Staphylococcus aureus, combination of MLST, staphylococcal protein A or SCCmec typing along with WGS could be the choice for epidemiological typing of hospital- and community-acquired strains. This review highlights that combinations of different typing methods should be used to get complete information since no one standalone method is sufficient to study the varying genome diversity.

In Silico Serotyping Based on Whole-Genome Sequencing Improves the Accuracy of Shigella Identification

Bacteria of the genus Shigella, consisting of 4 species and >50 serotypes, cause shigellosis, a foodborne disease of significant morbidity, mortality, and economic loss worldwide. Classical Shigella identification based on selective media and serology is tedious, time-consuming, expensive, and not always accurate. A molecular diagnostic assay does not distinguish Shigella at the species level or from enteroinvasive Escherichia coli (EIEC). We inspected genomic sequences from 221 Shigella isolates and observed low concordance rates between conventional designation and molecular serotyping: 86.4% and 80.5% at the species and serotype levels, respectively. Serotype determinants for 6 additional serotypes were identified. Examination of differentiation gene markers commonly perceived as characteristic hallmarks in Shigella showed high variability among different serotypes. Using this information, we developed ShigaTyper, an automated workflow that utilizes limited computational resources to accurately and rapidly determine 59 Shigella serotypes using Illumina paired-end whole-genome sequencing (WGS) reads. Shigella serotype determinants and species-specific diagnostic markers were first identified through read alignment to an in-house curated reference sequence database. Relying on sequence hits that passed a threshold level of coverage and accuracy, serotype could be unambiguously predicted within 1 min for an average-size WGS sample of ∼500 MB. Validation with WGS data from 380 isolates showed an accuracy rate of 98.2%. This pipeline is the first step toward building a comprehensive WGS-based analysis pipeline of Shigella spp. in a field laboratory setting, where speed is essential and resources need to be more cost-effectively dedicated.IMPORTANCE Shigella causes diarrheal disease with serious public health implications. However, conventional Shigella identification methods are laborious and time-consuming and can be erroneous due to the high similarity between Shigella and enteroinvasive Escherichia coli (EIEC) and cross-reactivity between serotyping antisera. Further, serotype interpretation is complicated for inexperienced users. To develop an easier method with higher accuracy based on whole-genome sequencing (WGS) for Shigella serotyping, we systematically examined genomic information of Shigella isolates from 53 serotypes to define rules for differentiation and serotyping. We created ShigaTyper, an automated pipeline that accurately and rapidly excludes non-Shigella isolates and identifies 59 Shigella serotypes using Illumina paired-end WGS reads. A serotype can be unambiguously predicted at a data processing speed of 538 MB/min with 98.2% accuracy from a regular laptop. Once it is installed, training in bioinformatics analysis and Shigella genetics is not required. This pipeline is particularly useful to general microbiologists in field laboratories.

Temporal Variability of Escherichia coli Diversity in the Gastrointestinal Tracts of Tanzanian Children with and without Exposure to Antibiotics

This study increases the number of resident Escherichia coli genome sequences, and explores E. coli diversity through longitudinal sampling. We investigate the genomes of E. coli isolated from human gastrointestinal tracts as part of an antibiotic treatment program among rural Tanzanian children. Phylogenomics demonstrates that resident E. coli are diverse, even within a single host. Though the E. coli isolates of the gastrointestinal community tend to be phylogenomically similar at a given time, they differed across the interrogated time points, demonstrating the variability of the members of the E. coli community in these subjects. Exposure to antibiotic treatment did not have an apparent impact on the E. coli community or the presence of resistance and virulence genes within E. coli genomes. The findings of this study highlight the variable nature of specific bacterial members of the human gastrointestinal tract.

The stability of the Escherichia coli populations in the human gastrointestinal tract is not fully appreciated, and represents a significant knowledge gap regarding gastrointestinal community structure, as well as resistance to incoming pathogenic bacterial species and antibiotic treatment. The current study examines the genomic content of 240 Escherichia coli isolates from 30 children, aged 2 to 35 months old, in Tanzania. The E. coli strains were isolated from three time points spanning a six-month time period, with and without antibiotic treatment. The resulting isolates were sequenced, and the genomes compared. The findings in this study highlight the transient nature of E. coli strains in the gastrointestinal tract of these children, as during a six-month interval, no one individual contained phylogenomically related isolates at all three time points. While the majority of the isolates at any one time point were phylogenomically similar, most individuals did not contain phylogenomically similar isolates at more than two time points. Examination of global genome content, canonical E. coli virulence factors, multilocus sequence type, serotype, and antimicrobial resistance genes identified diversity even among phylogenomically similar strains. There was no apparent increase in the antimicrobial resistance gene content after antibiotic treatment. The examination of the E. coli from longitudinal samples from multiple children in Tanzania provides insight into the genomic diversity and population variability of resident E. coli within the rapidly changing environment of the gastrointestinal tract of these children.

IMPORTANCE This study increases the number of resident Escherichia coli genome sequences, and explores E. coli diversity through longitudinal sampling. We investigate the genomes of E. coli isolated from human gastrointestinal tracts as part of an antibiotic treatment program among rural Tanzanian children. Phylogenomics demonstrates that resident E. coli are diverse, even within a single host. Though the E. coli isolates of the gastrointestinal community tend to be phylogenomically similar at a given time, they differed across the interrogated time points, demonstrating the variability of the members of the E. coli community in these subjects. Exposure to antibiotic treatment did not have an apparent impact on the E. coli community or the presence of resistance and virulence genes within E. coli genomes. The findings of this study highlight the variable nature of specific bacterial members of the human gastrointestinal tract.

Evaluation of a Culture-Dependent Algorithm and a Molecular Algorithm for Identification of Shigella spp., Escherichia coli, and Enteroinvasive E. coli

Identification of Shigella spp., Escherichia coli, and enteroinvasive E. coli (EIEC) is challenging because of their close relatedness. Distinction is vital, as infections with Shigella spp. are under surveillance of health authorities, in contrast to EIEC infections. In this study, a culture-dependent identification algorithm and a molecular identification algorithm were evaluated. Discrepancies between the two algorithms and original identification were assessed using whole-genome sequencing (WGS). After discrepancy analysis with the molecular algorithm, 100% of the evaluated isolates were identified in concordance with the original identification. However, the resolution for certain serotypes was lower than that of previously described methods and lower than that of the culture-dependent algorithm. Although the resolution of the culture-dependent algorithm is high, 100% of noninvasive E. coli, Shigella sonnei, and Shigella dysenteriae, 93% of Shigella boydii and EIEC, and 85% of Shigella flexneri isolates were identified in concordance with the original identification. Discrepancy analysis using WGS was able to confirm one of the used algorithms in four discrepant results. However, it failed to clarify three other discrepant results, as it added yet another identification. Both proposed algorithms performed well for the identification of Shigella spp. and EIEC isolates and are applicable in low-resource settings, in contrast to previously described methods that require WGS for daily diagnostics. Evaluation of the algorithms showed that both algorithms are capable of identifying Shigella species and EIEC isolates. The molecular algorithm is more applicable in clinical diagnostics for fast and accurate screening, while the culture-dependent algorithm is more suitable for reference laboratories to identify Shigella spp. and EIEC up to the serotype level.

Closed Genome and Comparative Phylogenetic Analysis of the Clinical Multidrug Resistant Shigella sonnei Strain 866

Shigella sonnei is responsible for the majority of shigellosis infections in the US with over 500,000 cases reported annually. Here, we present the complete genome of the clinical multidrug resistant (MDR) strain 866, which is highly susceptible to bacteriophage infections. The strain has a circular chromosome of 4.85 Mb and carries a 113 kb MDR plasmid. This IncB/O/K/Z-type plasmid, termed p866, confers resistance to five different classes of antibiotics including ß-lactamase, sulfonamide, tetracycline, aminoglycoside, and trimethoprim. Comparative analysis of the plasmid architecture and gene inventory revealed that p866 shares its plasmid backbone with previously described IncB/O/K/Z-type Shigella spp. and Escherichia coli plasmids, but is differentiated by the insertion of antibiotic resistance cassettes, which we found associated with mobile genetic elements such as Tn3, Tn7, and Tn10. A whole genome-derived phylogenetic reconstruction showed the evolutionary relationships of S. sonnei strain 866 and the four established Shigella species, highlighting the clonal nature of S. sonnei.

A Novel Shigella Proteome Microarray Discriminates Targets of Human Antibody Reactivity following Oral Vaccination and Experimental Challenge

Shigella spp. are a major cause of diarrhea and dysentery in children under 5 years old in the developing world. The development of an effective vaccine remains a public health priority, necessitating improved understanding of immune responses to Shigella and identification of protective antigens. We report the development of a core Shigella proteome microarray consisting of 2,133 antigen targets common to all Shigella species. We evaluated the microarray with serum samples from volunteers immunized with either an inactivated whole-cell S. flexneri serotype 2a (Sf2aWC) vaccine or a live attenuated S. flexneri 2a vaccine strain (CVD 1204) or challenged with wild-type S. flexneri 2a (Sf2a challenge). Baseline reactivities to most antigens were detected postintervention in all three groups. Similar immune profiles were observed after CVD 1204 vaccination and Sf2a challenge. Antigens with the largest increases in mean reactivity postintervention were members of the type three secretion system (T3SS), some of which are regarded as promising vaccine targets: these are the invasion plasmid antigens (Ipas) IpaB, IpaC, and IpaD. In addition, new immunogenic targets (IpaA, IpaH, and SepA) were identified. Importantly, immunoreactivities to antigens in the microarray correlated well with antibody titers determined by enzyme-linked immunosorbent assay (ELISA), validating the use of the microarray platform. Finally, our analysis uncovered an immune signature consisting of three conserved proteins (IpaA, IpaB, and IpaC) that was predictive of protection against shigellosis. In conclusion, the Shigella proteome microarray is a robust platform for interrogating serological reactivity to multiple antigens at once and identifying novel targets for the development of broadly protective vaccines.IMPORTANCE Each year, more than 180 million cases of severe diarrhea caused by Shigella occur globally. Those affected (mostly children in poor regions) experience long-term sequelae that severely impair quality of life. Without a licensed vaccine, the burden of disease represents a daunting challenge. An improved understanding of immune responses to Shigella is necessary to support ongoing efforts to identify a safe and effective vaccine. We developed a microarray containing >2,000 proteins common to all Shigella species. Using sera from human adults who received a killed whole-cell or live attenuated vaccine or were experimentally challenged with virulent organisms, we identified new immune-reactive antigens and defined a T3SS protein signature associated with clinical protection.

Responses of the Human Gut Escherichia coli Population to Pathogen and Antibiotic Disturbances

Research on human-associated E. coli tends to focus on pathogens, such as enterotoxigenic E. coli (ETEC) strains, which are a leading cause of diarrhea in developing countries. However, the severity of disease caused by these pathogens is thought to be influenced by the microbiome. The nonpathogenic E. coli community that resides in the human gastrointestinal tract may play a role in pathogen colonization and disease severity and may become a reservoir for virulence and antibiotic resistance genes. Our study used whole-genome sequencing of E. coli before, during, and after challenge with an archetype ETEC isolate, H10407, and antibiotic treatment to explore the diversity and resiliency of the resident E. coli population in response to the ecological disturbances caused by pathogen invasion and antibiotic treatment.

Studies of Escherichia coli in the human gastrointestinal tract have focused on pathogens, such as diarrhea-causing enterotoxigenic E. coli (ETEC), while overlooking the resident, nonpathogenic E. coli community. Relatively few genomes of nonpathogenic E. coli strains are available for comparative genomic analysis, and the ecology of these strains is poorly understood. This study examined the diversity and dynamics of resident human gastrointestinal E. coli communities in the face of the ecological challenges presented by pathogen (ETEC) challenge, as well as of antibiotic treatment. Whole-genome sequences obtained from E. coli isolates from before, during, and after ETEC challenge were used in phylogenomic and comparative genomic analyses to examine the diversity of the resident E. coli communities, as well as the dynamics of the challenge strain, H10407, a well-studied ETEC strain (serotype O78:H11) that produces both heat-labile and heat-stable enterotoxins. ETEC failed to become the dominant E. coli clone in two of the six challenge subjects, each of whom exhibited limited or no clinical presentation of diarrhea. The E. coli communities of the remaining four subjects became ETEC dominant during the challenge but reverted to their original, subject-specific populations following antibiotic treatment, suggesting resiliency of the resident E. coli population following major ecological disruptions. This resiliency is likely due in part to the abundance of antibiotic-resistant ST131 E. coli strains in the resident populations. This report provides valuable insights into the potential interactions of members of the gastrointestinal microbiome and its responses to challenge by an external pathogen and by antibiotic exposure.

IMPORTANCE Research on human-associated E. coli tends to focus on pathogens, such as enterotoxigenic E. coli (ETEC) strains, which are a leading cause of diarrhea in developing countries. However, the severity of disease caused by these pathogens is thought to be influenced by the microbiome. The nonpathogenic E. coli community that resides in the human gastrointestinal tract may play a role in pathogen colonization and disease severity and may become a reservoir for virulence and antibiotic resistance genes. Our study used whole-genome sequencing of E. coli before, during, and after challenge with an archetype ETEC isolate, H10407, and antibiotic treatment to explore the diversity and resiliency of the resident E. coli population in response to the ecological disturbances caused by pathogen invasion and antibiotic treatment.

Novel multiplex PCR assay for identification and subtyping of enteroinvasive Escherichia coli and differentiation from Shigella based on target genes selected by comparative genomics

PURPOSE

Both Shigella and enteroinvasive Escherichia coli (EIEC) can cause enterocolitis, but they have a distinct epidemiology and public health relevance. Current culture-independent testing (CIT) methods to identify Shigella in faecal samples rely on the ipaH gene as the target, which is also found in EIEC genomes. The aim of this study was to design an assay that can identify EIEC in cultures from CIT ipaH-positive samples.

METHODOLOGY

Shigella and EIEC genomes were screened to find unique regions present in EIEC genomes using a comparative genomics approach and differentiating genetic loci that are suitable PCR targets were identified. The primers for these loci were designed and tested in 6501 and 104 genomes of Shigella and EIEC, respectively.

RESULTS

An assay with two sets of multiplex PCR reactions that differentiates Shigella and EIEC based on the presence/absence of at least two out of six loci was developed and evaluated. The majority of Shigella genomes lacked all six loci, while at least two loci were present in most EIEC genomes. This assay successfully differentiated clinical EIEC from Shigella with a limit of detection of 10⁵ cells ml^-1. The sensitivity and specificity were over 95 and 99%, respectively. The assay can further subtype EIEC genomes into their genetic lineages.

CONCLUSION

This new highly specific assay can assist in the identification of EIEC in ipaH PCR-positive samples and augment the public health laboratory surveillance of EIEC and shigellosis.

Out of Asia: the independent rise and global spread of fluoroquinolone-resistant Shigella

Shigella are ranked among the most prevalent aetiologies of diarrhoeal disease worldwide, disproportionately affecting young children in developing countries and high-risk communities in developed settings. Antimicrobial treatment, most commonly with fluoroquinolones, is currently recommended for Shigella infections to alleviate symptoms and control disease transmission. Resistance to fluoroquinolones has emerged in differing Shigella species (S. dysenteriae, flexneri and sonnei) since the turn of the 21st century, originating in endemic areas, and latterly spreading into non-endemic regions. Despite occurring independently, the emergence of fluoroquinolone resistance in these different Shigella species shares striking similarities regarding their epidemiology and resistance mechanisms. Here, we review and discuss the current epidemiology of fluoroquinolone-resistant Shigella species, particularly in the light of recent genomic insights.

Direct Detection of Shigella in Stool Specimens by Use of a Metagenomic Approach

The underestimation of Shigella species as a cause of childhood diarrhea disease has become increasingly apparent with quantitative PCR (qPCR)-based diagnostic methods versus culture. We sought to confirm qPCR-based detection of Shigella via a metagenomics approach. Three groups of samples were selected from diarrheal cases from the Global Enteric Multicenter Study: nine Shigella culture-positive and qPCR-positive (culture⁺ qPCR⁺) samples, nine culture-negative but qPCR-positive (culture⁻ qPCR⁺) samples, and nine culture-negative and qPCR-negative (culture⁻ qPCR⁻) samples. Fecal DNA was sequenced using paired-end Illumina HiSeq, whereby 3.26 × 10⁸ ± 5.6 × 10⁷ high-quality reads were generated for each sample. We used Kraken software to compare the read counts specific to “Shigella” among the three groups. The proportions of Shigella-specific nonhuman sequence reads between culture⁺ qPCR⁺ (0.65 ± 0.42%) and culture⁻ qPCR⁺ (0.55 ± 0.31%) samples were similar (Mann-Whitney U test, P = 0.627) and distinct from the culture⁻ qPCR⁻ group (0.17 ± 0.15%, P < 0.05). The read counts of sequences previously targeted by Shigella/enteroinvasive Escherichia coli (EIEC) qPCR assays, namely, ipaH, virA, virG, ial, ShET2, and ipaH3, were also similar between the culture⁺ qPCR⁺ and culture⁻ qPCR⁺ groups and distinct from the culture⁻ qPCR⁻ groups (P < 0.001). Kraken performed well versus other methods: its precision and recall of Shigella were excellent at the genus level but variable at the species level. In summary, metagenomic sequencing indicates that Shigella/EIEC qPCR-positive samples are similar to those of Shigella culture-positive samples in Shigella sequence composition, thus supporting qPCR as an accurate method for detecting Shigella.

The Intriguing Evolutionary Journey of Enteroinvasive E. coli (EIEC) toward Pathogenicity

Among the intestinal pathogenic Escherichia coli, enteroinvasive E. coli (EIEC) are a group of intracellular pathogens able to enter epithelial cells of colon, multiplicate within them, and move between adjacent cells with a mechanism similar to Shigella, the ethiological agent of bacillary dysentery. Despite EIEC belong to the same pathotype of Shigella, they neither have the full set of traits that define Shigella nor have undergone the extensive gene decay observed in Shigella. Molecular analysis confirms that EIEC are widely distributed among E. coli phylogenetic groups and correspond to bioserotypes found in many E. coli serogroups. Like Shigella, also in EIEC the critical event toward a pathogenic life-style consisted in the acquisition by horizontal gene transfer of a large F-type plasmid (pINV) containing the genes required for invasion, intracellular survival, and spreading through the intestinal mucosa. In Shigella, the ample gain in virulence determinants has been counteracted by a substantial loss of functions that, although important for the survival in the environment, are redundant or deleterious for the life inside the host. The pathoadaptation process that has led Shigella to modify its metabolic profile and increase its pathogenic potential is still in infancy in EIEC, although maintenance of some features typical of E. coli might favor their emerging relevance as intestinal pathogens worldwide, as documented by recent outbreaks in industrialized countries. In this review, we will discuss the evolution of EIEC toward Shigella-like invasive forms going through the epidemiology, including the emergence of new virulent strains, their genome organization, and the complex interactions they establish with the host.

Comparative genomics and transcriptomics of Escherichia coli isolates carrying virulence factors of both enteropathogenic and enterotoxigenic E. coli

Escherichia coli that are capable of causing human disease are often classified into pathogenic variants (pathovars) based on their virulence gene content. However, disease-associated hybrid E. coli, containing unique combinations of multiple canonical virulence factors have also been described. Such was the case of the E. coli O104:H4 outbreak in 2011, which caused significant morbidity and mortality. Among the pathovars of diarrheagenic E. coli that cause significant human disease are the enteropathogenic E. coli (EPEC) and enterotoxigenic E. coli (ETEC). In the current study we use comparative genomics, transcriptomics, and functional studies to characterize isolates that contain virulence factors of both EPEC and ETEC. Based on phylogenomic analysis, these hybrid isolates are more genomically-related to EPEC, but appear to have acquired ETEC virulence genes. Global transcriptional analysis using RNA sequencing, demonstrated that the EPEC and ETEC virulence genes of these hybrid isolates were differentially-expressed under virulence-inducing laboratory conditions, similar to reference isolates. Immunoblot assays further verified that the virulence gene products were produced and that the T3SS effector EspB of EPEC, and heat-labile toxin of ETEC were secreted. These findings document the existence and virulence potential of an E. coli pathovar hybrid that blurs the distinction between E. coli pathovars.

Insights into enterotoxigenic Escherichia coli diversity in Bangladesh utilizing genomic epidemiology

Enterotoxigenic Escherichia coli (ETEC) cause more than 500,000 deaths each year in the developing world and are characterized on a molecular level by the presence of genes that encode the heat-stable (ST) and/or heat-labile (LT) enterotoxins, as well as surface structures, known as colonization factors (CFs). Genome sequencing and comparative genomic analyses of 94 previously uncharacterized ETEC isolates demonstrated remarkable genomic diversity, with 28 distinct sequence types identified in three phylogenomic groups. Interestingly, there is a correlation between the genomic sequence type and virulence factor profiles based on prevalence of the isolate, suggesting that there is an optimal combination of genetic factors required for survival, virulence and transmission in the most successful clones. A large-scale BLAST score ratio (LS-BSR) analysis was further applied to identify ETEC-specific genomic regions when compared to non-ETEC genomes, as well as genes that are more associated with clinical presentations or other genotypic markers. Of the strains examined, 21 of 94 ETEC isolates lacked any previously identified CF. Homology searches with the structural subunits of known CFs identified 6 new putative CF variants. These studies provide a roadmap to exploit genomic analyses by directing investigations of pathogenesis, virulence regulation and vaccine development.

Identification of Escherichia coli and Shigella Species from Whole-Genome Sequences

Escherichia coli and Shigella species are closely related and genetically constitute the same species. Differentiating between these two pathogens and accurately identifying the four species of Shigella are therefore challenging. The organism-specific bioinformatics whole-genome sequencing (WGS) typing pipelines at Public Health England are dependent on the initial identification of the bacterial species by use of a kmer-based approach. Of the 1,982 Escherichia coli and Shigella sp. isolates analyzed in this study, 1,957 (98.4%) had concordant results by both traditional biochemistry and serology (TB&S) and the kmer identification (ID) derived from the WGS data. Of the 25 mismatches identified, 10 were enteroinvasive E. coli isolates that were misidentified as Shigella flexneri or S. boydii by the kmer ID, and 8 were S. flexneri isolates misidentified by TB&S as S. boydii due to nonfunctional S. flexneri O antigen biosynthesis genes. Analysis of the population structure based on multilocus sequence typing (MLST) data derived from the WGS data showed that the remaining discrepant results belonged to clonal complex 288 (CC288), comprising both S. boydii and S. dysenteriae strains. Mismatches between the TB&S and kmer ID results were explained by the close phylogenetic relationship between the two species and were resolved with reference to the MLST data. Shigella can be differentiated from E. coli and accurately identified to the species level by use of kmer comparisons and MLST. Analysis of the WGS data provided explanations for the discordant results between TB&S and WGS data, revealed the true phylogenetic relationships between different species of Shigella, and identified emerging pathoadapted lineages.

The Effects of Signal Erosion and Core Genome Reduction on the Identification of Diagnostic Markers

Whole-genome sequence (WGS) data are commonly used to design diagnostic targets for the identification of bacterial pathogens. To do this effectively, genomics databases must be comprehensive to identify the strict core genome that is specific to the target pathogen. As additional genomes are analyzed, the core genome size is reduced and there is erosion of the target-specific regions due to commonality with related species, potentially resulting in the identification of false positives and/or false negatives.

A comparative analysis of 1,130 Burkholderia genomes identified unique markers for many named species, including the human pathogens B. pseudomallei and B. mallei. Due to core genome reduction and signature erosion, only 38 targets specific to B. pseudomallei/mallei were identified. By using only public genomes, a larger number of markers were identified, due to undersampling, and this larger number represents the potential for false positives. This analysis has implications for the design of diagnostics for other species where the genomic space of the target and/or closely related species is not well defined.

NASP: an accurate, rapid method for the identification of SNPs in WGS datasets that supports flexible input and output formats

Whole-genome sequencing (WGS) of bacterial isolates has become standard practice in many laboratories. Applications for WGS analysis include phylogeography and molecular epidemiology, using single nucleotide polymorphisms (SNPs) as the unit of evolution. NASP was developed as a reproducible method that scales well with the hundreds to thousands of WGS data typically used in comparative genomics applications. In this study, we demonstrate how NASP compares with other tools in the analysis of two real bacterial genomics datasets and one simulated dataset. Our results demonstrate that NASP produces similar, and often better, results in comparison with other pipelines, but is much more flexible in terms of data input types, job management systems, diversity of supported tools and output formats. We also demonstrate differences in results based on the choice of the reference genome and choice of inferring phylogenies from concatenated SNPs or alignments including monomorphic positions. NASP represents a source-available, version-controlled, unit-tested method and can be obtained from tgennorth.github.io/NASP.

Molecular Differentiation of Shigella Spp. from Enteroinvasive E. Coli

A real-time polymerase chain reaction (PCR) assay, amplifying the genes encoding lactose permease (lacY) and invasion plasmid antigen H (ipaH), was run on 121 isolates phenotypically classified as Shigella spp., enteroinvasive Escherichia coli (EIEC), or EIEC O nontypable (ONT). The results were compared with data from a generic E. coli multiple-locus variable-number of tandem repeat analysis (MLVA) and a Shigella MLVA.

The real-time PCR verified all Shigella spp. (n = 53) as Shigella (lacY negative) and all EIEC O121 (n = 15) and EIEC O124 (n = 2) as EIEC (lacY positive). However, the real-time PCR typed EIEC O164 as either EIEC (n = 2) or Shigella (n = 2) and, thus, was not suited for classifying this group of isolates. Interestingly, the majority (42/47, 89.4%) of the EIEC ONT were classified as Shigella (lacY negative) by the real-time PCR, and in nearly all cases, (92.9%, 39/42) data from both MLVA assays supported these findings. Overall, in 94.7% (114/121) of the isolates, the results from the real-time PCR were substantiated by the results from the MLVA assays.

In conclusion, the real-time PCR assay was fast and accurate in differentiating Shigella spp. from EIEC, with the exception of the EIEC O164 group. This molecular assay was particularly pragmatic for the challenging EIEC ONT group.

Investigating the Relatedness of Enteroinvasive Escherichia coli to Other E. coli and Shigella Isolates by Using Comparative Genomics

Enteroinvasive Escherichia coli (EIEC) is a unique pathovar that has a pathogenic mechanism nearly indistinguishable from that of Shigella species. In contrast to isolates of the four Shigella species, which are widespread and can be frequent causes of human illness, EIEC causes far fewer reported illnesses each year. In this study, we analyzed the genome sequences of 20 EIEC isolates, including 14 first described in this study. Phylogenomic analysis of the EIEC genomes demonstrated that 17 of the isolates are present in three distinct lineages that contained only EIEC genomes, compared to reference genomes from each of the E. coli pathovars and Shigella species. Comparative genomic analysis identified genes that were unique to each of the three identified EIEC lineages. While many of the EIEC lineage-specific genes have unknown functions, those with predicted functions included a colicin and putative proteins involved in transcriptional regulation or carbohydrate metabolism. In silico detection of the Shigella virulence plasmid (pINV), which is essential for the invasion of host cells, demonstrated that a form of pINV was present in nearly all EIEC genomes, but the Mxi-Spa-Ipa region of the plasmid that encodes the invasion-associated proteins was absent from several of the EIEC isolates. The comparative genomic findings in this study support the hypothesis that multiple EIEC lineages have evolved independently from multiple distinct lineages of E. coli via the acquisition of the Shigella virulence plasmid and, in some cases, the Shigella pathogenicity islands.