Mobile Genetic Element Flexibility as an Underlying Principle to Bacterial Evolution.

Mobile genetic elements are key to the evolution of bacteria and traits that affect host and ecosystem health. Here, we use a framework of a hierarchical and modular system that scales from genes to populations to synthesize recent findings on mobile genetic elements (MGEs) of bacteria. Doing so highlights the role that emergent properties of flexibility, robustness, and genetic capacitance of MGEs have on the evolution of bacteria. Some of their traits can be stored, shared, and diversified across different MGEs, taxa of bacteria, and time. Collectively, these properties contribute to maintaining functionality against perturbations while allowing changes to accumulate in order to diversify and give rise to new traits. These properties of MGEs have long challenged our abilities to study them. Implementation of new technologies and strategies allows for MGEs to be analyzed in new and powerful ways.

Natural variation of immune epitopes reveals intrabacterial antagonism.

Plants and animals detect biomolecules termed microbe-associated molecular patterns (MAMPs) and induce immunity. Agricultural production is severely impacted by pathogens which can be controlled by transferring immune receptors. However, most studies use a single MAMP epitope and the impact of diverse multicopy MAMPs on immune induction is unknown. Here, we characterized the epitope landscape from five proteinaceous MAMPs across 4,228 plant-associated bacterial genomes. Despite the diversity sampled, natural variation was constrained and experimentally testable. Immune perception in both Arabidopsis and tomato depended on both epitope sequence and copy number variation. For example, Elongation Factor Tu is predominantly single copy, and 92% of its epitopes are immunogenic. Conversely, 99.9% of bacterial genomes contain multiple cold shock proteins, and 46% carry a nonimmunogenic form. We uncovered a mechanism for immune evasion, intrabacterial antagonism, where a nonimmunogenic cold shock protein blocks perception of immunogenic forms encoded in the same genome. These data will lay the foundation for immune receptor deployment and engineering based on natural variation.

Biological activities and metabolomic profiles of extracts from the marine sediment bacterium Nocardiopsis alba DP1B cultivated in different media.

The soil bacterium DP1B was isolated from a marine sediment collected off the coast of Randayan Island, Kalimantan Barat, Indonesia and identified based on 16S rDNA as Nocardiopsis alba. The bacterium was cultivated in seven different media (A1, ISP1, ISP2, ISP4, PDB, PC-1, and SCB) with three different solvents [distilled water, 5 % NaCl solution, artificial seawater (ASW)] combinations, shaken at 200 rpm, 30 °C, for 7 days. The culture broths were extracted with ethyl acetate and each extract was tested for its antimicrobial activity and brine shrimp lethality, and the chemical diversity was assessed using thin-layer chromatography (TLC), gas chromatography (GC), and liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS). The result showed that almost all extracts showed antibacterial but not antifungal activity, whereas their brine shrimp toxicity levels vary from high to low. The best medium/solvent combinations for antibacterial activity and toxicity were PC-1 (in either distilled water, 5% NaCl solution, or ASW) and SCB in ASW. Different chemical diversity profiles were observed on TLC, GC-MS, and LC-MS/MS. Extracts from the PC-1 cultures seem to contain a significant number of cyclic dipeptides, whereas those from the SCB cultures contain sesquiterpenes, indicating that media and solvent compositions can affect the secondary metabolite profiles of DP1B. In addition, untargeted metabolomic analyses using LC-MS/MS showed many molecular ions that did not match with those in the Global Natural Products Social Molecular Networking (GNPS) database, suggesting that DP1B has great potential as a source of new natural products.

Dual adhesive unipolar polysaccharides synthesized by overlapping biosynthetic pathways in Agrobacterium tumefaciens.

Agrobacterium tumefaciens is a member of the Alphaproteobacteria that pathogenises plants and associates with biotic and abiotic surfaces via a single cellular pole. A. tumefaciens produces the unipolar polysaccharide (UPP) at the site of surface contact. UPP production is normally surface-contact inducible, but elevated levels of the second messenger cyclic diguanylate monophosphate (cdGMP) bypass this requirement. Multiple lines of evidence suggest that the UPP has a central polysaccharide component. Using an A. tumefaciens derivative with elevated cdGMP and mutationally disabled for other dispensable polysaccharides, a series of related genetic screens have identified a large number of genes involved in UPP biosynthesis, most of which are Wzx-Wzy-type polysaccharide biosynthetic components. Extensive analyses of UPP production in these mutants have revealed that the UPP is composed of two genetically, chemically, and spatially discrete forms of polysaccharide, and that each requires a specific Wzy-type polymerase. Other important biosynthetic, processing, and regulatory functions for UPP production are also revealed, some of which are common to both polysaccharides, and a subset of which are specific to each type. Many of the UPP genes identified are conserved among diverse rhizobia, whereas others are more lineage specific.

Genome features and secondary metabolite potential of the marine symbiont Streptomyces sp. RS2.

Streptomyces sp. RS2 was isolated from an unidentified sponge collected around Randayan Island, Indonesia. The genome of Streptomyces sp. RS2 consists of a linear chromosome of 9,391,717 base pairs with 71.9% of G + C content, 8270 protein-coding genes, as well as 18 rRNA and 85 tRNA loci. Twenty-eight putative secondary metabolites biosynthetic gene clusters (BGCs) were identified in the genome sequence. Nine of them have 100% similarity to BGCs for albaflavenone, α-lipomycin, coelibactin, coelichelin, ectoine, geosmin, germicidin, hopene, and lanthionine (SapB). The remaining 19 BGCs have low (< 50%) or moderate (50-80%) similarity to other known secondary metabolite BGCs. Biological activity assays of extracts from 21 different cultures of the RS2 strain showed that SCB ASW was the best medium for the production of antimicrobial and cytotoxic compounds. Streptomyces sp. RS2 has great potential to be a producer of novel secondary metabolites, particularly those with antimicrobial and antitumor activities.

Whole Genome Sequencing-Based Tracing of a 2022 Introduction and Outbreak of Xanthomonas hortorum pv. pelargonii.

Globalization has made agricultural commodities more accessible, available, and affordable. However, their global movement increases the potential for invasion by pathogens and necessitates development and implementation of sensitive, rapid, and scalable surveillance methods. Here, we used 35 strains, isolated by multiple diagnostic laboratories, as a case study for using whole genome sequence data in a plant disease diagnostic setting. Twenty-seven of the strains were isolated in 2022 and identified as Xanthomonas hortorum pv. pelargonii. Eighteen of these strains originated from material sold by a plant breeding company that had notified clients following a release of infected geranium cuttings. Analyses of whole genome sequences revealed epidemiological links among the 27 strains from different growers that confirmed a common source of the outbreak and uncovered likely secondary spread events within facilities that housed plants originating from different plant breeding companies. Whole genome sequencing data were also analyzed to reveal how preparatory and analytical methods can impact conclusions on outbreaks of clonal pathogenic strains. The results demonstrate the potential power of using whole genome sequencing among a network of diagnostic labs and highlight how sharing such data can help shorten response times to mitigate outbreaks more expediently and precisely than standard methods.

Modular evolution of secretion systems and virulence plasmids in a bacterial species complex.

BACKGROUND: Many named species as defined in current bacterial taxonomy correspond to species complexes. Uncertainties regarding the organization of their genetic diversity challenge research efforts. We utilized the Agrobacterium tumefaciens species complex (a.k.a. Agrobacterium biovar 1), a taxon known for its phytopathogenicity and applications in transformation, as a study system and devised strategies for investigating genome diversity and evolution of species complexes. RESULTS: We utilized 35 genome assemblies, including 14 newly generated ones, to achieve a phylogenetically balanced sampling of A. tumefaciens. Our genomic analysis suggested that the 10 genomospecies described previously are distinct biological species and supported a quantitative guideline for species delineation. Furthermore, our inference of gene content and core-genome phylogeny allowed for investigations of genes critical in fitness and ecology. For the type VI secretion system (T6SS) involved in interbacterial competition and thought to be conserved, we detected multiple losses and one horizontal gene transfer. For the tumor-inducing plasmids (pTi) and pTi-encoded type IV secretion system (T4SS) that are essential for agrobacterial phytopathogenicity, we uncovered novel diversity and hypothesized their involvement in shaping this species complex. Intriguingly, for both T6SS and T4SS, genes encoding structural components are highly conserved, whereas extensive diversity exists for genes encoding effectors and other proteins. CONCLUSIONS: We demonstrate that the combination of a phylogeny-guided sampling scheme and an emphasis on high-quality assemblies provides a cost-effective approach for robust analysis in evolutionary genomics. We show that the T6SS VgrG proteins involved in specific effector binding and delivery can be classified into distinct types based on domain organization. The co-occurrence patterns of VgrG-associated domains and the neighboring genes that encode different chaperones/effectors can be used to infer possible interacting partners. Similarly, the associations between plant host preference and the pTi type among these strains can be used to infer phenotype-genotype correspondence. Our strategies for multi-level investigations at scales that range from whole genomes to intragenic domains and phylogenetic depths from between- to within-species are applicable to other bacteria. Furthermore, modularity observed in the molecular evolution of genes and domains is useful for inferring functional constraints and informing experimental works.

Streptomyces caniscabiei sp. nov., which causes potato common scab and is distributed across the world.

Fourteen strains of Streptomyces isolated from scab lesions on potato are described as members of a novel species based on genetic distance, morphological observation and biochemical analyses. Morphological and biochemical characteristics of these strains are distinct from other described phytopathogenic species. Strain NE06-02DT has white aerial mycelium and grey, cylindrical, smooth spores on rectus-flexibilis spore chains. Members of this species group can utilize most of the International Streptomyces Project sugars, utilize melibiose and trehalose, produce melanin, grow on 6-7 % NaCl and pH 5-5.5 media, and are susceptible to oleandomycin (100 µg ml-1), streptomycin (20 µg ml-1) and penicillin G (30 µg ml-1). Though the 16S rRNA gene sequences from several members of this novel species are identical to the Streptomyces bottropensis 16S rRNA gene sequence, whole-genome average nucleotide identity and multi-locus sequence analysis confirm that the strains are members of a novel species. Strains belonging to this novel species have been isolated from the United States, Egypt and China with the earliest known members being isolated in 1961 from common scab lesions of potato in both California, USA, and Maine, USA. The name Streptomyces caniscabiei sp. nov. is proposed for strain NE06-02DT (=DSM111602T=ATCC TSD-236T) and the other members of this novel species group.

A Novel Species-Level Group of Streptomyces Exhibits Variation in Phytopathogenicity Despite Conservation of Virulence Loci.

The genus Streptomyces includes several phytopathogenic species that cause common scab, a devastating disease of tuber and root crops, in particular potato. The diversity of species that cause common scab is unknown. Likewise, the genomic context necessary for bacteria to incite common scab symptom development is not fully characterized. Here, we phenotyped and sequenced the genomes of five strains from a poorly studied Streptomyces lineage. These strains form a new species-level group. When genome sequences within just these five strains are compared, there are no polymorphisms of loci implicated in virulence. Each genome contains the pathogenicity island that encodes for the production of thaxtomin A, a phytotoxin necessary for common scab. Yet, not all sequenced strains produced thaxtomin A. Strains varied from nonpathogenic to highly virulent on two hosts. Unexpectedly, one strain that produced thaxtomin A and was pathogenic on radish was not aggressively pathogenic on potato. Therefore, while thaxtomin A biosynthetic genes and production of thaxtomin A are necessary, they are not sufficient for causing common scab of potato. Additionally, results show that even within a species-level group of Streptomyces strains, there can be aggressively pathogenic and nonpathogenic strains despite conservation of virulence genes.

Experimental evolution can enhance benefits of rhizobia to novel legume hosts.

Legumes preferentially associate with and reward beneficial rhizobia in root nodules, but the processes by which rhizobia evolve to provide benefits to novel hosts remain poorly understood. Using cycles of in planta and in vitro evolution, we experimentally simulated lifestyles where rhizobia repeatedly interact with novel plant genotypes with which they initially provide negligible benefits. Using a full-factorial replicated design, we independently evolved two rhizobia strains in associations with each of two Lotus japonicus genotypes that vary in regulation of nodule formation. We evaluated phenotypic evolution of rhizobia by quantifying fitness, growth effects and histological features on hosts, and molecular evolution via genome resequencing. Rhizobia evolved enhanced host benefits and caused changes in nodule development in one of the four host-symbiont combinations, that appeared to be driven by reduced costs during symbiosis, rather than increased nitrogen fixation. Descendant populations included genetic changes that could alter rhizobial infection or proliferation in host tissues, but lack of evidence for fixation of these mutations weakens the results. Evolution of enhanced rhizobial benefits occurred only in a subset of experiments, suggesting a role for host-symbiont genotype interactions in mediating the evolution of enhanced benefits from symbionts.

Association mapping, transcriptomics, and transient expression identify candidate genes mediating plant-pathogen interactions in a tree.

Invasive microbes causing diseases such as sudden oak death negatively affect ecosystems and economies around the world. The deployment of resistant genotypes for combating introduced diseases typically relies on breeding programs that can take decades to complete. To demonstrate how this process can be accelerated, we employed a genome-wide association mapping of ca 1,000 resequenced Populus trichocarpa trees individually challenged with Sphaerulina musiva, an invasive fungal pathogen. Among significant associations, three loci associated with resistance were identified and predicted to encode one putative membrane-bound L-type receptor-like kinase and two receptor-like proteins. A susceptibility-associated locus was predicted to encode a putative G-type D-mannose-binding receptor-like kinase. Multiple lines of evidence, including allele analysis, transcriptomics, binding assays, and overexpression, support the hypothesized function of these candidate genes in the P. trichocarpa response to S. musiva.

A DeoR-Type Transcription Regulator Is Required for Sugar-Induced Expression of Type III Secretion-Encoding Genes in Pseudomonas syringae pv. tomato DC3000.

The type III secretion system (T3SS) of plant-pathogenic Pseudomonas syringae is essential for virulence. Genes encoding the T3SS are not constitutively expressed and must be induced upon infection. Plant-derived metabolites, including sugars such as fructose and sucrose, are inducers of T3SS-encoding genes, yet the molecular mechanisms underlying perception of these host signals by P. syringae are unknown. Here, we report that sugar-induced expression of type III secretion A (setA), predicted to encode a DeoR-type transcription factor, is required for maximal sugar-induced expression of T3SS-associated genes in P. syringae DC3000. From a Tn5 transposon mutagenesis screen, we identified two independent mutants with insertions in setA. When both setA::Tn5 mutants were cultured in minimal medium containing fructose, genes encoding the T3SS master regulator HrpL and effector AvrRpm1 were expressed at lower levels relative to that of a wild-type strain. Decreased hrpL and avrRpm1 expression also occurred in a setA::Tn5 mutant in response to glucose, sucrose, galactose, and mannitol, demonstrating that setA is genetically required for T3SS induction by many different sugars. Expression of upstream regulators hrpR/S and rpoN was not altered in setA::Tn5, indicating that SetA positively regulates hrpL expression independently of increased transcription of these genes. In addition to decreased response to defined sugar signals, a setA::Tn5 mutant had decreased T3SS deployment during infection and was compromised in its ability to grow in planta and cause disease. These data suggest that SetA is necessary for P. syringae to effectively respond to T3SS-inducing sugar signals encountered during infection.

Recurrent mutualism breakdown events in a legume rhizobia metapopulation.

Bacterial mutualists generate major fitness benefits for eukaryotes, reshaping the host phenotype and its interactions with the environment. Yet, microbial mutualist populations are predicted to generate mutants that defect from providing costly services to hosts while maintaining the capacity to exploit host resources. Here, we examined the mutualist service of symbiotic nitrogen fixation in a metapopulation of root-nodulating Bradyrhizobium spp. that associate with the native legume Acmispon strigosus. We quantified mutualism traits of 85 Bradyrhizobium isolates gathered from a 700 km transect in California spanning 10 sampled A. strigosus populations. We clonally inoculated each Bradyrhizobium isolate onto A. strigosus hosts and quantified nodulation capacity and net effects of infection, including host growth and isotopic nitrogen concentration. Six Bradyrhizobium isolates from five populations were categorized as ineffective because they formed nodules but did not enhance host growth via nitrogen fixation. Six additional isolates from three populations failed to form root nodules. Phylogenetic reconstruction inferred two types of mutualism breakdown, including three to four independent losses of effectiveness and five losses of nodulation capacity on A. strigosus. The evolutionary and genomic drivers of these mutualism breakdown events remain poorly understood.

Tropical soils are a reservoir for fluorescent Pseudomonas spp. biodiversity.

Fluorescent Pseudomonas spp. are widely studied for their beneficial activities to plants. To explore the genetic diversity of Pseudomonas spp. in tropical regions, we collected 76 isolates from a Brazilian soil. Genomes were sequenced and compared to known strains, mostly collected from temperate regions. Phylogenetic analyses classified the isolates in the P. fluorescens (57) and P. putida (19) groups. Among the isolates in the P. fluorescens group, most (37) were classified in the P. koreensis subgroup and two in the P. jessenii subgroup. The remaining 18 isolates fell into two phylogenetic subclades distinct from currently recognized P. fluorescens subgroups, and probably represent new subgroups. Consistent with their phylogenetic distance from described subgroups, the genome sequences of strains in these subclades are asyntenous to the genome sequences of members of their neighbour subgroups. The tropical isolates have several functional genes also present in known fluorescent Pseudomonas spp. strains. However, members of the new subclades share exclusive genes not detected in other subgroups, pointing to the potential for novel functions. Additionally, we identified 12 potential new species among the 76 isolates from the tropical soil. The unexplored diversity found in the tropical soil is possibly related to biogeographical patterns.

Genome variations between rhizosphere and bulk soil ecotypes of a Pseudomonas koreensis population.

Bulk soil and rhizosphere are soil compartments selecting different microbial communities. However, it is unknown whether this selection also can change the genome content of specific bacterial taxa, splitting a population in distinct ecotypes. To answer this question we compared the genome sequences of 53 isolates obtained from sugarcane rhizosphere (28) and bulk soil (25). These isolates were previously classified in the Pseudomonas koreensis subgroup of the P. fluorescens complex. Phylogenomics showed a trend of separation between bulk soil and rhizosphere isolates. Discriminant analysis of principal components (DAPC) identified differences in the accessory genome of rhizosphere and bulk soil sub-populations. We found significant changes in gene frequencies distinguishing rhizosphere from bulk soil ecotypes, for example, enrichment of phosphatases and xylose utilization (xut) genes, respectively. Phenotypic assays and deletion of xutA gene indicated that accumulation of xut genes in the bulk soil sub-population provided a higher growth capacity in a d-xylose medium, supporting the corresponding genomic differences. Despite the clear differences distinguishing the two ecotypes, all 53 isolates were classified in a single 16S rRNA gene OTU. Collectively, our results revealed that the gene pool and ecological behavior of a bacterial population can be different for ecotypes living in neighbouring soil habitats.

Phylobiochemical characterization of class-Ib aspartate/prephenate aminotransferases reveals evolution of the plant arogenate phenylalanine pathway.

The aromatic amino acid Phe is required for protein synthesis and serves as the precursor of abundant phenylpropanoid plant natural products. While Phe is synthesized from prephenate exclusively via a phenylpyruvate intermediate in model microbes, the alternative pathway via arogenate is predominant in plant Phe biosynthesis. However, the molecular and biochemical evolution of the plant arogenate pathway is currently unknown. Here, we conducted phylogenetically informed biochemical characterization of prephenate aminotransferases (PPA-ATs) that belong to class-Ib aspartate aminotransferases (AspAT Ibs) and catalyze the first committed step of the arogenate pathway in plants. Plant PPA-ATs and succeeding arogenate dehydratases (ADTs) were found to be most closely related to homologs from Chlorobi/Bacteroidetes bacteria. The Chlorobium tepidum PPA-AT and ADT homologs indeed efficiently converted prephenate and arogenate into arogenate and Phe, respectively. A subset of AspAT Ib enzymes exhibiting PPA-AT activity was further identified from both Plantae and prokaryotes and, together with site-directed mutagenesis, showed that Thr-84 and Lys-169 play key roles in specific recognition of dicarboxylic keto (prephenate) and amino (aspartate) acid substrates. The results suggest that, along with ADT, a gene encoding prephenate-specific PPA-AT was transferred from a Chlorobi/Bacteroidetes ancestor to a eukaryotic ancestor of Plantae, allowing efficient Phe and phenylpropanoid production via arogenate in plants today.

A Proposal for a Genome Similarity-Based Taxonomy for Plant-Pathogenic Bacteria that Is Sufficiently Precise to Reflect Phylogeny, Host Range, and Outbreak Affiliation Applied to Pseudomonas syringae sensu lato as a Proof of Concept.

Taxonomy of plant pathogenic bacteria is challenging because pathogens of different crops often belong to the same named species but current taxonomy does not provide names for bacteria below the subspecies level. The introduction of the host range-based pathovar system in the 1980s provided a temporary solution to this problem but has many limitations. The affordability of genome sequencing now provides the opportunity for developing a new genome-based taxonomic framework. We already proposed to name individual bacterial isolates based on pairwise genome similarity. Here, we expand on this idea and propose to use genome similarity-based codes, which we now call life identification numbers (LINs), to describe and name bacterial taxa. Using 93 genomes of Pseudomonas syringae sensu lato, LINs were compared with a P. syringae genome tree whereby the assigned LINs were found to be informative of a majority of phylogenetic relationships. LINs also reflected host range and outbreak association for strains of P. syringae pathovar actinidiae, a pathovar for which many genome sequences are available. We conclude that LINs could provide the basis for a new taxonomic framework to address the shortcomings of the current pathovar system and to complement the current taxonomic system of bacteria in general.

Isothermal Amplification and Lateral-Flow Assay for Detecting Crown-Gall-Causing Agrobacterium spp.

Agrobacterium is a genus of soilborne gram-negative bacteria. Members carrying oncogenic plasmids can cause crown gall disease, which has significant economic costs, especially for the orchard and nursery industries. Early and rapid detection of pathogenic Agrobacterium spp. is key to the management of crown gall disease. To this end, we designed oligonucleotide primers and probes to target virD2 for use in a molecular diagnostic tool that relies on isothermal amplification and lateral-flow-based detection. The oligonucleotide tools were tested in the assay and evaluated for detection limit and specificity in detecting alleles of virD2. One set of primers that successfully amplified virD2 when used with an isothermal recombinase was selected. Both tested probes had detection limits in picogram amounts of DNA. Probe 1 could detect all tested pathogenic isolates that represented most of the diversity of virD2. Finally, the coupling of lateral-flow detection to the use of these oligonucleotide primers in isothermal amplification helped to reduce the onerousness of the process, and alleviated reliance on specialized tools necessary for molecular diagnostics. The assay is an advancement for the rapid molecular detection of pathogenic Agrobacterium spp.

Natural variation of immune epitopes reveals intrabacterial antagonism.

Plants and animals detect biomolecules termed Microbe-Associated Molecular Patterns (MAMPs) and induce immunity. Agricultural production is severely impacted by pathogens which can be controlled by transferring immune receptors. However, most studies use a single MAMP epitope and the impact of diverse multi-copy MAMPs on immune induction is unknown. Here we characterized the epitope landscape from five proteinaceous MAMPs across 4,228 plant-associated bacterial genomes. Despite the diversity sampled, natural variation was constrained and experimentally testable. Immune perception in both Arabidopsis and tomato depended on both epitope sequence and copy number variation. For example, Elongation Factor Tu is predominantly single copy and 92% of its epitopes are immunogenic. Conversely, 99.9% of bacterial genomes contain multiple Cold Shock Proteins and 46% carry a non-immunogenic form. We uncovered a new mechanism for immune evasion, intrabacterial antagonism, where a non-immunogenic Cold Shock Protein blocks perception of immunogenic forms encoded in the same genome. These data will lay the foundation for immune receptor deployment and engineering based on natural variation.

Genomic characterization of Listeria monocytogenes recovered from dairy facilities in British Columbia, Canada from 2007 to 2017.

Listeria monocytogenes is a foodborne pathogen of concern in dairy processing facilities, with the potential to cause human illness and trigger regulatory actions if found in the product. Monitoring for Listeria spp. through environmental sampling is recommended to prevent establishment of these microorganisms in dairy processing environments, thereby reducing the risk of product contamination. To inform on L. monocytogenes diversity and transmission, we analyzed genome sequences of L. monocytogenes strains (n = 88) obtained through the British Columbia Dairy Inspection Program. Strains were recovered from five different dairy processing facilities over a 10 year period (2007-2017). Analysis of whole genome sequences (WGS) grouped the isolates into nine sequence types and 11 cgMLST types (CT). The majority of isolates (93%) belonged to lineage II. Within each CT, single nucleotide polymorphism (SNP) differences ranged from 0 to 237 between isolates. A highly similar (0-16 SNPs) cluster of over 60 isolates, collected over 9 years within one facility (#71), was identified suggesting a possible persistent population. Analyses of genome content revealed a low frequency of genes associated with stress tolerance, with the exception of widely disseminated cadmium resistance genes cadA1 and cadA2. The distribution of virulence genes and mutations within internalin genes varied across the isolates and facilities. Further studies are needed to elucidate their phenotypic effect on pathogenicity and stress response. These findings demonstrate the diversity of L. monocytogenes isolates across dairy facilities in the same region. Findings also showed the utility of using WGS to discern potential persistence events within a single facility over time.