Strengthening Bordetella pertussis genomic surveillance by direct sequencing of residual positive specimens.

Whole-genome sequencing (WGS) of microbial pathogens recovered from patients with infectious disease facilitates high-resolution strain characterization and molecular epidemiology. However, increasing reliance on culture-independent methods to diagnose infectious diseases has resulted in few isolates available for WGS. Here, we report a novel culture-independent approach to genome characterization of Bordetella pertussis, the causative agent of pertussis and a paradigm for insufficient genomic surveillance due to limited culture of clinical isolates. Sequencing libraries constructed directly from residual pertussis-positive diagnostic nasopharyngeal specimens were hybridized with biotinylated RNA "baits" targeting B. pertussis fragments within complex mixtures that contained high concentrations of host and microbial background DNA. Recovery of B. pertussis genome sequence data was evaluated with mock and pooled negative clinical specimens spiked with reducing concentrations of either purified DNA or inactivated cells. Targeted enrichment increased the yield of B. pertussis sequencing reads up to 90% while simultaneously decreasing host reads to less than 10%. Filtered sequencing reads provided sufficient genome coverage to perform characterization via whole-genome single nucleotide polymorphisms and whole-genome multilocus sequencing typing. Moreover, these data were concordant with sequenced isolates recovered from the same specimens such that phylogenetic reconstructions from either consistently clustered the same putatively linked cases. The optimized protocol is suitable for nasopharyngeal specimens with diagnostic IS481 Ct < 35 and >10 ng DNA. Routine implementation of these methods could strengthen surveillance and study of pertussis resurgence by capturing additional cases with genomic characterization.

Exploring Antimicrobial Resistance Profiles of E. coli Isolates in Dairy Cattle: A Baseline Study across Dairy Farms with Varied Husbandry Practices in Puerto Rico.

Antimicrobial treatment in livestock can contribute to the emergence and spread of antimicrobial-resistant (AMR) microorganisms. Despite substantial surveillance of AMR bacteria in the continental United States, the prevalence of these AMR organisms in U.S. territories, such as Puerto Rico, remains understudied. The goals of this research included obtaining baseline data on the antimicrobial profile of E. coli isolates from Puerto Rico dairy farms with different husbandry practices. Seventy-nine fecal samples were collected from two types of conventional dairy farms: those that fed calves with tank milk and those that fed calves with waste milk. These samples were collected from the animals' rectums, culture, and subsequently confirmed through biochemical tests. Out of these samples, 32 isolates were analyzed phenotypically and genotypically to elucidate their AMR profiles. The results underscore a discrepancy in the occurrence of antimicrobial resistance genes between calves and adult cattle. Notably, waste milk-fed calves exhibited a significantly higher prevalence of antibiotic-resistant E. coli when compared to their tank milk-fed counterparts. These disparities emphasize the need for more comprehensive investigations to determine causative factors. These results underscore the urgency of comprehensive strategies to raise awareness about how management practices influence antimicrobial resistance, shifting the focus from treatment to prevention.

Genomic characterization of cocirculating Corynebacterium diphtheriae and non-diphtheritic Corynebacterium species among forcibly displaced Myanmar nationals, 2017-2019.

Respiratory diphtheria is a serious infection caused by toxigenic Corynebacterium diphtheriae, and disease transmission mainly occurs through respiratory droplets. Between 2017 and 2019, a large diphtheria outbreak among forcibly displaced Myanmar nationals densely settled in Bangladesh was investigated. Here we utilized whole-genome sequencing (WGS) to characterize recovered isolates of C. diphtheriae and two co-circulating non-diphtheritic Corynebacterium (NDC) species - C. pseudodiphtheriticum and C. propinquum. C. diphtheriae isolates recovered from all 53 positive cases in this study were identified as toxigenic biovar mitis, exhibiting intermediate resistance to penicillin, and formed four phylogenetic clusters circulating among multiple refugee camps. Additional sequenced isolates collected from two patients showed co-colonization with non-toxigenic C. diphtheriae biovar gravis, one of which exhibited decreased susceptibility to the first-line antibiotics and harboured a novel 23-kb multidrug resistance plasmid. Results of phylogenetic reconstruction and virulence-related gene contents of the recovered NDC isolates indicated they were likely commensal organisms, though 80.4 %(45/56) were not susceptible to erythromycin, and most showed high minimum inhibition concentrations against azithromycin. These results demonstrate the high resolution with which WGS can aid molecular investigation of diphtheria outbreaks, through the quantification of bacterial genetic relatedness, as well as the detection of virulence factors and antibiotic resistance markers among case isolates.

Virulence Potential of Escherichia coli Isolates from a Beef Farm in Puerto Rico.

Sixteen Escherichia coli isolates were obtained from fecal matter from a beef farm in Puerto Rico. Isolates were whole-genome sequenced for in silico characterization, including pathotype characterization, virulence, and plasmid identification. The results of the draft genomes identified potential pathogenic E. coli strains from beef cattle in Puerto Rico.

Benchmark datasets for SARS-CoV-2 surveillance bioinformatics.

Background: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the cause of coronavirus disease 2019 (COVID-19), has spread globally and is being surveilled with an international genome sequencing effort. Surveillance consists of sample acquisition, library preparation, and whole genome sequencing. This has necessitated a classification scheme detailing Variants of Concern (VOC) and Variants of Interest (VOI), and the rapid expansion of bioinformatics tools for sequence analysis. These bioinformatic tools are means for major actionable results: maintaining quality assurance and checks, defining population structure, performing genomic epidemiology, and inferring lineage to allow reliable and actionable identification and classification. Additionally, the pandemic has required public health laboratories to reach high throughput proficiency in sequencing library preparation and downstream data analysis rapidly. However, both processes can be limited by a lack of a standardized sequence dataset. Methods: We identified six SARS-CoV-2 sequence datasets from recent publications, public databases and internal resources. In addition, we created a method to mine public databases to identify representative genomes for these datasets. Using this novel method, we identified several genomes as either VOI/VOC representatives or non-VOI/VOC representatives. To describe each dataset, we utilized a previously published datasets format, which describes accession information and whole dataset information. Additionally, a script from the same publication has been enhanced to download and verify all data from this study. Results: The benchmark datasets focus on the two most widely used sequencing platforms: long read sequencing data from the Oxford Nanopore Technologies platform and short read sequencing data from the Illumina platform. There are six datasets: three were derived from recent publications; two were derived from data mining public databases to answer common questions not covered by published datasets; one unique dataset representing common sequence failures was obtained by rigorously scrutinizing data that did not pass quality checks. The dataset summary table, data mining script and quality control (QC) values for all sequence data are publicly available on GitHub: https://github.com/CDCgov/datasets-sars-cov-2. Discussion: The datasets presented here were generated to help public health laboratories build sequencing and bioinformatics capacity, benchmark different workflows and pipelines, and calibrate QC thresholds to ensure sequencing quality. Together, improvements in these areas support accurate and timely outbreak investigation and surveillance, providing actionable data for pandemic management. Furthermore, these publicly available and standardized benchmark data will facilitate the development and adjudication of new pipelines.

Complete Genome Sequences of Four Macrolide-Resistant Nondiphtheritic Corynebacterium Isolates.

This report describes the complete genome sequences of four isolates of the nondiphtheritic Corynebacterium (NDC) species Corynebacterium pseudodiphtheriticum and Corynebacterium propinquum, recovered during investigation of a large diphtheria outbreak in Bangladesh. These data will assist in better delineating the boundary between these related species and understanding their virulence potential.

Toxigenic Corynebacterium diphtheriae Infection in Cat, Texas, USA.

We report a toxigenic strain of Corynebacterium diphtheriae isolated from an oozing dermal wound in a pet cat in Texas, USA. We also describe the epidemiologic public health efforts conducted to identify potential sources of infection and mitigate its spread and the molecular and genetic studies performed to identify the bacterium.

Genomic epidemiology of nontoxigenic Corynebacterium diphtheriae from King County, Washington State, USA between July 2018 and May 2019.

Between July 2018 and May 2019, Corynebacterium diphtheriae was isolated from eight patients with non-respiratory infections, seven of whom experienced homelessness and had stayed at shelters in King County, WA, USA. All isolates were microbiologically identified as nontoxigenic C. diphtheriae biovar mitis. Whole-genome sequencing confirmed that all case isolates were genetically related, associated with sequence type 445 and differing by fewer than 24 single-nucleotide polymorphisms (SNPs). Compared to publicly available C. diphtheriae genomic data, these WA isolates formed a discrete cluster with SNP variation consistent with previously reported outbreaks. Virulence-related gene content variation within the highly related WA cluster isolates was also observed. These results indicated that genome characterization can readily support epidemiology of nontoxigenic C. diphtheriae.

Novel Effective Bacillus cereus Group Species "Bacillus clarus" Is Represented by Antibiotic-Producing Strain ATCC 21929 Isolated from Soil.

Gram-positive, spore-forming members of the Bacillus cereus group species complex are widespread in natural environments and display various degrees of pathogenicity. Recently, B. cereus group strain Bacillus mycoides Flugge ATCC 21929 was found to represent a novel lineage within the species complex, sharing a relatively low degree of genomic similarity with all B. cereus group genomes (average nucleotide identity [ANI] < 88). ATCC 21929 has been previously associated with the production of a patented antibiotic, antibiotic 60-6 (i.e., cerexin A); however, the virulence potential and growth characteristics of this lineage have never been assessed. Here, we provide an extensive genomic and phenotypic characterization of ATCC 21929, and we assess its pathogenic potential in vitro. ATCC 21929 most closely resembles Bacillus paramycoides NH24A2T (ANI and in silico DNA-DNA hybridization values of 86.70 and 34.10%, respectively). Phenotypically, ATCC 21929 does not possess cytochrome c oxidase activity and is able to grow at a range of temperatures between 15 and 43°C and a range of pH between 6 and 9. At 32°C, ATCC 21929 shows weak production of diarrheal enterotoxin hemolysin BL (Hbl) but no production of nonhemolytic enterotoxin (Nhe); at 37°C, neither Hbl nor Nhe is produced. Additionally, at 37°C, ATCC 21929 does not exhibit cytotoxic effects toward HeLa cells. With regard to fatty acid composition, ATCC 21929 has iso-C17:0 present in highest abundance. Based on the characterization provided here, ATCC 21929T (= PS00077AT = PS00077BT = PSU-0922T = BHPT) represents a novel effective B. cereus group species, which we propose as effective species "Bacillus clarus"IMPORTANCE The B. cereus group comprises numerous closely related lineages with various degrees of pathogenic potential and industrial relevance. Species-level taxonomic classification of B. cereus group strains is important for risk evaluation and communication but remains challenging. Biochemical and phenotypic assays are often used to assign B. cereus group strains to species but are insufficient for accurate taxonomic classification on a genomic scale. Here, we show that antibiotic-producing ATCC 21929 represents a novel lineage within the B. cereus group that, by all metrics used to delineate prokaryotic species, exemplifies a novel effective species. Furthermore, we show that ATCC 21929 is incapable of producing enterotoxins Hbl and Nhe or exhibiting cytotoxic effects on HeLa cells at human body temperature in vitro These results provide greater insight into the genomic and phenotypic diversity of the B. cereus group and may be leveraged to inform future public health and food safety efforts.

Detection and Characterization of Diphtheria Toxin Gene-Bearing Corynebacterium Species through a New Real-Time PCR Assay.

Respiratory diphtheria, characterized by a firmly adherent pseudomembrane, is caused by toxin-producing strains of Corynebacterium diphtheriae, with similar illness produced occasionally by toxigenic Corynebacterium ulcerans or, rarely, Corynebacterium pseudotuberculosis While diphtheria laboratory confirmation requires culture methods to determine toxigenicity, real-time PCR (RT-PCR) provides a faster method to detect the toxin gene (tox). Nontoxigenic tox-bearing (NTTB) Corynebacterium isolates have been described, but impact of these isolates on the accuracy of molecular diagnostics is not well characterized. Here, we describe a new triplex RT-PCR assay to detect tox and distinguish C. diphtheriae from the closely related species C. ulcerans and C. pseudotuberculosis Analytical sensitivity and specificity of the assay were assessed in comparison to culture using 690 previously characterized microbial isolates. The new triplex assay characterized Corynebacterium isolates accurately, with 100% analytical sensitivity for all targets. Analytical specificity with isolates was 94.1%, 100%, and 99.5% for tox, Diph_rpoB, and CUP_rpoB targets, respectively. Twenty-nine NTTB Corynebacterium isolates, representing 5.9% of 494 nontoxigenic isolates tested, were detected by RT-PCR. Whole-genome sequencing of NTTB isolates revealed varied mutations putatively underlying their lack of toxin production, as well as eight isolates with no mutation in tox or the promoter region. This new Corynebacterium RT-PCR method provides a rapid tool to screen isolates and identify probable diphtheria cases directly from specimens. However, the sporadic occurrence of NTTB isolates reinforces the viewpoint that diphtheria culture diagnostics continue to provide the most accurate case confirmation.

Whole-genome sequencing reveals resistome of highly drug-resistant retail meat and human Salmonella Dublin.

Non-typhoidal Salmonella (NTS) are a significant source of foodborne illness worldwide, with disease symptoms most often presenting as self-limiting gastroenteritis; however, occasionally the infection spreads and becomes invasive, frequently requiring anti-microbial treatment. The cattle-adapted Dublin serovar of NTS has commonly been associated with invasive illness and anti-microbial resistance (AMR). Here, the enhanced resolution conferred by whole-genome sequencing was utilized to elucidate and compare the resistome and genetic relatedness of 14 multidrug-resistant (MDR) and one pan-susceptible S. Dublin, isolated primarily in Pennsylvania, from fresh retail meat (one isolate) and humans (14 isolates). Twelve different genetic AMR determinants, including both acquired and chromosomal, were identified. Furthermore, comparative plasmid analysis indicated that AMR was primarily conferred by a putative IncA/C2 plasmid. A single pan-susceptible S. Dublin isolate, collected from the same timeframe and geographical region as the MDR isolates, did not carry an IncA/C2 replicon sequence within its genome. Moreover, the pan-susceptible isolate was genetically distinct from its MDR counterparts, as it was separated by ≥267 single nucleotide polymorphisms (SNPs), whereas there was a ≤38 SNP distance between the MDR isolates. Collectively, this data set advances our understanding of the genetic basis of the highly drug-resistant nature of S. Dublin, a serovar with significant public health implications.

A Putative Microcin Amplifies Shiga Toxin 2a Production of Escherichia coli O157:H7.

Escherichia coli O157:H7 is a foodborne pathogen implicated in various multistate outbreaks. It encodes Shiga toxin on a prophage, and Shiga toxin production is linked to phage induction. An E. coli strain, designated 0.1229, that amplified Stx2a production when cocultured with E. coli O157:H7 strain PA2 was identified. Growth of PA2 in 0.1229 cell-free supernatants had a similar effect, even when supernatants were heated to 100°C for 10 min, but not after treatment with proteinase K. The secreted molecule was shown to use TolC for export and the TonB system for import. The genes sufficient for production of this molecule were localized to a 5.2-kb region of a 12.8-kb plasmid. This region was annotated, identifying hypothetical proteins, a predicted ABC transporter, and a cupin superfamily protein. These genes were identified and shown to be functional in two other E. coli strains, and bioinformatic analyses identified related gene clusters in similar and distinct bacterial species. These data collectively suggest that E. coli 0.1229 and other E. coli strains produce a microcin that induces the SOS response in target bacteria. Besides adding to the limited number of microcins known to be produced by E. coli, this study provides an additional mechanism by which stx2a expression is increased in response to the gut microflora.IMPORTANCE How the gut microflora influences the progression of bacterial infections is only beginning to be understood. Antibiotics are counterindicated for E. coli O157:H7 infections, limiting treatment options. An increased understanding of how the gut microflora directs O157:H7 virulence gene expression may lead to additional treatment options. This work identified E. coli strains that enhance the production of Shiga toxin by O157:H7 through the secretion of a proposed microcin. Microcins are natural antimicrobial peptides that target specific species, can act as alternatives to antibiotics, and mediate microbial competition. This work demonstrates another mechanism by which non-O157 E. coli strains may increase Shiga toxin production and adds to our understanding of microcins, a group of antimicrobials less well understood than colicins.

Detection of CTX-M-1 extended-spectrum beta-lactamase among ceftiofur-resistant Salmonella enterica clinical isolates of poultry.

Salmonella enterica resistance to extended-spectrum cephalosporins (ESC) conferred by cefotaximases (blaCTX-M) is a growing concern in the United States. Among food-producing animals, poultry are a major reservoir of ESC-resistant Salmonella. A retrospective study was carried out to further characterize 38 ceftiofur-resistant clinical Salmonella enterica isolates obtained from poultry during 2007-2018. Of the isolates tested, 31 displayed resistance to ceftriaxone and harbored blaCMY-2, whereas 7 isolates demonstrated resistance or reduced susceptibility to cefepime in addition to ceftriaxone resistance. These 7 isolates displayed extended-spectrum ß-lactamase activity, harbored blaCTX-M-1, and were recovered only from recent poultry diagnostic submissions made in 2011-2018 as opposed to the 31 isolates that were recovered in 2007-2018. Further characterization of the blaCTX-M-1 gene determined that it was located on conjugative IncN/ST1 and IncI1/ST87 plasmids in the isolates from commercial turkeys and broilers, respectively. These plasmids have been responsible for extensive spread of blaCTX-M-1 in livestock, poultry, and humans in Europe. Potential transfer of IncN and IncI1 plasmids and/or nontyphoidal Salmonella carrying these plasmids through the food chain, or by other means to humans, may result in treatment failures. Our study demonstrates the importance of further characterization of ceftiofur-resistant S. enterica isolates detected by veterinary diagnostic laboratories to identify the sources of blaCTX-M-1 and to mitigate the spread of ESC-resistant Salmonella in the poultry production pyramid.

Draft Genome Sequence of Enterococcus durans OSY-EGY, a Multiple-Antimicrobial-Peptide Producer Isolated from Egyptian Hard Cheese.

Enterococcus durans OSY-EGY was isolated recently from cheese. The strain produces potent antimicrobial agents. Here, we present the draft genome sequence of the strain, with a genome size of 3,230,625 bp and an average G+C content of 37.69%. Draft genome mining identified several biosynthetic gene clusters encoding multiple antimicrobial peptides.

Retrospective whole-genome sequencing analysis distinguished PFGE and drug-resistance-matched retail meat and clinical Salmonella isolates.

Non-typhoidal Salmonella is a leading cause of outbreak and sporadic-associated foodborne illnesses in the United States. These infections have been associated with a range of foods, including retail meats. Traditionally, pulsed-field gel electrophoresis (PFGE) and antibiotic susceptibility testing (AST) have been used to facilitate public health investigations of Salmonella infections. However, whole-genome sequencing (WGS) has emerged as an alternative tool that can be routinely implemented. To assess its potential in enhancing integrated surveillance in Pennsylvania, USA, WGS was used to directly compare the genetic characteristics of 7 retail meat and 43 clinical historic Salmonella isolates, subdivided into 3 subsets based on PFGE and AST results, to retrospectively resolve their genetic relatedness and identify antimicrobial resistance (AMR) determinants. Single nucleotide polymorphism (SNP) analyses revealed that the retail meat isolates within S. Heidelberg, S. Typhimurium var. O5- subset 1 and S. Typhimurium var. O5- subset 2 were separated from each primary PFGE pattern-matched clinical isolate by 6-12, 41-96 and 21-81 SNPs, respectively. Fifteen resistance genes were identified across all isolates, including fosA7, a gene only recently found in a limited number of Salmonella and a ≥95 % phenotype to genotype correlation was observed for all tested antimicrobials. Moreover, AMR was primarily plasmid-mediated in S. Heidelberg and S. Typhimurium var. O5- subset 2, whereas AMR was chromosomally carried in S. Typhimurium var. O5- subset 1. Similar plasmids were identified in both the retail meat and clinical isolates. Collectively, these data highlight the utility of WGS in retrospective analyses and enhancing integrated surveillance for Salmonella from multiple sources.

Whole-Genome Sequences of Escherichia coli Isolates from Cocoa Beans Imported from Bolivia.

Four Escherichia coli isolates were obtained from a batch of cocoa beans imported from Bolivia. The cocoa beans were rejected by a U.S. chocolate manufacturer due to poor microbiological quality. The four isolates were whole-genome sequenced and the sequences analyzed to identify genotypes, serotypes, and virulence and antimicrobial resistance genes.

Non-pathogenic Escherichia coli Enhance Stx2a Production of E. coli O157:H7 Through Both bamA-Dependent and Independent Mechanisms.

Intestinal colonization by the foodborne pathogen Escherichia coli O157:H7 leads to serious disease symptoms, including hemolytic uremic syndrome (HUS) and hemorrhagic colitis (HC). Synthesis of one or more Shiga toxins (Stx) is essential for HUS and HC development. The genes encoding Stx, including Stx2a, are found within a lambdoid prophage integrated in the E. coli O157:H7 chromosome. Enhanced Stx2a expression was reported when specific non-pathogenic E. coli strains were co-cultured with E. coli O157:H7, and it was hypothesized that this phenotype required the non-pathogenic E. coli to be sensitive to stx-converting phage infection. We tested this hypothesis by generating phage resistant non-pathogenic E. coli strains where bamA (an essential gene and Stx phage receptor) was replaced with an ortholog from other species. Such heterologous gene replacement abolished the ability of the laboratory strain E. coli C600 to enhance toxin production when co-cultured with E. coli O157:H7 strain PA2, which belongs to the hypervirulent clade 8. The extracellular loops of BamA (loop 4, 6, 7) were further shown to be important for infection by stx2a-converting phages. However, similar gene replacement in another commensal E. coli, designated 1.1954, revealed a bamA-independent mechanism for toxin amplification. Toxin enhancement by 1.1954 was not the result of phage infection through an alternative receptor (LamB or FadL), lysogen formation by stx2a-converting phages, or the production of a secreted molecule. Collectively, these data suggest that non-pathogenic E. coli can enhance toxin production through at least two mechanisms.

Draft Genome Sequence of Bacillus velezensis GF610, a Producer of Potent Anti-Listeria Agents.

Bacillus velezensis GF610 was isolated from soil in Illinois, USA, and found to produce amyloliquecidin GF610, a potent two-component antimicrobial peptide. We report here the GF610 strain draft genome sequence, which contains 4.29 Mb and an overall GC content of 45.91%.

Draft Genome Sequence of Brevibacillus laterosporus OSY-I1, a Strain That Produces Brevibacillin, Which Combats Drug-Resistant Gram-Positive Bacteria.

Brevibacillus laterosporus OSY-I1 is a Gram-positive spore-forming bacterium isolated from soil. The bacterium produces brevibacillin, an antimicrobial lipopeptide effective against several drug-resistant Gram-positive bacteria. Here, we present the draft genome sequence of the strain OSY-I1 and the gene cluster responsible for the biosynthesis of brevibacillin.

Escherichia coli O157:H7 strains harbor at least three distinct sequence types of Shiga toxin 2a-converting phages.

BACKGROUND: Shiga toxin-producing Escherichia coli O157:H7 is a foodborne pathogen that causes severe human diseases including hemolytic uremic syndrome (HUS). The virulence factor that mediates HUS, Shiga toxin (Stx), is encoded within the genome of a lambdoid prophage. Although draft sequences are publicly available for a large number of E. coli O157:H7 strains, the high sequence similarity of stx-converting bacteriophages with other lambdoid prophages poses challenges to accurately assess the organization and plasticity among stx-converting phages due to assembly difficulties. METHODS: To further explore genome plasticity of stx-converting prophages, we enriched phage DNA from 45 ciprofloxacin-induced cultures for subsequent 454 pyrosequencing to facilitate assembly of the complete phage genomes. In total, 22 stx2a-converting phage genomes were closed. RESULTS: Comparison of the genomes distinguished nine distinct phage sequence types (PSTs) delineated by variation in obtained sequences, such as single nucleotide polymorphisms (SNPs) and insertion sequence element prevalence and location. These nine PSTs formed three distinct clusters, designated as PST1, PST2 and PST3. The PST2 cluster, identified in two clade 8 strains, was related to stx2a-converting phages previously identified in non-O157 Shiga-toxin producing E. coli (STEC) strains associated with a high incidence of HUS. The PST1 cluster contained phages related to those from E. coli O157:H7 strain Sakai (lineage I, clade 1), and PST3 contained a single phage that was distinct from the rest but most related to the phage from E. coli O157:H7 strain EC4115 (lineage I/II, clade 8). Five strains carried identical stx2a-converting phages (PST1-1) integrated at the same chromosomal locus, but these strains produced different levels of Stx2. CONCLUSION: The stx2a-converting phages of E. coli O157:H7 can be categorized into at least three phage types. Diversification within a phage type is mainly driven by IS629 and by a small number of SNPs. Polymorphisms between phage genomes may help explain differences in Stx2a production between strains, however our data indicates that genes encoded external to the phage affect toxin production as well.