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

Shigella and childhood stunting: Evidence, gaps, and future research directions

Early childhood growth deficits have been shown to have lifelong health and economic impacts, yet their connection to one of their underlying causes, diarrheal diseases, has remained difficult to characterize. Identifying the processes and mechanisms that underlie this link has remained a challenge due to the complexity of the relationship and limitations in access to more advanced laboratory methods. In recent years, however, several large-scale, multisite studies have extensively investigated and reported the prevalence, etiology, and impacts of diarrheal diseases in children under 5 years (CU5) in low- to middle-income countries (LMICs). These studies, in combination with several single-site studies, have applied more advanced laboratory methods to uncover the etiology, true prevalence, infection mechanisms, and inflammation biomarkers of diarrheal disease. Of the multiple pathogens that have been shown to be strongly associated with diarrheal disease in CU5, Shigella is one of the more prevalent and impactful of these pathogens. In this narrative review, we highlight key insights from these studies and identify knowledge gaps and directions for future research. According to these studies, Shigella is most commonly detected in toddlers and young children; however, it can cause more severe disease and has a greater impact on linear growth for infants. Shigella often has a stronger relationship to linear growth faltering (LGF) than other enteropathogens, with higher Shigella loads resulting in greater growth deficits. Future studies should employ more Shigella-specific molecular assays and identify diarrheal etiologies using standardized diagnostics to improve child anthropometric and Shigella surveillance. Also, they should focus on uncovering the mechanisms of the relationship underlying Shigella and growth faltering to better characterize the role of asymptomatic infections and intestinal inflammation in this relationship.

Cloning, improved expression and purification of invasion plasmid antigen D (IpaD): an effector protein of enteroinvasive Escherichia coli (EIEC)

The widespread increase in broad-spectrum antimicrobial resistance is making it more difficult to treat gastrointestinal infections. Enteroinvasive Escherichia coli is a prominent etiological agent of bacillary dysentery, invading via the fecal-oral route and exerting virulence on the host via the type III secretion system. IpaD, a surface-exposed protein on the T3SS tip that is conserved among EIEC and Shigella, may serve as a broad immunogen for bacillary dysentery protection. For the first time, we present an effective framework for improving the expression level and yield of IpaD in the soluble fraction for easy recovery, as well as ideal storage conditions, which may aid in the development of new protein therapies for gastrointestinal infections in the future. To achieve this, uncharacterized full length IpaD gene from EIEC was cloned into pHis-TEV vector and induction parameters were optimized for enhanced expression in the soluble fraction. After affinity-chromatography based purification, 61% pure protein with a yield of 0.33 mg per litre of culture was obtained. The purified IpaD was retained its secondary structure with a prominent α-helical structure as well as functional activity during storage, at 4°C, -20°C and -80°C using 5% sucrose as cryoprotectants, which is a critical criterion for protein-based treatments.

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.

A 16th century Escherichia coli draft genome associated with an opportunistic bile infection

Escherichia coli – one of the most characterized bacteria and a major public health concern – remains invisible across the temporal landscape. Here, we present the meticulous reconstruction of the first ancient E. coli genome from a 16^th century gallstone from an Italian mummy with chronic cholecystitis. We isolated ancient DNA and reconstructed the ancient E. coli genome. It consisted of one chromosome of 4446 genes and two putative plasmids with 52 genes. The E. coli strain belonged to the phylogroup A and an exceptionally rare sequence type 4995. The type VI secretion system component genes appears to be horizontally acquired from Klebsiella aerogenes, however we could not identify any pathovar specific genes nor any acquired antibiotic resistances. A sepsis mouse assay showed that a closely related contemporary E. coli strain was avirulent. Our reconstruction of this ancient E. coli helps paint a more complete picture of the burden of opportunistic infections of the past.

Ancient DNA from an Italian mummy’s gallstone provides insight into opportunistic E. coli infection.

Impact of biannual mass azithromycin treatment on enteropathogen carriage in children younger than 5 years in Niger

We analyzed samples obtained at baseline and 24 months in a mass azithromycin administration trial in Niger using quantitative polymerase chain reaction. In villages randomized to azithromycin, Shigella was the only pathogen reduced at 24 months (prevalence ratio, 0.36 [95% confidence interval: .17–.79]; difference in log quantity, −.42 [−.75 to −.10]).

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.

Shigellosis in young children in low-income and middle-income countries: insights from molecular diagnostics

PURPOSE OF REVIEW

To describe the impact of molecular diagnostics on our understanding of the burden and epidemiology of shigellosis in children in low-income and middle-income countries.

RECENT FINDINGS

The incorporation of molecular diagnostics has led to a substantial increase in estimates of the burden of shigellosis and have allowed for further resolution of other aspects of Shigella epidemiology, including the clinical characteristics of shigellosis, the association between clinical and subclinical Shigella infection and linear growth shortfalls, protection after natural infection, duration of convalescent shedding, and host determinants of susceptibility.

SUMMARY

The increased sensitivity and precision afforded by molecular approaches has represented a major advance in our understanding of the epidemiology and burden of shigellosis in the settings of highest importance.

The E phylogroup of Escherichia coli is highly diverse and mimics the whole E. coli species population structure

To get a global picture of the population structure of the Escherichia coli phylogroup E, encompassing the O157:H7 EHEC lineage, we analysed the whole genome of 144 strains isolated from various continents, hosts and lifestyles and representative of the phylogroup diversity. The strains possess 4331 to 5440 genes with a core genome of 2771 genes and a pangenome of 33 722 genes. The distribution of these genes among the strains shows an asymmetric U-shaped distribution. E phylogenetic strains have the largest genomes of the species, partly explained by the presence of mobile genetic elements. Sixty-eight lineages were delineated, some of them exhibiting extra-intestinal virulence genes and being virulent in the mouse sepsis model. Except for the EHEC lineages and the reference EPEC, EIEC and ETEC strains, very few strains possess intestinal virulence genes. Most of the strains were devoid of acquired resistance genes, but eight strains possessed extended-spectrum beta-lactamase genes. Human strains belong to specific lineages, some of them being virulent and antibiotic-resistant [sequence type complexes (STcs) 350 and 2064]. The E phylogroup mimics all the features of the species as a whole, a phenomenon already observed at the STc level, arguing for a fractal population structure of E. coli.

Draft Genome Sequences of Two Enteroinvasive Escherichia coli Strains Representative of Major Enteroinvasive E. coli Clades

There are six described pathotypes of Escherichia coli that cause significant clinical illness in humans. Enteroinvasive E. coli (EIEC) strains have been shown to be separated into three phylogenomic clades. To add to a limited body of EIEC genomic data, we report two high-quality draft genome sequences representing different EIEC phylogenomic clades.

Molecular determinants of surface colonisation in diarrhoeagenic Escherichia coli (DEC): from bacterial adhesion to biofilm formation

Escherichia coli is primarily known as a commensal colonising the gastrointestinal tract of infants very early in life but some strains being responsible for diarrhoea, which can be especially severe in young children. Intestinal pathogenic E. coli include six pathotypes of diarrhoeagenic E. coli (DEC), namely, the (i) enterotoxigenic E. coli, (ii) enteroaggregative E. coli, (iii) enteropathogenic E. coli, (iv) enterohemorragic E. coli, (v) enteroinvasive E. coli and (vi) diffusely adherent E. coli. Prior to human infection, DEC can be found in natural environments, animal reservoirs, food processing environments and contaminated food matrices. From an ecophysiological point of view, DEC thus deal with very different biotopes and biocoenoses all along the food chain. In this context, this review focuses on the wide range of surface molecular determinants acting as surface colonisation factors (SCFs) in DEC. In the first instance, SCFs can be broadly discriminated into (i) extracellular polysaccharides, (ii) extracellular DNA and (iii) surface proteins. Surface proteins constitute the most diverse group of SCFs broadly discriminated into (i) monomeric SCFs, such as autotransporter (AT) adhesins, inverted ATs, heat-resistant agglutinins or some moonlighting proteins, (ii) oligomeric SCFs, namely, the trimeric ATs and (iii) supramolecular SCFs, including flagella and numerous pili, e.g. the injectisome, type 4 pili, curli chaperone-usher pili or conjugative pili. This review also details the gene regulatory network of these numerous SCFs at the various stages as it occurs from pre-transcriptional to post-translocational levels, which remains to be fully elucidated in many cases.

A One Health Perspective for Defining and Deciphering Escherichia coli Pathogenic Potential in Multiple Hosts

E. coli is one of the most common species of bacteria colonizing humans and animals. The singularity of E. coli 's genus and species underestimates its multifaceted nature, which is represented by different strains, each with different combinations of distinct virulence factors. In fact, several E. coli pathotypes, or hybrid strains, may be associated with both subclinical infection and a range of clinical conditions, including enteric, urinary, and systemic infections. E. coli may also express DNA-damaging toxins that could impact cancer development. This review summarizes the different E. coli pathotypes in the context of their history, hosts, clinical signs, epidemiology, and control. The pathotypic characterization of E. coli in the context of disease in different animals, including humans, provides comparative and One Health perspectives that will guide future clinical and research investigations of E. coli infections.

Virulence-Related Gene Distribution Among Shigella Isolates in Anhui, China: The Association with Antimicrobial Resistance

Objective

The aim of this study was to investigate the antimicrobial resistance profiles and distribution of virulence-related genes (VRGs) among Shigella isolates in Anhui, China, and to identify the correlation between the VRGs and antimicrobial resistance.

Materials and Methods

A total of 525 non-duplicate Shigella isolates (449 S. flexneri, 68 S. sonnei, 3 S. boydii, and 5 S. dysenteriae) were collected in Anhui Province, China between September 2011 and September 2015. The antimicrobial resistance of the strains was determined by the agar dilution method according to CLSI guidelines. The presence of 16 VRGs, including ipaH, ipaA-D, ial, virB, virF, set, sen, icsA, icsB, sigA, sat, pic, and sepA, was evaluated using PCR amplification and sequencing.

Results

Shigella flexneri was the most abundant (85.5%), followed by S. sonnei (13.0%). The proportion of males with S. flexneri was higher than that of females (57% vs 43%; P<0.0001). The most common resistance pattern was the combination of ampicillin, nalidixic acid, and tetracycline for S. flexneri (90.2%) and S. sonnei (94.1%). Resistance to ciprofloxacin and levofloxacin was more common among S. flexneri than among S. sonnei (49.7% vs.19.1%, P<0.0001; 30.5% vs 10.3%, P=0.001, respectively). All the isolates were positive for the ipaH gene, while the set, sat, pic, and sepA genes were not detected among the S. sonnei isolates. Except for sigA and sen, resistance to chloramphenicol and ciprofloxacin was more common among VRG-positive S. flexneri than among VRG-negative S. flexneri (P<0.05). Furthermore, resistance to ceftriaxone and ceftazidime was more frequently detected among sat- and set-positive S. flexneri than among sat- and set-negative S. flexneri (P<0.05). However, gentamicin resistance was more prevalent among VRG-negative (ial, virF, set, sat, pic, and sepA) S. flexneri than among VRG-positive S. flexneri (P<0.05).

Conclusion

Shigella flexneri remains the predominant species in Anhui, China, and the resistance to fluoroquinolones was more widespread among S. flexneri than among S. sonnei. Shigella flexneri strains harboring specific VRGs were associated with antimicrobial resistance. To the best of our knowledge, this is the first report of the correlation between the VRGs and antimicrobial resistance in Anhui, China.

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.

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.

Phylogenetic background and habitat drive the genetic diversification of Escherichia coli

Escherichia coli is mostly a commensal of birds and mammals, including humans, where it can act as an opportunistic pathogen. It is also found in water and sediments. We investigated the phylogeny, genetic diversification, and habitat-association of 1,294 isolates representative of the phylogenetic diversity of more than 5,000 isolates from the Australian continent. Since many previous studies focused on clinical isolates, we investigated mostly other isolates originating from humans, poultry, wild animals and water. These strains represent the species genetic diversity and reveal widespread associations between phylogroups and isolation sources. The analysis of strains from the same sequence types revealed very rapid change of gene repertoires in the very early stages of divergence, driven by the acquisition of many different types of mobile genetic elements. These elements also lead to rapid variations in genome size, even if few of their genes rise to high frequency in the species. Variations in genome size are associated with phylogroup and isolation sources, but the latter determine the number of MGEs, a marker of recent transfer, suggesting that gene flow reinforces the association of certain genetic backgrounds with specific habitats. After a while, the divergence of gene repertoires becomes linear with phylogenetic distance, presumably reflecting the continuous turnover of mobile element and the occasional acquisition of adaptive genes. Surprisingly, the phylogroups with smallest genomes have the highest rates of gene repertoire diversification and fewer but more diverse mobile genetic elements. This suggests that smaller genomes are associated with higher, not lower, turnover of genetic information. Many of these genomes are from freshwater isolates and have peculiar traits, including a specific capsule, suggesting adaptation to this environment. Altogether, these data contribute to explain why epidemiological clones tend to emerge from specific phylogenetic groups in the presence of pervasive horizontal gene transfer across the species.

Previous large scale studies on the evolution of E. coli focused on clinical isolates emphasizing virulence and antibiotic resistance in medically important lineages. Yet, most E. coli strains are either human commensals or not associated with humans at all. Here, we analyzed a large collection of non-clinical isolates of the species to assess the mechanisms of gene repertoire diversification in the light of isolation sources and phylogeny. We show that gene repertoires evolve so rapidly by the high turnover of mobile genetic elements that epidemiologically indistinguishable strains can be phenotypically extremely heterogeneous, illustrating the velocity of bacterial adaptation and the importance of accounting for the information on the whole genome at the epidemiological scale. Phylogeny and habitat shape the genetic diversification of E. coli to similar extents. Surprisingly, freshwater strains seem specifically adapted to this environment, breaking the paradigm that E. coli environmental isolates are systematically fecal contaminations. As a consequence, the evolution of this species is also shaped by environmental habitats, and it may diversify by acquiring genes and mobile elements from environmental bacteria (and not just from gut bacteria). This may facilitate the acquisition of virulence factors and antibiotic resistance in the strains that become pathogenic.

Challenges in identification of enteroinvasive Escherichia coli and Shigella spp. in Lebanon

This study aimed to evaluate the routine identification tools available in Lebanon for differentiation of Escherichia coli and Shigella spp. The identification of 43 isolates defined as Shigella spp. by Api 20E was accessed using MALDI-TOF, serological testing, duplex PCR targeting ipaH (present in Shigella spp. and enteroinvasive E. coli "EIEC") and lacY (found in E. coli including EIEC but not Shigella spp.) as well as gyrB gene sequencing. Antibiotic susceptibility was investigated as well as Shiga-toxin production. All isolates were identified as E. coli by MALDI-TOF while the PCR showed a disparate group of 26 EIEC, 11 Shigella spp., 5 E. coli and 1 inactive E. coli. However, the sequencing of gyrB gene, which was recently described as a suitable marker for distinguishing E. coli and Shigella spp., identified all isolates as E. coli. Antibiotic resistance was noticeable against ß-lactams, rifampicin, trimethoprim-sulfamethoxazole, gentamicin, and ciprofloxacin. The most common variants of beta-lactamase genes were blaTEM-1, blaCTX-M-15, and blaCTX-M-3. A great discordance between the used methods in identification was revealed herein. An accurate identification technique able to distinguish E. coli from Shigella spp. in routine laboratories is a pressing need in order to select the appropriate treatment and assess the epidemiology of these bacteria.

Laboratory Cultivation and Storage of Shigella

Shigella species, which are closely related to Escherichia coli, can easily be maintained and stored in the laboratory. This article includes protocols for preparation of routine growth conditions and media, for storage of the bacteria, and for monitoring of the presence of the virulence plasmid. © 2019 by John Wiley & Sons, Inc. Basic Protocol 1: Growth of S. flexneri from frozen stocks or agar stabs Basic Protocol 2: Growth of S. flexneri in rich liquid medium Alternate Protocol 1: Growth of S. flexneri in rich defined medium Alternate Protocol 2: Growth of S. flexneri in minimal medium Basic Protocol 3: Storage of S. flexneri in frozen stocks Alternate Protocol 3: Storage of S. flexneri in agar stabs.

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.

The Complete Genome of the Atypical Enteropathogenic Escherichia coli Archetype Isolate E110019 Highlights a Role for Plasmids in Dissemination of the Type III Secreted Effector EspT

Enteropathogenic Escherichia coli (EPEC) is a leading cause of moderate to severe diarrhea among young children in developing countries, and EPEC isolates can be subdivided into two groups. Typical EPEC (tEPEC) bacteria are characterized by the presence of both the locus of enterocyte effacement (LEE) and the plasmid-encoded bundle-forming pilus (BFP), which are involved in adherence and translocation of type III effectors into the host cells. Atypical EPEC (aEPEC) bacteria also contain the LEE but lack the BFP. In the current report, we describe the complete genome of outbreak-associated aEPEC isolate E110019, which carries four plasmids. Comparative genomic analysis demonstrated that the type III secreted effector EspT gene, an autotransporter gene, a hemolysin gene, and putative fimbrial genes are all carried on plasmids. Further investigation of 65 espT-containing E. coli genomes demonstrated that different espT alleles are associated with multiple plasmids that differ in their overall gene content from the E110019 espT-containing plasmid. EspT has been previously described with respect to its role in the ability of E110019 to invade host cells. While other type III secreted effectors of E. coli have been identified on insertion elements and prophages of the chromosome, we demonstrated in the current study that the espT gene is located on multiple unique plasmids. These findings highlight a role of plasmids in dissemination of a unique E. coli type III secreted effector that is involved in host invasion and severe diarrheal illness.

Existence of virulence genes in clinical Shigella sonnei isolates from Jiangsu Province of China: a multicenter study

Background

The ability of Shigella to invade, colonizes, and eventually kill host cells is influenced by many virulence factors. The aims of this study were to assess the presence of 11 virulence genes of S. sonnei strains isolated in this country.

Methods

A total of 166 S. sonnei was collected from 13 cities of Jiangsu province through the provincial Centers for Disease Control (CDC) from 2010 to 2015 and then the distribution of virulence genes was detected by polymerase chain reaction (PCR) technology.

Results

Invasive virulence genes included ipaH and ial, in which the positive rate of ipaH was 100% while the positive rate of ial was 15.1% in S. sonnei. The classic pathway of regulating expression of Shigella virulence gene involved virF and virB gene, which positive rates were 33.7% and 24.1% respectively. The most common serine protease autotransporters of Enterobacteriaceae among S. sonnei were sigA (100%), followed by sepA (3.0%), sat (3.0%), pic (1.2%). Shigella enterotoxin genes include sen, set1A, set1B were found in 16.3%, 6.0% and 1.8% of the isolates, respectively.

Conclusions

This study provides baseline information on the distribution of virulence genes in clinical S. sonnei trains in Jiangsu province in China, which will be important for implementation of effective control strategies.

The Autotransporter IcsA Promotes Shigella flexneri Biofilm Formation in the Presence of Bile Salts

Shigella flexneri is an intracellular bacterial pathogen that invades epithelial cells in the colonic mucosa, leading to bloody diarrhea. A previous study showed that S. flexneri forms biofilms in the presence of bile salts, through an unknown mechanism. Here, we investigated the potential role of adhesin-like autotransporter proteins in S. flexneri biofilm formation. BLAST search analysis revealed that the S. flexneri 2457T genome harbors 4 genes, S1242, S1289, S2406, and icsA, encoding adhesin-like autotransporter proteins. Deletion mutants of the S1242, S1289, S2406 and icsA genes were generated and tested for biofilm formation. Phenotypic analysis of the mutant strains revealed that disruption of icsA abolished bile salt-induced biofilm formation. IcsA is an outer membrane protein secreted at the bacterial pole that is required for S. flexneri actin-based motility during intracellular infection. In extracellular biofilms, IcsA was also secreted at the bacterial pole and mediated bacterial cell-cell contacts and aggregative growth in the presence of bile salts. Dissecting individual roles of bile salts showed that deoxycholate is a robust biofilm inducer compared to cholate. The release of the extracellular domain of IcsA through IcsP-mediated cleavage was greater in the presence of cholate, suggesting that the robustness of biofilm formation was inversely correlated with IcsA processing. Accordingly, deletion of icsP abrogated IcsA processing in biofilms and enhanced biofilm formation.

Shifting national surveillance of Shigella infections toward geno-serotyping by the development of a tailored Luminex assay and NGS workflow

The phylogenetically closely related Shigella species and enteroinvasive Escherichia coli (EIEC) are responsible for millions of episodes of bacterial dysenteriae worldwide. Given its distinct epidemiology and public health relevance, only Shigellae are subject to mandatory reporting and follow‐up by public health authorities. However, many clinical laboratories struggle to differentiate non‐EIEC, EIEC, and Shigella in their current workflows, leading to inaccuracies in surveillance and rising numbers of misidentified E. coli samples at the National Reference Centre (NRC). In this paper, we describe two novel tools to enhance Shigella surveillance. First, we developed a low‐cost Luminex‐based multiplex assay combining five genetic markers for species identification with 11 markers for serotype prediction for S. sonnei and S. flexneri isolates. Using a test panel of 254 clinical samples, this assay has a sensitivity of 100% in differentiation of EIEC/Shigella pathotype from non‐EIEC strains, and 68.7% success rate in distinction of Shigella and EIEC. A novel, and particularly successful marker was a Shigella‐specific deletion in the spermidine acetyltransferase gene speG, reflecting its metabolic decay. For Shigella serotype prediction, the multiplex assay scored a sensitivity and specificity of 96.6% and 98.4%, respectively. All discrepancies were analyzed with whole‐genome sequencing and shown to be related to causative mutations (stop codons, indels, and promoter mutations) in glycosyltransferase genes. This observation spurred the development of an in silico workflow which extracts the Shigella serotype from Next‐Generation Sequencing (NGS) data, taking into account gene functionality. Both tools will be implemented in the workflow of the NRC, and will play a major role in the shift from phenotypic to genotyping‐based surveillance of shigellosis in Belgium.

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.

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.

Phylogenetic comparison of enteroinvasive Escherichia coli isolated from cases of diarrhoeal disease in England, 2005-2016

We compared the genomes of 60 isolates of enteroinvasive Escherichia coli (EIEC) in order to better understand the step-wise evolutionary process from non-pathogenic to host-adapted pathogenic E. coli. All isolates belonged to either sequence type (ST) 6, ST99 or ST270. Each ST was located on different branches of the E. coli phylogeny and had invasion plasmids (pINVs) belonging to FII-21 (ST99, ST270), FII-27 (ST270) or FII-28 (ST6, ST270) incompatibility groups. A higher number of insertion sequence (IS) elements were identified in ST6 and ST270 than in ST99, and appeared to be driving the loss of functional genes. Comparison of the pINV from each ST revealed different degrees of gene loss, with pINV from ST270 being most similar to that found in Shigella species. We captured three EIEC STs at different stages of patho-adaptation, with ST270 being the most 'shigella-like' and the most divergent from non-pathogenic E. coli, and ST99 being the least divergent.

Cooperative Immune Suppression by Escherichia coli and Shigella Effector Proteins

The enteric attaching and effacing (A/E) pathogens enterohemorrhagic Escherichia coli (EHEC) and enteropathogenic E. coli (EPEC) and the invasive pathogens enteroinvasive E. coli (EIEC) and Shigella encode type III secretion systems (T3SS) used to inject effector proteins into human host cells during infection. Among these are a group of effectors required for NF-κB-mediated host immune evasion. Recent studies have identified several effector proteins from A/E pathogens and EIEC/Shigella that are involved in suppression of NF-κB and have uncovered their cellular and molecular functions. A novel mechanism among these effectors from both groups of pathogens is to coordinate effector function during infection. This cooperativity among effector proteins explains how bacterial pathogens are able to effectively suppress innate immune defense mechanisms in response to diverse classes of immune receptor signaling complexes (RSCs) stimulated during infection.

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.

Characterization of a Large Antibiotic Resistance Plasmid Found in Enteropathogenic Escherichia coli Strain B171 and Its Relatedness to Plasmids of Diverse E. coli and Shigella Strains

Enteropathogenic Escherichia coli (EPEC) is a leading cause of severe infantile diarrhea in developing countries. Previous research has focused on the diversity of the EPEC virulence plasmid, whereas less is known regarding the genetic content and distribution of antibiotic resistance plasmids carried by EPEC. A previous study demonstrated that in addition to the virulence plasmid, reference EPEC strain B171 harbors a second, larger plasmid that confers antibiotic resistance. To further understand the genetic diversity and dissemination of antibiotic resistance plasmids among EPEC strains, we describe the complete sequence of an antibiotic resistance plasmid from EPEC strain B171. The resistance plasmid, pB171_90, has a completed sequence length of 90,229 bp, a GC content of 54.55%, and carries protein-encoding genes involved in conjugative transfer, resistance to tetracycline (tetA), sulfonamides (sulI), and mercury, as well as several virulence-associated genes, including the transcriptional regulator hha and the putative calcium sequestration inhibitor (csi). In silico detection of the pB171_90 genes among 4,798 publicly available E. coli genome assemblies indicates that the unique genes of pB171_90 (csi and traI) are primarily restricted to genomes identified as EPEC or enterotoxigenic E. coli However, conserved regions of the pB171_90 plasmid containing genes involved in replication, stability, and antibiotic resistance were identified among diverse E. coli pathotypes. Interestingly, pB171_90 also exhibited significant similarity with a sequenced plasmid from Shigella dysenteriae type I. Our findings demonstrate the mosaic nature of EPEC antibiotic resistance plasmids and highlight the need for additional sequence-based characterization of antibiotic resistance plasmids harbored by pathogenic E. coli.

Uropathogenic Escherichia coli pathogenicity islands and other ExPEC virulence genes may contribute to the genome variability of enteroinvasive E. coli

BACKGROUND

Enteroinvasive Escherichia coli (EIEC) may be the causative agent of part of those million cases of diarrhea illness reported worldwide every year and attributable to Shigella. That is because both enteropathogens have many common characteristics that difficult their identification either by traditional microbiological methods or by molecular tools used in the clinical laboratory settings. While Shigella has been extensively studied, EIEC remains barely characterized at the molecular level. Recent EIEC important outbreaks, apparently generating more life-threatening cases, have prompted us to screen EIEC for virulence traits usually related to extraintestinal pathogenic E. coli (ExPEC). That could explain the appearance of EIEC strains presenting higher virulence potential.

RESULTS

EIEC strains were distributed mainly in three phylogroups in a serogroup-dependent manner. Serogroups O124, O136, O144, and O152 were exclusively classified in phylogroup A; O143 in group E; and O28ac and O29 in group B1. Only two serogroups showed diverse phylogenetic origin as follows: O164 was assigned to groups A, B1, C, and B2 (one strain each), and O167 in groups E (five strains), and A (one strain) (Table 1). Eleven of 20 virulence genes (VGs) searched were detected, and the majority of the 19 different VGs combinations found were serogroup-specific. Uropathogenic E. coli (UPEC) PAI genetic markers were detected in all EIEC strains. PAIs I_J96 and II_CFT073 were the most frequent (92.1 and 80.4%, respectively). PAI IV₅₃₆ was restricted to some serogroups from phylogroups A, B1 and E. PAI I_CFT073 was uniquely detected in phylogroups B2 and E. A total of 45 (88%) strains presented multiple PAI markers (two to four). PAIs I_J96 and II_CFT073 were found together in 80% of strains.

CONCLUSIONS

EIEC is a DEC pathovar that presents VGs and pathogenicity island genetic markers typically associated with ExPEC, especially UPEC. These features are distributed in a phylogenetic and serogroup-dependent manner suggesting the existence of stable EIEC subclones. The presence of phylogroups B2 and E strains allied to the presence of UPEC virulence-associated genes may underscore the ongoing evolution of EIEC towards a hypervirulent pathotype.

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.