Search

Abstract As observed in cancers, individual mutagens and defects in DNA repair create distinctive mutational signatures that combine to form context-specific spectra within cells. We reasoned that similar processes must occur in bacterial lineages, potentially allowing decomposition analysis to detect both disruption of DNA repair processes and exposure to niche-specific mutagens. Here we reconstruct mutational spectra for 84 clades from 31 diverse bacterial species and find distinct mutational patterns. We extract signatures driven by specific DNA repair defects using hypermutator lineages, and further deconvolute the spectra into multiple signatures operating within different clades. We show that these signatures are explained by both bacterial phylogeny and replication niche. By comparing mutational spectra of clades from different environmental and biological locations, we identify niche-associated mutational signatures, and then employ these signatures to infer the predominant replication niches for several clades where this was previously obscure. Our results show that mutational spectra may be associated with sites of bacterial replication when mutagen exposures differ, and can be used in these cases to infer transmission routes for established and emergent human bacterial pathogens.

Salmonella enterica (S. enterica) is a causative agent of multiple outbreaks of foodborne illness associated with fresh produce, including pre-cut melon and leafy vegetables. Current industrial antimicrobial interventions have been shown to reduce microbial populations by 3 log CFU/g and S. Newport S2 by 1 log CFU/g on both lettuce and cantaloupe tissues at all sampling times. In contrast, populations of strains S. Thompson S193 and S194 were reduced by 2 log CFU/g on day 0 on lettuce, but were not significantly different (P > 0.05) from the controls thereafter, S. Newport S195 populations were reduced on lettuce by 1 log CFU/g on day 0 and no reductions were found on cantaloupe tissue. Both aerobic populations and water activity were higher on cantaloupe than on lettuce. The water activity of lettuce decreased significantly (P

The rapid detection of foodborne microbial pathogens contaminating fresh fruits and vegetables during the intervening period between harvest and consumption could revolutionize microbial quality assurance of food usually consumed raw and those with a limited shelf life. We have developed a sensitive, shotgun whole genome sequencing protocol capable of detecting as few as 1 colony forming unit (cfu) of Salmonella enterica serovar Typhimurium spiked on 25 g of lettuce. The Ion Torrent sequencing platform was used to generate reads of globally amplified DNA from microbes recovered from the surface of lettuce followed by bioinformatic analyses of the nucleotide sequences to detect the presence of Salmonella. The test is rapid and sensitive, and appropriate for testing perishable foods, and those consumed raw, for Salmonella contamination. The test has the potential to be universally applicable to any microbial contaminant on lettuce as long as a suitable bioinformatics pipeline is available and validated. A universal test is expected to pave the way for preventive and precision food safety and the re-shaping of the entire spectrum of food safety investigations from the current disease-limiting, reactive procedure to a proactive, disease prevention process.

The goal of this study was to isolate, screen, and characterize Arctic microbial isolates from Expedition Fjord, Axel Heiberg Island, Nunavut, Canada capable of inhibiting the growth of foodborne and clinically relevant pathogens. Arctic bacteria were isolated from twelve different high Arctic habitats pertaining to active layer permafrost soil, saline spring sediments, lake sediments, and endoliths. This was achieved using (1) the cryo-iPlate, an innovative in situ cultivation device within active layer permafrost soil and (2) bulk plating of Arctic samples by undergraduate students that applied standard culturing methods. To mitigate the possibility of identifying isolates with already-known antibacterial activities, a cell-based dereplication platform was used. Ten out of the twelve Arctic habitats tested were found to yield cold-adapted isolates with antibacterial activity. Eight cold-adapted Arctic isolates were identified with the ability to inhibit the entire dereplication platform, suggesting the possibility of new mechanisms of action. Two promising isolates, initially cultured from perennial saline spring sediments and from active layer permafrost soil (Paenibacillus sp. GHS.8.NWYW.5 and Pseudomonas sp. AALPS.10.MNAAK.13, respectively), displayed antibacterial activity against foodborne and clinically relevant pathogens. Paenibacillus sp. GHS.8.NWYW.5 was capable of inhibiting methicillin resistant and susceptible Staphylococcus aureus (MRSA and MSSA), Listeria monocytogenes, Salmonella enterica and Escherichia coli O157:H7. Pseudomonas sp. AALPS.10.MNAAK.13 was observed to have antagonistic activity against MRSA, MSSA, Acinetobacter baumanii, Enterococcus faecium, and Enterococcus faecalis. After whole genome sequencing and mining, the genome of Paenibacillus sp. GHS.8.NWYW.5 was found to contain seven putative secondary metabolite biosynthetic gene clusters that displayed low homology ( 0) to clusters identified within the genome of the type strain pertaining to the same species. These findings suggest that cold-adapted Arctic microbes may be a promising source of novel secondary metabolites for potential use in both industrial and medical settings.

ABSTRACT The O-specific antigen (OSA) in Pseudomonas aeruginosa lipopolysaccharide is highly varied by sugar identity, side chains, and bond between O-repeats. These differences classified P. aeruginosa into 20 distinct serotypes. In the past few decades, O12 has emerged as the predominant serotype in clinical settings and outbreaks. These serotype O12 isolates exhibit high levels of resistance to various classes of antibiotics. Here, we explore how the P. aeruginosa OSA biosynthesis gene clusters evolve in the population by investigating the association between the phylogenetic relationships among 83 P. aeruginosa strains and their serotypes. While most serotypes were closely linked to the core genome phylogeny, we observed horizontal exchange of OSA biosynthesis genes among phylogenetically distinct P. aeruginosa strains. Specifically, we identified a “serotype island” ranging from 62 kb to 185 kb containing the P. aeruginosa O12 OSA gene cluster, an antibiotic resistance determinant (gyrAC248T), and other genes that have been transferred between P. aeruginosa strains with distinct core genome architectures. We showed that these genes were likely acquired from an O12 serotype strain that is closely related to P. aeruginosa PA7. Acquisition and recombination of the “serotype island” resulted in displacement of the native OSA gene cluster and expression of the O12 serotype in the recipients. Serotype switching by recombination has apparently occurred multiple times involving bacteria of various genomic backgrounds. In conclusion, serotype switching in combination with acquisition of an antibiotic resistance determinant most likely contributed to the dissemination of the O12 serotype in clinical settings. IMPORTANCE Infection rates in hospital settings by multidrug-resistant (MDR) Pseudomonas aeruginosa clones have increased during the past decades, and serotype O12 is predominant among these epidemic strains. It is not known why the MDR phenotype is associated with serotype O12 and how this clone type has emerged. This study shows that evolution of MDR O12 strains involved a switch from an ancestral O4 serotype to O12. Serotype switching was the result of horizontal transfer and genetic recombination of lipopolysaccharide (LPS) biosynthesis genes originating from an MDR taxonomic outlier P. aeruginosa strain. Moreover, the recombination event also resulted in acquisition of antibiotic resistance genes. These results impact on our understanding of MDR outbreak strain and serotype evolution and can potentially assist in better monitoring and prevention.

Abstract Papaver genus, commonly known as popies, is a valuable source of alkaloids used in medicine, including papaverine, morphine, codeine, and thebaine. We isolated six endophytic fungal isolates producing morphinan alkaloids from four Papaver species growing in Kurdistan Province, Iran. To do this, a 1:1 mixture of methanol and chloroform was used to extract fungal cultures. The contents of morphinan alkaloids in the extracts were subsequently determined using phase high-performance liquid chromatography. Among the morphinan alkaloid-producing fungal isolates, IRAN 4653C had the highest yield giving 23.06 (mg/g) morphine and 2.03 (mg/g) codeine when grown in potato dextrose liquid medium. The identity of this isolate was examined and recognized as a new fungal species named as Pithoascus kurdistanesis sp. nov. based on multi-gene phylogenetic analyses of ITS, TEF-1α, and TUB2 sequences data and morphological features. The morphinan-producing endophytic fungus and the isolated Pithoascus species from Papaver are being reported for the first time. Accordingly, this fungus shows promise as a new source of valuable compounds which is illustrated and introduced here as a new Microascaceae member belonging to Pithoascus from Kurdistan Province, Iran. Moreover, the morphinan productivity of P. kurdistanesis was further validated by gas chromatography-mass spectrometry (GC–MS).

Abstract Background Flight can drastically enhance dispersal capacity and is a key trait defining the potential of exotic insect species to spread and invade new habitats. The phytophagous European spongy moths (ESM, Lymantria dispar dispar) and Asian spongy moths (ASM; a multi–species group represented here by L. d. asiatica and L. d. japonica), are globally invasive species that vary in adult female flight capability—female ASM are typically flight capable, whereas female ESM are typically flightless. Genetic markers of flight capability would supply a powerful tool for flight profiling of these species at any intercepted life stage. To assess the functional complexity of spongy moth flight and to identify potential markers of flight capability, we used multiple genetic approaches aimed at capturing complementary signals of putative flight–relevant genetic divergence between ESM and ASM: reduced representation genome–wide association studies, whole genome sequence comparisons, and developmental transcriptomics. We then judged the candidacy of flight–associated genes through functional analyses aimed at addressing the proximate demands of flight and salient features of the ecological context of spongy moth flight evolution. Results Candidate gene sets were typically non–overlapping across different genetic approaches, with only nine gene annotations shared between any pair of approaches. We detected an array of flight–relevant functional themes across gene sets that collectively suggest divergence in flight capability between European and Asian spongy moth lineages has coincided with evolutionary differentiation in multiple aspects of flight development, execution, and surrounding life history. Overall, our results indicate that spongy moth flight evolution has shaped or been influenced by a large and functionally broad network of traits. Conclusions Our study identified a suite of flight–associated genes in spongy moths suited to exploration of the genetic architecture and evolution of flight, or validation for flight profiling purposes. This work illustrates how complementary genetic approaches combined with phenotypically targeted functional analyses can help to characterize genetically complex traits.

IntroductionThere is little information on evolutionarily ancient eukaryotes, which are often referred to as basal eukaryotes, in Arctic waters. Despite earlier studies being conducted in the Russian White Sea, only few have been reported.MethodsFollowing a shotgun sequence survey of diatom cultures from Sugluk Inlet off the Hudson Strait in Northern Québec, we obtained the complete mitochondrial genome and the operon of nuclear ribosomal RNA genes from a strain that matches that of Ancyromonas sigmoides (Kent, 1881).ResultsThe sequence of the mitogenome retrieved was 41,889 bp in length and encoded 38 protein-coding genes, 5 non-conserved open-reading frames, and 2 rRNA and 24 tRNA genes. The mitogenome has retained sdh2 and sdh3, two genes of the succinate dehydrogenase complex, which are sometimes found among basal eukaryotes but seemingly missing among the Malawimonadidae, a lineage sister to Ancyromonadida in some phylogenies. The phylogeny inferred from the 18S rRNA gene associated A. sigmoides from Sugluk Inlet with several other strains originating from the Arctic. The study also unveiled the presence of a metagenomic sequence ascribed to bacteria in GenBank, but it was clearly a mitochondrial genome with a gene content highly similar to that of A. sigmoides, including the non-conserved open-reading frames.DiscussionAfter re-annotation, a phylogeny was inferred from mitochondrial protein sequences, and it strongly associated A. sigmoides with the misidentified organism, with the two being possibly conspecific or sibling species as they are more similar to one another than to species of the genus Malawimonas. Overall our phylogeny showed that the ice associated ancryomonads were clearly distinct from more southerly strains.

With the increasing occurrence and severity of cyanobacterial harmful algal blooms (cHAB) at the global scale, there is an urgent need for rapid, accurate, accessible, and cost-effective detection tools. Here, we detail the RosHAB workflow, an innovative, in-the-field applicable genomics approach for real-time, early detection of cHAB outbreaks. We present how the proposed workflow offers consistent taxonomic identification of water samples in comparison to traditional microscopic analyses in a few hours and discuss how the generated data can be used to deepen our understanding on cyanobacteria ecology and forecast HABs events. In parallel, processed water samples will be used to iteratively build the International cyanobacterial toxin database (ICYATOX; http://icyatox.ibis.ulaval.ca) containing the analysis of novel cyanobacterial genomes, including phenomics and genomics metadata. Ultimately, RosHAB will (1) improve the accuracy of on-site rapid diagnostics, (2) standardize genomic procedures in the field, (3) facilitate these genomics procedures for non-scientific personnel, and (4) identify prognostic markers for evidence-based decisions in HABs surveillance.

Bacteriophages (phages) are potential alternatives to chemical antimicrobials against pathogens of public health significance. Understanding the diversity and host specificity of phages is important for developing effective phage biocontrol approaches. Here, we assessed the host range, morphology, and genetic diversity of eight Salmonella enterica phages isolated from a wastewater treatment plant. The host range analysis revealed that six out of eight phages lysed more than 81% of the 43 Salmonella enterica isolates tested. The genomic sequences of all phages were determined. Whole-genome sequencing (WGS) data revealed that phage genome sizes ranged from 41 to 114 kb, with GC contents between 39.9 and 50.0%. Two of the phages SB13 and SB28 represent new species, Epseptimavirus SB13 and genera Macdonaldcampvirus, respectively, as designated by the International Committee for the Taxonomy of Viruses (ICTV) using genome-based taxonomic classification. One phage (SB18) belonged to the Myoviridae morphotype while the remaining phages belonged to the Siphoviridae morphotype. The gene content analyses showed that none of the phages possessed virulence, toxin, antibiotic resistance, type I–VI toxin–antitoxin modules, or lysogeny genes. Three (SB3, SB15, and SB18) out of the eight phages possessed tailspike proteins. Whole-genome-based phylogeny of the eight phages with their 113 homologs revealed three clusters A, B, and C and seven subclusters (A1, A2, A3, B1, B2, C1, and C2). While cluster C1 phages were predominantly isolated from animal sources, cluster B contained phages from both wastewater and animal sources. The broad host range of these phages highlights their potential use for controlling the presence of S. enterica in foods.

ABSTRACT Epidemic strains of Pseudomonas aeruginosa often have increased levels of antibiotic resistance, and these resistance phenotypes are difficult to predict. We used quantitative proteomics to compare the responses of laboratory strain PAO1 and two isolates of the Liverpool epidemic strain (LES) to four antibiotics and hydrogen peroxide. The majority of proteome changes were unique to either one isolate or treatment, but smaller groups of proteins were differentially abundant in more than one sample. Proteins in these shared adaptive responses represent promising avenues for further investigation to understand the contribution of these proteins to resistance and their potential as targets for novel treatments. We observed the largest overlap between the proteome profiles of PAO1 challenged with tobramycin and hydrogen peroxide, and our data support previous work suggesting a role for heat shock protein IbpA in response to tobramycin. Our proteome profiling uncovered extensive changes in LESB58 in response to carbenicillin, with more than 1,000 proteins significantly changed in abundance. This included unique changes in the abundance of proteins involved in cell wall synthesis and division. We used phase contrast microscopy to check for corresponding changes in cell morphology and show that LESB58 maintained shorter cell lengths under treatment with β-lactams. We propose that the ability to maintain the processes of cell wall synthesis and division through changes in protein abundances contributes to the high levels of β-lactam resistance in LESB58. IMPORTANCE Pseudomonas aeruginosa is an important pathogen often associated with hospital-acquired infections and chronic lung infections in people with cystic fibrosis. P. aeruginosa possesses a wide array of intrinsic and adaptive mechanisms of antibiotic resistance, and the regulation of these mechanisms is complex. Label-free quantitative proteomics is a powerful tool to compare susceptible and resistant strains of bacteria and their responses to antibiotic treatments. Here we compare the proteomes of three isolates of P. aeruginosa with different antibiotic resistance profiles in response to five challenge conditions. We uncover unique and shared proteome changes for the widely used laboratory strain PAO1 and two isolates of the Liverpool epidemic strain of P. aeruginosa, LESlike1 and LESB58. Our data set provides insight into antibiotic resistance in clinically relevant Pseudomonas isolates and highlights proteins, including those with uncharacterized functions, which can be further investigated for their role in adaptive responses to antibiotic treatments.

IntroductionPersistent symptoms after COVID-19 infection (“long COVID”) negatively affects almost half of COVID-19 survivors. Despite its prevalence, its pathophysiology is poorly understood, with multiple host systems likely affected. Here, we followed patients from hospital to discharge and used a systems-biology approach to identify mechanisms of long COVID.MethodsRNA-seq was performed on whole blood collected early in hospital and 4-12 weeks after discharge from 24 adult COVID-19 patients (10 reported post-COVID symptoms after discharge). Differential gene expression analysis, pathway enrichment, and machine learning methods were used to identify underlying mechanisms for post-COVID symptom development.ResultsCompared to patients with post-COVID symptoms, patients without post-COVID symptoms had larger temporal gene expression changes associated with downregulation of inflammatory and coagulation genes over time. Patients could also be separated into three patient endotypes with differing mechanistic trajectories, which was validated in another published patient cohort. The “Resolved” endotype (lowest rate of post-COVID symptoms) had robust inflammatory and hemostatic responses in hospital that resolved after discharge. Conversely, the inflammatory/hemostatic responses of “Suppressive” and “Unresolved” endotypes (higher rates of patients with post-COVID symptoms) were persistently dampened and activated, respectively. These endotypes were accurately defined by specific blood gene expression signatures (6-7 genes) for potential clinical stratification.DiscussionThis study allowed analysis of long COVID whole blood transcriptomics trajectories while accounting for the issue of patient heterogeneity. Two of the three identified and externally validated endotypes (“Unresolved” and “Suppressive”) were associated with higher rates of post-COVID symptoms and either persistently activated or suppressed inflammation and coagulation processes. Gene biomarkers in blood could potentially be used clinically to stratify patients into different endotypes, paving the way for personalized long COVID treatment.

Abstract Severely-afflicted COVID-19 patients can exhibit disease manifestations representative of sepsis, including acute respiratory distress syndrome and multiple organ failure. We hypothesized that diagnostic tools used in managing all-cause sepsis, such as clinical criteria, biomarkers, and gene expression signatures, should extend to COVID-19 patients. Here we analyzed the whole blood transcriptome of 124 early (1–5 days post-hospital admission) and late (6–20 days post-admission) sampled patients with confirmed COVID-19 infections from hospitals in Quebec, Canada. Mechanisms associated with COVID-19 severity were identified between severity groups (ranging from mild disease to the requirement for mechanical ventilation and mortality), and established sepsis signatures were assessed for dysregulation. Specifically, gene expression signatures representing pathophysiological events, namely cellular reprogramming, organ dysfunction, and mortality, were significantly enriched and predictive of severity and lethality in COVID-19 patients. Mechanistic endotypes reflective of distinct sepsis aetiologies and therapeutic opportunities were also identified in subsets of patients, enabling prediction of potentially-effective repurposed drugs. The expression of sepsis gene expression signatures in severely-afflicted COVID-19 patients indicates that these patients should be classified as having severe sepsis. Accordingly, in severe COVID-19 patients, these signatures should be strongly considered for the mechanistic characterization, diagnosis, and guidance of treatment using repurposed drugs.

Abstract Background Structural variants (SVs), including deletions, insertions, duplications, and inversions, are relatively long genomic variations implicated in a diverse range of processes from human disease to ecology and evolution. Given their complex signatures, tendency to occur in repeated regions, and large size, discovering SVs based on short reads is challenging compared to single-nucleotide variants. The increasing availability of long-read technologies has greatly facilitated SV discovery; however, these technologies remain too costly to apply routinely to population-level studies. Here, we combined short-read and long-read sequencing technologies to provide a comprehensive population-scale assessment of structural variation in a panel of Canadian soybean cultivars. Results We used Oxford Nanopore long-read sequencing data (~12× mean coverage) for 17 samples to both benchmark SV calls made from Illumina short-read data and predict SVs that were subsequently genotyped in a population of 102 samples using Illumina data. Benchmarking results show that variants discovered using Oxford Nanopore can be accurately genotyped from the Illumina data. We first use the genotyped deletions and insertions for population genetics analyses and show that results are comparable to those based on single-nucleotide variants. We observe that the population frequency and distribution within the genome of deletions and insertions are constrained by the location of genes. Gene Ontology and PFAM domain enrichment analyses also confirm previous reports that genes harboring high-frequency deletions and insertions are enriched for functions in defense response. Finally, we discover polymorphic transposable elements from the deletions and insertions and report evidence of the recent activity of a Stowaway MITE. Conclusions We show that structural variants discovered using Oxford Nanopore data can be genotyped with high accuracy from Illumina data. Our results demonstrate that long-read and short-read sequencing technologies can be efficiently combined to enhance SV analysis in large populations, providing a reusable framework for their study in a wider range of samples and non-model species.

Driven in part by its metabolic versatility, high intrinsic antibiotic resistance, and a large repertoire of virulence factors, Pseudomonas aeruginosa is expertly adapted to thrive in a wide variety of environments, and in the process, making it a notorious opportunistic pathogen. Apart from the extensively studied chronic infection in the lungs of people with cystic fibrosis (CF), P. aeruginosa also causes multiple serious infections encompassing essentially all organs of the human body, among others, lung infection in patients with chronic obstructive pulmonary disease, primary ciliary dyskinesia and ventilator-associated pneumonia; bacteremia and sepsis; soft tissue infection in burns, open wounds and postsurgery patients; urinary tract infection; diabetic foot ulcers; chronic suppurative otitis media and otitis externa; and keratitis associated with extended contact lens use. Although well characterized in the context of CF, pathogenic processes mediated by various P. aeruginosa virulence factors in other organ systems remain poorly understood. In this review, we use an organ system-based approach to provide a synopsis of disease mechanisms exerted by P. aeruginosa virulence determinants that contribute to its success as a versatile pathogen.

ABSTRACT Eight isolates of the Liverpool epidemic strain (LES) of Pseudomonas aeruginosa have previously been characterized using comparative genomics and preliminary phenotypic assays. Here, we extend the characterization of these clinically relevant P. aeruginosa isolates with planktonic and biofilm growth assays and analysis of antibiotic susceptibility for both planktonic and biofilm cultures. Laboratory strains PAO1 and PA14 were included as comparator strains. Antibiotic susceptibility to eight classes of antibiotics was determined. MICs were determined to measure susceptibility of planktonic cultures, and minimum biofilm eradication concentration (MBEC) assays were used to estimate levels of resistance during the production of biofilm. LES isolates had high levels of resistance compared with laboratory reference strains when grown planktonically (up to nine 2-fold dilutions higher), and resistance was increased in the biofilm mode of growth. Measurements of biofilm biomass in the MBEC assays showed that certain isolates often show increased biofilm biomass in the presence of antibiotics. IMPORTANCE Pseudomonas aeruginosa is an opportunistic pathogen with high intrinsic antibiotic resistance. This resistance is typically increased in clinical isolates through adaptations to the host and production of small-colony variants (SCVs) and when P. aeruginosa forms biofilms, which are surface-attached communities that are protected by a self-produced matrix. Understanding the combination of SCVs, biofilm production, and the diversity of drug resistance phenotypes in clinical isolates can lead to improved treatments for P. aeruginosa infections.

Abstract Background Salmonella enterica is a leading cause of foodborne illness worldwide resulting in considerable public health and economic costs. Testing for the presence of this pathogen in food is often hampered by the presence of background microflora that may present as Salmonella (false positives). False positive isolates belonging to the genus Citrobacter can be difficult to distinguish from Salmonella due to similarities in their genetics, cell surface antigens, and other phenotypes. In order to understand the genetic basis of these similarities, a comparative genomic approach was used to define the pan-, core, accessory, and unique coding sequences of a representative population of Salmonella and Citrobacter strains. Results Analysis of the genomic content of 58 S. enterica strains and 37 Citrobacter strains revealed the presence of 31,130 and 1540 coding sequences within the pan- and core genome of this population. Amino acid sequences unique to either Salmonella (n = 1112) or Citrobacter (n = 195) were identified and revealed potential niche-specific adaptations. Phylogenetic network analysis of the protein families encoded by the pan-genome indicated that genetic exchange between Salmonella and Citrobacter may have led to the acquisition of similar traits and also diversification within the genera. Conclusions Core genome analysis suggests that the Salmonella enterica and Citrobacter populations investigated here share a common evolutionary history. Comparative analysis of the core and pan-genomes was able to define the genetic features that distinguish Salmonella from Citrobacter and highlight niche specific adaptations.

The emergence of multidrug resistance (MDR) in bacteria represents a global threat to human health. Here, Cazares et al. identify a family of MDR megaplasmids carrying large arrays of antibiotic resistance genes in Pseudomonas strains from various sources, including P. aeruginosa clinical isolates.

Abstract The spruce budworm, Choristoneura fumiferana, is presumed to be panmictic across vast regions of North America. We examined the extent of panmixia by genotyping 3,650 single nucleotide polymorphism (SNP) loci in 1975 individuals from 128 collections across the continent. We found three spatially structured subpopulations: Western (Alaska, Yukon), Central (southeastern Yukon to the Manitoba–Ontario border), and Eastern (Manitoba–Ontario border to the Atlantic). Additionally, the most diagnostic genetic differentiation between the Central and Eastern subpopulations was chromosomally restricted to a single block of SNPs that may constitute an island of differentiation within the species. Geographic differentiation in the spruce budworm parallels that of its principal larval host, white spruce (Picea glauca), providing evidence that spruce budworm and spruce trees survived in the Beringian refugium through the Last Glacial Maximum and that at least two isolated spruce budworm populations diverged with spruce/fir south of the ice sheets. Gene flow in the spruce budworm may also be affected by mountains in western North America, habitat isolation in West Virginia, regional adaptations, factors related to dispersal, and proximity of other species in the spruce budworm species complex. The central and eastern geographic regions contain individuals that assign to Eastern and Central subpopulations, respectively, indicating that these barriers are not complete. Our discovery of previously undetected geographic and genomic structure in the spruce budworm suggests that further population modelling of this ecologically important insect should consider regional differentiation, potentially co‐adapted blocks of genes, and gene flow between subpopulations.

Abstract Spatial synchrony is a common characteristic of spatio‐temporal population dynamics across many taxa. While it is known that both dispersal and spatially autocorrelated environmental variation (i.e., the Moran effect) can synchronize populations, the relative contributions of each, and how they interact, are generally unknown. Distinguishing these mechanisms and their effects on synchrony can help us to better understand spatial population dynamics, design conservation and management strategies, and predict climate change impacts. Population genetic data can be used to tease apart these two processes as the spatio‐temporal genetic patterns they create are expected to be different. A challenge, however, is that genetic data are often collected at a single point in time, which may introduce context‐specific bias. Spatio‐temporal sampling strategies can be used to reduce bias and to improve our characterization of the drivers of spatial synchrony. Using spatio‐temporal analyses of genotypic data, our objective was to identify the relative support for these two mechanisms to the spatial synchrony in population dynamics of the irruptive forest insect pest, the spruce budworm (Choristoneura fumiferana), in Quebec (Canada). AMOVA, cluster analysis, isolation by distance, and sPCA were used to characterize spatio‐temporal genomic variation using 1,370 SBW larvae sampled over four years (2012–2015) and genotyped at 3,562 SNP loci. We found evidence of overall weak spatial genetic structure that decreased from 2012 to 2015 and a genetic diversity homogenization among the sites. We also found genetic evidence of a long‐distance dispersal event over >140 km. These results indicate that dispersal is the key mechanism involved in driving population synchrony of the outbreak. Early intervention management strategies that aim to control source populations have the potential to be effective through limiting dispersal. However, the timing of such interventions relative to outbreak progression is likely to influence their probability of success.

Aeromonas salmonicida subsp. salmonicida is a Gram-negative bacterium that causes furunculosis, a fish disease claiming substantial economic losses in the aquaculture industry. Major challenges exist in monitoring and controlling fish infections in aquaculture farms. Development of management practices to improve the sustainability of fish farming with disease prevention necessitates studies using well-defined systems and well-characterized bacterial isolates. Even though several A. salmonicida subsp. salmonicida genomes have been completely assembled and thoroughly annotated, in vivo pathogenicity data are lacking. Here we present A. salmonicida subsp. salmonicida 890054 as a prototype strain for standardized furunculosis challenges with survival data. Computational analysis of sequencing results provided a complete circular genome with annotations of plasmids carrying virulence factors, antimicrobial resistance, and secondary metabolite coding genes. The analysis also revealed the presence of an IncU plasmid distinct from other IncU plasmids previously associated with Aeromonas.

We have developed a targeted, amplicon-based next-generation sequencing method to detect and analyze 227 virulence genes (VG) of Salmonella (AmpliSeqSalm_227VG) for assessing the pathogenicity potential of Salmonella. The procedure was developed using 80 reference genomes representing 75 epidemiologically-relevant serovars associated with human salmonellosis. We applied the AmpliSeqSalm_227VG assay to (a) 35 previously characterized field strains of Salmonella consisting of serovars commonly incriminated in foodborne illnesses and (b) 34 Salmonella strains with undisclosed serological or virulence attributes, and were able to divide Salmonella VGs into two groups: core VGs and variable VGs. The commonest serovars causing foodborne illnesses such as Enteritidis, Typhimurium, Heidelberg and Newport had a high number of VGs (217–227). In contrast, serovars of subspecies not commonly associated with human illnesses, such as houtenae, arizonae and salame, tended to have fewer VGs (177–195). Variable VGs were not only infrequent but, when present, displayed considerable sequence variation: safC, sseL, sseD, sseE, ssaK and stdB showed the highest variation and were linked to strain pathogenicity. In a chicken infection model, VGs belonging to rfb and sse operons showed differences and were linked with pathogenicity. The high-throughput, targeted NGS-based AmpliSeqSalm_227VG procedure provided previously unknown information about variation in select virulence genes that can now be applied to a much larger population of Salmonella for evaluating pathogenicity of various serovars of Salmonella and for risk assessment of foodborne salmonellosis.

ABSTRACT Salmonella comprises more than 2,600 serovars. Very few environmental and uncommon serovars have been characterized for their potential role in virulence and human infections. A complementary in vitro and in vivo systematic high-throughput analysis of virulence was used to elucidate the association between genetic and phenotypic variations across Salmonella isolates. The goal was to develop a strategy for the classification of isolates as a benchmark and predict virulence levels of isolates. Thirty-five phylogenetically distant strains of unknown virulence were selected from the Salmonella Foodborne Syst-OMICS (SalFoS) collection, representing 34 different serovars isolated from various sources. Isolates were evaluated for virulence in 4 complementary models of infection to compare virulence traits with the genomics data, including interactions with human intestinal epithelial cells, human macrophages, and amoeba. In vivo testing was conducted using the mouse model of Salmonella systemic infection. Significant correlations were identified between the different models. We identified a collection of novel hypothetical and conserved proteins associated with isolates that generate a high burden. We also showed that blind prediction of virulence of 33 additional strains based on the pan-genome was high in the mouse model of systemic infection (82% agreement) and in the human epithelial cell model (74% agreement). These complementary approaches enabled us to define virulence potential in different isolates and present a novel strategy for risk assessment of specific strains and for better monitoring and source tracking during outbreaks. IMPORTANCE Salmonella species are bacteria that are a major source of foodborne disease through contamination of a diversity of foods, including meat, eggs, fruits, nuts, and vegetables. More than 2,600 different Salmonella enterica serovars have been identified, and only a few of them are associated with illness in humans. Despite the fact that they are genetically closely related, there is enormous variation in the virulence of different isolates of Salmonella enterica. Identification of foodborne pathogens is a lengthy process based on microbiological, biochemical, and immunological methods. Here, we worked toward new ways of integrating whole-genome sequencing (WGS) approaches into food safety practices. We used WGS to build associations between virulence and genetic diversity within 83 Salmonella isolates representing 77 different Salmonella serovars. Our work demonstrates the potential of combining a genomics approach and virulence tests to improve the diagnostics and assess risk of human illness associated with specific Salmonella isolates.

Current treatments for Pseudomonas aeruginosa infections are becoming less effective because of the increasing rates of multi-antibiotic resistance. Pharmacological targeting of virulence through inhibition of quorum sensing (QS) dependent virulence gene regulation has considerable therapeutic potential. In P. aeruginosa, the pqs QS system regulates the production of multiple virulence factors as well as biofilm maturation and is a promising approach for developing antimicrobial adjuvants for combatting drug resistance. In this work, we report the hit optimisation for a series of potent novel inhibitors of PqsR, a key regulator of the pqs system, bearing a 2-((5-methyl-5H-[1,2,4]triazino[5,6-b]indol-3-yl)thio) acetamide scaffold. The initial hit compound 7 (PAO1-L IC50 0.98 ± 0.02 μM, PA14 inactive at 10 μM) was obtained through a virtual screening campaign performed on the PqsR ligand binding domain using the University of Nottingham Managed Chemical Compound Collection. Hit optimisation gave compounds with enhanced potency against strains PAO1-L and PA14, evaluated using P. aeruginosa pqs-based QS bioreporter assays. Compound 40 (PAO1-L IC50 0.25 ± 0.12 μM, PA14 IC50 0.34 ± 0.03 μM) is one of the most potent PqsR antagonists reported showing significant inhibition of P. aeruginosa pyocyanin production and pqs system signaling in both planktonic cultures and biofilms. The co-crystal structure of 40 with the PqsR ligand binding domain revealed the specific binding interactions occurring between inhibitor and this key regulatory protein.

Abstract Objective Pseudomonas aeruginosa is an opportunistic bacterial pathogen well known to cause chronic lung infections in individuals with cystic fibrosis (CF). Some strains adapted to this particular niche show distinct phenotypes, such as biofilm hyperproduction. It is necessary to study CF clinical P. aeruginosa isolates, such as Liverpool Epidemic Strains (LES), to acquire a better understanding of the key genes essential for in vivo maintenance and the major virulence mechanisms involved in CF lung infections. Previously, a library of 9216 mutants of the LESB58 strain were generated by signature-tagged mutagenesis (STM) and screened in the rat model of chronic lung infection, allowing the identification of 163 STM mutants showing defects in in vivo maintenance. Results In the present study, these 163 mutants were successively screened in two additional surrogate host models (the amoeba and the fruit fly). The STM PALES_11731 mutant was the unique non-virulent in the three hosts. A competitive index study in rat lungs confirmed that the mutant was 20-fold less virulent than the wild-type strain. This study demonstrated the pertinence to use a multi-host approach to study the genetic determinants of P. aeruginosa strains infecting CF patients.

The emergence of multidrug-resistant bacterial strains worldwide has become a serious problem for public health over recent decades. The increase in antimicrobial resistance has been expanding via plasmids as mobile genetic elements encoding antimicrobial resistance (AMR) genes that are transferred vertically and horizontally. This study focuses on Salmonella enterica, one of the leading foodborne pathogens in industrialized countries. S. enterica is known to carry several plasmids involved not only in virulence but also in AMR. In the current paper, we present an integrated strategy to detect plasmid scaffolds in whole genome sequencing (WGS) assemblies. We developed a two-step procedure to predict plasmids based on i) the presence of essential elements for plasmid replication and mobility, as well as ii) sequence similarity to a reference plasmid. Next, to confirm the accuracy of the prediction in 1750 S. enterica short-read sequencing data, we combined Oxford Nanopore MinION long-read sequencing with Illumina MiSeq short-read sequencing in hybrid assemblies for 84 isolates to evaluate the proportion of plasmid that has been detected. At least one scaffold with an origin of replication (ORI) was predicted in 61.3% of the Salmonella isolates tested. The results indicated that IncFII and IncI1 ORIs were distributed in many S. enterica serotypes and were the most prevalent AMR genes carrier, whereas IncHI2A/IncHI2 and IncA/C2 were more serotype restricted but bore several AMR genes. Comparison between hybrid and short-read assemblies revealed that 81.1% of plasmids were found in the short-read sequencing using our pipeline. Through this process, we established that plasmids are prevalent in S. enterica and we also substantially expand the AMR genes in the resistome of this species.

Bacteriophages have shown promise as therapeutic alternatives to antibiotics for the control of infectious bacteria, including the human pathogen Salmonella. However, the development of effective phage-based applications requires the elucidation of key interactions between phages and target hosts, particularly since host resistance to phage is inevitable. Little is known about the alteration of host phenotypes following the development of resistance to phage. The aim of this study is to evaluate the antibiotic susceptibility and virulence of a Salmonella isolate following the development of resistance to bacteriophage SI1. We observed enhanced susceptibility to tetracycline and decreased invasion capacity in a differentiated Caco-2 intestinal cell line. Whole genome sequence analysis revealed an array of mutations, most notably, truncations in vgrG1_2, a core gene involved in Type VI secretion and mutations in the lipopolysaccharide, thereby indicating the plausible attachment site of phage SI1. These findings shed light on understanding the underlying mechanism for phage immunity within the host. Importantly, we reveal an associated genetic cost to the bacterial host with developing resistance to phages. Taken together, these results will aid in advancing strategies to delay or eliminate the development of host resistance when designing informed phage-based antimicrobials.

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.

Alfalfa sprouts have been linked to numerous North American outbreaks of Salmonella in recent years. Conventionally, treatments involving chlorine, heat, and irradiation are used for alfalfa seed sanitation. However, such treatments may be highly variable in their efficacy for pathogen control and/or detrimental to sprout quality, therefore negatively perceived by consumers advocating for natural alternatives. The usage of bacteriophages for pathogen control in sprouts has been previously explored, although with conflicting and inconsistent results. Lytic phages, viral predators of bacteria, represent an attractive approach as they provide several advantages compared to conventional treatments, such as their high specificity for bacterial targets and their ubiquity in nature. In this study, four Salmonella phages were isolated from British Columbia, Canada and characterized with respect to host range, burst size, latent period, and environmental stability to assess their potential to control Salmonella. Phage isolate SI1 showed the greatest host range, highest burst size and shortest latent period, greatest stability across all pH and temperatures and was the most effective in control of S. Enteritidis in vitro. Therefore, SI1 was chosen for treatment of sprouting alfalfa seeds artificially contaminated with S. Enteritidis with a multiplicity of infection (MOI) of ∼110 PFU/CFU. A significant (p < 0.05) reduction of 38.3 ± 3.0% of viable Salmonella cells was observed following two h of phage treatment. On days two to six of the sprouting process, reductions of Salmonella were also observed, but were not significant compared to the control (p > 0.05). It was further demonstrated that the sprout yield was not significantly (p > 0.05) affected by phage treatment. These results highlight the potential of phages recovered from the British Columbia environment for use as biocontrol agents against Salmonella, although differing efficacies in vitro was observed. Moreover, the effectiveness of SI1 to significantly (p < 0.05) control Salmonella on sprouting alfalfa seeds on day 1 of treatment was demonstrated. Although promising, future work should aim to optimize this treatment to achieve more effective, and longer lasting, biocontrol of Salmonella in sprouting alfalfa seeds.

Salmonella enterica is a bacterial species that is a major cause of illness in humans and food-producing animals. S. enterica exhibits considerable inter-serovar diversity, as evidenced by the large number of host adapted serovars that have been identified. The development of methods to assess genome diversity in S. enterica will help to further define the limits of diversity in this foodborne pathogen. Thus, we evaluated a PCR assay, which targets prophage integrase genes, as a rapid method to investigate S. enterica genome diversity. To evaluate the PCR prophage integrase assay, 49 isolates of S. enterica were selected, including 19 clinical isolates from clonal serovars (Enteritidis and Heidelberg) that commonly cause human illness, and 30 isolates from food-associated Salmonella serovars that rarely cause human illness. The number of integrase genes identified by the PCR assay was compared to the number of integrase genes within intact prophages identified by whole genome sequencing and phage finding program PHASTER. The PCR assay identified a total of 147 prophage integrase genes within the 49 S. enterica genomes (79 integrase genes in the food-associated Salmonella isolates, 50 integrase genes in S. Enteritidis, and 18 integrase genes in S. Heidelberg). In comparison, whole genome sequencing and PHASTER identified a total of 75 prophage integrase genes within 102 intact prophages in the 49 S. enterica genomes (44 integrase genes in the food-associated Salmonella isolates, 21 integrase genes in S. Enteritidis, and 9 integrase genes in S. Heidelberg). Collectively, both the PCR assay and PHASTER identified the presence of a large diversity of prophage integrase genes in the food-associated isolates compared to the clinical isolates, thus indicating a high degree of diversity in the food-associated isolates, and confirming the clonal nature of S. Enteritidis and S. Heidelberg. Moreover, PHASTER revealed a diversity of 29 different types of prophages and 23 different integrase genes within the food-associated isolates, but only identified four different phages and integrase genes within clonal isolates of S. Enteritidis and S. Heidelberg. These results demonstrate the potential usefulness of PCR based detection of prophage integrase genes as a rapid indicator of genome diversity in S. enterica.

The Salmonella Syst-OMICS consortium is sequencing 4,500 Salmonella genomes and building an analysis pipeline for the study of Salmonella genome evolution, antibiotic resistance and virulence genes. Metadata, including phenotypic as well as genomic data, for isolates of the collection are provided through the Salmonella Foodborne Syst-OMICS database (SalFoS), at https://salfos.ibis.ulaval.ca/. Here, we present our strategy and the analysis of the first 3,377 genomes. Our data will be used to draw potential links between strains found in fresh produce, humans, animals and the environment. The ultimate goals are to understand how Salmonella evolves over time, improve the accuracy of diagnostic methods, develop control methods in the field, and identify prognostic markers for evidence-based decisions in epidemiology and surveillance.

Phages infecting members of the opportunistic human pathogen, Salmonella enterica, are widespread in natural environments and offer a potential source of agents that could be used for controlling populations of this bacterium; yet, relatively little is known about these phages. Here we describe the isolation and characterization of 45 phages of Salmonella enterica from disparate geographic locations within British Columbia, Canada. Host-range profiling revealed host-specific patterns of susceptibility and resistance, with several phages identified that have a broad-host range (i.e., able to lyse >40% of bacterial hosts tested). One phage in particular, SE13, is able to lyse 51 out of the 61 Salmonella strains tested. Comparative genomic analyses also revealed an abundance of sequence diversity in the sequenced phages. Alignment of the genomes grouped the phages into 12 clusters with three singletons. Phages within certain clusters exhibited extraordinarily high genome homology (>98% nucleotide identity), yet between clusters, genomes exhibited a span of diversity (

Here we report the construction of three different vectors for the identification of bacterial genes induced in vitro and/or in vivo. These plasmids contain kanamycin, gentamicin, or tetracycline resistance genes as selectable markers. A promoterless cat and an improved GFP (mut3-gfp) can be used to follow the induction of gene expression by measuring chloramphenicol resistance and fluorescence, respectively.

We describe a fast and simple method for signature-tagged mutagenesis (STM) using defined oligonucleotides for tag construction into mini-Tn5 and PCR instead of hybridization for rapid screening of bacterial mutants in vivo. A collection of 12 unique 21-mers were synthesized as complementary DNA strands to tag bacterial mutants constructed by insertional mutagenesis using pUTmini-Tn5Km2 plasmids. Tags were tested in a combination of assays by PCR and compared to hybridization for specificity and for large-scale screening. Each defined tag has the same melting temperature, an invariable region to optimize PCRs and a variable region for specific amplification by PCR. A series of “suicide” plasmids carrying mini-Tn5s, each with a specific tag, were transferred into Pseudomonas aeruginosa, giving 12 libraries of mutants; groups of 12 mutants were pooled and arrayed into 96-well microplates, representing approximately one-sixth of the P. aeruginosa 5.9-Mb genome. This simple STM method can be adapted to any bacterial system and used for genome scanning in various growth conditions.