Comparison of advanced whole genome sequence-based methods to distinguish strains of Salmonella enterica serovar Heidelberg involved in foodborne outbreaks in Québec.

Salmonella enterica serovar Heidelberg (S. Heidelberg) is one of the top serovars causing human salmonellosis. This serovar ranks second and third among serovars that cause human infections in Québec and Canada, respectively, and has been associated with severe infections. Traditional typing methods such as PFGE do not display adequate discrimination required to resolve outbreak investigations due to the low level of genetic diversity of isolates belonging to this serovar. This study evaluates the ability of four whole genome sequence (WGS)-based typing methods to differentiate among 145 S. Heidelberg strains involved in four distinct outbreak events and sporadic cases of salmonellosis that occurred in Québec between 2007 and 2016. Isolates from all outbreaks were indistinguishable by PFGE. The core genome single nucleotide variant (SNV), core genome multilocus sequence typing (MLST) and whole genome MLST approaches were highly discriminatory and separated outbreak strains into four distinct phylogenetic clusters that were concordant with the epidemiological data. The clustered regularly interspaced short palindromic repeats (CRISPR) typing method was less discriminatory. However, CRISPR typing may be used as a secondary method to differentiate isolates of S. Heidelberg that are genetically similar but epidemiologically unrelated to outbreak events. WGS-based typing methods provide a highly discriminatory alternative to PFGE for the laboratory investigation of foodborne outbreaks.

Roles of type 1A topoisomerases in genome maintenance in Escherichia coli.

In eukaryotes, type 1A topoisomerases (topos) act with RecQ-like helicases to maintain the stability of the genome. Despite having been the first type 1A enzymes to be discovered, much less is known about the involvement of the E. coli topo I (topA) and III (topB) enzymes in genome maintenance. These enzymes are thought to have distinct cellular functions: topo I regulates supercoiling and R-loop formation, and topo III is involved in chromosome segregation. To better characterize their roles in genome maintenance, we have used genetic approaches including suppressor screens, combined with microscopy for the examination of cell morphology and nucleoid shape. We show that topA mutants can suffer from growth-inhibitory and supercoiling-dependent chromosome segregation defects. These problems are corrected by deleting recA or recQ but not by deleting recJ or recO, indicating that the RecF pathway is not involved. Rather, our data suggest that RecQ acts with a type 1A topo on RecA-generated recombination intermediates because: 1-topo III overproduction corrects the defects and 2-recQ deletion and topo IIII overproduction are epistatic to recA deletion. The segregation defects are also linked to over-replication, as they are significantly alleviated by an oriC::aph suppressor mutation which is oriC-competent in topA null but not in isogenic topA+ cells. When both topo I and topo III are missing, excess supercoiling triggers growth inhibition that correlates with the formation of extremely long filaments fully packed with unsegregated and diffuse DNA. These phenotypes are likely related to replication from R-loops as they are corrected by overproducing RNase HI or by genetic suppressors of double topA rnhA mutants affecting constitutive stable DNA replication, dnaT::aph and rne::aph, which initiates from R-loops. Thus, bacterial type 1A topos maintain the stability of the genome (i) by preventing over-replication originating from oriC (topo I alone) and R-loops and (ii) by acting with RecQ.

Interplay between type 1A topoisomerases and gyrase in chromosome segregation in Escherichia coli.

Escherichia coli possesses two type 1A topoisomerases, Topo I (topA) and Topo III (topB). Topo I relaxes excess negative supercoiling, and topA mutants can grow only in the presence of compensatory mechanisms, such as gyrase mutations. topB mutants grow as well as wild-type cells. In vitro, Topo III, but not Topo I, can efficiently decatenate DNA during replication. However, in vivo, a chromosome segregation defect is seen only when both type 1A topoisomerases are absent. Here we present experimental evidence for an interplay between gyrase and type 1A topoisomerases in chromosome segregation. We found that both the growth defect and the Par(-) phenotypes of a gyrB(Ts) mutant at nonpermissive temperatures were significantly corrected by deleting topA, but only when topB was present. Overproducing Topo IV, the major cellular decatenase, could not substitute for topB. We also show that overproducing Topo III at a very high level could suppress the Par(-) phenotype. We previously found that the growth and chromosome segregation defects of a triple topA rnhA gyrB(Ts) mutant in which gyrase supercoiling activity was strongly inhibited could be corrected by overproducing Topo III (V. Usongo, F. Nolent, P. Sanscartier, C. Tanguay, S. Broccoli, I. Baaklini, K. Drlica, and M. Drolet, Mol. Microbiol. 69:968-981, 2008). We show here that this overproduction could be bypassed by substituting the gyrB(Ts) allele for a gyrB(+) one or by growing cells in a minimal medium, conditions that reduced both topA- and rnhA-dependent unregulated replication. Altogether, our data point to a role for Topo III in chromosome segregation when gyrase is inefficient and suggest that Topo I plays an indirect role via supercoiling regulation.

Analysis of an IncR Plasmid Carrying blaNDM-1 Linked to an Azithromycin Resistance Region in Enterobacter hormaechei Involved in an Outbreak in Quebec.

In the context of a recent rise in prevalence of NDM-encoding carbapenemase-producing Enterobacterales (CPE) in the province of QC, Canada, the genetic environment of blaNDM-1 was investigated. Three NDM-producing clinical isolates of Enterobacter hormaechei recovered from hospitalized patients involved in a putative outbreak were further characterized by whole-genome sequencing (WGS). Two isolates were confirmed by pulsed-field gel electrophoresis and WGS to be closely related. In addition to a ∼128 kb IncFII conjugative multidrug-resistance (MDR) plasmid, these isolates possessed a ∼45 kb mobilizable IncR MDR plasmid containing 2 MDR regions: a complex class 1 integron harboring blaNDM-1 and 7 other AMR genes, and the IS26-mph(A)-mrx-mphR(A)-IS6100 azithromycin resistance unit. The predicted antimicrobial resistance (AMR) genes correlated with the antimicrobial susceptibility testing results. The multidrug-resistant phenotype in addition to the presence of two important mobile genetic elements, suggest a potent role as a reservoir of antibiotic resistance for such a small IncR plasmid. IMPORTANCE Analyzing the genetic environment of clinically relevant MDR genes can provide information on the way in which such genes are maintained and disseminated. Understanding this phenomenon is of interest for clinicians as it can also provide insight on where these genes might have been sourced, possibly supporting outbreak investigations.

Prevalence and characteristics of Livestock-Associated Methicillin-Resistant Staphylococcus aureus (LA-MRSA) isolated from chicken meat in the province of Quebec, Canada.

This study was conducted to estimate the prevalence of Livestock-Associated Methicillin-Resistant Staphylococcus aureus (LA-MRSA) in retail chicken meat and broiler chickens from the Province of Quebec, Canada, and to characterize LA-MRSA isolates. A total of 309 chicken drumsticks and thighs were randomly selected in 2013 from 43 retail stores in the Monteregie. In addition, nasal swabs and caeca samples were collected in 2013-2014 from 200 broiler chickens of 38 different flocks. LA-MRSA was not detected in broiler chickens. Fifteen LA-MRSA isolates were recovered from four (1.3%) of the 309 chicken meat samples. Multi-Locus Sequence Typing (MLST) and SCCmec typing revealed two profiles (ST398-MRSA-V and ST8-MRSA-IVa), which were distinct using pulse-field gel electrophoresis (PFGE) and microarray (antimicrobial resistance and virulence genes) analyses. In addition to beta-lactam resistance, tetracycline and spectinomycin resistance was detected in all isolates from the 3 positive samples of the ST398 profile. Southern blot hybridization revealed that the resistance genes aad(D) and lnu(A), encoding resistances to aminoglycosides and lincosamides respectively, were located on plasmid. All isolates were able to produce biofilms, but biofilm production was not correlated with hld gene expression. Our results show the presence of two separate lineages of MRSA in retail chicken meat in Quebec, one of which is likely of human origin.

Source Tracking Based on Core Genome SNV and CRISPR Typing of Salmonella enterica Serovar Heidelberg Isolates Involved in Foodborne Outbreaks in Québec, 2012.

Whole-genome sequencing (WGS) is the method of choice for bacterial subtyping and it is rapidly replacing the more traditional methods such as pulsed-field gel electrophoresis (PFGE). Here we used the high-resolution core genome single nucleotide variant (cgSNV) typing method to characterize clinical and food from Salmonella enterica serovar Heidelberg isolates in the context of source attribution. Additionally, clustered regularly interspaced short palindromic repeats (CRISPR) analysis was included to further support this method. Our results revealed that cgSNV was highly discriminatory and separated the outbreak isolates into distinct clusters (0-4 SNVs). CRISPR analysis was also able to distinguish outbreak strains from epidemiologically unrelated isolates. Specifically, our data clearly demonstrated the strength of these two methods to determine the probable source(s) of a 2012 epidemiologically characterized outbreak of S. Heidelberg. Using molecular cut-off of 0-10 SNVs, the cgSNV analysis of 246 clinical and food isolates of S. Heidelberg collected in Québec, in the same year of the outbreak event, revealed that retail and abattoir chicken isolates likely represent an important source of human infection to S. Heidelberg. Interestingly, the isolates genetically related by cgSNV also harbored the same CRISPR as outbreak isolates and clusters. This indicates that CRISPR profiles can be useful as a complementary approach to determine source attribution in foodborne outbreaks. Use of the genomic analysis also allowed to identify a large number of cases that were missed by PFGE, indicating that most outbreaks are probably underestimated. Although epidemiological information must still support WGS-based results, cgSNV method is a highly discriminatory method for the resolution of outbreak events and the attribution of these events to their respective sources. CRISPR typing can serve as a complimentary tool to this analysis during source tracking.

Depletion of RNase HI activity in Escherichia coli lacking DNA topoisomerase I leads to defects in DNA supercoiling and segregation.

Gyrase-mediated hypernegative supercoiling is one manifestation of R-loop formation, a phenomenon that is normally suppressed by topoisomerase I (topA) in Escherichia coli. Overproduction of RNase HI (rnhA), an enzyme that removes the RNA moiety of R-loops, prevents hypernegative supercoiling and allows growth of topA null mutants. We previously showed that topA and rnhA null mutations are incompatible. We now report that such mutants were viable when RNase HI or topoisomerase III was expressed from a plasmid-borne gene. Surprisingly, DNA of topA null mutants became relaxed rather than hypernegatively supercoiled following depletion of RNase HI activity. This result failed to correlate with the cellular concentration of gyrase or topoisomerase IV (the other relaxing enzyme in the cell) or with transcription-induced supercoiling. Rather, intracellular DNA relaxation in the absence of RNase HI was related to inhibition of gyrase activity both in vivo and in extracts. Cells lacking topA and rnhA also exhibited properties consistent with segregation defects. Overproduction of topoisomerase III, an enzyme that can carry out DNA decatenation, corrected the segregation defects without restoring supercoiling activity. Collectively these data reveal (i) the existence of a cellular response to loss of RNase HI that counters the supercoiling activity of gyrase, and (ii) supercoiling-independent segregation defects due to loss of RNase HI from topA null mutants. Thus RNase HI plays a more central role in DNA topology than previously thought.

Hypernegative supercoiling inhibits growth by causing RNA degradation.

Transcription-induced hypernegative supercoiling is a hallmark of Escherichia coli topoisomerase I (topA) mutants. However, its physiological significance has remained unclear. Temperature downshift of a mutant yielded transient growth arrest and a parallel increase in hypernegative supercoiling that was more severe with lower temperature. Both properties were alleviated by overexpression of RNase HI. While ribosomes in extracts showed normal activity when obtained during growth arrest, mRNA on ribosomes was reduced for fis and shorter for crp, polysomes were much less abundant relative to monosomes, and protein synthesis rate dropped, as did the ratio of large to small proteins. Altered processing and degradation of lacA and fis mRNA was also observed. These data are consistent with truncation of mRNA during growth arrest. These effects were not affected by a mutation in the gene encoding RNase E, indicating that this endonuclease is not involved in the abnormal mRNA processing. They were also unaffected by spectinomycin, an inhibitor of protein synthesis, which argued against induction of RNase activity. In vitro transcription revealed that R-loop formation is more extensive on hypernegatively supercoiled templates. These results allow us, for the first time, to present a model by which hypernegative supercoiling inhibits growth. In this model, the introduction of hypernegative supercoiling by gyrase facilitates degradation of nascent RNA; overproduction of RNase HI limits the accumulation of hypernegative supercoiling, thereby preventing extensive RNA degradation.

Colistin Induces S. aureus Susceptibility to Bacitracin.

Bacitracin has been used in topical preparations with polymyxin B for bacterial infections. Colistin belongs to the polymyxin group of antibiotics and is effective against most Gram-negative bacilli. This study investigated whether colistin could affect the susceptibility of S. aureus to bacitracin. S. aureus isolates were first incubated with colistin and the susceptibility of S. aureus to bacitracin was increased. The effect of the combination of colistin and bacitracin on S. aureus was then confirmed by the checkerboard assay and the time-kill kinetics. The Triton X-100-induced autolysis was significantly increased after S. aureus was exposed to colistin. Exposure to colistin also led to a less positive charge on the cell surface and a significant leakage of Na+, Mg2, K+, Ca2+, Mn2+, Cu2+, and Zn2+. Finally, disruptions on the cell surface and an irregular morphology were observed when the bacteria were exposed to colistin and bacitracin. Bacitracin had a stronger antibacterial activity against S. aureus in the presence of colistin. This could be due to the fact that colistin damaged the bacterial membrane. This study suggests that combination of colistin with bacitracin has a potential for treating clinical S. aureus infections.

Impact of the choice of reference genome on the ability of the core genome SNV methodology to distinguish strains of Salmonella enterica serovar Heidelberg.

Salmonella enterica serovar Heidelberg (S. Heidelberg) is one of the top serovars causing human salmonellosis. The core genome single nucleotide variant pipeline (cgSNV) is one of several whole genome based sequence typing methods used for the laboratory investigation of foodborne pathogens. SNV detection using this method requires a reference genome. The purpose of this study was to investigate the impact of the choice of the reference genome on the cgSNV-informed phylogenetic clustering and inferred isolate relationships. We found that using a draft or closed genome of S. Heidelberg as reference did not impact the ability of the cgSNV methodology to differentiate among 145 S. Heidelberg isolates involved in foodborne outbreaks. We also found that using a distantly related genome such as S. Dublin as choice of reference led to a loss in resolution since some sporadic isolates were found to cluster together with outbreak isolates. In addition, the genetic distances between outbreak isolates as well as between outbreak and sporadic isolates were overall reduced when S. Dublin was used as the reference genome as opposed to S. Heidelberg.

Salmonella enterica Prophage Sequence Profiles Reflect Genome Diversity and Can Be Used for High Discrimination Subtyping.

Non-typhoidal Salmonella is a leading cause of foodborne illness worldwide. Prompt and accurate identification of the sources of Salmonella responsible for disease outbreaks is crucial to minimize infections and eliminate ongoing sources of contamination. Current subtyping tools including single nucleotide polymorphism (SNP) typing may be inadequate, in some instances, to provide the required discrimination among epidemiologically unrelated Salmonella strains. Prophage genes represent the majority of the accessory genes in bacteria genomes and have potential to be used as high discrimination markers in Salmonella. In this study, the prophage sequence diversity in different Salmonella serovars and genetically related strains was investigated. Using whole genome sequences of 1,760 isolates of S. enterica representing 151 Salmonella serovars and 66 closely related bacteria, prophage sequences were identified from assembled contigs using PHASTER. We detected 154 different prophages in S. enterica genomes. Prophage sequences were highly variable among S. enterica serovars with a median ± interquartile range (IQR) of 5 ± 3 prophage regions per genome. While some prophage sequences were highly conserved among the strains of specific serovars, few regions were lineage specific. Therefore, strains belonging to each serovar could be clustered separately based on their prophage content. Analysis of S. Enteritidis isolates from seven outbreaks generated distinct prophage profiles for each outbreak. Taken altogether, the diversity of the prophage sequences correlates with genome diversity. Prophage repertoires provide an additional marker for differentiating S. enterica subtypes during foodborne outbreaks.

Mutations reducing replication from R-loops suppress the defects of growth, chromosome segregation and DNA supercoiling in cells lacking topoisomerase I and RNase HI activity.

R-loop formation occurs when the nascent RNA hybridizes with the template DNA strand behind the RNA polymerase. R-loops affect a wide range of cellular processes and their use as origins of replication was the first function attributed to them. In Escherichia coli, R-loop formation is promoted by the ATP-dependent negative supercoiling activity of gyrase (gyrA and gyrB) and is inhibited by topoisomerase (topo) I (topA) relaxing transcription-induced negative supercoiling. RNase HI (rnhA) degrades the RNA moiety of R-loops. The depletion of RNase HI activity in topA null mutants was previously shown to lead to extensive DNA relaxation, due to DNA gyrase inhibition, and to severe growth and chromosome segregation defects that were partially corrected by overproducing topo III (topB). Here, DNA gyrase assays in crude cell extracts showed that the ATP-dependent activity (supercoiling) of gyrase but not its ATP-independent activity (relaxation) was inhibited in topA null cells lacking RNase HI. To characterize the cellular event(s) triggered by the absence of RNase HI, we performed a genetic screen for suppressors of the growth defect of topA rnhA null cells. Suppressors affecting genes in replication (holC2::aph and dnaT18::aph) nucleotide metabolism (dcd49::aph), RNA degradation (rne59::aph) and fimbriae synthesis (fimD22::aph) were found to reduce replication from R-loops and to restore supercoiling, thus pointing to a correlation between R-loop-dependent replication in topA rnhA mutants and the inhibition of gyrase activity and growth. Interestingly, the position of fimD on the E. coli chromosome corresponds to the site of one of the five main putative origins of replication from R-loops in rnhA null cells recently identified by next-generation sequencing, thus suggesting that the fimD22::aph mutation inactivated one of these origins. Furthermore, we show that topo III overproduction is unable to complement the growth defect of topA rnhA null mutants at low temperatures that stabilizes hyper-negatively supercoiled DNA.