Pangenomic analyses of antibiotic-resistant Campylobacter jejuni reveal unique lineage distributions and epidemiological associations

Application of whole-genome sequencing (WGS) to characterize foodborne pathogens has advanced our understanding of circulating genotypes and evolutionary relationships. Herein, we used WGS to investigate the genomic epidemiology of Campylobacter jejuni, a leading cause of foodborne disease. Among the 214 strains recovered from patients with gastroenteritis in Michigan, USA, 85 multilocus sequence types (STs) were represented and 135 (63.1 %) were phenotypically resistant to at least one antibiotic. Horizontally acquired antibiotic resistance genes were detected in 128 (59.8 %) strains and the genotypic resistance profiles were mostly consistent with the phenotypes. Core-gene phylogenetic reconstruction identified three sequence clusters that varied in frequency, while a neighbour-net tree detected significant recombination among the genotypes (pairwise homoplasy index P<0.01). Epidemiological analyses revealed that travel was a significant contributor to pangenomic and ST diversity of C. jejuni, while some lineages were unique to rural counties and more commonly possessed clinically important resistance determinants. Variation was also observed in the frequency of lineages over the 4 year period with chicken and cattle specialists predominating. Altogether, these findings highlight the importance of geographically specific factors, recombination and horizontal gene transfer in shaping the population structure of C. jejuni. They also illustrate the usefulness of WGS data for predicting antibiotic susceptibilities and surveillance, which are important for guiding treatment and prevention strategies.

Whole-Genome Sequencing of Shiga Toxin-Producing Escherichia coli for Characterization and Outbreak Investigation

Shiga toxin-producing Escherichia coli (STEC) causes high frequencies of foodborne infections worldwide and has been linked to numerous outbreaks each year. Pulsed-field gel electrophoresis (PFGE) has been the gold standard for surveillance until the recent transition to whole-genome sequencing (WGS). To further understand the genetic diversity and relatedness of outbreak isolates, a retrospective analysis of 510 clinical STEC isolates was conducted. Among the 34 STEC serogroups represented, most (59.6%) belonged to the predominant six non-O157 serogroups. Core genome single nucleotide polymorphism (SNP) analysis differentiated clusters of isolates with similar PFGE patterns and multilocus sequence types (STs). One serogroup O26 outbreak strain and another non-typeable (NT) strain, for instance, were identical by PFGE and clustered together by MLST; however, both were distantly related in the SNP analysis. In contrast, six outbreak-associated serogroup O5 strains clustered with five ST-175 serogroup O5 isolates, which were not part of the same outbreak as determined by PFGE. The use of high-quality SNP analyses enhanced the discrimination of these O5 outbreak strains into a single cluster. In all, this study demonstrates how public health laboratories can more rapidly use WGS and phylogenetics to identify related strains during outbreak investigations while simultaneously uncovering important genetic attributes that can inform treatment practices.

Long read genome assemblers struggle with small plasmids

Whole-genome sequencing has become a preferred method for studying bacterial plasmids, as it is generally assumed to capture the entire genome. However, long-read genome assemblers have been shown to sometimes miss plasmid sequences - an issue that has been associated with plasmid size. The purpose of this study was to investigate the relationship between plasmid size and plasmid recovery by the long-read-only assemblers, Flye, Raven, Miniasm, and Canu. This was accomplished by determining the number of times each assembler successfully recovered 33 plasmids, ranging from 1919 to 194 062 bp in size and belonging to 14 bacterial isolates from six bacterial genera, using Oxford Nanopore long reads. These results were additionally compared to plasmid recovery rates by the short-read-first assembler, Unicycler, using both Oxford Nanopore long reads and Illumina short reads. Results from this study indicate that Canu, Flye, Miniasm, and Raven are prone to missing plasmid sequences, whereas Unicycler was successful at recovering 100 % of plasmid sequences. Excluding Canu, most plasmid loss by long-read-only assemblers was due to failure to recover plasmids smaller than 10 kb. As such, it is recommended that Unicycler be used to increase the likelihood of plasmid recovery during bacterial genome assembly.

Evaluation of a laboratory-developed multiplex real-time PCR assay for diagnosis of syphilis, herpes and chancroid genital ulcers in four public health laboratories in the USA

OBJECTIVE

To evaluate the field performance of a multiplex PCR (M-PCR) assay for detection of herpes simplex virus (HSV)-1 and HSV-2, Treponema pallidum (T. pallidum) and Haemophilus ducreyi (H. ducreyi) in genital ulcer disease (GUD) specimens.

METHODS

GUD M-PCR was performed on 186 remnant specimens, previously collected for HSV testing, by four public health laboratories (PHLs) and the Laboratory Reference and Research Branch (LRRB) at the Centers for Disease Control and Prevention. The results from the PHLs were compared with those of LRRB, which served as the reference testing method, and percentage agreement was calculated.

RESULTS

HSV was detected in 31 of 52 (59.6%), 20 of 40 (50%), 43 of 44 (97.7%) and 19 of 50 (38.0%) specimens from PHL1, PHL2, PHL3 and PHL4, respectively. There were seven discrepant results for HSV, and the overall percent agreement between the PHLs and the LRRB was 94%-100%, with a kappa value of 0.922, which demonstrates high agreement. T. pallidum was identified in 7 of 51 (13.7%) specimens from PHL1 with 94.1% agreement and in 2 of 40 (5.0%) specimens from PHL2 with 100% agreement. The LRRB identified three additional T. pallidum-positive specimens from PHL1. The kappa value (0.849) for T. pallidum testing suggests good agreement. Consistent with the LRRB results, no T. pallidum was detected in specimens from PHL3 and PHL4, and H. ducreyi was not detected at any of the study sites.

CONCLUSIONS

The GUD M-PCR assay performed well in four independent PHLs and 12 suspected syphilis cases were identified in this study. The M-PCR assay could provide improved diagnostic options for GUD infections in state and local PHLs.

Population structure and genetic diversity of non-O157 Shiga toxin-producing Escherichia coli (STEC) clinical isolates from Michigan

Non-O157 STEC are increasingly linked to foodborne infections, yet little is known about the diversity and molecular epidemiology across locations. Herein, we used whole genome sequencing to examine genetic variation in 894 isolates collected from Michigan patients between 2001 and 2018. In all, 67 serotypes representing 69 multilocus sequence types were identified. Serotype diversity increased from an average of four (2001-2006) to 17 (2008-2018) serotypes per year. The top six serogroups reported nationally caused > 60% of infections in 16 of the 18 years; serogroups O111 and O45 were associated with hospitalization as were age ≥ 65 years, diarrhea with blood and female sex. Phylogenetic analyses of seven multilocus sequence typing (MLST) loci identified three clades as well as evidence of parallel evolution and recombination. Most (95.5%) isolates belonged to one clade, which could be further differentiated into seven subclades comprising isolates with varying virulence gene profiles and serotypes. No association was observed between specific clades and the epidemiological data, suggesting that serogroup- and serotype-specific associations are more important predictors of disease outcomes than lineages defined by MLST. Molecular epidemiological studies of non-O157 STEC are important to enhance understanding of circulating strain distributions and traits, genetic variation, and factors that may impact disease risk and severity.

Shigella sonnei Outbreak Investigation During a Municipal Water Crisis-Genesee and Saginaw Counties, Michigan, 2016

Objectives. To investigate a shigellosis outbreak in Genesee County, Michigan (including the City of Flint), and Saginaw County, Michigan, in 2016 and address community concerns about the role of the Flint water system.Methods. We met frequently with community members to understand concerns and develop the investigation. We surveyed households affected by the outbreak, analyzed Shigella isolate data, examined the geospatial distribution of cases, and reviewed available water quality data.Results. We surveyed 83 households containing 158 cases; median age was 10 years. Index case-patients from 55 of 83 households (66%) reported contact with a person outside their household who wore diapers or who had diarrhea in the week before becoming ill; results were similar regardless of household drinking water source. Genomic diversity was not consistent with a point source. In Flint, no space-time clustering was identified, and average free chlorine residual values remained above recommended levels throughout the outbreak period.Conclusions. The outbreak was most likely caused by person-to-person contact and not by the Flint water system. Consistent community engagement was essential to the design and implementation of the investigation.

Onsite investigation at a mail-order hatchery following a multistate Salmonella illness outbreak linked to live poultry-United States, 2018

Centers for Disease Control and Prevention (CDC), health departments, and other state and federal partners have linked contact with live poultry to 70 human Salmonella outbreaks in the United States from 2000 to 2017, which resulted in a total of 4,794 illnesses, 894 hospitalizations, and 7 deaths. During human salmonellosis outbreaks environmental sampling is rarely conducted as part of the outbreak investigation. CDC was contacted by state health officials on June 12, 2018, to provide support during an investigation of risk factors for Salmonella infections linked to live poultry originating at a mail-order hatchery. From January 1, 2018, to June 15, 2018, 13 human Salmonella infections in multiple states were attributed to exposure to live poultry from a single hatchery. Two serotypes of Salmonella were associated with these infections, Salmonella Enteritidis and Salmonella Litchfield. Molecular subtyping of the S. Enteritidis clinical isolates revealed they were closely related genetically (within 0 to 9 alleles) by core genome multi-locus sequence typing (cgMLST) to isolates obtained from environmental samples taken from hatchery shipping containers received at retail outlets. Environmental sampling and onsite investigation of practices was conducted at the mail-order hatchery during an investigation on June 19, 2018. A total of 45 environmental samples were collected, and 4 (9%) grew Salmonella. A chick box liner from a box in the pre-shipping area yielded an isolate closely related to the S. Enteritidis outbreak strain (within 1 to 9 alleles by cgMLST). The onsite investigation revealed lapses in biosecurity, sanitation, quality assurance, and education of consumers. Review of Salmonella serotype testing performed by the hatchery revealed that the number of samples and type of samples collected monthly varied. Also, S. Enteritidis was identified at the hatchery every year since testing began in 2016. Recommendations to the hatchery for biosecurity, testing, and sanitation measures were made to help reduce burden of Salmonella in the hatchery and breeding flocks, thereby reducing the occurrence of human illness.

Insights from the draft genome into the pathogenicity of a clinical isolate of Elizabethkingia meningoseptica Em3

Elizabethkingia meningoseptica is an emerging, healthcare-associated pathogen causing a high mortality rate in immunocompromised patients. We report the draft genome sequence of E. meningoseptica Em3, isolated from sputum from a patient with multiple underlying diseases. The genome has a length of 4,037,922 bp, a GC-content 36.4%, and 3673 predicted protein-coding sequences. Average nucleotide identity analysis (>95%) assigned the bacterium to the species E. meningoseptica. Genome analysis showed presence of the curli formation and assembly operon and a gene encoding hemagglutinins, indicating ability to form biofilm. In vitro biofilm assays demonstrated that E. meningoseptica Em3 formed more biofilm than E. anophelis Ag1 and E. miricola Emi3, both lacking the curli operon. A gene encoding thiol-activated cholesterol-dependent cytolysin in E. meningoseptica Em3 (potentially involved in lysing host immune cells) was also absent in E. anophelis Ag1 and E. miricola Emi3. Strain Em3 showed α-hemolysin activity on blood agar medium, congruent with presence of hemolysin and cytolysin genes. Furthermore, presence of heme uptake and utilization genes demonstrated adaptations for bloodstream infections. Strain Em3 contained 12 genes conferring resistance to β-lactams, including β-lactamases class A, class B, and metallo-β-lactamases. Results of comparative genomic analysis here provide insights into the evolution of E. meningoseptica Em3 as a pathogen.