Morag Graham

High-throughput genome sequencing of two Listeria monocytogenes clinical isolates during a large foodborne outbreak.

BackgroundA large, multi-province outbreak of listeriosis associated with ready-to-eat meat products contaminated with Listeria monocytogenes serotype 1/2a occurred in Canada in 2008. Subtyping of outbreak-associated isolates using pulsed-field gel electrophoresis (PFGE) revealed two similar but distinct Asc I PFGE patterns. High-throughput pyrosequencing of two L. monocytogenes isolates was used to rapidly provide the genome sequence of the primary outbreak strain and to investigate the extent of genetic diversity associated with a change of a single restriction enzyme fragment during PFGE.ResultsThe chromosomes were collinear, but differences included 28 single nucleotide polymorphisms (SNPs) and three indels, including a 33 kbp prophage that accounted for the observed difference in Asc I PFGE patterns. The distribution of these traits was assessed within further clinical, environmental and food isolates associated with the outbreak, and this comparison indicated that three distinct, but highly related strains may have been involved in this nationwide outbreak. Notably, these two isolates were found to harbor a 50 kbp putative mobile genomic island encoding translocation and efflux functions that has not been observed in other Listeria genomes.ConclusionsHigh-throughput genome sequencing provided a more detailed real-time assessment of genetic traits characteristic of the outbreak strains than could be achieved with routine subtyping methods. This study confirms that the latest generation of DNA sequencing technologies can be applied during high priority public health events, and laboratories need to prepare for this inevitability and assess how to properly analyze and interpret whole genome sequences in the context of molecular epidemiology.

Comparative Genomics of Vibrio cholerae from Haiti, Asia, and Africa

A strain from Haiti shares genetic ancestry with those from Asia and Africa.

Evolutionary dynamics of Vibrio cholerae O1 following a single-source introduction to Haiti

ABSTRACT Prior to the epidemic that emerged in Haiti in October of 2010, cholera had not been documented in this country. After its introduction, a strain of Vibrio cholerae O1 spread rapidly throughout Haiti, where it caused over 600,000 cases of disease and >7,500 deaths in the first two years of the epidemic. We applied whole-genome sequencing to a temporal series of V. cholerae isolates from Haiti to gain insight into the mode and tempo of evolution in this isolated population of V. cholerae O1. Phylogenetic and Bayesian analyses supported the hypothesis that all isolates in the sample set diverged from a common ancestor within a time frame that is consistent with epidemiological observations. A pangenome analysis showed nearly homogeneous genomic content, with no evidence of gene acquisition among Haiti isolates. Nine nearly closed genomes assembled from continuous-long-read data showed evidence of genome rearrangements and supported the observation of no gene acquisition among isolates. Thus, intrinsic mutational processes can account for virtually all of the observed genetic polymorphism, with no demonstrable contribution from horizontal gene transfer (HGT). Consistent with this, the 12 Haiti isolates tested by laboratory HGT assays were severely impaired for transformation, although unlike previously characterized noncompetent V. cholerae isolates, each expressed hapR and possessed a functional quorum-sensing system. Continued monitoring of V. cholerae in Haiti will illuminate the processes influencing the origin and fate of genome variants, which will facilitate interpretation of genetic variation in future epidemics. IMPORTANCE Vibrio cholerae is the cause of substantial morbidity and mortality worldwide, with over three million cases of disease each year. An understanding of the mode and rate of evolutionary change is critical for proper interpretation of genome sequence data and attribution of outbreak sources. The Haiti epidemic provides an unprecedented opportunity to study an isolated, single-source outbreak of Vibrio cholerae O1 over an established time frame. By using multiple approaches to assay genetic variation, we found no evidence that the Haiti strain has acquired any genes by horizontal gene transfer, an observation that led us to discover that it is also poorly transformable. We have found no evidence that environmental strains have played a role in the evolution of the outbreak strain. Vibrio cholerae is the cause of substantial morbidity and mortality worldwide, with over three million cases of disease each year. An understanding of the mode and rate of evolutionary change is critical for proper interpretation of genome sequence data and attribution of outbreak sources. The Haiti epidemic provides an unprecedented opportunity to study an isolated, single-source outbreak of Vibrio cholerae O1 over an established time frame. By using multiple approaches to assay genetic variation, we found no evidence that the Haiti strain has acquired any genes by horizontal gene transfer, an observation that led us to discover that it is also poorly transformable. We have found no evidence that environmental strains have played a role in the evolution of the outbreak strain.

Livestock-associated Methicillin- Resistant Staphylococcus aureus Sequence Type 398 in Humans, Canada

Recent emergence of infections resulting from this strain is of public health concern.

Sequence Typing Confirms that a Predominant Listeria monocytogenes Clone Caused Human Listeriosis Cases and Outbreaks in Canada from 1988 to 2010

ABSTRACT Human listeriosis outbreaks in Canada have been predominantly caused by serotype 1/2a isolates with highly similar pulsed-field gel electrophoresis (PFGE) patterns. Multilocus sequence typing (MLST) and multi-virulence-locus sequence typing (MVLST) each identified a diverse population of Listeria monocytogenes isolates, and within that, both methods had congruent subtypes that substantiated a predominant clone (clonal complex 8; virulence type 59; proposed epidemic clone 5 [ECV]) that has been causing human illness across Canada for more than 2 decades.

Molecular definition of vaginal microbiota in East African commercial sex workers.

ABSTRACT Resistance to HIV infection in a cohort of commercial sex workers living in Nairobi, Kenya, is linked to mucosal and antiinflammatory factors that may be influenced by the vaginal microbiota. Since bacterial vaginosis (BV), a polymicrobial dysbiosis characterized by low levels of protective Lactobacillus organisms, is an established risk factor for HIV infection, we investigated whether vaginal microbiology was associated with HIV-exposed seronegative (HESN) or HIV-seropositive (HIV+) status in this cohort. A subset of 44 individuals was selected for deep-sequencing analysis based on the chaperonin 60 (cpn60) universal target (UT), including HESN individuals (n = 16), other HIV-seronegative controls (HIV-N, n = 16), and HIV+ individuals (n = 12). Our findings indicate exceptionally high phylogenetic resolution of the cpn60 UT using reads as short as 200 bp, with 54 species in 29 genera detected in this group. Contrary to our initial hypothesis, few differences between HESN and HIV-N women were observed. Several HIV+ women had distinct profiles dominated by Escherichia coli. The deep-sequencing phylogenetic profile of the vaginal microbiota corresponds closely to BV+ and BV− diagnoses by microscopy, elucidating BV at the molecular level. A cluster of samples with intermediate abundance of Lactobacillus and dominant Gardnerella was identified, defining a distinct BV phenotype that may represent a transitional stage between BV+ and BV−. Several alpha- and betaproteobacteria, including the recently described species Variovorax paradoxus, were found to correlate positively with increased Lactobacillus levels that define the BV− (“normal”) phenotype. We conclude that cpn60 UT is ideally suited to next-generation sequencing technologies for further investigation of microbial community dynamics and mucosal immunity underlying HIV resistance in this cohort.

Whole-Genome Phylogenomic Heterogeneity of Neisseria gonorrhoeae Isolates with Decreased Cephalosporin Susceptibility Collected in Canada between 1989 and 2013

ABSTRACT A large-scale, whole-genome comparison of Canadian Neisseria gonorrhoeae isolates with high-level cephalosporin MICs was used to demonstrate a genomic epidemiology approach to investigate strain relatedness and dynamics. Although current typing methods have been very successful in tracing short-chain transmission of gonorrheal disease, investigating the temporal evolutionary relationships and geographical dissemination of highly clonal lineages requires enhanced resolution only available through whole-genome sequencing (WGS). Phylogenomic cluster analysis grouped 169 Canadian strains into 12 distinct clades. While some N. gonorrhoeae multiantigen sequence types (NG-MAST) agreed with specific phylogenomic clades or subclades, other sequence types (ST) and closely related groups of ST were widely distributed among clades. Decreased susceptibility to extended-spectrum cephalosporins (ESC-DS) emerged among a group of diverse strains in Canada during the 1990s with a variety of nonmosaic penA alleles, followed in 2000/2001 with the penA mosaic X allele and then in 2007 with ST1407 strains with the penA mosaic XXXIV allele. Five genetically distinct ESC-DS lineages were associated with penA mosaic X, XXXV, and XXXIV alleles and nonmosaic XII and XIII alleles. ESC-DS with coresistance to azithromycin was observed in 5 strains with 23S rRNA C2599T or A2143G mutations. As the costs associated with WGS decline and analysis tools are streamlined, WGS can provide a more thorough understanding of strain dynamics, facilitate epidemiological studies to better resolve social networks, and improve surveillance to optimize treatment for gonorrheal infections.

Characterization of a stable, metronidazole-resistant Clostridium difficile clinical isolate.

Background Clostridium difficile are Gram-positive, spore forming anaerobic bacteria that are the leading cause of healthcare-associated diarrhea, usually associated with antibiotic usage. Metronidazole is currently the first-line treatment for mild to moderate C. difficile diarrhea however recurrence occurs at rates of 15–35%. There are few reports of C. difficile metronidazole resistance in the literature, and when observed, the phenotype has been transient and lost after storage or exposure of the bacteria to freeze/thaw cycles. Owing to the unstable nature of the resistance phenotype in the laboratory, clinical significance and understanding of the resistance mechanisms is lacking. Methodology/Principal Findings Genotypic and phenotypic characterization was performed on a metronidazole resistant clinical isolate of C. difficile. Whole-genome sequencing was used to identify potential genetic contributions to the phenotypic variation observed with molecular and bacteriological techniques. Phenotypic observations of the metronidazole resistant strain revealed aberrant growth in broth and elongated cell morphology relative to a metronidazole-susceptible, wild type NAP1 strain. Comparative genomic analysis revealed single nucleotide polymorphism (SNP) level variation within genes affecting core metabolic pathways such as electron transport, iron utilization and energy production. Conclusions/Significance This is the first characterization of stable, metronidazole resistance in a C. difficile isolate. The study provides an in-depth genomic and phenotypic analysis of this strain and provides a foundation for future studies to elucidate mechanisms conferring metronidazole resistance in C. difficile that have not been previously described.

Genomic Epidemiology and Molecular Resistance Mechanisms of Azithromycin-Resistant Neisseria gonorrhoeae in Canada from 1997 to 2014

ABSTRACT The emergence of Neisseria gonorrhoeae strains with decreased susceptibility to cephalosporins and azithromycin (AZM) resistance (AZMr) represents a public health threat of untreatable gonorrhea infections. Genomic epidemiology through whole-genome sequencing was used to describe the emergence, dissemination, and spread of AZMr strains. The genomes of 213 AZMr and 23 AZM-susceptible N. gonorrhoeae isolates collected in Canada from 1989 to 2014 were sequenced. Core single nucleotide polymorphism (SNP) phylogenomic analysis resolved 246 isolates into 13 lineages. High-level AZMr (MICs ≥ 256 μg/ml) was found in 5 phylogenetically diverse isolates, all of which possessed the A2059G mutation (Escherichia coli numbering) in all four 23S rRNA alleles. One isolate with high-level AZMr collected in 2009 concurrently had decreased susceptibility to ceftriaxone (MIC = 0.125 μg/ml). An increase in the number of 23S rRNA alleles with the C2611T mutations (E. coli numbering) conferred low to moderate levels of AZMr (MICs = 2 to 4 and 8 to 32 μg/ml, respectively). Low-level AZMr was also associated with mtrR promoter mutations, including the −35A deletion and the presence of Neisseria meningitidis-like sequences. Geographic and temporal phylogenetic clustering indicates that emergent AZMr strains arise independently and can then rapidly expand clonally in a region through local sexual networks.

Nationwide Molecular Surveillance of Pandemic H1N1 Influenza A Virus Genomes: Canada, 2009

Background In April 2009, a novel triple-reassortant swine influenza A H1N1 virus (“A/H1N1pdm”; also known as SOIV) was detected and spread globally as the first influenza pandemic of the 21st century. Sequencing has since been conducted at an unprecedented rate globally in order to monitor the diversification of this emergent virus and to track mutations that may affect virus behavior. Methodology/Principal Findings By Sanger sequencing, we determined consensus whole-genome sequences for A/H1N1pdm viruses sampled nationwide in Canada over 33 weeks during the 2009 first and second pandemic waves. A total of 235 virus genomes sampled from unique subjects were analyzed, providing insight into the temporal and spatial trajectory of A/H1N1pdm lineages within Canada. Three clades (2, 3, and 7) were identifiable within the first two weeks of A/H1N1pdm appearance, with clades 5 and 6 appearing thereafter; further diversification was not apparent. Only two viral sites displayed evidence of adaptive evolution, located in hemagglutinin (HA) corresponding to D222 in the HA receptor-binding site, and to E374 at HA2-subunit position 47. Among the Canadian sampled viruses, we observed notable genetic diversity (1.47×10−3 amino acid substitutions per site) in the gene encoding PB1, particularly within the viral genomic RNA (vRNA)-binding domain (residues 493–757). This genome data set supports the conclusion that A/H1N1pdm is evolving but not excessively relative to other H1N1 influenza A viruses. Entropy analysis was used to investigate whether any mutated A/H1N1pdm protein residues were associated with infection severity; however no virus genotypes were observed to trend with infection severity. One virus that harboured heterozygote coding mutations, including PB2 D567D/G, was attributed to a severe and potentially mixed infection; yet the functional significance of this PB2 mutation remains unknown. Conclusions/Significance These findings contribute to enhanced understanding of Influenza A/H1N1pdm viral dynamics.