James E. Thomas

Pan-genome sequence analysis using Panseq: an online tool for the rapid analysis of core and accessory genomic regions

BackgroundThe pan-genome of a bacterial species consists of a core and an accessory gene pool. The accessory genome is thought to be an important source of genetic variability in bacterial populations and is gained through lateral gene transfer, allowing subpopulations of bacteria to better adapt to specific niches. Low-cost and high-throughput sequencing platforms have created an exponential increase in genome sequence data and an opportunity to study the pan-genomes of many bacterial species. In this study, we describe a new online pan-genome sequence analysis program, Panseq.ResultsPanseq was used to identify Escherichia coli O157:H7 and E. coli K-12 genomic islands. Within a population of 60 E. coli O157:H7 strains, the existence of 65 accessory genomic regions identified by Panseq analysis was confirmed by PCR. The accessory genome and binary presence/absence data, and core genome and single nucleotide polymorphisms (SNPs) of six L. monocytogenes strains were extracted with Panseq and hierarchically clustered and visualized. The nucleotide core and binary accessory data were also used to construct maximum parsimony (MP) trees, which were compared to the MP tree generated by multi-locus sequence typing (MLST). The topology of the accessory and core trees was identical but differed from the tree produced using seven MLST loci. The Loci Selector module found the most variable and discriminatory combinations of four loci within a 100 loci set among 10 strains in 1 s, compared to the 449 s required to exhaustively search for all possible combinations; it also found the most discriminatory 20 loci from a 96 loci E. coli O157:H7 SNP dataset.ConclusionPanseq determines the core and accessory regions among a collection of genomic sequences based on user-defined parameters. It readily extracts regions unique to a genome or group of genomes, identifies SNPs within shared core genomic regions, constructs files for use in phylogeny programs based on both the presence/absence of accessory regions and SNPs within core regions and produces a graphical overview of the output. Panseq also includes a loci selector that calculates the most variable and discriminatory loci among sets of accessory loci or core gene SNPs.AvailabilityPanseq is freely available online at http://76.70.11.198/panseq. Panseq is written in Perl.

Molecular subtypes of Campylobacter spp., Salmonella enterica, and Escherichia coli O157:H7 isolated from faecal and surface water samples in the Oldman River watershed, Alberta, Canada

Campylobacter spp., Salmonella enterica, and Escherichia coli O157:H7 isolated from 898 faecal, 43 sewage, and 342 surface water samples from the Oldman River were characterized using bacterial subtyping methods in order to investigate potential sources of contamination of the watershed. Among these pathogens, Campylobacter spp. were the most frequently isolated from faecal, sewage, and surface water samples (266/895, 11/43, and 91/342, respectively), followed by Salmonella (67/898, 8/43, and 29/342, respectively), and E. coli O157:H7 (16/898, 2/43, and 8/342, respectively). Salmonella Rubislaw was the most common serovar isolated from water. This serovar was also isolated from two wild bird species. Most other serovars isolated from water were either not isolated from animals or were isolated from multiple species. E. coli O157:H7 was predominantly isolated from cattle. The most common phage-types of this pathogen from cattle were also the most common among water isolates, and there were exact pulsed field gel electrophoresis and comparative genomic fingerprint matches between cattle, sewage, and water isolates. Campylobacters were commonly isolated from surface waters and faeces from most animal species. Restriction fragment length polymorphism of the Campylobacter flaA gene identified several location and host species-specific (cattle, goose, pig) fingerprints. Molecular subtyping of these bacterial pathogens shows considerable promise as a tool for determining the sources of faecal pollution of water.

Development and validation of a comparative genomic fingerprinting method for high-resolution genotyping of Campylobacter jejuni.

ABSTRACT Campylobacter spp. are a leading cause of bacterial gastroenteritis worldwide. The need for molecular subtyping methods with enhanced discrimination in the context of surveillance- and outbreak-based epidemiologic investigations of Campylobacter spp. is critical to our understanding of sources and routes of transmission and the development of mitigation strategies to reduce the incidence of campylobacteriosis. We describe the development and validation of a rapid and high-resolution comparative genomic fingerprinting (CGF) method for C. jejuni. A total of 412 isolates from agricultural, environmental, retail, and human clinical sources obtained from the Canadian national integrated enteric pathogen surveillance program (C-EnterNet) were analyzed using a 40-gene assay (CGF40) and multilocus sequence typing (MLST). The significantly higher Simpsons index of diversity (ID) obtained with CGF40 (ID = 0.994) suggests that it has a higher discriminatory power than MLST at both the level of clonal complex (ID = 0.873) and sequence type (ID = 0.935). High Wallace coefficients obtained when CGF40 was used as the primary typing method suggest that CGF and MLST are highly concordant, and we show that isolates with identical MLST profiles are comprised of isolates with distinct but highly similar CGF profiles. The high concordance with MLST coupled with the ability to discriminate between closely related isolates suggests that CFG40 is useful in differentiating highly prevalent sequence types, such as ST21 and ST45. CGF40 is a high-resolution comparative genomics-based method for C. jejuni subtyping with high discriminatory power that is also rapid, low cost, and easily deployable for routine epidemiologic surveillance and outbreak investigations.

BACTERIAL PATHOGENS IN RURAL WATER SUPPLIES IN SOUTHERN ALBERTA, CANADA

Raw river and irrigation water in the Oldman River Basin in southern Alberta was tested for the presence of two bacterial pathogens, Escherichia coli O157:H7 and Salmonella spp., over the last 2 yr (2000–2001). The number of E. coli O157:H7 and Salmonella spp. isolated from raw water peaked during the summer months. While E. coli O157:H7 was only isolated from 11/802 (1.3%) of raw water samples over the entire sampling season in 2000 and from 16/806 (2.0%) of the samples in 2001, the pathogen was isolated one or more times from 10/35 (28.5%) sampling sites in 2000 and from 13/40 (32.5%) sampling sites in 2001. Salmonella was isolated from 44/802 (5.5%) of raw water samples in 2000 and from 122/822 (14.9%) of the samples in 2001; the pathogen was isolated one or more times from 25/35 (71.4%) sampling sites in 2000 and from 29/40 (72.5%) sampling sites in 2001. Certain sites had multiple pathogen isolations in the same year and from year to year. Salmonella Rublislaw was the most common Salmonella serovar isolated in both years, accounting for 52.4% of isolates.

The occurrence and sources of Campylobacter spp., Salmonella enterica and Escherichia coli O157:H7 in the Salmon River, British Columbia, Canada

In this study, we wished to assess the prevalence and determine the sources of three zoonotic bacterial pathogens (Salmonella, Campylobacter, and Escherichia coli O157:H7) in the Salmon River watershed in southwestern British Columbia. Surface water, sewage, and animal faecal samples were collected from the watershed. Selective bacterial culture and PCR techniques were used to isolate these three pathogens and indicator bacteria from these samples and characterize them. Campylobacter was the most prevalent pathogen in all samples, followed by Salmonella, and E. coli O157:H7. E. coli O157:H7 and Salmonella isolation rates from water, as well as faecal coliform densities correlated positively with precipitation, while Campylobacter isolation rates correlated negatively with precipitation. Analysis of DNA extracted from water samples for the presence of Bacteroides host-species markers, and comparisons of C. jejuni flaA-RFLP types and Salmonella serovars from faecal and water samples provided evidence that human sewage and specific domestic and wild animal species were sources of these pathogens; however, in most cases the source could not be determined or more than one source was possible. The frequent isolation of these zoonotic pathogens in the Salmon River highlights the risks to human health associated with intentional and unintentional consumption of untreated surface waters.

Lineage and Host Source Are Both Correlated with Levels of Shiga Toxin 2 Production by Escherichia coli O157:H7 Strains

ABSTRACT Escherichia coli O157:H7 strains fall into three major genetic lineages that differ in their distribution among humans and cattle. Several recent studies have reported differences in the expression of virulence factors between E. coli O157:H7 strains from these two host species. In this study, we wished to determine if important virulence-associated “mobile genetic elements” such as Shiga toxin 2 (Stx2)-encoding prophage are lineage restricted or are host source related and acquired independently of the pathogen genotype. DNA sequencing of the stx2 flanking region from a lineage II (LII) strain, EC970520, revealed that the transcriptional activator gene Q in LI strain EDL933 (upstream of stx2) is replaced by a pphA (serine/threonine phosphatase) homologue and an altered Q gene in this and all other LII strains tested. In addition, nearly all LI strains carried stx2, whereas all LII strains carried variant stx2c and 4 of 14 LI/II strains had copies of both stx2 and variant stx2c. Real-time PCR (RT-PCR) and enzyme-linked immunosorbent assay (ELISA) demonstrated that LI and LI/II strains produce significantly more stx2 mRNA and Stx2 than LII strains. However, among LI strains significantly more Stx2 is also produced by strains from humans than from cattle. Therefore, lineage-associated differences among E. coli O157:H7 strains such as prophage content, toxin type, and toxin expression may contribute to host isolation bias. However, the level of Stx2 production alone may also play an important role in the within-lineage association of E. coli O157:H7 strains with human clinical disease.

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.

Spatial and temporal drivers of zoonotic pathogen contamination of an agricultural watershed.

In regions where animal agriculture is prominent, such as southern Alberta, higher rates of gastrointestinal illness have been reported when compared with nonagricultural regions. This difference in the rate of illness is thought to be a result of increased zoonotic pathogen exposure through environmental sources such as water. In this study, temporal and spatial factors associated with bacterial pathogen contamination of the Oldman River, which transverses this region, were analyzed using classification and regression tree analysis. Significantly higher levels of fecal indicators; more frequent isolations of Campylobacter spp., Escherichia coli O157:H7, and Salmonella enterica spp.; and higher rates of detection of pig-specific Bacteroides markers occurred at downstream sites than at upstream sites, suggesting additive stream inputs. Fecal indicator densities were also significantly higher when any one of these three bacterial pathogens was present and where there were higher total animal manure units; however, occasionally pathogens were present when fecal indicator levels were low or undetectable. Overall, Salmonella spp., Campylobacter spp., and E. coli O157:H7 presence was associated with season, animal manure units, and total rainfall on the day of sampling and 3 d in advance of sampling. Several of the environmental variables analyzed in this study appear to influence pathogen prevalence and therefore may be useful in predicting water quality and safety and in the improvement of watershed management practices in this and other agricultural regions.

Escherichia coli O157:H7 Strain Origin, Lineage, and Shiga Toxin 2 Expression Affect Colonization of Cattle†

ABSTRACT Enterohemorrhagic Escherichia coli O157:H7 has evolved into an important human pathogen with cattle as the main reservoir. The recent discovery of E. coli O157:H7-induced pathologies in challenged cattle has suggested that previously discounted bacterial virulence factors may contribute to the colonization of cattle. The objective of the present study was to examine the impact of lineage type, cytotoxin activity, and cytotoxin expression on the amount of E. coli O157:H7 colonization of cattle tissue and cells in vitro. Using selected bovine- and human-origin strains, we determined that lineage type predicted the amount of E. coli O157:H7 strain colonization: lineage I > intermediate lineages > lineage II. All E. coli O157:H7 strain colonization was dose dependent, with threshold colonization at 103 to 105 CFU and maximum colonization at 107 CFU. We also determined that an as-yet-unknown factor of strain origin was the most dominant predictor of the amount of strain colonization in vitro. The amount of E. coli O157:H7 colonization was also influenced by strain cytotoxin activity and the inclusion of cytotoxins from lineage I or intermediate lineage strains increased colonization of a lineage II strain. There was a higher level of expression of the Shiga toxin 1 gene (stx1) in human-origin strains than in bovine-origin strains. In addition, lineage I strains expressed higher levels of the Shiga toxin 2 gene (stx2). The present study supports a role for strain origin, lineage type, cytotoxin activity, and stx2 expression in modulating the amount of E. coli O157:H7 colonization of cattle.

Further development of sample preparation and detection methods for O157 and the top 6 non-O157 STEC serogroups in cattle feces.

Shiga toxin-producing Escherichia coli (STEC) are food-borne pathogens responsible for outbreaks of human infections worldwide. Ruminant livestock harbor STEC in their intestinal tract, and through fecal contamination possess the potential to compromise the safety of food and water. As a human health safety risk, STEC detection methods on beef carcasses and trim are needed as mandated by the USDA-FSIS. In order to monitor STEC prior to harvest and human consumption, our goal was to evaluate and/or improve detection of seven STEC serogroups in cattle feces. In comparison to traditional approaches, sample processing methods in bovine feces were evaluated using a multi-factorial Latin square design which involved freezing or freeze drying feces. Autoclaved versus non-autoclaved feces were spiked with O26:H11 or O157:H7 serotypes in various dilutions and enriched for up to 6h. Each hour, enriched aliquots were compared using traditional culture methods and quantitative polymerase chain reaction (qPCR). Furthermore, a 7-serogroup multiplex PCR (mPCR) was developed to detect O26, O45, O103, O111, O121, O145 and O157 serogroups simultaneously. The diagnostic sensitivity of our mPCR assay following 6h enrichment was superior (10CFU/g across all serogroups) compared to a previously established PCR assay (10CFU/g for O26, and O103; ≥10(4)CFU/g for all other serogroups). Obtaining viable isolates appeared to be limited by the efficiency of current immunomagnetic separation (IMS) methods, which ranged from 20 to 100% effectiveness at retrieving colonies depending on serogroup. After IMS, 70 putative STEC isolates were screened for Shiga toxin and attachment genes by mPCR. Sixty-five isolates contained one or both Shiga toxin genes.