Genome sequence of Shigella flexneri 2a: insights into pathogenicity through comparison with genomes of Escherichia coli K12 and O157

Pseudogenes in plasmid genomes reveal past transitions in plasmid mobility

Evidence for gene non-functionalization due to mutational processes is found in genomes in the form of pseudogenes. Pseudogenes are known to be rare in prokaryote chromosomes, with the exception of lineages that underwent an extreme genome reduction (e.g. obligatory symbionts). Much less is known about the frequency of pseudogenes in prokaryotic plasmids; those are genetic elements that can transfer between cells and may encode beneficial traits for their host. Non-functionalization of plasmid-encoded genes may alter the plasmid characteristics, e.g. mobility, or their effect on the host. Analyzing 10 832 prokaryotic genomes, we find that plasmid genomes are characterized by threefold-higher pseudogene density compared to chromosomes. The majority of plasmid pseudogenes correspond to deteriorated transposable elements. A detailed analysis of enterobacterial plasmids furthermore reveals frequent gene non-functionalization events associated with the loss of plasmid self-transmissibility. Reconstructing the evolution of closely related plasmids reveals that non-functionalization of the conjugation machinery led to the emergence of non-mobilizable plasmid types. Examples are virulence plasmids in Escherichia and Salmonella. Our study highlights non-functionalization of core plasmid mobility functions as one route for the evolution of domesticated plasmids. Pseudogenes in plasmids supply insights into past transitions in plasmid mobility that are akin to transitions in bacterial lifestyle.

Novel inhibitors that target bacterial virulence identified via HTS against intra-macrophage survival of Shigella flexneri

Shigella flexneri is a facultative intracellular pathogen that causes shigellosis, a human diarrheal disease characterized by the destruction of the colonic epithelium. Novel antimicrobial compounds to treat infections are urgently needed due to the proliferation of bacterial antibiotic resistance and lack of new effective antimicrobials in the market. Our approach to find compounds that block the Shigella virulence pathway has three potential advantages: (i) resistance development should be minimized due to the lack of growth selection pressure, (ii) no resistance due to environmental antibiotic exposure should be developed since the virulence pathways are not activated outside of host infection, and (iii) the normal intestinal microbiota, which do not have the targeted virulence pathways, should be unharmed. We chose to utilize two phenotypic assays, inhibition of Shigella survival in macrophages and Shigella growth inhibition (minimum inhibitory concentration), to interrogate the 1.7 M compound screening collection subset of the GlaxoSmithKline drug discovery chemical library. A number of secondary assays on the hit compounds resulting from the primary screens were conducted, which, in combination with chemical, structural, and physical property analyses, narrowed the final hit list to 44 promising compounds for further drug discovery efforts. The rapid development of antibiotic resistance is a critical problem that has the potential of returning the world to a "pre-antibiotic" type of environment, where millions of people will die from previously treatable infections. One relatively newer approach to minimize the selection pressures for the development of resistance is to target virulence pathways. This is anticipated to eliminate any resistance selection pressure in environmental exposure to virulence-targeted antibiotics and will have the added benefit of not affecting the non-virulent microbiome. This paper describes the development and application of a simple, reproducible, and sensitive assay to interrogate an extensive chemical library in high-throughput screening format for activity against the survival of Shigella flexneri 2457T-nl in THP-1 macrophages. The ability to screen very large numbers of compounds in a reasonable time frame (~1.7 M compounds in ~8 months) distinguishes this assay as a powerful tool in further exploring new compounds with intracellular effect on S. flexneri or other pathogens with similar pathways of pathogenesis. The assay utilizes a luciferase reporter which is extremely rapid, simple, relatively inexpensive, and sensitive and possesses a broad linear range. The assay also utilized THP-1 cells that resemble primary monocytes and macrophages in morphology and differentiation properties. THP-1 cells have advantages over human primary monocytes or macrophages because they are highly plastic and their homogeneous genetic background minimizes the degree of variability in the cell phenotype (1). The intracellular and virulence-targeted selectivity of our methodology, determined via secondary screening, is an enormous advantage. Our main interest focuses on hits that are targeting virulence, and the most promising compounds with adequate physicochemical and drug metabolism and pharmacokinetic (DMPK) properties will be progressed to a suitable in vivo shigellosis model to evaluate the therapeutic potential of this approach. Additionally, compounds that act via a host-directed mechanism could be a promising source for further research given that it would allow a whole new, specific, and controlled approach to the treatment of diseases caused by some pathogenic bacteria.

Bacteriophage P22 SieA mediated superinfection exclusion

Many temperate phages encode prophage-expressed functions that interfere with superinfection of the host bacterium by external phages. Salmonella phage P22 has four such systems that are expressed from the prophage in a lysogen that are encoded by the c2 (repressor), gtrABC, sieA, and sieB genes. Here we report that the P22-encoded SieA protein is the only phage protein required for exclusion by the SieA system, and that it is an inner membrane protein that blocks DNA injection by P22 and its relatives, but has no effect on infection by other tailed phage types. The P22 virion injects its DNA through the host cell membranes and periplasm via a conduit assembled from three "ejection proteins" after their release from the virion. Phage P22 mutants were isolated that overcome the SieA block, and they have amino acid changes in the C-terminal regions of the gene 16 and 20 encoded ejection proteins. Three different single amino acid changes in these proteins are required to obtain nearly full resistance to SieA. Hybrid P22 phages that have phage HK620 ejection protein genes are also partially resistant to SieA. There are three sequence types of extant phage-encoded SieA proteins that are less than 30% identical to one another, yet comparison of two of these types found no differences in target specificity. Our data are consistent with a model in which the inner membrane protein SieA interferes with the assembly or function of the periplasmic gp20 and membrane-bound gp16 DNA delivery conduit.

Toward a Shigella Vaccine: Opportunities and Challenges to Fight an Antimicrobial-Resistant Pathogen

Shigellosis causes more than 200,000 deaths worldwide and most of this burden falls on Low- and Middle-Income Countries (LMICs), with a particular incidence in children under 5 years of age. In the last decades, Shigella has become even more worrisome because of the onset of antimicrobial-resistant strains (AMR). Indeed, the WHO has listed Shigella as one of the priority pathogens for the development of new interventions. To date, there are no broadly available vaccines against shigellosis, but several candidates are being evaluated in preclinical and clinical studies, bringing to light very important data and information. With the aim to facilitate the understanding of the state-of-the-art of Shigella vaccine development, here we report what is known about Shigella epidemiology and pathogenesis with a focus on virulence factors and potential antigens for vaccine development. We discuss immunity after natural infection and immunization. In addition, we highlight the main characteristics of the different technologies that have been applied for the development of a vaccine with broad protection against Shigella.

Prenylation: A Critical Step for Biomanufacturing of Prenylated Aromatic Natural Products

Prenylated aromatic natural products (PANPs) have received much attention due to their biomedical benefits for human health. The prenylation of aromatic natural products (ANPs), which is mainly catalyzed by aromatic prenyltransferases (aPTs), contributes significantly to their structural and functional diversity by providing higher lipophilicity and enhanced bioactivity. aPTs are widely distributed in bacteria, fungi, animals, and plants and play a key role in the regiospecific prenylation of ANPs. Recent studies have greatly advanced our understanding of the characteristics and application of aPTs. In this review, we comment on research progress regarding sources, evolutionary relationships, structural features, reaction mechanism, engineering modification, and application of aPTs. Particular emphasis is also placed on recent advances, challenges, and prospects about applications of aPTs in microbial cell factories for producing PANPs. Generally, this review could provide guidance for using aPTs as robust biocatalytic tools to produce various PANPs with high efficiency.

The Diversity of Escherichia coli Pathotypes and Vaccination Strategies against This Versatile Bacterial Pathogen

Escherichia coli (E. coli) is a gram-negative bacillus and resident of the normal intestinal microbiota. However, some E. coli strains can cause diseases in humans, other mammals and birds ranging from intestinal infections, for example, diarrhea and dysentery, to extraintestinal infections, such as urinary tract infections, respiratory tract infections, meningitis, and sepsis. In terms of morbidity and mortality, pathogenic E. coli has a great impact on public health, with an economic cost of several billion dollars annually worldwide. Antibiotics are not usually used as first-line treatment for diarrheal illness caused by E. coli and in the case of bloody diarrhea, antibiotics are avoided due to the increased risk of hemolytic uremic syndrome. On the other hand, extraintestinal infections are treated with various antibiotics depending on the site of infection and susceptibility testing. Several alarming papers concerning the rising antibiotic resistance rates in E. coli strains have been published. The silent pandemic of multidrug-resistant bacteria including pathogenic E. coli that have become more difficult to treat favored prophylactic approaches such as E. coli vaccines. This review provides an overview of the pathogenesis of different pathotypes of E. coli, the virulence factors involved and updates on the major aspects of vaccine development against different E. coli pathotypes.

Transient Glycolytic Complexation of Arsenate Enhances Resistance in the Enteropathogen Vibrio cholerae

The ubiquitous presence of toxic arsenate (As^V) in the environment has raised mechanisms of resistance in all living organisms. Generally, bacterial detoxification of As^V relies on its reduction to arsenite (As^III) by ArsC, followed by the export of As^III by ArsB. However, how pathogenic species resist this metalloid remains largely unknown. Here, we found that Vibrio cholerae, the etiologic agent of the diarrheal disease cholera, outcompetes other enteropathogens when grown on millimolar concentrations of As^V. To do so, V. cholerae uses, instead of ArsCB, the As^V-inducible vc1068-1071 operon (renamed var for vibrio arsenate resistance), which encodes the arsenate repressor ArsR, an alternative glyceraldehyde-3-phosphate dehydrogenase, a putative phosphatase, and the As^V transporter ArsJ. In addition to Var, V. cholerae induces oxidative stress-related systems to counter reactive oxygen species (ROS) production caused by intracellular As^V. Characterization of the var mutants suggested that these proteins function independently from one another and play critical roles in preventing deleterious effects on the cell membrane potential and growth derived from the accumulation As^V. Mechanistically, we demonstrate that V. cholerae complexes As^V with the glycolytic intermediate 3-phosphoglycerate into 1-arseno-3-phosphoglycerate (1As3PG). We further show that 1As3PG is not transported outside the cell; instead, it is subsequently dissociated to enable extrusion of free As^V through ArsJ. Collectively, we propose the formation of 1As3PG as a transient metabolic storage of As^V to curb the noxious effect of free As^V. This study advances our understanding of As^V resistance in bacteria and underscores new points of vulnerability that might be an attractive target for antimicrobial interventions. IMPORTANCE Even though resistance to arsenate has been extensively investigated in environmental bacteria, how enteric pathogens tolerate this toxic compound remains unknown. Here, we found that the cholera pathogen V. cholerae exhibits increased resistance to arsenate compared to closely related enteric pathogens. Such resistance is promoted not by ArsC-dependent reduction of arsenate to arsenite but by an operon encoding an arsenate transporter (ArsJ), an alternative glyceraldehyde 3-phosphate dehydrogenase (VarG), and a putative, uncharacterized phosphatase (VarH). Mechanistically, we demonstrate that V. cholerae detoxifies arsenate by complexing it with the glycolytic intermediate 3-phosphoglycerate into 1-arseno-3-phosphoglycerate (1As3PG). 1As3PG is not transported outside the cell; instead, it is subsequently dissociated by VarH to enable extrusion of free arsenate through ArsJ. Collectively, this study proposes a novel mechanism for arsenate detoxification, entirely independent of arsenate reduction and arsenite extrusion, that enhances V. cholerae resistance to this metalloid compared to other enteric pathogens.

Genome Analysis of Enterobacter asburiae and Lelliottia spp. Proliferating in Oligotrophic Drinking Water Reservoirs and Lakes

Surface waters are one of the main sources for drinking water production, and thus microbial contamination should be as minimal as possible. However, high concentrations of coliform bacteria were detected in reservoirs and lakes used for drinking water production during summer months due to autochthonous proliferation processes. Here, we present the genomic analyses of 17 strains of Enterobacter asburiae and Lelliottia spp. proliferating in reservoirs and lakes with special focus on the hygienic relevance, antibiotic resistance, and adaptations to the oligotrophic environments. The genomes contain neither genes for the type III secretion system nor cytotoxins or hemolysins, which are considered typical virulence factors. Examination of antibiotic resistance genes revealed mainly efflux pumps and β-lactamase class C (ampC) genes. Phenotypically, single isolates of Enterobacter asburiae showed resistance to fosfomycin and ceftazidime. The genome analyses further suggest adaptations to oligotrophic and changing environmental conditions in reservoirs and lakes, e.g., genes to cope with low nitrate and phosphate levels and the ability to utilize substances released by algae, like amino acids, chitin, alginate, rhamnose, and fucose. This leads to the hypothesis that the proliferation of the coliform bacteria could occur at the end of summer due to algae die-off. IMPORTANCE Certain strains of coliform bacteria have been shown to proliferate in the oligotrophic water of drinking water reservoirs and lakes, reaching values above 10⁴ per 100 mL. Such high concentrations challenge drinking water treatment, and occasionally the respective coliform bacteria have been detected in the treated drinking water. Thus, the question of their hygienic relevance is of high importance for water suppliers and authorities. Our genomic analyses suggest that the strains are not hygienically relevant, as typical virulence factors are absent and antibiotic resistance genes in the genomes most likely are of natural origin. Furthermore, their presence in the water is not related to fecal contamination. The proliferation in reservoirs and lakes during stable summer stratification is an autochthonic process of certain E. asburiae and Lelliottia strains that are well adapted to the surrounding oligotrophic environment.

The NEL Family of Bacterial E3 Ubiquitin Ligases

Some pathogenic or symbiotic Gram-negative bacteria can manipulate the ubiquitination system of the eukaryotic host cell using a variety of strategies. Members of the genera Salmonella, Shigella, Sinorhizobium, and Ralstonia, among others, express E3 ubiquitin ligases that belong to the NEL family. These bacteria use type III secretion systems to translocate these proteins into host cells, where they will find their targets. In this review, we first introduce type III secretion systems and the ubiquitination process and consider the various ways bacteria use to alter the ubiquitin ligation machinery. We then focus on the members of the NEL family, their expression, translocation, and subcellular localization in the host cell, and we review what is known about the structure of these proteins, their function in virulence or symbiosis, and their specific targets.

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

Background

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

Results

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

Conclusion

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

Supplementary Information

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

Distinct Potentially Adaptive Accumulation of Truncation Mutations in Salmonella enterica serovar Typhi and Salmonella enterica serovar Paratyphi A

Gene inactivation through the accumulation of truncation (or premature stop codon) mutations is a common mode of evolution in bacteria. It is frequently believed to result from reductive evolutionary processes allowing purging of superfluous traits. However, several works have demonstrated that, similar to the occurrences of inactivating nonsynonymous (i.e., amino acid replacement) mutations under positive selection pressures, truncation mutations can also be adaptive where specific traits deleterious in particular environmental conditions need to be inactivated for survival. Here, we performed a comparative analysis of genome-wide accumulation of truncation mutations in Salmonella enterica serovar Typhi and Salmonella enterica serovar Paratyphi A. Considering the known convergent evolutionary trajectories in these two serovars, we expected a strong overlap of truncated genes in S. Typhi and S. Paratyphi A, emerging through either reductive or adaptive dynamics. However, we detected a distinct set of core truncated genes encoding different overrepresented functional clusters in each serovar. In 54% and 28% truncated genes in S. Typhi and S. Paratyphi A, respectively, inactivating mutations were acquired by only different subsets of isolates, instead of all isolates analyzed for that serovar. Importantly, 62% truncated genes (P < 0.001) in S. Typhi and S. Paratyphi A were also targeted by convergent amino acid mutations in different serovars, suggesting those genes to be under selection pressures. Our findings indicate significant presence of potentially adaptive truncation mutations in conjunction with the ones emerging due to reductive evolution. Further experimental and large-scale bioinformatic studies are necessary to better explore the impact of such adaptive footprints of truncation mutations in the evolution of bacterial virulence. IMPORTANCE Detecting the adaptive mutations leading to gene inactivation or loss of function is crucial for understanding their contribution in the evolution of bacterial virulence and antibiotic resistance. Such inactivating mutations, apart from being of nonsynonymous (i.e., amino acid replacement) nature, can also be truncation mutations, abruptly trimming the length of encoded proteins. Importantly, the notion of reductive evolutionary dynamics is primarily accepted toward the accumulation of truncation mutations. However, our case study on S. Typhi and S. Paratyphi A, two human-restricted systemically invasive pathogens exerting similar clinical manifestations, indicated that a significant proportion of truncation mutations emerge from positive selection pressures. The candidate genes from our study will enable directed functional assays for deciphering the adaptive role of truncation mutations in S. Typhi and S. Paratyphi A pathogenesis. Also, our genome-level analytical approach will pave the way to understand the contribution of truncation mutations in the adaptive evolution of other bacterial pathogens.

Impact of horizontal gene transfer on emergence and stability of cooperative virulence in Salmonella Typhimurium

Intestinal inflammation fuels the transmission of Salmonella Typhimurium (S.Tm). However, a substantial fitness cost is associated with virulence expression. Mutations inactivating transcriptional virulence regulators generate attenuated variants profiting from inflammation without enduring virulence cost. Such variants interfere with the transmission of fully virulent clones. Horizontal transfer of functional regulatory genes (HGT) into attenuated variants could nevertheless favor virulence evolution. To address this hypothesis, we cloned hilD, coding for the master regulator of virulence, into a conjugative plasmid that is highly transferrable during intestinal colonization. The resulting mobile hilD allele allows virulence to emerge from avirulent populations, and to be restored in attenuated mutants competing against virulent clones within-host. However, mutations inactivating the mobile hilD allele quickly arise. The stability of virulence mediated by HGT is strongly limited by its cost, which depends on the hilD expression level, and by the timing of transmission. We conclude that robust evolution of costly virulence expression requires additional selective forces such as narrow population bottlenecks during transmission.

Salmonella Typhimurium virulence is costly and can be lost by mutation during infection. Bakkeren et al. show that virulence restoration via horizontal gene transfer is only transient while transmission bottlenecks promote long-term virulence stability.

Effect of Different Initial Fermentation pH on Exopolysaccharides Produced by Pseudoalteromonas agarivorans Hao 2018 and Identification of Key Genes Involved in Exopolysaccharide Synthesis via Transcriptome Analysis

Exopolysaccharides (EPSs) are carbohydrate polymers produced and secreted by microorganisms. In a changing marine environment, EPS secretion can reduce damage from external environmental disturbances to microorganisms. Meanwhile, EPSs have promising application prospects in the fields of food, cosmetics, and pharmaceuticals. Changes in external environmental pH have been shown to affect the synthesis of EPSs in microorganisms. In this study, we analyzed the effects of different initial fermentation pHs on the production, monosaccharide composition, and antioxidant activity of the EPSs of Pseudoalteromonas agarivorans Hao 2018. In addition, the transcriptome sequence of P. agarivorans Hao 2018 under different initial fermentation pH levels was determined. GO and KEGG analyses showed that the differentially expressed genes were concentrated in the two-component regulatory system and bacterial chemotaxis pathways. We further identified the expression of key genes involved in EPS synthesis during pH changes. In particular, the expression of genes encoding the glucose/galactose MFS transporter, phosphomannomutase, and mannose-1-phosphate guanylyltransferase was upregulated when the environmental pH increased, thus promoting EPS synthesis. This study not only contributes to elucidating the environmental adaptation mechanisms of P. agarivorans, but also provides important theoretical guidance for the directed development of new products using biologically active polysaccharides.

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.

Prevalence of Methicillin-Resistant Staphylococcus Species Among Filarial Lymphedema Patients in Ahanta West District of Ghana

Background

Filarial pathologies such as lymphedema may be associated with complications such as chronic non-healing wounds. Nonetheless, the role of bacterial population colonizing the lymphedematous legs has been posited to worsen the conditions of those living with the infection. These bacteria are usually composed of staphylococcal species partly because they are commensals. Thus, this present study sought to type the methicillin-resistant Staphylococcus aureus (MRSA) prevalence among individuals presenting with filarial lymphedema, particularly as MRSA tends to affect treatments options.

Methods

We recruited individuals (n = 321) with stages I–VII of lymphedema in a cross-sectional study in the Ahanta West district of the Western Region of Ghana. Swabs from lymphedematous limb ulcers, pus, and cutaneous surfaces were cultured using standard culture-based techniques. The culture isolates were later identified using Matrix-assisted Laser Desorption/Ionization Time of Flight (MALDI-TOF) mass spectrometry.

Results

A total of 192 Staphylococci species were isolated, with an overall prevalence of 39.7% (95% CI: 35%–44%; N = 483). S. hominis was the most prevalent species (23.95%), followed by S. haemolyticus (20.83%), S. epidermidis (15.10%), S. aureus (10.41%), and S. saprophyticus (9.32%). The remaining 20.34% were distributed among S. wanneri, S. sciuri, S. pasteuri, S. xylosus, S. simulans, S. cohnii, S. caprae, S. lugdunensis, and S. capitis. MRSA, containing mecA gene, was detected in 21 out of 31 Staphylococci isolates tested, with an overall prevalence of 68% (95% CI: 51%–84%). In addition, a virulent gene, Panton–Valentine leukocidin (PVL), which is usually associated with S. aureus, was detected in 20/31 (64.5%) S. aureus in the study.

Conclusion

These results suggest that MRSA species may pose a challenge to the treatment of filarial lymphedema with antibiotics particularly, as doxycycline is currently being piloted in some endemic areas to treat the infection. Thus, intensive antimicrobial resistance surveillance should be conducted in endemic areas by health authorities to forestall the dilemma of multidrug resistance not only against lymphatic filariasis (LF) infection but other diseases.

Shigella flexneri Diguanylate Cyclases Regulate Virulence

Shigella flexneri is an intracellular human pathogen that invades colonic cells and causes bloody diarrhea. S. flexneri evolved from commensal Escherichia coli, and genome comparisons reveal that S. flexneri has lost approximately 20% of its genes through the process of pathoadaptation, including a disproportionate number of genes associated with the turnover of the nucleotide-based second messenger cyclic di-GMP (c-di-GMP); however, the remaining c-di-GMP turnover enzymes are highly conserved. c-di-GMP regulates many behavioral changes in other bacteria in response to changing environmental conditions, including biofilm formation, but this signaling system has not been examined in S. flexneri. In this study, we expressed VCA0956, a constitutively active c-di-GMP synthesizing diguanylate cyclase (DGC) from Vibrio cholerae, in S. flexneri to determine if virulence phenotypes were regulated by c-di-GMP. We found that expressing VCA0956 in S. flexneri increased c-di-GMP levels, and this corresponds with increased biofilm formation and reduced acid resistance, host cell invasion, and plaque size. We examined the impact of VCA0956 expression on the S. flexneri transcriptome and found that genes related to acid resistance were repressed, and this corresponded with decreased survival to acid shock. We also found that individual S. flexneri DGC mutants exhibit reduced biofilm formation and reduced host cell invasion and plaque size, as well as increased resistance to acid shock. This study highlights the importance of c-di-GMP signaling in regulating S. flexneri virulence phenotypes. IMPORTANCE The intracellular human pathogen Shigella causes dysentery, resulting in as many as one million deaths per year. Currently, there is no approved vaccine for the prevention of shigellosis, and the incidence of antimicrobial resistance among Shigella species is on the rise. Here, we explored how the widely conserved c-di-GMP bacterial signaling system alters Shigella behaviors associated with pathogenesis. We found that expressing or removing enzymes associated with c-di-GMP synthesis results in changes in Shigella's ability to form biofilms, invade host cells, form lesions in host cell monolayers, and resist acid stress.

Evaluation of whole-genome sequencing-based subtyping methods for the surveillance of Shigella spp. and the confounding effect of mobile genetic elements in long-term outbreaks

Many public health laboratories across the world have implemented whole-genome sequencing (WGS) for the surveillance and outbreak detection of foodborne pathogens. PulseNet-affiliated laboratories have determined that most single-strain foodborne outbreaks are contained within 0–10 multi-locus sequence typing (MLST)-based allele differences and/or core genome single-nucleotide variants (SNVs). In addition to being a food- and travel-associated outbreak pathogen, most Shigella spp. cases occur through continuous person-to-person transmission, predominantly involving men who have sex with men (MSM), leading to long-term and recurrent outbreaks. Continuous transmission patterns coupled to genetic evolution under antibiotic treatment pressure require an assessment of existing WGS-based subtyping methods and interpretation criteria for cluster inclusion/exclusion. An evaluation of 4 WGS-based subtyping methods [SNVPhyl, coreMLST, core genome MLST (cgMLST) and whole-genome MLST (wgMLST)] was performed on 9 foodborne-, travel- and MSM-related retrospective outbreaks from a collection of 91 Shigella flexneri and 232  Shigella sonnei isolates to determine the methods’ epidemiological concordance, discriminatory power, robustness and ability to generate stable interpretation criteria. The discriminatory powers were ranked as follows: coreMLST<SNVPhyl<cgMLST<wgMLST (range: 0.970–1.000). The genetic differences observed for non-MSM-related Shigella spp. outbreaks respect the standard 0–10 allele/SNV guideline; however, mobile genetic element (MGE)-encoded loci caused inflated genetic variation and discrepant phylogenies for prolonged MSM-related S. sonnei outbreaks via wgMLST. The S. sonnei correlation coefficients of wgMLST were also the lowest at 0.680, 0.703 and 0.712 for SNVPhyl, coreMLST and cgMLST, respectively. Plasmid maintenance, mobilization and conjugation-associated genes were found to be the main source of genetic distance inflation in addition to prophage-related genes. Duplicated alleles arising from the repeated nature of IS elements were also responsible for many false cg/wgMLST differences. The coreMLST approach was shown to be the most robust, followed by SNVPhyl and wgMLST for inter-laboratory comparability. Our results highlight the need for validating species-specific subtyping methods based on microbial genome plasticity and outbreak dynamics in addition to the importance of filtering confounding MGEs for cluster detection.

VirB, a key transcriptional regulator of virulence plasmid genes in Shigella flexneri, forms DNA-binding site dependent foci in the bacterial cytoplasm

VirB is a key regulator of genes located on the large virulence plasmid (pINV) in the bacterial pathogen Shigella flexneri VirB is unusual; it is not related to other transcriptional regulators, instead, it belongs to a family of proteins that primarily function in plasmid and chromosome partitioning; exemplified by ParB. Despite this, VirB does not function to segregate DNA, but rather counters transcriptional silencing mediated by the nucleoid structuring protein, H-NS. Since ParB localizes subcellularly as discrete foci in the bacterial cytoplasm, we chose to investigate the subcellular localization of VirB to gain novel insight into how VirB functions as a transcriptional anti-silencer. To do this, a GFP-VirB fusion that retains the regulatory activity of VirB and yet, does not undergo significant protein degradation in S. flexneri, was used. Surprisingly, discrete fluorescent foci were observed in live wild-type S. flexneri cells and an isogenic virB mutant using fluorescence microscopy. In contrast, foci were rarely observed (<10%) in pINV-cured cells or in cells expressing a GFP-VirB fusion carrying amino acid substitutions in the VirB DNA binding domain. Finally, the 25 bp VirB-binding site was demonstrated to be sufficient and necessary for GFP-VirB focus formation using a set of small surrogate plasmids. Combined, these data demonstrate that the VirB:DNA interactions required for the transcriptional anti-silencing activity of VirB on pINV are a prerequisite for the subcellular localization of VirB in the bacterial cytoplasm. The significance of these findings, in light of the anti-silencing activity of VirB, is discussed.ImportanceThis study reveals the subcellular localization of VirB, a key transcriptional regulator of virulence genes found on the large virulence plasmid (pINV) in Shigella. Fluorescent signals generated by an active GFP-VirB fusion form 2, 3, or 4 discrete foci in the bacterial cytoplasm, predominantly at the quarter cell position. These signals are completely dependent upon VirB interacting with its DNA binding site found either on the virulence plasmid or an engineered surrogate. Our findings: 1) provide novel insight into VirB:pINV interactions, 2) suggest that VirB may have utility as a DNA marker, and 3) raise questions about how and why this anti-silencing protein that controls virulence gene expression on pINV of Shigella spp. forms discrete foci/hubs within the bacterial cytoplasm.

Binding of Campylobacter jejuni FliW Adjacent to the CsrA RNA-Binding Pockets Modulates CsrA Regulatory Activity

Campylobacter jejuni CsrA is an mRNA-binding, post-transcriptional regulator that controls many metabolic- and virulence-related characteristics of this important pathogen. In contrast to E. coli CsrA, whose activity is modulated by binding to small non-coding RNAs (sRNAs), C. jejuni CsrA activity is controlled by binding to the CsrA antagonist FliW. In this study, we identified the FliW binding site on CsrA. Deletion of the C-terminus of C. jejuni CsrA, which is extended relative to sRNA-binding CsrA proteins, abrogated FliW binding. Bacterial two-hybrid experiments were used to assess the interaction of FliW with wild-type CsrA and mutants thereof, in which every amino acid was individually mutated. Two CsrA mutations (V51A and N55A) resulted in a significant decrease in FliW binding. The V51A and N55A mutants also showed a decrease in CsrA-FliW complex formation, as assessed by size-exclusion chromatography and surface plasmon resonance. These residues were highly conserved in bacterial species containing CsrA orthologs whose activities are predicted to be regulated by FliW. The location of FliW binding was immediately adjacent to the two RNA-binding sites of the CsrA homodimer, suggesting the model that FliW binding to CsrA modulates its ability to bind to its mRNA targets either by steric hindrance, electrostatic repulsion, or by altering the overall structure of the RNA-binding sites.

Genetic Dissection of the Fermentative and Respiratory Contributions Supporting Vibrio cholerae Hypoxic Growth

Both fermentative and respiratory processes contribute to bacterial metabolic adaptations to low oxygen tension (hypoxia). In the absence of O2 as a respiratory electron sink, many bacteria utilize alternative electron acceptors, such as nitrate (NO3-). During canonical NO3- respiration, NO3- is reduced in a stepwise manner to N2 by a dedicated set of reductases. Vibrio cholerae, the etiological agent of cholera, requires only a single periplasmic NO3- reductase (NapA) to undergo NO3- respiration, suggesting that the pathogen possesses a noncanonical NO3- respiratory chain. In this study, we used complementary transposon-based screens to identify genetic determinants of general hypoxic growth and NO3- respiration in V. cholerae We found that while the V. cholerae NO3- respiratory chain is primarily composed of homologues of established NO3- respiratory genes, it also includes components previously unlinked to this process, such as the Na+-NADH dehydrogenase Nqr. The ethanol-generating enzyme AdhE was shown to be the principal fermentative branch required during hypoxic growth in V. cholerae Relative to single adhE or napA mutant strains, a V. cholerae strain lacking both genes exhibited severely impaired hypoxic growth in vitro and in vivo Our findings reveal the genetic basis of a specific interaction between disparate energy production pathways that supports pathogen fitness under shifting conditions. Such metabolic specializations in V. cholerae and other pathogens are potential targets for antimicrobial interventions.IMPORTANCE Bacteria reprogram their metabolism in environments with low oxygen levels (hypoxia). Typically, this occurs via regulation of two major, but largely independent, metabolic pathways: fermentation and respiration. In this study, we found that the diarrheal pathogen Vibrio cholerae has a respiratory chain for NO3- that consists largely of components found in other NO3- respiratory systems but also contains several proteins not previously linked to this process. Both AdhE-dependent fermentation and NO3- respiration were required for efficient pathogen growth under both laboratory conditions and in an animal infection model. These observations provide a specific example of fermentative respiratory interactions and identify metabolic vulnerabilities that may be targetable for new antimicrobial agents in V. cholerae and related pathogens.

Multiple PCR assay based on the cigR gene for detection of Salmonella spp. and Salmonella Pullorum/Gallinarum identification

Salmonella spp. are important zoonotic pathogens that are responsible for severe diseases in both animals and humans. Salmonella enterica subsp. enterica serovar Gallinarum biovar Gallinarum (S. Gallinarum) and biovar Pullorum (S. Pullorum) are typical infectious pathogens detected in the chicken industry that have caused great economic losses. To facilitate their detection and prevent contamination, we developed a rapid multiple PCR method, which can simultaneously detect Salmonella spp. and further identify the biovars S. Pullorum/Gallinarum. This PCR detection method is based on the cigR gene, which is conserved among Salmonella spp. but has a 42-bp deletion in S. Pullorum/Gallinarum. The specificity and sensitivity of the PCR assay was evaluated with 41 different strains: 34 Salmonella strains, including 5 S. Pullorum/Gallinarum strains, and 7 non-Salmonella strains. The lower limit of detection was 8.15 pg of S. Pullorum (S06004) genomic DNA and 20 cfu in PCR, which shows a great sensitivity. In addition, this method was applied to detect or identify Salmonella from processing chicken liver and egg samples, and the results corresponded to those obtained from serotype analysis using the conventional slide agglutination test. Overall, the new cigR-based PCR assay is efficient and practical for Salmonella detection and S. Pullorum/Gallinarum identification and will greatly reduce the workload of epidemiologic investigation.

Phylogenetic Characterization Reveals Prevalent Shigella flexneri ST100 Clone in Beijing, China, 2005 to 2018

Beijing is the largest transportation hub in China, with a highly mobile population. Shigella flexneri is a major cause of bacillary dysentery in Beijing. However, little is known about the genetic features and population structure of locally circulating S. flexneri clones. Whole-genome sequencing of 93 S. flexneri isolates revealed that S. flexneri epidemics in Beijing were predominantly caused by an ST100 clone. Interregional spread, rapid local expansion, and acquirement of antimicrobial resistance determinants have cocontributed to the epidemics of this clone. Another ST18 clone was also identified and showed long-term colonization in Beijing. Our study provides comprehensive insights into the population structure and evolutionary history of S. flexneri in Beijing.

Shigella flexneri is a major cause of bacillary dysentery in Beijing, China. The genetic features and population structure of locally circulating clones remained unclear. In this study, we sequenced the genomes of 93 S. flexneri isolates from patients in Beijing from 2005 to 2018. Phylogenetic analysis revealed a predominant lineage comprised of ST100 isolates that had acquired an extensive repertoire of antimicrobial resistance determinants. A rapid local expansion of the largest clade of this lineage began in 2008 and gradually resulted in the dominance of serotype 2a. Other clades showed substantial evidence of interregional spread from other areas of China. Another lineage consisting of ST18 isolates was also identified and appeared to have persisted locally for nearly 6 decades. These findings suggest that S. flexneri epidemics in Beijing were caused by both local expansion and interregional transmission.

IMPORTANCE Beijing is the largest transportation hub in China, with a highly mobile population. Shigella flexneri is a major cause of bacillary dysentery in Beijing. However, little is known about the genetic features and population structure of locally circulating S. flexneri clones. Whole-genome sequencing of 93 S. flexneri isolates revealed that S. flexneri epidemics in Beijing were predominantly caused by an ST100 clone. Interregional spread, rapid local expansion, and acquirement of antimicrobial resistance determinants have cocontributed to the epidemics of this clone. Another ST18 clone was also identified and showed long-term colonization in Beijing. Our study provides comprehensive insights into the population structure and evolutionary history of S. flexneri in Beijing.

Complete genome sequence and annotation of the laboratory reference strain Shigella flexneri serotype 5a M90T and genome-wide transcriptional start site determination

Background

Shigella is a Gram-negative facultative intracellular bacterium that causes bacillary dysentery in humans. Shigella invades cells of the colonic mucosa owing to its virulence plasmid-encoded Type 3 Secretion System (T3SS), and multiplies in the target cell cytosol. Although the laboratory reference strain S. flexneri serotype 5a M90T has been extensively used to understand the molecular mechanisms of pathogenesis, its complete genome sequence is not available, thereby greatly limiting studies employing high-throughput sequencing and systems biology approaches.

Results

We have sequenced, assembled, annotated and manually curated the full genome of S. flexneri 5a M90T. This yielded two complete circular contigs, the chromosome and the virulence plasmid (pWR100). To obtain the genome sequence, we have employed long-read PacBio DNA sequencing followed by polishing with Illumina RNA-seq data. This provides a new hybrid strategy to prepare gapless, highly accurate genome sequences, which also cover AT-rich tracks or repetitive sequences that are transcribed. Furthermore, we have performed genome-wide analysis of transcriptional start sites (TSS) and determined the length of 5′ untranslated regions (5′-UTRs) at typical culture conditions for the inoculum of in vitro infection experiments. We identified 6723 primary TSS (pTSS) and 7328 secondary TSS (sTSS). The S. flexneri 5a M90T annotated genome sequence and the transcriptional start sites are integrated into RegulonDB (http://regulondb.ccg.unam.mx) and RSAT (http://embnet.ccg.unam.mx/rsat/) databases to use their analysis tools in the S. flexneri 5a M90T genome.

Conclusions

We provide the first complete genome for S. flexneri serotype 5a, specifically the laboratory reference strain M90T. Our work opens the possibility of employing S. flexneri M90T in high-quality systems biology studies such as transcriptomic and differential expression analyses or in genome evolution studies. Moreover, the catalogue of TSS that we report here can be used in molecular pathogenesis studies as a resource to know which genes are transcribed before infection of host cells. The genome sequence, together with the analysis of transcriptional start sites, is also a valuable tool for precise genetic manipulation of S. flexneri 5a M90T. Further, we present a new hybrid strategy to prepare gapless, highly accurate genome sequences. Unlike currently used hybrid strategies combining long- and short-read DNA sequencing technologies to maximize accuracy, our workflow using long-read DNA sequencing and short-read RNA sequencing provides the added value of using non-redundant technologies, which yield distinct, exploitable datasets.

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.

Selection-Driven Gene Inactivation in Salmonella

Bacterial genes are sometimes found to be inactivated by mutation. This inactivation may be observable simply because selection for function is intermittent or too weak to eliminate inactive alleles quickly. Here, I investigate cases in Salmonella enterica where inactivation is instead positively selected. These are identified by a rate of introduction of premature stop codons to a gene that is higher than expected under selective neutrality, as assessed by comparison to the rate of synonymous changes. I identify 84 genes that meet this criterion at a 10% false discovery rate. Many of these genes are involved in virulence, motility and chemotaxis, biofilm formation, and resistance to antibiotics or other toxic substances. It is hypothesized that most of these genes are subject to an ongoing process in which inactivation is favored under rare conditions, but the inactivated allele is deleterious under most other conditions and is subsequently driven to extinction by purifying selection.

Investigating host-bacterial interactions among enteric pathogens

Background

In 2017, World Health Organization (WHO) published a catalogue of 12 families of antibiotic-resistant “priority pathogens” that are posing the greatest threats to human health. Six of these dreaded pathogens are known to infect the human gastrointestinal system. In addition to causing gastrointestinal and systemic infections, these pathogens can also affect the composition of other microbes constituting the healthy gut microbiome. Such aberrations in gut microbiome can significantly affect human physiology and immunity. Identifying the virulence mechanisms of these enteric pathogens are likely to help in developing newer therapeutic strategies to counter them.

Results

Using our previously published in silico approach, we have evaluated (and compared) Host-Pathogen Protein-Protein Interaction (HPI) profiles of four groups of enteric pathogens, namely, different species of Escherichia, Shigella, Salmonella and Vibrio. Results indicate that in spite of genus/ species specific variations, most enteric pathogens possess a common repertoire of HPIs. This core set of HPIs are probably responsible for the survival of these pathogen in the harsh nutrient-limiting environment within the gut. Certain genus/ species specific HPIs were also observed.

Conslusions

The identified bacterial proteins involved in the core set of HPIs are expected to be helpful in understanding the pathogenesis of these dreaded gut pathogens in greater detail. Possible role of genus/ species specific variations in the HPI profiles in the virulence of these pathogens are also discussed. The obtained results are likely to provide an opportunity for development of novel therapeutic strategies against the most dreaded gut pathogens.

Human Enteric Defensin 5 Promotes Shigella Infection of Macrophages

Human α-defensins are 3- to 5-kDa disulfide-bridged peptides with a multitude of antimicrobial activities and immunomodulatory functions. Recent studies show that human enteric α-defensin 5 (HD5), a host defense peptide important for intestinal homeostasis and innate immunity, aids the highly infectious enteropathogen Shigella in breaching the intestinal epithelium in vitro and in vivo Whether and how HD5 influences Shigella infection of resident macrophages following its invasion of the intestinal epithelium remain poorly understood. Here, we report that HD5 greatly promoted phagocytosis of Shigella by macrophages by targeting the bacteria to enhance bacterium-to-cell contacts in a structure- and sequence-dependent fashion. Subsequent intracellular multiplication of phagocytosed Shigella led to massive necrotic cell death and release of the bacteria. HD5-promoted phagocytosis of Shigella was independent of the status of the type 3 secretion system. Furthermore, HD5 neither inhibited nor enhanced phagosomal escape of Shigella Collectively, these findings confirm a potential pathogenic role of HD5 in Shigella infection of not only epithelial cells but also macrophages, illuminating how an enteropathogen exploits a host protective factor for virulence and infection.

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.

Structural analysis of Shigella flexneri bi-functional enzyme HisIE in histidine biosynthesis

Histidine biosynthesis, which is absent in animals, was shown to be highly conserved among gram-negative bacteria, thus making it an attractive target for antibiotic design. There are many fusion forms of enzymes in the histidine biosynthetic pathway and people still have limited knowledge about their domain organizations and catalytic mechanisms, due to the lack of structural information. Here we report the first crystal structure of Shigella flexneri bi-functional enzyme HisIE (SfHisIE) that functions in the 2nd and 3rd steps in the histidine biosynthetic pathway. This structure shows that HisIE exists as dimers with two loops (fusion loop) connecting the individual dimer of HisE and HisI in its N-terminus and C-terminus respectively. Our mutagenesis study shows mutations in this fusion loop are lethal for bacteria indicating the advantage of gene fusion in Histidine biosynthesis. Structural analysis revealed several highly conserved residues in the putative ligand binding grooves of HisE and HisI, showing an evolutionarily conserved catalytic mechanism shared among gram negative-bacteria.

Critical determinants of human neutrophil peptide 1 for enhancing host epithelial adhesion of Shigella flexneri

Human neutrophil peptides (HNPs), also known as human myeloid α-defensins degranulated by infiltrating neutrophils at bacterial infection loci, exhibit broad antomicrobial activities against bacteria, fungi, and viruses. We have made a surprising recent finding that Shigella, a highly contagious, yet poorly adhesive enteric pathogen, exploits human α-defensins including HNP1 to enhance its adhesion to and invasion of host epithelial cells. However, the critical molecular determinants responsible for HNP1-enhanced Shigella adhesion and invasion have yet to be investigated. Using cultured epithelial cells and polarised Caco2 cells as an in vitro infection model, we demonstrated that HNP1 promoted Shigella infection in a structure- and sequence-dependent manner, with two bulky hydrophobic residues, Trp26 and Phe28 important for HNP1 self-assembly, being most critical. The functional importance of hydrophobicity for HNP1-enhanced Shigella infection was further verified by substitutions for Trp26 of a series of unnatural amino acids with straight aliphatic side chains of different lengths. Dissection of the Shigella infection process revealed that bacteria-rather than host cells-bound HNP1 contributed most to the enhancement. Further, mutagenesis analysis of bacterial surface components, while precluding the involvement of lipopolysaccharides (LPS) in the interaction with HNP1, identified outer membrane proteins and the Type 3 secretion apparatus as putative binding targets of HNP1 involved in enhanced Shigella adhesion and invasion. Our findings provide molecular and mechanistic insights into the mode of action of HNP1 in promoting Shigella infection, thus showcasing another example of how innate immune factors may serve as a double-edged sword in health and disease.

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.

RNA-Seq analysis of the T3SA regulon in Shigella flexneri reveals two new chromosomal genes upregulated in the on-state

Shigella spp. are enterobacteria that invade human colonic mucosal cells using their Type Three Secretion Apparatus (T3SA). Shigella spp. possess a large plasmid that encodes most of its virulence factors and has been the focus of seminal work that defined the T3SA regulon. Thus, a global assessment of the transcriptional response regulated by the T3SA has been lacking. Herein we used RNA-Seq to identify genes that are differentially expressed when the T3SA is active (on-state) versus inactive (off-state). The quality of the RNA-Seq dataset was validated by its correlation with a prior microarray study. Using novel insights about the expression of non-coding regions, bioinformatic tools and experimentations, we demonstrated the existence of six operons and evidence that ipaH2.5 is a pseudogene. In addition, 86 chromosomal genes were downregulated in the on-state including several non-coding transcripts corresponding to short antisense RNA embedded in the 16S and 23S RNA genes, and 40 coding transcripts, whose cognate proteins were highly connected at the genetic and biochemical levels. Finally, we identified two novel chromosomal genes dubbed gem1 and gem3, which were upregulated in the on-state similarly to genes belonging to the T3SA regulon. The latter findings were validated on biological triplicates by droplet digital PCR. To our knowledge gem1 and gem3 are the first chromosomal members of the T3SA regulon that have no homologs on the plasmid. Our approach provides a path to optimizing RNA-Seq studies in case of bacterial models that had previously been the subject of medium to large scale studies.

In silicoAnalysis of L-Glutaminase from Extremophiles

Background:

L-glutaminase enzyme belongs to the family of hydrolases, those acting on carbon-nitrogen bonds other than peptide bonds, specifically in linear amides. Protein L-glutaminase, which converts amino acid glutamine to a glutamate residue, is useful as antileukemic agent, antiretroviral agent and a new food-processing enzyme.

Objective:

The sequences representing L-glutaminase from extremophiles were analyzed for different physico-chemical properties and to relate these observed differences to their extremophilic properties, phylogenetic tree construction and the evolutionary relationship among them.

Methods:

In this work, in silico analysis of amino acid sequences of extremophilic (thermophile, halophile and psychrophiles) proteins has been done. The physiochemical properties of these four groups of proteins for L-glutaminase also differ in number of amino acids, aliphatic index and grand average of hydropathicity (GRAVY).

Result:

The GRAVY was found to be significantly high in thermophilic (2.29 fold) and psychrophilic bacteria (3.3 fold) as compare to mesophilic bacteria. The amino acid Cys (C) was found to be statistically significant in mesophilic bacteria (approximately or more than 3 fold) as compared to the abundance of this amino acid in extremophilic bacteria.

Conclusion:

Multiple sequence alignment revealed the domain/motif for glutaminase that consists of Ser-74, Lys-77, Asn-126, Lys-268, and Ser-269, which is highly conserved in all microorganisms.

Investigating the DNA-Binding Site for VirB, a Key Transcriptional Regulator of Shigella Virulence Genes, Using an In Vivo Binding Tool

The transcriptional anti-silencing and DNA-binding protein, VirB, is essential for the virulence of Shigella species and, yet, sequences required for VirB-DNA binding are poorly understood. While a 7-8 bp VirB-binding site has been proposed, it was derived from studies at a single VirB-dependent promoter, icsB. Our previous in vivo studies at a different VirB-dependent promoter, icsP, found that the proposed VirB-binding site was insufficient for regulation. Instead, the required site was found to be organized as a near-perfect inverted repeat separated by a single nucleotide spacer. Thus, the proposed 7-8 bp VirB-binding site needed to be re-evaluated. Here, we engineer and validate a molecular tool to capture protein-DNA binding interactions in vivo. Our data show that a sequence organized as a near-perfect inverted repeat is required for VirB-DNA binding interactions in vivo at both the icsB and icsP promoters. Furthermore, the previously proposed VirB-binding site and multiple sites found as a result of its description (i.e., sites located at the virB, virF, spa15, and virA promoters) are not sufficient for VirB to bind in vivo using this tool. The implications of these findings are discussed.

The species-spanning family of LPX-motif harbouring effector proteins

The delivery of effector proteins into infected eukaryotic cells represents a key virulence feature of many microbial pathogens in order to derail essential cellular processes and effectively counter the host defence system. Although bacterial effectors are truly numerous and exhibit a wide range of biochemical activities, commonalities in terms of protein structure and function shared by many bacterial pathogens exist. Recent progress has shed light on a species-spanning family of bacterial effectors containing an LPX repeat motif as a subtype of the leucine-rich repeat superfamily, partially combined with a novel E3 ubiquitin ligase domain. This review highlights the immunomodulatory effects of LPX effector proteins, with particular emphasis on the exploitation of the host ubiquitin system.

A low-complexity metabolic network model for the respiratory and fermentative metabolism of Escherichia coli

Over the last decades, predictive microbiology has made significant advances in the mathematical description of microbial spoiler and pathogen dynamics in or on food products. Recently, the focus of predictive microbiology has shifted from a (semi-)empirical population-level approach towards mechanistic models including information about the intracellular metabolism in order to increase model accuracy and genericness. However, incorporation of this subpopulation-level information increases model complexity and, consequently, the required run time to simulate microbial cell and population dynamics. In this paper, results of metabolic flux balance analyses (FBA) with a genome-scale model are used to calibrate a low-complexity linear model describing the microbial growth and metabolite secretion rates of Escherichia coli as a function of the nutrient and oxygen uptake rate. Hence, the required information about the cellular metabolism (i.e., biomass growth and secretion of cell products) is selected and included in the linear model without incorporating the complete intracellular reaction network. However, the applied FBAs are only representative for microbial dynamics under specific extracellular conditions, viz., a neutral medium without weak acids at a temperature of 37℃. Deviations from these reference conditions lead to metabolic shifts and adjustments of the cellular nutrient uptake or maintenance requirements. This metabolic dependency on extracellular conditions has been taken into account in our low-complex metabolic model. In this way, a novel approach is developed to take the synergistic effects of temperature, pH, and undissociated acids on the cell metabolism into account. Consequently, the developed model is deployable as a tool to describe, predict and control E. coli dynamics in and on food products under various combinations of environmental conditions. To emphasize this point,three specific scenarios are elaborated: (i) aerobic respiration without production of weak acid extracellular metabolites, (ii) anaerobic fermentation with secretion of mixed acid fermentation products into the food environment, and (iii) respiro-fermentative metabolic regimes in between the behaviors at aerobic and anaerobic conditions.

Rapid Classification and Identification of Multiple Microorganisms with Accurate Statistical Significance via High-Resolution Tandem Mass Spectrometry

Rapid and accurate identification and classification of microorganisms is of paramount importance to public health and safety. With the advance of mass spectrometry (MS) technology, the speed of identification can be greatly improved. However, the increasing number of microbes sequenced is complicating correct microbial identification even in a simple sample due to the large number of candidates present. To properly untwine candidate microbes in samples containing one or more microbes, one needs to go beyond apparent morphology or simple "fingerprinting"; to correctly prioritize the candidate microbes, one needs to have accurate statistical significance in microbial identification. We meet these challenges by using peptide-centric representations of microbes to better separate them and by augmenting our earlier analysis method that yields accurate statistical significance. Here, we present an updated analysis workflow that uses tandem MS (MS/MS) spectra for microbial identification or classification. We have demonstrated, using 226 MS/MS publicly available data files (each containing from 2500 to nearly 100,000 MS/MS spectra) and 4000 additional MS/MS data files, that the updated workflow can correctly identify multiple microbes at the genus and often the species level for samples containing more than one microbe. We have also shown that the proposed workflow computes accurate statistical significances, i.e., E values for identified peptides and unified E values for identified microbes. Our updated analysis workflow MiCId, a freely available software for Microorganism Classification and Identification, is available for download at https://www.ncbi.nlm.nih.gov/CBBresearch/Yu/downloads.html . Graphical Abstract ᅟ.

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.

Human Enteric α-Defensin 5 Promotes Shigella Infection by Enhancing Bacterial Adhesion and Invasion

Shigella is a Gram-negative bacterium that causes bacillary dysentery worldwide. It invades the intestinal epithelium to elicit intense inflammation and tissue damage, yet the underlying mechanisms of its host selectivity and low infectious inoculum remain perplexing. Here, we report that Shigella co-opts human α-defensin 5 (HD5), a host defense peptide important for intestinal homeostasis and innate immunity, to enhance its adhesion to and invasion of mucosal tissues. HD5 promoted Shigella infection in vitro in a structure-dependent manner. Shigella, commonly devoid of an effective host-adhesion apparatus, preferentially targeted HD5 to augment its ability to colonize the intestinal epithelium through interactions with multiple bacterial membrane proteins. HD5 exacerbated infectivity and Shigella-induced pathology in a culture of human colorectal tissues and three animal models. Our findings illuminate how Shigella exploits innate immunity by turning HD5 into a virulence factor for infection, unveiling a mechanism of action for this highly proficient human pathogen.

Identification of potential drug targets and inhibitor of the pathogenic bacteria Shigella flexneri 2a through the subtractive genomic approach

Shigella flexneri 2a is one of the most pathogenic bacteria among the Shigella spp., which is responsible for dysentery and causes masses of deaths throughout the world per year. A proper identification of the potential drug targets and inhibitors is crucial for the treatment of the shigellosis due to their emerging multidrug resistance (MDR) patterns. In this study, a systematic subtractive approach was implemented for the identification of novel therapeutic targets of S. flexneri 2a (301) through genome-wide metabolic pathway analysis of the essential genes and proteins. Ligand-based virtual screening and ADMET analyses were also made for the identification of potential inhibitors as well. Initially, we found 70 essential unique proteins as novel targets. After subsequent prioritization, finally we got six unique targets as the potential therapeutic targets and their three-dimensional models were built thereafter. Aspartate-β-semialdehyde dehydrogenase (ASD), was the most potent target among them and used for docking analysis through ligand-based virtual screening. The compound 3 (PubChem CID: 11319750) suited well as the best inhibitor of the ASD through ADMET and enzyme inhibition capacity analysis. To end, we hope that our proposed therapeutic targets and its inhibitors might give some breakthrough to treat shigellosis efficiently in in vitro.

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.

Plasmid-Mediated Quinolone Resistance in Shigella flexneri Isolated From Macaques

Non-human primates (NHPs) for biomedical research are commonly infected with Shigella spp. that can cause acute dysentery or chronic episodic diarrhea. These animals are often prophylactically and clinically treated with quinolone antibiotics to eradicate these possible infections. However, chromosomally- and plasmid-mediated antibiotic resistance has become an emerging concern for species in the family Enterobacteriaceae. In this study, five individual isolates of multi-drug resistant Shigella flexneri were isolated from the feces of three macaques. Antibiotic susceptibility testing confirmed resistance or decreased susceptibility to ampicillin, amoxicillin-clavulanic acid, cephalosporins, gentamicin, tetracycline, ciprofloxacin, enrofloxacin, levofloxacin, and nalidixic acid. S. flexneri isolates were susceptible to trimethoprim-sulfamethoxazole, and this drug was used to eradicate infection in two of the macaques. Plasmid DNA from all isolates was positive for the plasmid-encoded quinolone resistance gene qnrS, but not qnrA and qnrB. Conjugation and transformation of plasmid DNA from several S. flexneri isolates into antibiotic-susceptible Escherichia coli strains conferred the recipients with resistance or decreased susceptibility to quinolones and beta-lactams. Genome sequencing of two representative S. flexneri isolates identified the qnrS gene on a plasmid-like contig. These contigs showed >99% homology to plasmid sequences previously characterized from quinolone-resistant Shigella flexneri 2a and Salmonella enterica strains. Other antibiotic resistance genes and virulence factor genes were also identified in chromosome and plasmid sequences in these genomes. The findings from this study indicate macaques harbor pathogenic S. flexneri strains with chromosomally- and plasmid-encoded antibiotic resistance genes. To our knowledge, this is the first report of plasmid-mediated quinolone resistance in S. flexneri isolated from NHPs and warrants isolation and antibiotic testing of enteric pathogens before treating macaques with quinolones prophylactically or therapeutically.

Virulence and Stress Responses of Shigella flexneri Regulated by PhoP/PhoQ

The two-component signal transduction system PhoP/PhoQ is an important regulator for stress responses and virulence in most Gram-negative bacteria, but characterization of PhoP/PhoQ in Shigella has not been thoroughly investigated. In the present study, we found that deletion of phoPQ (ΔphoPQ) from Shigella flexneri 2a 301 (Sf301) resulted in a significant decline (reduced by more than 15-fold) in invasion of HeLa cells and Caco-2 cells, and less inflammation (− or +) compared to Sf301 (+++) in the guinea pig Sereny test. In low Mg²⁺ (10 μM) medium or pH 5 medium, the ΔphoPQ strain exhibited a growth deficiency compared to Sf301. The ΔphoPQ strain was more sensitive than Sf301 to polymyxin B, an important antimicrobial agent for treating multi-resistant Gram-negative infections. By comparing the transcriptional profiles of ΔphoPQ and Sf301 using DNA microarrays, 117 differentially expressed genes (DEGs) were identified, which were involved in Mg²⁺ transport, lipopolysaccharide modification, acid resistance, bacterial virulence, respiratory, and energy metabolism. Based on the reported PhoP box motif [(T/G) GTTTA-5nt-(T/G) GTTTA], we screened 38 suspected PhoP target operons in S. flexneri, and 11 of them (phoPQ, mgtA, slyB, yoaE, yrbL, icsA, yhiWX, rstA, hdeAB, pagP, and shf–rfbU-virK-msbB2) were demonstrated to be PhoP-regulated genes based on electrophoretic mobility shift assays and β-galactosidase assays. One of these PhoP-regulated genes, icsA, is a well-known virulence factor in S. flexneri. In conclusion, our data suggest that the PhoP/PhoQ system modulates S. flexneri virulence (in an icsA-dependent manner) and stress responses of Mg²⁺, pH and antibacterial peptides.

Shigella and Enteroinvasive Escherichia Coli

Shigella and enteroinvasive Escherichia coli (EIEC) are gram-negative bacteria responsible for bacillary dysentery (shigellosis) in humans, which is characterized by invasion and inflammatory destruction of the human colonic epithelium. Different EIEC and Shigella subgroups rose independently from commensal E. coli through patho-adaptive evolution that included loss of functional genes interfering with the virulence and/or with the intracellular lifestyle of the bacteria, as well as acquisition of genetic elements harboring virulence genes. Among the latter is the large virulence plasmid encoding for a type three secretion system (T3SS), which enables translocation of virulence proteins (effectors) from the bacterium directly into the host cell cytoplasm. These effectors enable the pathogen to subvert epithelial cell functions, promoting its own uptake, replication in the host cytosol, and dissemination to adjacent cells while concomitantly inhibiting pro-inflammatory cell death. Furthermore, T3SS effectors are directly involved in Shigella manipulation of immune cells causing their dysfunction and promoting cell death. In the current chapter, we first describe the evolution of the enteroinvasive pathovars and then summarize the overall knowledge concerning the pathogenesis of these bacteria, with a particular focus on Shigella flexneri. Subversion of host cell functions in the human gut, both epithelial and immune cells, by different virulence factors is especially highlighted.

An H-NS Family Protein, Sfh, Regulates Acid Resistance by Inhibition of Glutamate Decarboxylase Expression in Shigella flexneri 2457T

The glutamate-dependent acid-resistance system is the most effective acid tolerance pathway in Shigella, allowing survival in extremely acidic environments. However, the regulation of this system in Shigella remains elusive. In the current study, we identified significant differences in the levels of glutamate decarboxylase between three Shigella flexneri strains with different levels of acid resistance using blue native-polyacrylamide gel electrophoresis (PAGE) and isoelectric focusing (IEF)/sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) analysis. The results showed that the degree of acid resistance and the levels of GadA/B were significantly lower in strain 2457T compared with two other S. flexneri strains. It has been reported that plasmid pSf-R27 is expressed in strain 2457T but not in the other 142 sequenced S. flexneri isolates. pSf-R27 encodes protein Sfh, which belongs to a family of histone-like nucleoid-structuring (H-NS) proteins that participate in the transcriptional control of glutamate-dependent acid resistance, implicating pSf-R27 in the lower acid resistance of strain 2457T. Transformation of pSf-R27 or sfh alone into strain 301 resulted in decreased expression of GadA/B in the recombinant strains. Thus, we confirmed that H-NS family protein Sfh, bound to the gadA/B regulatory region and regulates the expression of glutamate decarboxylase at the transcriptional level. We also examined the acid tolerance of the wild-type and recombinant strains using flow cytometry and determined that the acid tolerance of S. flexneri is closely related to the expression of GadA/B. These findings further our understanding of the acid tolerance of S. flexneri, especially via the glutamate-dependent pathway.

Genetic plasticity of the Shigella virulence plasmid is mediated by intra- and inter-molecular events between insertion sequences

Acquisition of a single copy, large virulence plasmid, pINV, led to the emergence of Shigella spp. from Escherichia coli. The plasmid encodes a Type III secretion system (T3SS) on a 30 kb pathogenicity island (PAI), and is maintained in a bacterial population through a series of toxin:antitoxin (TA) systems which mediate post-segregational killing (PSK). The T3SS imposes a significant cost on the bacterium, and strains which have lost the plasmid and/or genes encoding the T3SS grow faster than wild-type strains in the laboratory, and fail to bind the indicator dye Congo Red (CR). Our aim was to define the molecular events in Shigella flexneri that cause loss of Type III secretion (T3S), and to examine whether TA systems exert positional effects on pINV. During growth at 37°C, we found that deletions of regions of the plasmid including the PAI lead to the emergence of CR-negative colonies; deletions occur through intra-molecular recombination events between insertion sequences (ISs) flanking the PAI. Furthermore, by repositioning MvpAT (which belongs to the VapBC family of TA systems) near the PAI, we demonstrate that the location of this TA system alters the rearrangements that lead to loss of T3S, indicating that MvpAT acts both globally (by reducing loss of pINV through PSK) as well as locally (by preventing loss of adjacent sequences). During growth at environmental temperatures, we show for the first time that pINV spontaneously integrates into different sites in the chromosome, and this is mediated by inter-molecular events involving IS1294. Integration leads to reduced PAI gene expression and impaired secretion through the T3SS, while excision of pINV from the chromosome restores T3SS function. Therefore, pINV integration provides a reversible mechanism for Shigella to circumvent the metabolic burden imposed by pINV. Intra- and inter-molecular events between ISs, which are abundant in Shigella spp., mediate plasticity of S. flexneri pINV.