Cloning and analysis of the glucosyl transferase gene encoding type I antigen in Shigella flexneri

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.

The genomic signatures of Shigella evolution, adaptation and geographical spread

Shigella spp. are some of the key pathogens responsible for the global burden of diarrhoeal disease. These facultative intracellular bacteria belong to the family Enterobacteriaceae, together with other intestinal pathogens, such as Escherichia coli and Salmonella spp. The genus Shigella comprises four different species, each consisting of several serogroups, all of which show phenotypic similarity, including invasive pathogenicity. DNA sequencing suggests that this similarity results from the convergent evolution of different Shigella spp. founders. Here, we review the evolutionary relationships between Shigella spp. and E . coli, and we highlight how the genomic plasticity of these bacteria and their acquisition of a distinctive virulence plasmid have enabled the development of such highly specialized pathogens. Furthermore, we discuss the insights that genotyping and whole-genome sequencing have provided into the phylogenetics and intercontinental spread of Shigella spp.

Isolation and genomic characterization of SfI, a serotype-converting bacteriophage of Shigella flexneri

Background

All Shigella flexneri serotypes except serotype 6 share a common O-antigen tetrasaccharide backbone and nearly all variations between serotypes are due to glucosyl and/or O-acetyl modifications of the common O unit mediated by glycosyltransferases encoded by serotype-converting bacteriophages. Several S. flexneri serotype-converting phages including SfV, SfX, Sf6 and SfII have been isolated and characterized. However, S. flexneri serotype-converting phage SfI which encodes a type I modification of serotype 1 (1a, 1b, 1c and 1d) had not yet been characterized.

Results

The SfI phage was induced and purified from a S. flexneri serotype 1a clinical strain 019. Electron microscopy showed that the SfI phage has a hexagonal head and a long contractile tail, characteristic of the members of Myoviridae family. SfI can convert serotype Y to serotype 1a and serotype X to serotype 1d, but cannot convert 10 other S. flexneri serotypes (1a, 1b, 2a, 2b, 3a, 3b, 4a, 4b, 5a, Xv) tested, suggesting that SfI has a narrow host range. Similar to other S. flexneri serotype-converting phages, SfI integrates into the tRNA-thrW gene adjacent to proA of the host chromosome when lysogenized. The complete sequence of the SfI genome was 38,389 bp, encoding 66 open reading frames and two tRNA genes. Phage SfI shares significant homology with S. flexneri phage SfV, Escherichia coli prophage e14 and lambda, and is classified into the lambdoid phage family. SfI was found to use a cos mechanism for DNA packaging similar to that of phage SfV.

Conclusions

SfI contains features of lambdoid phages and is closely related to S. flexneri phage SfV, E. coli prophage e14 and lambda. The characterization of SfI enhances our understanding of serotype conversion of S. flexneri.

Comparative genomic analysis of mycobacteriophage Tweety: evolutionary insights and construction of compatible site-specific integration vectors for mycobacteria

Mycobacteriophage Tweety is a newly isolated phage of Mycobacterium smegmatis. It has a viral morphology with an isometric head and a long flexible tail, and forms turbid plaques from which stable lysogens can be isolated. The Tweety genome is 58 692 bp in length, contains 109 protein-coding genes, and shows significant but interrupted nucleotide sequence similarity with the previously described mycobacteriophages Llij, PMC and Che8. However, overall the genome possesses mosaic architecture, with gene products being related to other mycobacteriophages such as Che9d, Omega and Corndog. A gene encoding an integrase of the tyrosine-recombinase family is located close to the centre of the genome, and a putative attP site has been identified within a short intergenic region immediately upstream of int. This Tweety attP–int cassette was used to construct a new set of integration-proficient plasmid vectors that efficiently transform both fast- and slow-growing mycobacteria through plasmid integration at a chromosomal locus containing a tRNA^Lys gene. These vectors are maintained well in the absence of selection and are completely compatible with integration vectors derived from mycobacteriophage L5, enabling the simple construction of complex recombinants with genes integrated simultaneously at different chromosomal positions.

Coexistence of two distinct versions of O-antigen polymerase, Wzy-alpha and Wzy-beta, in Pseudomonas aeruginosa serogroup O2 and their contributions to cell surface diversity

Assembly of B-band lipopolysaccharide (LPS) in Pseudomonas aeruginosa follows a Wzy-dependent pathway, requiring the O-antigen polymerase Wzy and other proteins. The peptide sequences of the wzy(alpha) product from strains of serotypes O2, O5, and O16 are identical, but the O units in O5 are alpha-glycosidically linked, while those in O2 and O16 are beta-linked. We hypothesized that a derivative of the D3 bacteriophage wzy(beta) is present in the chromosomes of O2 and O16 and that this gene is responsible for the beta-linkage. By a combination of PCR and primer walking, wzy(beta) genes of both serotypes have been amplified and cloned. They are identical but share only 87.42% sequence identity with their xenolog in D3. A chromosomal knockout mutant of O16 wzy(beta) was made, and it produces semirough LPS devoid of B-band O antigen. The cloned wzy(beta) is capable of complementing the O16 wzy(beta) mutant, as well as cross-complementing a wzy(alpha) knockout mutant. However, in the latter case, the restored O antigen was beta-linked. Using reverse transcription-PCR, we showed that wzy(alpha) was transcribed in O2 and O16 strains and was functional, since both of these genes could complement the wzy(alpha) mutant of O5. With the coexistence of wzy(alpha) and wzy(beta) in O2 and O16 and the B-band O polysaccharides in these being beta-linked, we hypothesized that iap, an inhibitor of the alpha-polymerase gene, must be present in these serotypes. Indeed, through PCR, TOPO-cloning, and nucleotide-sequencing results, we verified the presence of iap in both O2 and O16 serotypes.

Construction of a multivalent vaccine strain ofShigella flexneriand evaluation of serotype-specific immunity

Shigella flexneri causes more fatalities by shigellosis than any other Shigella species. There are 13 different serotypes of S. flexneri and their distribution varies between endemic geographical regions. The immune response against S. flexneri is serotype-specific, so current immunization strategies have required the administration of multiple vaccine strains to provide protection against multiple serotypes. In this study, we report the construction of a multivalent S. flexneri vaccine strain, SFL1425, expressing the O-antigen structure specific for serotypes 2a and 5a. This combination of type antigens has not previously been reported for S. flexneri. The multivalent vaccine strain, SFL1425 was able to induce a specific immune response against both serotypes 2a and 5a in a mouse pulmonary model.

Complete genomic sequence of SfV, a serotype-converting temperate bacteriophage of Shigella flexneri

Bacteriophage SfV is a temperate serotype-converting phage of Shigella flexneri. SfV encodes the factors involved in type V O-antigen modification, and the serotype conversion and integration-excision modules of the phage have been isolated and characterized. We now report on the complete sequence of the SfV genome (37,074 bp). A total of 53 open reading frames were predicted from the nucleotide sequence, and analysis of the corresponding proteins was used to construct a functional map. The general organization of the genes in the SfV genome is similar to that of bacteriophage lambda, and numerous features of the sequence are described. The superinfection immunity system of SfV includes a lambda-like repression system and a P4-like transcription termination mechanism. Sequence analysis also suggests that SfV encodes multiple DNA methylases, and experiments confirmed that orf-41 encodes a Dam methylase. Studies conducted to determine if the phage-encoded methylase confers host DNA methylation showed that the two S. flexneri strains analyzed encode their own Dam methylase. Restriction mapping and sequence analysis revealed that the phage genome has cos sites at the termini. The tail assembly and structural genes of SfV show homology to those of phage Mu and Mu-like prophages in the genome of Escherichia coli O157:H7 and Haemophilus influenzae. Significant homology (30% of the genome in total) between sections of the early, regulatory, and structural regions of the SfV genome and the e14 and KpLE1 prophages in the E. coli K-12 genome were noted, suggesting that these three phages have common evolutionary origins.