Whole genome sequencing for genomics-guided investigations of Escherichia coli O157:H7 outbreaks

Multi isolate whole genome sequencing (WGS) and typing for outbreak investigations has become a reality in the post-genomics era. We applied this technology to strains from Escherichia coli O157:H7 outbreaks. These include isolates from seven North America outbreaks, as well as multiple isolates from the same patient and from different infected individuals in the same household. Customized high-resolution bioinformatics sequence typing strategies were developed to assess the core genome and mobilome plasticity. Sequence typing was performed using an in-house single nucleotide polymorphism (SNP) discovery and validation pipeline. Discriminatory power becomes of particular importance for the investigation of isolates from outbreaks in which macrogenomic techniques such as pulse-field gel electrophoresis or multiple locus variable number tandem repeat analysis do not differentiate closely related organisms. We also characterized differences in the phage inventory, allowing us to identify plasticity among outbreak strains that is not detectable at the core genome level. Our comprehensive analysis of the mobilome identified multiple plasmids that have not previously been associated with this lineage. Applied phylogenomics approaches provide strong molecular evidence for exceptionally little heterogeneity of strains within outbreaks and demonstrate the value of intra-cluster comparisons, rather than basing the analysis on archetypal reference strains. Next generation sequencing and whole genome typing strategies provide the technological foundation for genomic epidemiology outbreak investigation utilizing its significantly higher sample throughput, cost efficiency, and phylogenetic relatedness accuracy. These phylogenomics approaches have major public health relevance in translating information from the sequence-based survey to support timely and informed countermeasures. Polymorphisms identified in this work offer robust phylogenetic signals that index both short- and long-term evolution and can complement currently employed typing schemes for outbreak ex- and inclusion, diagnostics, surveillance, and forensic studies

Detection and Identification of Common Food-Borne Viruses with a Tiling Microarray

Microarray hybridization based identification of viral genotypes is increasingly assuming importance due to outbreaks of multiple pathogenic viruses affecting humans causing wide-spread morbidity and mortality. Surprisingly, microarray based identification of food-borne viruses, one of the largest groups of pathogenic viruses, causing more than 1.5 billion infections world-wide every year, has lagged behind. Cell-culture techniques are either unavailable or time consuming for routine application. Consequently, current detection methods for these pathogens largely depend on polymerase chain reaction (PCR) based techniques, generally requiring an investigator to preselect the target virus of interest. Here we describe the first attempt to use the microarray as an identification tool for these viruses. We have developed methodology to synthesize targets for virus identification without using PCR, making the process genuinely sequence independent. We show here that a tiling microarray can simultaneously detect and identify the genotype and strain of common food-borne viruses in a single experiment

A Microarray Based Approach for the Identification of Common Foodborne Viruses

An oligonucleotide array (microarray) incorporating 13,000 elements representing selected strains of hepatitis A virus (HAV), human coxsackieviruses A and B (CVA and CVB), genogroups I and II of Norovirus (NV), and human rotavirus (RV) gene segments 3,4,10, and 11 was designed based on the principle of tiling. Each oligonucleotide was 29 bases long, starting at every 5th base of every sequence, resulting in an overlap of 24 bases in two consecutive oligonucleotides. The applicability of the array for virus identification was examined using PCR amplified products from multiple HAV and CV strains. PCR products labeled with biotin were hybridized to the array, and the biotin was detected using a brief reaction with Cy3-labeled streptavidin, the array subjected to laser scanning, and the hybridization data plotted as fluorescence intensity against each oligonucleotide in the array. The combined signal intensities of all probes representing a particular strain of virus were calculated and plotted against all virus strains identified on a linear representation of the array. The profile of the total signal intensity identified the strain that is most likely represented in the amplified cDNA target. The results obtained with HAV and CV indicated that the hybridization profile thus generated can be used to identify closely related viral strains. This represents a significant improvement over current methods for virus identification using PCR amplification and amplicon sequencing

Multiple antimicrobial resistance in plague: An emerging public health risk

Antimicrobial resistance in Yersinia pestis is rare, yet constitutes a significant international public health and biodefense threat. In 1995, the first multidrug resistant (MDR) isolate of Y. pestis (strain IP275) was identified, and was shown to contain a self-transmissible plasmid (pIP1202) that conferred resistance to many of the antimicrobials recommended for plague treatment and prophylaxis. Comparative analysis of the DNA sequence of Y. pestis plasmid pIP1202 revealed a near identical IncA/C plasmid backbone that is shared by MDR plasmids isolated from Salmonella enterica serotype Newport SL254 and the fish pathogen Yersinia ruckeri YR71. The high degree of sequence identity and gene synteny between the plasmid backbones suggests recent acquisition of these plasmids from a common ancestor. In addition, the Y. pestis pIP1202-like plasmid backbone was detected in numerous MDR enterobacterial pathogens isolated from retail meat samples collected between 2002 and 2005 in the United States. Plasmid-positive strains were isolated from beef, chicken, turkey and pork, and were found in samples from the following states: California, Colorado, Connecticut, Georgia, Maryland, Minnesota, New Mexico, New York and Oregon. Our studies reveal that this common plasmid backbone is broadly disseminated among MDR zoonotic pathogens associated with agriculture. This reservoir of mobile resistance determinants has the potential to disseminate to Y. pestis and other human and zoonotic bacterial pathogens and therefore represents a significant public health concern

Development of a targeted amplicon sequencing method for genotyping Cyclospora cayetanensis from fresh produce and clinical samples with enhanced genomic resolution and sensitivity

Outbreaks of cyclosporiasis, an enteric illness caused by the parasite Cyclospora cayetanensis, have been associated with consumption of various types of fresh produce. Although a method is in use for genotyping C. cayetanensis from clinical specimens, the very low abundance of C. cayetanensis in food and environmental samples presents a greater challenge. To complement epidemiological investigations, a molecular surveillance tool is needed for use in genetic linkage of food vehicles to cyclosporiasis illnesses, estimation of the scope of outbreaks or clusters of illness, and determination of geographical areas involved. We developed a targeted amplicon sequencing (TAS) assay that incorporates a further enrichment step to gain the requisite sensitivity for genotyping C. cayetanensis contaminating fresh produce samples. The TAS assay targets 52 loci, 49 of which are located in the nuclear genome, and encompasses 396 currently known SNP sites. The performance of the TAS assay was evaluated using lettuce, basil, cilantro, salad mix, and blackberries inoculated with C. cayetanensis oocysts. A minimum of 24 markers were haplotyped even at low contamination levels of 10 oocysts in 25 g leafy greens. The artificially contaminated fresh produce samples were included in a genetic distance analysis based on haplotype presence/absence with publicly available C. cayetanensis whole genome sequence assemblies. Oocysts from two different sources were used for inoculation, and samples receiving the same oocyst preparation clustered together, but separately from the other group, demonstrating the utility of the assay for genetically linking samples. Clinical fecal samples with low parasite loads were also successfully genotyped. This work represents a significant advance in the ability to genotype C. cayetanensis contaminating fresh produce along with greatly expanding the genomic diversity included for genetic clustering of clinical specimens

Phylogenomic Analysis of Salmonella enterica subsp. enterica Serovar Bovismorbificans from Clinical and Food Samples Using Whole Genome Wide Core Genes and kmer Binning Methods to Identify Two Distinct Polyphyletic Genome Pathotypes

Salmonella enterica subsp. enterica serovar Bovismorbificans has caused multiple outbreaks involving the consumption of produce, hummus, and processed meat products worldwide. To elucidate the intra-serovar genomic structure of S. Bovismorbificans, a core-genome analysis with 2690 loci (based on 150 complete genomes representing Salmonella enterica serovars developed as part of this study) and a k-mer-binning based strategy were carried out on 95 whole genome sequencing (WGS) assemblies from Swiss, Canadian, and USA collections of S. Bovismorbificans strains from foodborne infections. Data mining of a digital DNA tiling array of legacy SARA and SARB strains was conducted to identify near-neighbors of S. Bovismorbificans. The core genome analysis and the k-mer-binning methods identified two polyphyletic clusters, each with emerging evolutionary properties. Four STs (2640, 142, 1499, and 377), which constituted the majority of the publicly available WGS datasets from >260 strains analyzed by k-mer-binning based strategy, contained a conserved core genome backbone with a different evolutionary lineage as compared to strains comprising the other cluster (ST150). In addition, the assortment of genotypic features contributing to pathogenesis and persistence, such as antimicrobial resistance, prophage, plasmid, and virulence factor genes, were assessed to understand the emerging characteristics of this serovar that are relevant clinically and for food safety concerns. The phylogenomic profiling of polyphyletic S. Bovismorbificans in this study corresponds to intra-serovar variations observed in S. Napoli and S. Newport serovars using similar high-resolution genomic profiling approaches and contributes to the understanding of the evolution and sequence divergence of foodborne Salmonellae. These intra-serovar differences may have to be thoroughly understood for the accurate classification of foodborne Salmonella strains needed for the uniform development of future food safety mitigation strategies

High Density Microarray Analysis Reveals New Insights into Genetic Footprints of Listeria monocytogenes Strains Involved in Listeriosis Outbreaks

Listeria monocytogenes, a foodborne bacterial pathogen, causes invasive and febrile gastroenteritis forms of listeriosis in humans. Both invasive and febrile gastroenteritis listeriosis is caused mostly by serotypes 1/2a, 1/2b and 4b strains. The outbreak strains of serotype 1/2a and 4b could be further classified into several epidemic clones but the genetic bases for the diverse pathophysiology have been unsuccessful. DNA microarray provides an important tool to scan the entire genome for genetic signatures that may distinguish the L. monocytogenes strains belonging to different outbreaks. We have designed a pan-genomic microarray chip (Listeria GeneChip) containing sequences from 24 L. monocytogenes strains. The chip was designed to identify the presence/absence of genomic sequences, analyze transcription profiles and identify SNPs. Analysis of the genomic profiles of 38 outbreak strains representing 1/2a, 1/2b and 4b serotypes, revealed that the strains formed distinct genetic clusters adhering to their serotypes and epidemic clone types. Although serologically 1/2a and 1/b strains share common antigenic markers microarray analysis revealed that 1/2a strains are further apart from the closely related 1/2b and 4b strains. Within any given serotype and epidemic clone type the febrile gastroenteritis and invasive strains can be further distinguished based on several genetic markers including large numbers of phage genome, and intergenic sequences. Our results showed that the microarray-based data can be an important tool in characterization of L. monocytogenes strains involved in both invasive and gastroenteritis outbreaks. The results for the first time showed that the serotypes and epidemic clones are based on extensive pan-genomic variability and the 1/2b and 4bstrains are more closely related to each other than the 1/2a strains. The data also supported the hypothesis that the strains causing these two diverse outbreaks are genotypically different and this finding might be important in understanding the pathophysiology of this organism

Inter- and Intra-Host Nucleotide Variations in Hepatitis A Virus in Culture and Clinical Samples Detected by Next-Generation Sequencing

The accurate virus detection, strain discrimination, and source attribution of contaminated food items remains a persistent challenge because of the high mutation rates anticipated to occur in foodborne RNA viruses, such as hepatitis A virus (HAV). This has led to predictions of the existence of more than one sequence variant between the hosts (inter-host) or within an individual host (intra-host). However, there have been no reports of intra-host variants from an infected single individual, and little is known about the accuracy of the single nucleotide variations (SNVs) calling with various methods. In this study, the presence and identity of viral SNVs, either between HAV clinical specimens or among a series of samples derived from HAV clone1-infected FRhK4 cells, were determined following analyses of nucleotide sequences generated using next-generation sequencing (NGS) and pyrosequencing methods. The results demonstrate the co-existence of inter- and intra-host variants both in the clinical specimens and the cultured samples. The discovery and confirmation of multi-viral RNAs in an infected individual is dependent on the strain discrimination at the SNV level, and critical for successful outbreak traceback and source attribution investigations. The detection of SNVs in a time series of HAV infected FRhK4 cells improved our understanding on the mutation dynamics determined probably by different selective pressures. Additionally, it demonstrated that NGS could potentially provide a valuable investigative approach toward SNV detection and identification for other RNA viruses