Conversion of Shigella flexneri serotype 2a to serotype Y in a shigellosis patient due to a single amino acid substitution in the protein product of the bacterial glucosyltransferase gtrII gene

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.

Influence of Shigella flexneri 2a O Antigen Acetylation on Its Bacteriophage Sf6 Receptor Activity and Bacterial Interaction with Human Cells

Shigella flexneri is a major causative agent of bacillary dysentery in developing countries, where serotype 2a₂ is the prevalent strain. To date, approximately 30 serotypes have been identified for S. flexneri, and the major contribution to the emergence of new serotypes is chemical modifications of the lipopolysaccharide (LPS) component O antigen (Oag). Glucosylation, O-acetylation, and phosphoethanolamine (PEtN) modifications increase the Oag diversity, providing benefits to S. flexneri LPS Oag acts as a primary receptor for bacteriophage Sf6, which infects only a limited range of S. flexneri serotypes (Y and X). It uses its tailspike protein (Sf6TSP) to establish initial interaction with LPS Oags that it then hydrolyzes. Currently, there is a lack of comprehensive study on the parent and serotype variant strains from the same genetic background and an understanding of the importance of LPS Oag O-acetylations. Therefore, a set of isogenic strains (based on S. flexneri 2457T [2a₂]) with deletions of different Oag modification genes (oacB, oacD, and gtrII) that resemble different naturally occurring serotype Y and 2a strains was created. The impacts of these Oag modifications on S. flexneri sensitivity to Sf6 and the pathogenesis-related properties were then compared. We found that Sf6TSP can hydrolyze serotype 2a LPS Oag, identified that 3/4-O-acetylation is essential for resistance of serotype 2a strains to Sf6, and showed that serotype 2a strains have better invasion ability. Lastly, we revealed two new serotype conversions for S. flexneri, thereby contributing to understanding the evolution of this important human pathogen.IMPORTANCE The emergence of antibiotic-resistant strains and lack of efficient vaccines have made Shigella a priority organism for the World Health Organization (1). Therefore, bacteriophage therapy has received increasing attention as an alternative therapeutic approach. LPS Oag is the most variable part of LPS due to chemical modifications and is the target of bacteriophage Sf6 (S. flexneri specific). We dissected the evolution of S. flexneri serotype Y to 2a₂, which revealed a new role for a gene acquired during serotype conversion and furthermore identified new specific forms of LPS receptor for Sf6. Collectively, these results unfold the importance of the acquisition of those Oag modification genes and further our understanding of the relationship between Sf6 and S. flexneri.

Contamination of Currency Notes with Kanamycin ResistantShigella flexneri.. [DOI: 10.1101/2020.03.07.982017]

Shigella flexneriis the main causative agent of shigellosis commonly distributed in developing countries with high morbidity and mortality rates. This study aimed to examine the presence ofShigellaspecies in Sudanese currency notes using both traditional and molecular techniques. One hundred thirty five currency notes were collected and their contaminants were isolated and identified conventionally and genetically using 16S rRNA gene amplification and sequencing. Eight isolates were identified asShigellaspecies in different notes, and 3 of them were resistant to penicillin, kanamycin and nitrofurantoin. OneS. flexneriisolate has insertion mutation of guanine nucleotide at position 730 of life’s essential gene 16S rRNA which known evolutionarily to be stable gene. Banknotes are highly circulating items and therefore, appropriate measures such as regular replacement of the dirty notes with new papers are necessary to protect peoples from being infected with drug resistant pathogens.

Association of Anti-GT1a Antibodies with an Outbreak of Guillain-Barré Syndrome and Analysis of Ganglioside Mimicry in an Associated Campylobacter jejuni Strain

An outbreak of Guillain-Barré syndrome (GBS), subsequent to Campylobacter jejuni enteritis, occurred in China in 2007. Serum anti-ganglioside antibodies were measured in GBS patients and controls. Genome sequencing was used to determine the phylogenetic relationship among three C. jejuni strains from a patient with GBS (ICDCCJ07001), a patient with gastroenteritis (ICDCCJ07002) and a healthy carrier (ICDCCJ07004), which were all associated with the outbreak. The ganglioside-like structures of the lipo-oligosaccharides of these strains were determined by mass spectrometry. Seventeen (53%) of the GBS patients had anti-GT1a IgG antibodies. GT1a mimicry was found in the lipo-oligosaccharides of strain ICDCCJ07002 and ICDCCJ07004; but a combination of GM3/GD3 mimics was observed in ICDCCJ07001, although this patient had anti-GT1a IgG antibodies. A single-base deletion in a glycosyltransferase gene caused the absence of GT1a mimicry in ICDCCJ07001. The phylogenetic tree showed that ICDCCJ07002 and ICDCCJ07004 were genetically closer to each other than to ICDCCJ07001. C. jejuni, bearing a GT1a-like lipo-oligosaccharide, might have caused the GBS outbreak and the loss of GT1a mimicry may have helped ICDCCJ07001 to survive in the host.

Identification and characterization of a novel Shigella flexneri serotype Yv in China

Shigella flexneri is the major cause of bacterial shigellosis in developing countries. S. flexneri is divided into at least 19 serotypes, the majority of which are modifications of the same basic O-antigen by glucosylation and/or O-acetylation of its sugar residues by phage encoded serotype-converting genes. Recently, a plasmid encoded phosphoethanolamine (PEtN) modification of the O-antigen has been reported, which is responsible for the presence of the MASF IV-1 determinant and results in conversion of traditional serotypes X, 4a and Y to novel serotypes Xv, 4av and Yv, respectively. In this study, we characterized 19 serotype Yv strains isolated in China. A variant of the O-antigen phosphoethanolamine transferase gene opt (formerly called lpt-O) carried by a pSFxv_2-like plasmid was found in serotype Yv strains, which specifies the phosphorylation pattern on the O-antigen of this serotype. For the majority of the O-antigen units, the PEtN modification occurs on Rha^III, while for a minority, modifications occur on both Rha^II and Rha^III. Serotype-specific gene detection and PFGE analysis suggested that these serotype Yv isolates were originated from serotypes Y, Xv and 2a by acquisition of an opt-carrying plasmid and/or inactivation of serotype-specific gene gtrII or gtrX. These data, combined with those of serotypes Xv and 4av reported earlier, demonstrate that the plasmid-encoded PEtN modification is an important serotype conversion mechanism in S. flexneri, in addition to glucosylation and O-acetylation.

Development of a multiplex PCR assay targeting O-antigen modification genes for molecular serotyping of Shigella flexneri

Shigella flexneri is the major Shigella species that causes diarrheal disease in developing countries. It is further subdivided into 15 serotypes based on O-antigen structure. Serotyping of S. flexneri is important for epidemiological purposes. In this study, we developed a multiplex PCR assay targeting the O-antigen synthesis gene wzx and the O-antigen modification genes gtrI, gtrIC, gtrII, oac, gtrIV, gtrV, and gtrX for molecular serotyping of S. flexneri. The multiplex PCR assay contained eight sets of specific PCRs in a single tube and can identify 14 of the 15 serotypes (the exception being serotype Xv) of S. flexneri recognized thus far. A nearly perfect concordance (97.8%) between multiplex PCR assay and slide agglutination was observed when 358 S. flexneri strains of various serotypes were analyzed, except that 8 strains were carrying additional cryptic and/or defective serotype-specific genes. The multiplex PCR assay provides a rapid and specific method for the serotype identification of S. flexneri.

Multilocus variable-number tandem repeat analysis for molecular typing and phylogenetic analysis of Shigella flexneri

Background

Shigella flexneri is one of the causative agents of shigellosis, a major cause of childhood mortality in developing countries. Multilocus variable-number tandem repeat (VNTR) analysis (MLVA) is a prominent subtyping method to resolve closely related bacterial isolates for investigation of disease outbreaks and provide information for establishing phylogenetic patterns among isolates. The present study aimed to develop an MLVA method for S. flexneri and the VNTR loci identified were tested on 242 S. flexneri isolates to evaluate their variability in various serotypes. The isolates were also analyzed by pulsed-field gel electrophoresis (PFGE) to compare the discriminatory power and to evaluate the usefulness of MLVA as a tool for phylogenetic analysis of S. flexneri.

Results

Thirty-six VNTR loci were identified by exploring the repeat sequence loci in genomic sequences of Shigella species and by testing the loci on nine isolates of different subserotypes. The VNTR loci in different serotype groups differed greatly in their variability. The discriminatory power of an MLVA assay based on four most variable VNTR loci was higher, though not significantly, than PFGE for the total isolates, a panel of 2a isolates, which were relatively diverse, and a panel of 4a/Y isolates, which were closely-related. Phylogenetic groupings based on PFGE patterns and MLVA profiles were considerably concordant. The genetic relationships among the isolates were correlated with serotypes. The phylogenetic trees constructed using PFGE patterns and MLVA profiles presented two distinct clusters for the isolates of serotype 3 and one distinct cluster for each of the serotype groups, 1a/1b/NT, 2a/2b/X/NT, 4a/Y, and 6. Isolates that had different serotypes but had closer genetic relatedness than those with the same serotype were observed between serotype Y and subserotype 4a, serotype X and subserotype 2b, subserotype 1a and 1b, and subserotype 3a and 3b.

Conclusions

The 36 VNTR loci identified exhibited considerably different degrees of variability among S. flexneri serotype groups. VNTR locus could be highly variable in a serotype but invariable in others. MLVA assay based on four highly variable loci could display a comparable resolving power to PFGE in discriminating isolates. MLVA is also a prominent molecular tool for phylogenetic analysis of S. flexneri; the resulting data are beneficial to establish clear clonal patterns among different serotype groups and to discern clonal groups among isolates within the same serotype. As highly variable VNTR loci could be serotype-specific, a common MLVA protocol that consists of only a small set of loci, for example four to eight loci, and that provides high resolving power to all S. flexneri serotypes may not be obtainable.

Bacteriophage-encoded glucosyltransferase GtrII of Shigella flexneri: membrane topology and identification of critical residues

The Shigella flexneri serotypes differ in the nature of their O-antigens. The addition of glucosyl or O-acetyl groups to the common backbone repeat units gives rise to the different serotypes. GtrII glucosylates rhamnose III of the O-antigen repeat unit, thus converting serotype Y (which has no modifications to the basic O-antigen repeat unit) into serotype 2a, the most prevalent serotype. In the present study, the topology of GtrII has been determined. GtrII has nine transmembrane helices, a re-entrant loop and three large periplasmic regions. Four critical residues (Glu40, Phe414, Cys435 and Lys478) were identified in two of the periplasmic regions. Despite the lack of sequence similarity between GtrII and the Gtrs from other serotypes, three of the critical residues identified are conserved in the remaining Gtrs. This is consistent with some degree of mechanistic conservation in this functionally related group of proteins.

The Shigella flexneri serotype Y vaccine candidate SFL124 originated from a serotype 2a background

Shigella flexneri is endemic in most developing countries and responsible for the highest mortality rate among the Shigella species. The attenuated serotype Y S. flexneri strain SFL124 has been used as the parental strain for the development of recombinant vaccines expressing multiple O-antigen structures. During the development of one such multivalent vaccine, a region of gtrII homology was found in SFL124. Sequencing and analysis of this region revealed the presence of an insertion element interrupted serotype 2a serotype-conversion locus in the serotype Y vaccine strain SFL124. The data presented suggests that SFL124 has derived from a serotype 2a background.

Single-molecule approach to bacterial genomic comparisons via optical mapping

Modern comparative genomics has been established, in part, by the sequencing and annotation of a broad range of microbial species. To gain further insights, new sequencing efforts are now dealing with the variety of strains or isolates that gives a species definition and range; however, this number vastly outstrips our ability to sequence them. Given the availability of a large number of microbial species, new whole genome approaches must be developed to fully leverage this information at the level of strain diversity that maximize discovery. Here, we describe how optical mapping, a single-molecule system, was used to identify and annotate chromosomal alterations between bacterial strains represented by several species. Since whole-genome optical maps are ordered restriction maps, sequenced strains of Shigella flexneri serotype 2a (2457T and 301), Yersinia pestis (CO 92 and KIM), and Escherichia coli were aligned as maps to identify regions of homology and to further characterize them as possible insertions, deletions, inversions, or translocations. Importantly, an unsequenced Shigella flexneri strain (serotype Y strain AMC[328Y]) was optically mapped and aligned with two sequenced ones to reveal one novel locus implicated in serotype conversion and several other loci containing insertion sequence elements or phage-related gene insertions. Our results suggest that genomic rearrangements and chromosomal breakpoints are readily identified and annotated against a prototypic sequenced strain by using the tools of optical mapping.