The evolution of the pan-genome of Shiga-toxin (Stx) producing Escherichia coli and the Stx₂ bacteriophage

xvi, 206 leaves : ill. (some col.) ; 29 cmHuman infections with Shiga-toxin (Stx)-producing E. coli (STEC) vary in severity of illness. The pan-genome of a bacterial species contains a shared, essential core genome, and a variably distributed accessory genome. While single nucleotide changes likely influence virulence in STEC, horizontal gene transfer (HGT) on elements such as bacteriophage are thought to be most important. My thesis objectives were to: 1) develop tools for the pangenomic analyses of bacterial genomes; 2) describe the phylogeny of STEC and; 3) determine if the evolution of the Stx2-bacteriophage parallels that of its bacterial host. For this thesis, the software program Panseq was created and used to identify pan-genomic differences among STEC. Whole-genome phylogenies showed all serotypes as discrete clusters, with O157:H7 having three distinct lineages and grouping separately from all other STEC. Finally, the phylogenies of Stx2-bacteriophage and their bacterial hosts were largely concordant, with occasional instances of HGT having led to novel pathogen emergence

In silico genomic analyses reveal three distinct lineages of Escherichia coli O157:H7, one of which is associated with hyper-virulence

<p>Abstract</p> <p>Background</p> <p>Many approaches have been used to study the evolution, population structure and genetic diversity of <it>Escherichia coli </it>O157:H7; however, observations made with different genotyping systems are not easily relatable to each other. Three genetic lineages of <it>E. coli </it>O157:H7 designated I, II and I/II have been identified using octamer-based genome scanning and microarray comparative genomic hybridization (mCGH). Each lineage contains significant phenotypic differences, with lineage I strains being the most commonly associated with human infections. Similarly, a clade of hyper-virulent O157:H7 strains implicated in the 2006 spinach and lettuce outbreaks has been defined using single-nucleotide polymorphism (SNP) typing. In this study an <it>in silico </it>comparison of six different genotyping approaches was performed on 19 <it>E. coli </it>genome sequences from 17 O157:H7 strains and single O145:NM and K12 MG1655 strains to provide an overall picture of diversity of the <it>E. coli </it>O157:H7 population, and to compare genotyping methods for O157:H7 strains.</p> <p>Results</p> <p><it>In silico </it>determination of lineage, Shiga-toxin bacteriophage integration site, comparative genomic fingerprint, mCGH profile, novel region distribution profile, SNP type and multi-locus variable number tandem repeat analysis type was performed and a supernetwork based on the combination of these methods was produced. This supernetwork showed three distinct clusters of strains that were O157:H7 lineage-specific, with the SNP-based hyper-virulent clade 8 synonymous with O157:H7 lineage I/II. Lineage I/II/clade 8 strains clustered closest on the supernetwork to <it>E. coli </it>K12 and <it>E. coli </it>O55:H7, O145:NM and sorbitol-fermenting O157 strains.</p> <p>Conclusion</p> <p>The results of this study highlight the similarities in relationships derived from multi-locus genome sampling methods and suggest a "common genotyping language" may be devised for population genetics and epidemiological studies. Future genotyping methods should provide data that can be stored centrally and accessed locally in an easily transferable, informative and extensible format based on comparative genomic analyses.</p

Multi-Year Persistence of Verotoxigenic Escherichia coli (VTEC) in a Closed Canadian Beef Herd: A Cohort Study

In this study, fecal samples were collected from a closed beef herd in Alberta, Canada from 2012 to 2015. To limit serotype bias, which was observed in enrichment broth cultures, Verotoxigenic Escherichia coli (VTEC) were isolated directly from samples using a hydrophobic grid-membrane filter verotoxin immunoblot assay. Overall VTEC isolation rates were similar for three different cohorts of yearling heifers on both an annual (68.5 to 71.8%) and seasonal basis (67.3 to 76.0%). Across all three cohorts, O139:H19 (37.1% of VTEC-positive samples), O22:H8 (15.8%) and O?(O108):H8 (15.4%) were among the most prevalent serotypes. However, isolation rates for serotypes O139:H19, O130:H38, O6:H34, O91:H21, and O113:H21 differed significantly between cohort-years, as did isolation rates for some serotypes within a single heifer cohort. There was a high level of VTEC serotype diversity with an average of 4.3 serotypes isolated per heifer and 65.8% of the heifers classified as “persistent shedders” of VTEC based on the criteria of &gt;50% of samples positive and ≥4 consecutive samples positive. Only 26.8% (90/336) of the VTEC isolates from yearling heifers belonged to the human disease-associated seropathotypes A (O157:H7), B (O26:H11, O111:NM), and C (O22:H8, O91:H21, O113:H21, O137:H41, O2:H6). Conversely, seropathotypes B (O26:NM, O111:NM) and C (O91:H21, O2:H29) strains were dominant (76.0%, 19/25) among VTEC isolates from month-old calves from this herd. Among VTEC from heifers, carriage rates of vt1, vt2, vt1+vt2, eae, and hlyA were 10.7, 20.8, 68.5, 3.9, and 88.7%, respectively. The adhesin gene saa was present in 82.7% of heifer strains but absent from all of 13 eae+ve strains (from serotypes/intimin types O157:H7/γ1, O26:H11/β1, O111:NM/θ, O84:H2/ζ, and O182:H25/ζ). Phylogenetic relationships inferred from wgMLST and pan genome-derived core SNP analysis showed that strains clustered by phylotype and serotype. Further, VTEC strains of the same serotype usually shared the same suite of antibiotic resistance and virulence genes, suggesting the circulation of dominant clones within this distinct herd. This study provides insight into the diverse and dynamic nature of VTEC populations within groups of cattle and points to a broad spectrum of human health risks associated with these E. coli strains

Comparative genomic analysis of Escherichia coli 0157:H7 isolated from super-shedder and low-shedder cattle

Sherpa Romeo green journal; open accessCattle are the primary reservoir of the foodborne pathogen Escherichia coli O157:H7, with the concentration and frequency of E. coli O157:H7 shedding varying substantially among individual hosts. The term ‘‘super-shedder” has been applied to cattle that shed 104 cfu E. coli O157:H7/g of feces. Super-shedders have been reported to be responsible for the majority of E. coli O157:H7 shed into the environment. The objective of this study was to determine if there are phenotypic and/or genotypic differences between E. coli O157:H7 isolates obtained from super-shedder compared to low-shedder cattle. From a total of 784 isolates, four were selected from low-shedder steers and six isolates from super-shedder steers (4.01–8.45 log cfu/g feces) for whole genome sequencing. Isolates were phage and clade typed, screened for substrate utilization, pH sensitivity, virulence gene profiles and Stx bacteriophage insertion (SBI) sites. A range of 89–2473 total single nucleotide polymorphisms (SNPs) were identified when sequenced strains were compared to E. coli O157:H7 strain Sakai. More non-synonymous SNP mutations were observed in low-shedder isolates. Pan-genomic and SNPs comparisons did not identify genetic segregation between supershedder or low-shedder isolates. All super-shedder isolates and 3 of 4 of low-shedder isolates were typed as phage type 14a, SBI cluster 3 and SNP clade 2. Super-shedder isolates displayed increased utilization of galactitol, thymidine and 3-O-β-D-galactopyranosyl-Darabinose when compared to low-shedder isolates, but no differences in SNPs were observed in genes encoding for proteins involved in the metabolism of these substrates. While genetic traits specific to super-shedder isolates were not identified in this study, differences in the level of gene expression or genes of unknown function may still contribute to some strains of E. coli 0157:H7 reaching high densities within bovine feces.Ye

A Comparison of Shiga-Toxin 2 Bacteriophage from Classical Enterohemorrhagic Escherichia coli Serotypes and the German E. coli O104:H4 Outbreak Strain

Escherichia coli O104:H4 was associated with a severe foodborne disease outbreak originating in Germany in May 2011. More than 4000 illnesses and 50 deaths were reported. The outbreak strain was a typical enteroaggregative E. coli (EAEC) that acquired an antibiotic resistance plasmid and a Shiga-toxin 2 (Stx2)-encoding bacteriophage. Based on whole-genome phylogenies, the O104:H4 strain was most closely related to other EAEC strains; however, Stx2-bacteriophage are mobile, and do not necessarily share an evolutionary history with their bacterial host. In this study, we analyzed Stx2-bacteriophage from the E. coli O104:H4 outbreak isolates and compared them to all available Stx2-bacteriophage sequences. We also compared Stx2 production by an E. coli O104:H4 outbreak-associated isolate (ON-2011) to that of E. coli O157:H7 strains EDL933 and Sakai. Among the E. coli Stx2-phage sequences studied, that from O111:H- strain JB1-95 was most closely related phylogenetically to the Stx2-phage from the O104:H4 outbreak isolates. The phylogeny of most other Stx2-phage was largely concordant with their bacterial host genomes. Finally, O104:H4 strain ON-2011 produced less Stx2 than E. coli O157:H7 strains EDL933 and Sakai in culture; however, when mitomycin C was added, ON-2011 produced significantly more toxin than the E. coli O157:H7 strains. The Stx2-phage from the E. coli O104:H4 outbreak strain and the Stx2-phage from O111:H- strain JB1-95 likely share a common ancestor. Incongruence between the phylogenies of the Stx2-phage and their host genomes suggest the recent Stx2-phage acquisition by E. coli O104:H4. The increase in Stx2-production by ON-2011 following mitomycin C treatment may or may not be related to the high rates of hemolytic uremic syndrome associated with the German outbreak strain. Further studies are required to determine whether the elevated Stx2-production levels are due to bacteriophage or E. coli O104:H4 host related factors

Finishing the euchromatic sequence of the human genome

The sequence of the human genome encodes the genetic instructions for human physiology, as well as rich information about human evolution. In 2001, the International Human Genome Sequencing Consortium reported a draft sequence of the euchromatic portion of the human genome. Since then, the international collaboration has worked to convert this draft into a genome sequence with high accuracy and nearly complete coverage. Here, we report the result of this finishing process. The current genome sequence (Build 35) contains 2.85 billion nucleotides interrupted by only 341 gaps. It covers ∼99% of the euchromatic genome and is accurate to an error rate of ∼1 event per 100,000 bases. Many of the remaining euchromatic gaps are associated with segmental duplications and will require focused work with new methods. The near-complete sequence, the first for a vertebrate, greatly improves the precision of biological analyses of the human genome including studies of gene number, birth and death. Notably, the human enome seems to encode only 20,000-25,000 protein-coding genes. The genome sequence reported here should serve as a firm foundation for biomedical research in the decades ahead

PANC Study (Pancreatitis: A National Cohort Study): national cohort study examining the first 30 days from presentation of acute pancreatitis in the UK

Abstract Background Acute pancreatitis is a common, yet complex, emergency surgical presentation. Multiple guidelines exist and management can vary significantly. The aim of this first UK, multicentre, prospective cohort study was to assess the variation in management of acute pancreatitis to guide resource planning and optimize treatment. Methods All patients aged greater than or equal to 18 years presenting with acute pancreatitis, as per the Atlanta criteria, from March to April 2021 were eligible for inclusion and followed up for 30 days. Anonymized data were uploaded to a secure electronic database in line with local governance approvals. Results A total of 113 hospitals contributed data on 2580 patients, with an equal sex distribution and a mean age of 57 years. The aetiology was gallstones in 50.6 per cent, with idiopathic the next most common (22.4 per cent). In addition to the 7.6 per cent with a diagnosis of chronic pancreatitis, 20.1 per cent of patients had a previous episode of acute pancreatitis. One in 20 patients were classed as having severe pancreatitis, as per the Atlanta criteria. The overall mortality rate was 2.3 per cent at 30 days, but rose to one in three in the severe group. Predictors of death included male sex, increased age, and frailty; previous acute pancreatitis and gallstones as aetiologies were protective. Smoking status and body mass index did not affect death. Conclusion Most patients presenting with acute pancreatitis have a mild, self-limiting disease. Rates of patients with idiopathic pancreatitis are high. Recurrent attacks of pancreatitis are common, but are likely to have reduced risk of death on subsequent admissions. </jats:sec

Pan-genome Analyses of the Species Salmonella enterica, and Identification of Genomic Markers Predictive for Species, Subspecies, and Serovar

Food safety is a global concern, with upward of 2.2 million deaths due to enteric disease every year. Current whole-genome sequencing platforms allow routine sequencing of enteric pathogens for surveillance, and during outbreaks; however, a remaining challenge is the identification of genomic markers that are predictive of strain groups that pose the most significant health threats to humans, or that can persist in specific environments. We have previously developed the software program Panseq, which identifies the pan-genome among a group of sequences, and the SuperPhy platform, which utilizes this pan-genome information to identify biomarkers that are predictive of groups of bacterial strains. In this study, we examined the pan-genome of 4893 genomes of Salmonella enterica, an enteric pathogen responsible for the loss of more disability adjusted life years than any other enteric pathogen. We identified a pan-genome of 25.3 Mbp, a strict core of 1.5 Mbp present in all genomes, and a conserved core of 3.2 Mbp found in at least 96% of these genomes. We also identified 404 genomic regions of 1000 bp that were specific to the species S. enterica. These species-specific regions were found to encode mostly hypothetical proteins, effectors, and other proteins related to virulence. For each of the six S. enterica subspecies, markers unique to each were identified. No serovar had pan-genome regions that were present in all of its genomes and absent in all other serovars; however, each serovar did have genomic regions that were universally present among all constituent members, and statistically predictive of the serovar. The phylogeny based on SNPs within the conserved core genome was found to be highly concordant to that produced by a phylogeny using the presence/absence of 1000 bp regions of the entire pan-genome. Future studies could use these predictive regions as components of a vaccine to prevent salmonellosis, as well as in simple and rapid diagnostic tests for both in silico and wet-lab applications, with uses ranging from food safety to public health. Lastly, the tools and methods described in this study could be applied as a pan-genomics framework to other population genomic studies seeking to identify markers for other bacterial species and their sub-groups

Characteristics of Clinical Shiga Toxin-Producing Escherichia coli Isolated from British Columbia

Shiga toxin-producing Escherichia coli (STEC) are significant public health threats. Although STEC O157 are recognized foodborne pathogens, non-O157 STEC are also important causes of human disease. We characterized 10 O157:H7 and 15 non-O157 clinical STEC derived from British Columbia (BC). Eae, hlyA, and stx were more frequently observed in STEC O157, and 80 and 100% of isolates possessed stx1 and stx2, respectively. In contrast, stx1 and stx2 occurred in 80 and 40% of non-O157 STEC, respectively. Comparative genomic fingerprinting (CGF) revealed three distinct clusters (C). STEC O157 was identified as lineage I (LI; LSPA-6 111111) and clustered as a single group (C1). The cdi gene previously observed only in LII was seen in two LI O157 isolates. CGF C2 strains consisted of diverse non-O157 STEC while C3 included only O103:H25, O118, and O165 serogroup isolates. With the exception of O121 and O165 isolates which were similar in virulence gene complement to STEC O157, C1 O157 STEC produced more Stx2 than non-O157 STEC. Antimicrobial resistance (AMR) screening revealed resistance or reduced sensitivity in all strains, with higher levels occurring in non-O157 STEC. One STEC O157 isolate possessed a mobile blaCMY-2 gene transferrable across genre via conjugation