Precision metagenomics sequencing for food safety: hybrid assembly of Shiga toxin-producing Escherichia coli in enriched agricultural water

Culture-independent metagenomic sequencing of enriched agricultural water could expedite the detection and virulotyping of Shiga toxin-producing Escherichia coli (STEC). We previously determined the limits of a complete, closed metagenome-assembled genome (MAG) assembly and of a complete, fragmented MAG assembly for O157:H7 in enriched agricultural water using long reads (Oxford Nanopore Technologies, Oxford), which were 107 and 105 CFU/ml, respectively. However, the nanopore assemblies did not have enough accuracy to be used in Single Nucleotide Polymorphism (SNP) phylogenies and cannot be used for the precise identification of an outbreak STEC strain. The present study aimed to determine the limits of detection and assembly for STECs in enriched agricultural water by Illumina MiSeq sequencing technology alone, followed by establishing the limit of hybrid assembly with nanopore long-read sequencing using three different hybrid assemblers (SPAdes, Unicycler, and OPERA-MS). We also aimed to generate a genome with enough accuracy to be used in a SNP phylogeny. The classification of MiSeq and nanopore sequencing identified the same highly abundant species. Using the totality of the MiSeq output and a precision metagenomics approach in which the E. coli reads are binned before assembly, the limit of detection and assembly of STECs by MiSeq were determined to be 105 and 107 CFU/ml, respectively. While a complete, closed MAG could not be generated at any concentration, a complete, fragmented MAG was produced using the SPAdes assembler with an STEC concentration of at least 107 CFU/ml. At this concentration, hybrid assembled contigs aligned to the nanopore-assembled genome could be accurately placed in a neighbor-joining tree. The MiSeq limit of detection and assembly was less sensitive than nanopore sequencing, which was likely due to factors including the small starting material (50 vs. 1 μg) and the dilution of the library loaded on the cartridge. This pilot study demonstrates that MiSeq sequencing requires higher coverage in precision metagenomic samples; however, with sufficient concentration, STECs can be characterized and phylogeny can be accurately determined.

Whole-genome sequencing reveals changes in genomic diversity and distinctive repertoires of T3SS and T6SS effector candidates in Chilean clinical Campylobacter strains

Campylobacter is the leading cause of bacterial gastroenteritis worldwide and an emerging and neglected pathogen in South America. This zoonotic pathogen colonizes the gastrointestinal tract of a wide range of mammals and birds, with poultry as the most important reservoir for human infections. Apart from its high morbidity rates, the emergence of resistant strains is of global concern. The aims of this work were to determine genetic diversity, presence of antimicrobial resistance determinants and virulence potential of Campylobacter spp. isolated from patients with acute gastrointestinal disease at 'Clinica Alemana', Santiago de Chile. The study considered the isolation of Campylobacter spp., from stool samples during a 20-month period (January 2020 to September 2021). We sequenced (NextSeq, Illumina) and performed an in-depth analysis of the genome sequences of 88 Campylobacter jejuni and 2 Campylobacter coli strains isolated from clinical samples in Chile. We identified a high genetic diversity among C. jejuni strains and the emergence of prevalent clonal complexes, which were not identified in our previous reports. While ~40% of strains harbored a mutation in the gyrA gene associated with fluoroquinolone resistance, no macrolide-resistance determinants were detected. Interestingly, gene clusters encoding virulence factors such as the T6SS or genes associated with long-term sequelae such as Guillain-Barré syndrome showed lineage-relatedness. In addition, our analysis revealed a high degree of variability regarding the presence of fT3SS and T6SS effector proteins in comparison to type strains 81-176, F38011, and NCTC 11168 and 488. Our study provides important insights into the molecular epidemiology of this emerging foodborne pathogen. In addition, the differences observed regarding the repertoire of fT3SS and T6SS effector proteins could have an impact on the pathogenic potential and transmissibility of these Latin American isolates, posing another challenge in characterizing the infection dynamics of this emergent and neglected bacterial pathogen.

Multiple Displacement Amplification as a Solution for Low Copy Number Plasmid Sequencing

Plasmids play a major role in bacterial adaptation to environmental stress and often contribute to antibiotic resistance and disease virulence. Although the complete sequence of each plasmid is essential for studying plasmid biology, most antibiotic resistance and virulence plasmids in Salmonella are present only in a low copy number, making extraction and sequencing difficult. Long read sequencing technologies require higher concentrations of DNA to provide optimal results. To resolve this problem, we assessed the sufficiency of multiple displacement amplification (MDA) for replicating Salmonella plasmid DNA to a satisfactory concentration for accurate sequencing and multiplexing. Nine Salmonella enterica isolates, representing nine different serovars carrying plasmids for which sequence data are already available at NCBI, were cultured and their plasmids isolated using an alkaline lysis extraction protocol. We then used the Phi29 polymerase to perform MDA, thereby obtaining enough plasmid DNA for long read sequencing. These amplified plasmids were multiplexed and sequenced on one single molecule, real-time (SMRT) cell with the Pacific Biosciences (Pacbio) Sequel sequencer. We were able to close all Salmonella plasmids (sizes ranged from 38 to 166 Kb) with sequencing coverage from 24 to 2,582X. This protocol, consisting of plasmid isolation, MDA, and multiplex sequencing, is an effective and fast method for closing high-molecular weight and low-copy-number plasmids. This high throughput protocol reduces the time and cost of plasmid closure.

Whole-Genome Phylogenetic Analysis Reveals a Wide Diversity of Non-O157 STEC Isolated From Ground Beef and Cattle Feces

Shiga toxin-producing Escherichia coli (STEC) causes foodborne outbreaks that can lead to complications such as hemolytic uremic syndrome. Their main reservoir is cattle, and ground beef has been frequently associated with disease and outbreaks. In this study, we attempted to understand the genetic relationship among STEC isolated in Chile from different sources, their relationship to STEC from the rest of the world, and to identify molecular markers of Chilean STEC. We sequenced 62 STEC isolated in Chile using MiSeq Illumina. In silico typing was determined using tools of the Center Genomic Epidemiology, Denmark University (CGE/DTU). Genomes of our local STEC collection were compared with 113 STEC isolated worldwide through a core genome MLST (cgMLST) approach, and we also searched for distinct genes to be used as molecular markers of Chilean isolates. Genomes in our local collection were grouped based on serogroup and sequence type, and clusters were formed within local STEC. In the worldwide STEC analysis, Chilean STEC did not cluster with genomes of the rest of the world suggesting that they are not phylogenetically related to previously described STEC. The pangenome of our STEC collection was 11,650 genes, but we did not identify distinct molecular markers of local STEC. Our results showed that there may be local emerging STEC with unique features, nevertheless, no molecular markers were detected. Therefore, there might be elements such as a syntenic organization that might explain differential clustering detected between local and worldwide STEC.

Genomic analyses of multidrug-resistant Salmonella Indiana, Typhimurium, and Enteritidis isolates using MinION and MiSeq sequencing technologies

We sequenced 25 isolates of phenotypically multidrug-resistant Salmonella Indiana (n = 11), Typhimurium (n = 8), and Enteritidis (n = 6) using both MinION long-read [SQK-LSK109 and flow cell (R9.4.1)] and MiSeq short-read (Nextera XT and MiSeq Reagent Kit v2) sequencing technologies to determine the advantages of each approach in terms of the characteristics of genome structure, antimicrobial resistance (AMR), virulence potential, whole-genome phylogeny, and pan-genome. The MinION reads were base-called in real-time using MinKnow 3.4.8 integrated with Guppy 3.0.7. The long-read-only assembly, Illumina-only assembly, and hybrid assembly pipelines of Unicycler 0.4.8 were used to generate the MinION, MiSeq, and hybrid assemblies, respectively. The MinION assemblies were highly contiguous compared to the MiSeq assemblies but lacked accuracy, a deficiency that was mitigated by adding the MiSeq short reads through the Unicycler hybrid assembly which corrected erroneous single nucleotide polymorphisms (SNPs). The MinION assemblies provided similar predictions of AMR and virulence potential compared to the MiSeq and hybrid assemblies, although they produced more total false negatives of AMR genotypes, primarily due to failure in identifying tetracycline resistance genes in 11 of the 19 MinION assemblies of tetracycline-resistant isolates. The MinION assemblies displayed a large genetic distance from their corresponding MiSeq and hybrid assemblies on the whole-genome phylogenetic tree, indicating that the lower read accuracy of MinION sequencing caused incorrect clustering. The pan-genome of the MinION assemblies contained significantly more accessory genes and less core genes compared to the MiSeq and hybrid assemblies, suggesting that although these assemblies were more contiguous, their sequencing errors reduced accurate genome annotations. Our research demonstrates that MinION sequencing by itself provides an efficient assessment of the genome structure, antimicrobial resistance, and virulence potential of Salmonella; however, it is not sufficient for whole-genome phylogenetic and pan-genome analyses. MinION in combination with MiSeq facilitated the most accurate genomic analyses.

Nanopore sequencing for fast determination of plasmids, phages, virulence markers, and antimicrobial resistance genes in Shiga toxin-producing Escherichia coli

Whole genome sequencing can provide essential public health information. However, it is now known that widely used short-read methods have the potential to miss some randomly-distributed segments of genomes. This can prevent phages, plasmids, and virulence factors from being detected or properly identified. Here, we compared assemblies of three complete Shiga toxin-producing Escherichia coli (STEC) O26:H11/H- genomes from two different sequence types (ST21 and 29), each acquired using the Nextera XT MiSeq, MinION nanopore-based sequencing, and Pacific Biosciences (PacBio) sequencing. Each closed genome consisted of a single chromosome, approximately 5.7 Mb for CFSAN027343, 5.6 Mb for CFSAN027346, and 5.4 MB for CFSAN027350. However, short-read whole genome sequencing (WGS) using Nextera XT MiSeq failed to identify some virulence genes in plasmids and on the chromosome, both of which were detected using the long-read platforms. Results from long-read MinION and PacBio allowed us to identify differences in plasmid content: a single 88 kb plasmid in CFSAN027343; a 157kb plasmid in CFSAN027350; and two plasmids in CFSAN027346 (one 95 Kb, one 72 Kb). These data enabled rapid characterization of the virulome, detection of antimicrobial genes, and composition/location of Stx phages. Taken together, positive correlations between the two long-read methods for determining plasmids, virulome, antimicrobial resistance genes, and phage composition support MinION sequencing as one accurate and economical option for closing STEC genomes and identifying specific virulence markers.

Virulence gene profiles and phylogeny of Shiga toxin-positive Escherichia coli strains isolated from FDA regulated foods during 2010-2017

Illnesses caused by Shiga toxin-producing Escherichia coli (STECs) can be life threatening, such as hemolytic uremic syndrome (HUS). The STECs most frequently identified by USDA's Microbiological Data Program (MDP) carried toxin gene subtypes stx1a and/or stx2a. Here we described the genome sequences of 331 STECs isolated from foods regulated by the FDA 2010-2017, and determined their genomic identity, serotype, sequence type, virulence potential, and prevalence of antimicrobial resistance. Isolates were selected from the MDP archive, routine food testing by FDA field labs (ORA), and food testing by a contract company. Only 276 (83%) strains were confirmed as STECs by in silico analysis. Foods from which STECs were recovered included cilantro (6%), spinach (25%), lettuce (11%), and flour (9%). Phylogenetic analysis using core genome MLST revealed these STEC genomes were highly variable, with some clustering associated with ST types and serotypes. We detected 95 different sequence types (ST); several ST were previously associated with HUS: ST21 and ST29 (O26:H11), ST11 (O157:H7), ST33 (O91:H14), ST17 (O103:H2), and ST16 (O111:H-). in silico virulome analyses showed ~ 51% of these strains were potentially pathogenic [besides stx gene they also carried eae (25%) or 26% saa (26%)]. Virulence gene prevalence was also determined: stx1 only (19%); stx2 only (66%); and stx1/sxt2 (15%). Our data form a new WGS dataset that can be used to support food safety investigations and monitor the recurrence/emergence of E. coli in foods.

Genetic characterization of trh positive Vibrio spp. isolated from Norway

The thermostable direct hemolysin (TDH) and/or TDH-related hemolysin (TRH) genes are carried by most virulent Vibrio parahaemolyticus serovars. In Norway, trh+ V. parahaemolyticus constitute 4.4 and 4.5% of the total number of V. parahaemolyticus isolated from blue mussel (Mytilus edulis) and water, respectively. The trh gene is located in a region close to the gene cluster for urease production (ure). This region was characterized in V. parahaemolyticus strain TH3996 and it was found that a nickel transport operon (nik) was located between the first gene (ureR) and the rest of the ure cluster genes. The organization of the trh-ureR-nik-ure gene cluster in the Norwegian trh+ isolates was unknown. In this study, we explore the gene organization within the trh-ureR-nik-ure cluster for these isolates. PCR analyses revealed that the genes within the trh-ureR-nik-ure gene cluster of Norwegian trh+ isolates were organized in a similar fashion as reported previously for TH33996. Additionally, the phylogenetic relationship among these trh+ isolates was investigated using Multilocus Sequence Typing (MLST). Analysis by MLST or ureR-trh sequences generated two different phylogenetic trees for the same strains analyzed, suggesting that ureR-trh genes have been acquired at different times in Norwegian V. parahaemolyticus isolates. MLST results revealed that some pathogenic and non-pathogenic V. parahaemolyticus isolates in Norway appear to be highly genetically related.

Population structure of clinical and environmental Vibrio parahaemolyticus from the Pacific Northwest coast of the United States

Vibrio parahaemolyticus is a common marine bacterium and a leading cause of seafood-borne bacterial gastroenteritis worldwide. Although this bacterium has been the subject of much research, the population structure of cold-water populations remains largely undescribed. We present a broad phylogenetic analysis of clinical and environmental V. parahaemolyticus originating largely from the Pacific Northwest coast of the United States. Repetitive extragenic palindromic PCR (REP-PCR) separated 167 isolates into 39 groups and subsequent multilocus sequence typing (MLST) separated a subset of 77 isolates into 24 sequence types. The Pacific Northwest population exhibited a semi-clonal structure attributed to an environmental clade (ST3, N = 17 isolates) clonally related to the pandemic O3:K6 complex and a clinical clade (ST36, N = 20 isolates) genetically related to a regionally endemic O4:K12 complex. Further, the identification of at least five additional clinical sequence types (i.e., ST43, 50, 65, 135 and 417) demonstrates that V. parahaemolyticus gastroenteritis in the Pacific Northwest is polyphyletic in nature. Recombination was evident as a significant source of genetic diversity and in particular, the recA and dtdS alleles showed strong support for frequent recombination. Although pandemic-related illnesses were not documented during the study, the environmental occurrence of the pandemic clone may present a significant threat to human health and warrants continued monitoring. It is evident that V. parahaemolyticus population structure in the Pacific Northwest is semi-clonal and it would appear that multiple sequence types are contributing to the burden of disease in this region.

Phylogenetic evaluation of the 'Typhimurium' complex of Salmonella strains using a seven-gene multi-locus sequence analysis

Salmonella enterica comprises over 2500 serovars, many of which are significant foodborne pathogens in humans. The ability to subtype these microbes is difficult due to the highly clonal nature of many Salmonella strains and a lack of congruence among traditional typing approaches. This work examines the phylogenetic utility of a multi-locus sequence typing (MLST) approach to discriminate between members of a closely related collection of salmonellae, the Salmonella reference collection A (SARA). This 72 strain collection, referred to as the 'Typhimurium' complex, consists of S. Typhimurium and its four closest serological relatives. In this analysis, nucleotide sequences from seven housekeeping genes (aroC, dnaN, hemD, hisD, purE, sucA, and thrA) were PCR amplified, sequenced, and combined into a single concatenated character matrix providing 3360bp for cladistic analysis. The resultant most parsimonious tree yielded seven clades of Salmonella strains that partitioned largely along serovar divisions within the collection except for five 'Paratyphi B' strains, two 'Saintpaul' strains, and two 'Typhimurium' strains. Convergence in the SARA tree was approximately 20% indicating that the vast majority of sequence changes were phylogenetically informative. Despite a high consistency among nucleotide substitutions, analysis of congruence identified several SARA strains with recombinant alleles in the concatenated matrix. These findings point to important differences among phylogenetic contributions made by the individual genes comprising this MLST dataset.

Achromobacter marplatensis sp. nov., isolated from a pentachlorophenol-contaminated soil

A polyphasic taxonomic approach was applied to the study of a Gram-negative bacterium (B2(T)) isolated from soil by selective enrichment with pentachlorophenol. 16S rRNA gene sequence analysis of strain B2(T) showed that the strain belongs to the genus Achromobacter within the Betaproteobacteria. The 16S rRNA gene sequence displayed more than 99 % similarity to the sequences of the type strains of all species of Achromobacter, with the highest sequence similarity to those of Achromobacter spanius CCM 7183(T) and A. piechaudii CCM 2986(T) (99.8 %). On the basis of phylogenetic analysis, genomic DNA-DNA relatedness and phenotypic characteristics, including chemotaxonomic (cellular fatty acid profile) analysis, a novel species is proposed, Achromobacter marplatensis sp. nov., with the type strain B2(T) ( = CCM 7608(T)  = CCUG 56371(T)  = CECT 7342(T)).

Epidemiology of Vibrio parahaemolyticus outbreaks, southern Chile

One-sentence summary for table of contents: Outbreaks are decreasing and the O3:K6 pandemic strain is being replaced by a new serotype and new strains.

Disease outbreaks caused by Vibrio parahaemolyticus in Puerto Montt, Chile, began in 2004 and reached a peak in 2005 at 3,600 clinical cases. Until 2006, every analyzed case was caused by the serovar O3:K6 pandemic strain. In the summer of 2007, only 475 cases were reported; 73% corresponded to the pandemic strain. This decrease was associated with a change in serotype of many pandemic isolates to O3:K59 and the emergence of new clinical strains. One of these strains, associated with 11% of the cases, was genotypically different from the pandemic strain but contained genes that were identical to those found on its pathogenicity island. These findings suggest that pathogenicity-related genes were laterally transferred from the pandemic strain to one of the different V. parahaemolyticus groups comprising the diverse and shifting bacterial population in shellfish in this region.

Typing of Vibrio vulnificus strains by variability in their 16S-23S rRNA intergenic spacer regions

Amplification of the 16S-23S rDNA spacer region (ISR1) is a simple and rapid procedure for subtyping bacteria, especially those with several ribosomal operons including Vibrio vulnificus. V. vulnificus contains nine ribosomal operons with four or five ISR1 classes that differ in size and sequence. In the present study, 47 V. vulnificus strains of both shellfish and clinical origin were subtyped by their ISR1 patterns using "universal" primers, which target conserved sequences located in the 16S and the 23S rRNA genes. Sixteen different ISR1 patterns were observed that were grouped into two major clusters. Most (21/27, 77.8%) clinical isolates examined in this study grouped into a single cluster containing ISR1 patterns I, V, XI, and XII and these were highly similar (75%). This cluster was restricted to strains carrying the type B 16S rDNA (rrs) sequence which has been associated with human illness in previous studies. The remaining cluster consisted primarily of shellfish isolates. The highest variability in the ISR1 patterns was observed among shellfish isolates. Sequence analysis of the ISR1 region of selected strains demonstrated that all of them possess five ISR1 classes, with two "conserved sequence blocks" at the 5' and 3' end of the ISR1. All of these strains carried at least one tRNA gene and different classes differed in their tRNA gene composition. Some of the same ISR1 classes differed in size mainly due to an insertion of 35 bp in either of the conserved sequence blocks. These results demonstrate the feasibility of the ISR1 technique for V. vulnificus subtyping and suggest that ISR1 patterns appear to be linked to rrs sequence types and perhaps with virulence.

Comparison of direct genome restriction enzyme analysis and pulsed-field gel electrophoresis for typing of Vibrio vulnificus and their correspondence with multilocus sequence typing data

We compared the potential of direct genome restriction enzyme analysis (DGREA) and pulsed-field gel electrophoresis (PFGE) for discriminating Vibrio vulnificus isolates from clinical (23) and environmental (17) sources. The genotypes generated by both methodologies were compared to previous multilocus sequence typing (MLST) data. DGREA established clearer relationships among V. vulnificus strains and was more consistent with MLST than with PFGE. DGREA is a very promising tool for epidemiological and ecological studies of V. vulnificus.

Variation in the 16S-23S rRNA intergenic spacer regions in Vibrio parahaemolyticus strains are due to indels nearby their tRNAGlu

Vibrio parahaemolyticus contains 11 rRNA operons each including one of six 16S-23S rRNA gene intergenic spacer classes differing in size and nucleotide sequence. Some of the spacer classes may differ between isolates. We observed that the differences in the spacers between isolates are generally in two spacer classes present in single copies in the genome, one class containing tRNA(Ala) and tRNA(Glu) and the other tRNA(Glu) exclusively. Moreover, these differences are due to indels located nearby their tRNA(Glu). Comparison of the nucleotide sequence between spacer classes suggests that intragenomic nonreciprocal recombination causes the size variations observed in the spacer regions of V. parahaemolyticus strains.

Characterization of a Vibrio alginolyticus strain, isolated from Alaskan oysters, carrying a hemolysin gene similar to the thermostable direct hemolysin-related hemolysin gene (trh) of Vibrio parahaemolyticus

A Vibrio strain isolated from Alaskan oysters and classified by its biochemical characteristics as Vibrio alginolyticus possessed a thermostable direct hemolysin-related hemolysin (trh) gene previously reported only in Vibrio parahaemolyticus. This trh-like gene was cloned and sequenced and was 98% identical to the trh2 gene of V. parahaemolyticus. This gene seems to be functional since it was transcriptionally active in early-stationary-phase growing cells. To our knowledge, this is the first report of V. alginolyticus possessing a trh gene.

Quantitative reverse transcription polymerase chain reaction analysis of Vibrio cholerae cells entering the viable but non-culturable state and starvation in response to cold shock

We performed a comparative analysis of the Vibrio cholerae strain El Tor 3083 entering the viable but non-culturable (VBNC) state and starvation after incubation in artificial seawater (ASW) at 4 and 15 degrees C respectively. To this end, we determined bacterial culturability and membrane integrity, as well as the cellular levels of 16S rRNA and mRNA for the tuf, rpoS and relA genes, which were assessed by real-time quantitative reverse transcription polymerase chain reaction (Q-RT-PCR). Bacterial cells entering the VBNC state showed a 154, 5.1 x 10(3), 24- and 23-fold reduction in the number of copies of 16S rRNA and mRNA for tuf, rpoS and relA, in comparison to exponentially growing cells. The differences were less striking between cells in the VBNC and starvation states. The mRNA for relA was selectively increased in VBNC cells (3.2-folds), whereas a 3.9-fold reduction was observed for 16S rRNA. The obtained results confirmed that key activities of the cellular metabolism (i.e. tuf representing protein synthesis, and relA or rpoS stress response) were still detected in bacteria entering the VBNC state and starvation. These data suggest that the new Q-RT-PCR methodology, based on the selected RNA targets, could be successfully exploited for the identification (rRNA) of V. cholerae and assessment of its metabolic activity (tuf, rpoS, relA mRNA) in environmental samples.

Polymorphism and gene conversion of the 16S rRNA genes in the multiple rRNA operons of Vibrio parahaemolyticus

The genome sequence of a strain of Vibrio parahaemolyticus holds 11 copies of rRNA operons (rrn) with identical 16S rRNA genes (rrs). Conversely, the species type strain contains two rrs classes differing in 10 nucleotide sites within a short segment of 25 bp. Furthermore, we show here that the sequence of this particular segment largely differs between some strains of this species. We also show that of the eleven rrn operons in the species type strain, seven contain one rrs class and four the other, indicating gene conversion. Our results support the hypothesis that the rrs differences observed between strains of this species were caused by lateral transfer of an rrs segment and subsequent conversion.

Vibrio parahaemolyticus diarrhea, Chile, 1998 and 2004

Analysis of clinical isolates of Vibrio parahaemolyticus from outbreaks in Chile in the cities of Puerto Montt in 2004 and in Antofagasta in 1998 indicated that 23 of 24 isolates from Puerto Montt and 19 of 20 from Antofagasta belonged to the pandemic clonal complex that emerged in Southeast Asia in 1996.