Assessment of plasmids for relating the 2020 Salmonella enterica serovar Newport onion outbreak to farms implicated by the outbreak investigation.

BACKGROUND: The Salmonella enterica serovar Newport red onion outbreak of 2020 was the largest foodborne outbreak of Salmonella in over a decade. The epidemiological investigation suggested two farms as the likely source of contamination. However, single nucleotide polymorphism (SNP) analysis of the whole genome sequencing data showed that none of the Salmonella isolates collected from the farm regions were linked to the clinical isolates-preventing the use of phylogenetics in source identification. Here, we explored an alternative method for analyzing the whole genome sequencing data driven by the hypothesis that if the outbreak strain had come from the farm regions, then the clinical isolates would disproportionately contain plasmids found in isolates from the farm regions due to horizontal transfer. RESULTS: SNP analysis confirmed that the clinical isolates formed a single, nearly-clonal clade with evidence for ancestry in California going back a decade. The clinical clade had a large core genome (4,399 genes) and a large and sparsely distributed accessory genome (2,577 genes, at least 64% on plasmids). At least 20 plasmid types occurred in the clinical clade, more than were found in the literature for Salmonella Newport. A small number of plasmids, 14 from 13 clinical isolates and 17 from 8 farm isolates, were found to be highly similar (> 95% identical)-indicating they might be related by horizontal transfer. Phylogenetic analysis was unable to determine the geographic origin, isolation source, or time of transfer of the plasmids, likely due to their promiscuous and transient nature. However, our resampling analysis suggested that observing a similar number and combination of highly similar plasmids in random samples of environmental Salmonella enterica within the NCBI Pathogen Detection database was unlikely, supporting a connection between the outbreak strain and the farms implicated by the epidemiological investigation. CONCLUSION: Horizontally transferred plasmids provided evidence for a connection between clinical isolates and the farms implicated as the source of the outbreak. Our case study suggests that such analyses might add a new dimension to source tracking investigations, but highlights the need for detailed and accurate metadata, more extensive environmental sampling, and a better understanding of plasmid molecular evolution.

Genomic analysis of Listeria monocytogenes from US food processing environments reveals a high prevalence of QAC efflux genes but limited evidence of their contribution to environmental persistence.

BACKGROUND: Quaternary ammonium compound (QAC) efflux genes increase the minimum inhibitory concentration of Listeria monocytogenes (Lm) to benzalkonium chloride sanitizer, but the contribution of these genes to persistence in food processing environments is unclear. The goal of this study was to leverage genomic data and associated metadata for 4969 Lm isolates collected between 1999 and 2019 to: (1) evaluate the prevalence of QAC efflux genes among Lm isolates from diverse US food processors, (2) use comparative genomic analyses to assess confounding factors, such as clonal complex identity and stress tolerance genotypes, and (3) identify patterns in QAC efflux gene gain and loss among persistent clones within specific facilities over time. RESULTS: The QAC efflux gene cassette bcrABC was present in nearly half (46%) of all isolates. QAC efflux gene prevalence among isolates was associated with clonal complex (𝛘2 < 0.001) and clonal complex was associated with the facility type (𝛘2 < 0.001). Consequently, changes in the prevalence of QAC efflux genes within individual facilities were generally attributable to changes in the prevalence of specific clonal complexes. Additionally, a GWAS and targeted BLAST search revealed that clonal complexes with a high prevalence of QAC efflux genes commonly possessed other stress tolerance genes. For example, a high prevalence of bcrABC in a clonal complex was significantly associated with the presence of the SSI-1 gene cluster (p < 0.05). QAC efflux gene gain and loss were both observed among persistent populations of Lm in individual facilities, suggesting a limited direct role for QAC efflux genes as predictors of persistence. CONCLUSION: This study suggests that although there is evidence that QAC efflux genes are part of a suite of adaptations common among Lm isolated from some food production environments, these genes may be neither sufficient nor necessary to enhance persistence. This is a crucial distinction for decision making in the food industry. For example, changes to sanitizer regimen targeting QAC tolerance would not address other contributing genetic or non-genetic factors, such as equipment hygienic design which physically mediates sanitizer exposure.

Interpretative Labor and the Bane of Nonstandardized Metadata in Public Health Surveillance and Food Safety.

Open-source DNA sequence databases have long been touted as beneficial to public health, including the facilitation of earlier detection and response to infectious disease outbreaks. Of critical importance to harnessing these benefits is the metadata that describe general and other domain-specific attributes (eg, collection location, isolate type) of a sample. Unlike the sequence data, metadata are often incomplete and lack adherence to an international standard. Here, we describe the problem posed by such variable and incomplete metadata in terms of interpretative labor costs (the time and energy necessary to make sense of the signal in the genetic data) and the impact such metadata have on foodborne outbreak detection and response. Improving the quality of sequence-associated metadata would allow for earlier detection of emerging food safety hazards and allow faster response to foodborne outbreaks.

Phylogenetic and Biogeographic Patterns of Vibrio parahaemolyticus Strains from North America Inferred from Whole-Genome Sequence Data.

Vibrio parahaemolyticus is the most common cause of seafood-borne illness reported in the United States. The draft genomes of 132 North American clinical and oyster V. parahaemolyticus isolates were sequenced to investigate their phylogenetic and biogeographic relationships. The majority of oyster isolate sequence types (STs) were from a single harvest location; however, four were identified from multiple locations. There was population structure along the Gulf and Atlantic Coasts of North America, with what seemed to be a hub of genetic variability along the Gulf Coast, with some of the same STs occurring along the Atlantic Coast and one shared between the coastal waters of the Gulf and those of Washington State. Phylogenetic analyses found nine well-supported clades. Two clades were composed of isolates from both clinical and oyster sources. Four were composed of isolates entirely from clinical sources, and three were entirely from oyster sources. Each single-source clade consisted of one ST. Some human isolates lack tdh, trh, and some type III secretion system (T3SS) genes, which are established virulence genes of V. parahaemolyticus Thus, these genes are not essential for pathogenicity. However, isolates in the monophyletic groups from clinical sources were enriched in several categories of genes compared to those from monophyletic groups of oyster isolates. These functional categories include cell signaling, transport, and metabolism. The identification of genes in these functional categories provides a basis for future in-depth pathogenicity investigations of V. parahaemolyticusIMPORTANCEVibrio parahaemolyticus is the most common cause of seafood-borne illness reported in the United States and is frequently associated with shellfish consumption. This study contributes to our knowledge of the biogeography and functional genomics of this species around North America. STs shared between the Gulf Coast and the Atlantic seaboard as well as Pacific waters suggest possible transport via oceanic currents or large shipping vessels. STs frequently isolated from humans but rarely, if ever, isolated from the environment are likely more competitive in the human gut than other STs. This could be due to additional functional capabilities in areas such as cell signaling, transport, and metabolism, which may give these isolates an advantage in novel nutrient-replete environments such as the human gut.

Investigating the Distribution of Strains of Erwinia amylovora and Streptomycin Resistance in Apple Orchards in New York Using Clustered Regularly Interspaced Short Palindromic Repeat Profiles: A 6-Year Follow-Up.

Fire blight, caused by the bacterium Erwinia amylovora, is one of the most important diseases of apple. The antibiotic streptomycin is routinely used in the commercial apple industries of New York (NY) and New England to manage the disease. In 2002 and again, from 2011 to 2014, outbreaks of streptomycin resistance (SmR) were reported and investigated in NY. Motivated by new grower reports of control failures, we conducted a follow-up investigation of the distribution of SmR and E. amylovora strains for major apple production regions of NY over the last 6 years (2015 to 2020). Characterization of clustered regularly interspaced short palindromic repeat (CRISPR) profiles revealed that a few "cosmopolitan" strains were widely prevalent across regions, whereas many other "resident" strains were confined to one location. In addition, we uncovered novel CRISPR profile diversity in all investigated regions. SmR E. amylovora was detected only in a small area spanning two counties from 2017 to 2020 and was always associated with one CRISPR profile (41:23:38), which matched the profile of SmR E. amylovora, discovered in 2002. This suggests the original SmR E. amylovora was never fully eradicated and went undetected because of several seasons of low disease pressure in this region. Investigation of several representative isolates under controlled greenhouse conditions indicated significant differences in aggressiveness on 'Gala' apples. Potential implications of strain differences include the propensity of strains to become distributed across wide geographic regions and associated resistance management practices. Results from this work will directly influence sustainable fire blight management recommendations for commercial apple industries in NY state and other regions.

Genomic diversity of Salmonella enterica isolated from papaya samples collected during multiple outbreaks in 2017.

In 2017, the US Food and Drug Administration investigated the sources of multiple outbreaks of salmonellosis. Epidemiologic and traceback investigations identified Maradol papayas as the suspect vehicles. During the investigations, the genomes of 55 Salmonella enterica that were isolated from papaya samples were sequenced. Serovar assignments and phylogenetic analysis placed the 55 isolates into ten distinct groups, each representing a different serovar. Within-serovar SNP differences are generally between 0 and 20 SNPs, while the median between-serovar distance is 51 812 SNPs. We observed two groups with SNP distances between 21 and 100 SNPs. These relatively large within-serovar SNP distances may indicate that the isolates represent either diverse populations or multiple, genetically distinct subpopulations. Further inspection of these cases with traceback evidence allowed us to identify an 11th population. We observed that high levels of genomic diversity from individual firms is possible, with one firm yielding five of the ten serovars. Also, high levels of diversity are possible within small geographic regions, as five of the serovars were isolated from papayas that originated from farms located in Armería and Tecomán, Colima. In addition, we identified AMR genes that are present in three of the serovars studied here (aph(3')-lb, aph(6)-ld, tet(C), fosA7, and qnrB19) and we detected the presence of the plasmid IncHI2A among S. Urbana isolates.

Outbreak of Tattoo-associated Nontuberculous Mycobacterial Skin Infections.

BACKGROUND: On 29 April 2015, the Florida Department of Health in Miami-Dade County (DOH Miami-Dade) was notified by a local dermatologist of 3 patients with suspected nontuberculous mycobacterial (NTM) infection after receiving tattoos at a local tattoo studio. METHODS: DOH Miami-Dade conducted interviews and offered testing, described below, to tattoo studio clients reporting rashes. Culture of clinical isolates and identification were performed at the Florida Bureau of Public Health Laboratories. Characterization of NTM was performed by the Centers for Disease Control and Prevention and the US Food and Drug Administration (FDA), respectively. Whole-genome sequencing (WGS) and single-nucleotide polymorphism (SNP) analyses were used to construct a phylogeny among 21 Mycobacterium isolates at the FDA. RESULTS: Thirty-eight of 226 interviewed clients were identified as outbreak-associated cases. Multivariate logistic regression revealed that individuals who reported gray tattoo ink in their tattoos were 8.2 times as likely to report a rash (95% confidence interval, 3.1-22.1). Multiple NTM species were identified in clinical and environmental specimens. Phylogenetic results from environmental samples and skin biopsies indicated that 2 Mycobacterium fortuitum isolates (graywash ink and a skin biopsy) and 11 Mycobacterium abscessus isolates (5 from the implicated bottle of graywash tattoo ink, 2 from tap water, and 4 from skin biopsies) were indistinguishable. In addition, Mycobacterium chelonae was isolated from 5 unopened bottles of graywash ink provided by 2 other tattoo studios in Miami-Dade County. CONCLUSIONS: WGS and SNP analyses identified the tap water and the bottle of graywash tattoo ink as the sources of the NTM infections.

TreeToReads - a pipeline for simulating raw reads from phylogenies.

BACKGROUND: Using phylogenomic analysis tools for tracking pathogens has become standard practice in academia, public health agencies, and large industries. Using the same raw read genomic data as input, there are several different approaches being used to infer phylogenetic tree. These include many different SNP pipelines, wgMLST approaches, k-mer algorithms, whole genome alignment and others; each of these has advantages and disadvantages, some have been extensively validated, some are faster, some have higher resolution. A few of these analysis approaches are well-integrated into the regulatory process of US Federal agencies (e.g. the FDA's SNP pipeline for tracking foodborne pathogens). However, despite extensive validation on benchmark datasets and comparison with other pipelines, we lack methods for fully exploring the effects of multiple parameter values in each pipeline that can potentially have an effect on whether the correct phylogenetic tree is recovered. RESULTS: To resolve this problem, we offer a program, TreeToReads, which can generate raw read data from mutated genomes simulated under a known phylogeny. This simulation pipeline allows direct comparisons of simulated and observed data in a controlled environment. At each step of these simulations, researchers can vary parameters of interest (e.g., input tree topology, amount of sequence divergence, rate of indels, read coverage, distance of reference genome, etc) to assess the effects of various parameter values on correctly calling SNPs and reconstructing an accurate tree. CONCLUSIONS: Such critical assessments of the accuracy and robustness of analytical pipelines are essential to progress in both research and applied settings.

Singleton Sequence Type 382, an Emerging Clonal Group of Listeria monocytogenes Associated with Three Multistate Outbreaks Linked to Contaminated Stone Fruit, Caramel Apples, and Leafy Green Salad.

Three multistate outbreaks between 2014 and 2016, involving case patients in and outside the United States, were linked to stone fruit, caramel apples, and packaged leafy green salad contaminated with Listeria monocytogenes singleton sequence type 382 (ST382), a serotype IVb-v1 clone with limited genomic divergence. Isolates from these outbreaks and other ST382 isolates not associated with these outbreaks were analyzed by whole-genome sequencing (WGS) analysis. The primary differences among ST382 strains were single nucleotide polymorphisms (SNPs). WGS analysis differentiated ST382 from a clonal complex 1 outbreak strain co-contaminating the caramel apples. WGS clustered food, environmental, and clinical isolates within each outbreak, and also differentiated among the three outbreak strains and epidemiologically unrelated ST382 isolates, which were indistinguishable by pulsed-field gel electrophoresis. ST382 appeared to be an emerging clone that began to diverge from its ancestor approximately 32 years before 2016. We estimated that there was 1.29 nucleotide substitution per genome (2.94 Mbp) per year for this clone.

Mating system divergence affects the distribution of sequence diversity within and among populations of recently diverged subspecies of Clarkia xantiana (Onagraceae).

PREMISE OF THE STUDY: The population biology of outcrossing and self-fertilizing taxa is thought to differ because of the advantage that selfers have in colonizing unoccupied sites where mates and pollinators may be limiting (Baker's Law). This reduced tendency for outcrossers to colonize new sites, along with their greater dependence on pollinators to disperse pollen, has the potential to differently influence the genetic diversity and structure of outcrossing and selfing populations. METHODS: We conducted a comparative population genetic study of two sister outcrossing and selfing subspecies of Clarkia xantiana that have very recently diverged. We used DNA sequence variation (>40 kb from eight nuclear loci) from large samples of individuals from 14 populations to assess geographic patterns of genetic diversity and make inferences about the demographic and colonization histories of each subspecies. KEY RESULTS: We show that sequence variation is strongly reduced across all selfing populations. The demographic history of selfing populations exhibits recent colonization bottlenecks, whereas such bottlenecks are rarely observed for the outcrosser. The greater effect of genetic drift in the selfer has resulted in strong population genetic structure, but with no pattern of isolation by distance. By contrast, the stronger effect of gene flow in the outcrosser has resulted in considerably less structure, but a significant pattern of isolation by distance. CONCLUSIONS: Taken together, our results suggest that selfing populations are not at migration-drift equilibrium, are affected by strong episodes of genetic drift during colonization, and experience little or no subsequent gene flow from other populations after those founder events.

Characterization and evolution of Salmonella CRISPR-Cas systems.

Prokaryotic CRISPR-Cas (clustered regularly interspaced short palindromic repeats and CRISPR-associated genes) systems provide adaptive immunity from invasive genetic elements and encompass three essential features: (i) cas genes, (ii) a CRISPR array composed of spacers and direct repeats and (iii) an AT-rich leader sequence upstream of the array. We performed in-depth sequence analysis of the CRISPR-Cas systems in >600 Salmonella, representing four clinically prevalent serovars. Each CRISPR-Cas feature is extremely conserved in the Salmonella, and the CRISPR1 locus is more highly conserved than CRISPR2. Array composition is serovar-specific, although no convincing evidence of recent spacer acquisition against exogenous nucleic acids exists. Only 12 % of spacers match phage and plasmid sequences and self-targeting spacers are associated with direct repeat variants. High nucleotide identity (>99.9 %) exists across the cas operon among isolates of a single serovar and in some cases this conservation extends across divergent serovars. These observations reflect historical CRISPR-Cas immune activity, showing that this locus has ceased undergoing adaptive events. Intriguingly, the high level of conservation across divergent serovars shows that the genetic integrity of these inactive loci is maintained over time, contrasting with the canonical view that inactive CRISPR loci degenerate over time. This thorough characterization of Salmonella CRISPR-Cas systems presents new insights into Salmonella CRISPR evolution, particularly with respect to cas gene conservation, leader sequences, organization of direct repeats and protospacer matches. Collectively, our data suggest that Salmonella CRISPR-Cas systems are no longer immunogenic; rather, their impressive conservation indicates they may have an alternative function in Salmonella.

Characterization and evolution of Salmonella CRISPR-Cas systems.

Prokaryotic CRISPR-Cas (clustered regularly interspaced short palindromic repeats and CRISPR-associated genes) systems provide adaptive immunity from invasive genetic elements and encompass three essential features: (i) cas genes, (ii) a CRISPR array composed of spacers and direct repeats and (iii) an AT-rich leader sequence upstream of the array. We performed in-depth sequence analysis of the CRISPR-Cas systems in >600 Salmonella, representing four clinically prevalent serovars. Each CRISPR-Cas feature is extremely conserved in the Salmonella, and the CRISPR1 locus is more highly conserved than CRISPR2. Array composition is serovar-specific, although no convincing evidence of recent spacer acquisition against exogenous nucleic acids exists. Only 12% of spacers match phage and plasmid sequences and self-targeting spacers are associated with direct repeat variants. High nucleotide identity (>99.9%) exists across the cas operon among isolates of a single serovar and in some cases this conservation extends across divergent serovars. These observations reflect historical CRISPR-Cas immune activity, showing that this locus has ceased undergoing adaptive events. Intriguingly, the high level of conservation across divergent serovars shows that the genetic integrity of these inactive loci is maintained over time, contrasting with the canonical view that inactive CRISPR loci degenerate over time. This thorough characterization of Salmonella CRISPR-Cas systems presents new insights into Salmonella CRISPR evolution, particularly with respect to cas gene conservation, leader sequences, organization of direct repeats and protospacer matches. Collectively, our data suggest that Salmonella CRISPR-Cas systems are no longer immunogenic; rather, their impressive conservation indicates they may have an alternative function in Salmonella.

Cilantro microbiome before and after nonselective pre-enrichment for Salmonella using 16S rRNA and metagenomic sequencing.

BACKGROUND: Salmonella enterica is a common cause of foodborne gastroenteritis in the United States and is associated with outbreaks in fresh produce such as cilantro. Salmonella culture-based detection methods are complex and time consuming, and improvments to increase detection sensitivity will benefit consumers. In this study, we used 16S rRNA sequencing to determine the microbiome of cilantro. We also investigated changes to the microbial community prior to and after a 24-hour nonselective pre-enrichment culture step commonly used by laboratory analysts to resuscitate microorganisms in foods suspected of contamination with pathogens. Cilantro samples were processed for Salmonella detection according to the method in the United States Food and Drug Administration Bacteriological Analytical Manual. Genomic DNA was extracted from culture supernatants prior to and after a 24-hour nonselective pre-enrichment step and 454 pyrosequencing was performed on 16S rRNA amplicon libraries. A database of Enterobacteriaceae 16S rRNA sequences was created, and used to screen the libraries for Salmonella, as some samples were known to be culture positive. Additionally, culture positive cilantro samples were examined for the presence of Salmonella using shotgun metagenomics on the Illumina MiSeq. RESULTS: Time zero uncultured samples had an abundance of Proteobacteria while the 24-hour enriched samples were composed mostly of Gram-positive Firmicutes. Shotgun metagenomic sequencing of Salmonella culture positive cilantro samples revealed variable degrees of Salmonella contamination among the sequenced samples. CONCLUSIONS: Our cilantro study demonstrates the use of high-throughput sequencing to reveal the microbiome of cilantro, and how the microbiome changes during the culture-based protocols employed by food safety laboratories to detect foodborne pathogens. Finding that culturing the cilantro shifts the microbiome to a predominance of Firmicutes suggests that changing our culture-based methods will improve detection sensitivity for foodborne enteric pathogens.

Proficiency testing for bacterial whole genome sequencing: an end-user survey of current capabilities, requirements and priorities.

The advent of next-generation sequencing (NGS) has revolutionised public health microbiology. Given the potential impact of NGS, it is paramount to ensure standardisation of 'wet' laboratory and bioinformatic protocols and promote comparability of methods employed by different laboratories and their outputs. Therefore, one of the ambitious goals of the Global Microbial Identifier (GMI) initiative (http://www.globalmicrobialidentifier.org/) has been to establish a mechanism for inter-laboratory NGS proficiency testing (PT). This report presents findings from the survey recently conducted by Working Group 4 among GMI members in order to ascertain NGS end-use requirements and attitudes towards NGS PT. The survey identified the high professional diversity of laboratories engaged in NGS-based public health projects and the wide range of capabilities within institutions, at a notable range of costs. The priority pathogens reported by respondents reflected the key drivers for NGS use (high burden disease and 'high profile' pathogens). The performance of and participation in PT was perceived as important by most respondents. The wide range of sequencing and bioinformatics practices reported by end-users highlights the importance of standardisation and harmonisation of NGS in public health and underpins the use of PT as a means to assuring quality. The findings of this survey will guide the design of the GMI PT program in relation to the spectrum of pathogens included, testing frequency and volume as well as technical requirements. The PT program for external quality assurance will evolve and inform the introduction of NGS into clinical and public health microbiology practice in the post-genomic era.

The impact of systemic and copper pesticide applications on the phyllosphere microflora of tomatoes.

BACKGROUND: Contamination of tomatoes by Salmonella can occur in agricultural settings. Little is currently understood about how agricultural inputs such as pesticide applications may impact epiphytic crop microflora and potentially play a role in contamination events. We examined the impact of two materials commonly used in Virginia tomato agriculture: acibenzolar-S-methyl (crop protectant) and copper oxychloride (pesticide) to identify the effects these materials may exert on baseline tomato microflora and on the incidence of three specific genera; Salmonella, Xanthomonas and Paenibacillus. RESULTS: Approximately 186 441 16S rRNA gene and 39 381 18S rRNA gene sequences per independent replicate were used to analyze the impact of the pesticide applications on tomato microflora. An average of 3 346 677 (634 892 974 bases) shotgun sequences per replicate were used for metagenomic analyses. CONCLUSION: A significant decrease in the presence of Gammaproteobacteria was observed between controls and copper-treated plants, suggesting that copper is effective at suppressing growth of certain taxa in this class. A higher mean abundance of Salmonella and Paenibacillus in control samples compared to treatments may suggest that both systemic and copper applications diminish the presence of these genera in the phyllosphere; however, owing to the lack of statistical significance, this could also be due to other factors. The most distinctive separation of shared membership was observed in shotgun data between the two different sampling time-points (not between treatments), potentially supporting the hypothesis that environmental pressures may exert more selective pressures on epiphytic microflora than do certain agricultural management practices.

Genetic diversity and evolution of Salmonella enterica serovar Enteritidis strains with different phage types.

Phage typing has been used for the epidemiological surveillance of Salmonella enterica serovar Enteritidis for over 2 decades. However, knowledge of the genetic and evolutionary relationships between phage types is very limited, making differences difficult to interpret. Here, single nucleotide polymorphisms (SNPs) identified from whole-genome comparisons were used to determine the relationships between some S. Enteritidis phage types (PTs) commonly associated with food-borne outbreaks in the United States. Emphasis was placed on the predominant phage types PT8, PT13a, and PT13 in North America. With >89,400 bp surveyed across 98 S. Enteritidis isolates representing 14 distinct phage types, 55 informative SNPs were discovered within 23 chromosomally anchored loci. To maximize the discriminatory and evolutionary partitioning of these highly homogeneous strains, sequences comprising informative SNPs were concatenated into a single combined data matrix and subjected to phylogenetic analysis. The resultant phylogeny allocated most S. Enteritidis isolates into two distinct clades (clades I and II) and four subclades. Synapomorphic (shared and derived) sets of SNPs capable of distinguishing individual clades/subclades were identified. However, individual phage types appeared to be evolutionarily disjunct when mapped to this phylogeny, suggesting that phage typing may not be valid for making phylogenetic inferences. Furthermore, the set of SNPs identified here represents useful genetic markers for strain differentiation of more clonal S. Enteritidis strains and provides core genotypic markers for future development of a SNP typing scheme with S. Enteritidis.

SARS-CoV-2 wastewater variant surveillance: pandemic response leveraging FDA's GenomeTrakr network.

Wastewater surveillance has emerged as a crucial public health tool for population-level pathogen surveillance. Supported by funding from the American Rescue Plan Act of 2021, the FDA's genomic epidemiology program, GenomeTrakr, was leveraged to sequence SARS-CoV-2 from wastewater sites across the United States. This initiative required the evaluation, optimization, development, and publication of new methods and analytical tools spanning sample collection through variant analyses. Version-controlled protocols for each step of the process were developed and published on protocols.io. A custom data analysis tool and a publicly accessible dashboard were built to facilitate real-time visualization of the collected data, focusing on the relative abundance of SARS-CoV-2 variants and sub-lineages across different samples and sites throughout the project. From September 2021 through June 2023, a total of 3,389 wastewater samples were collected, with 2,517 undergoing sequencing and submission to NCBI under the umbrella BioProject, PRJNA757291. Sequence data were released with explicit quality control (QC) tags on all sequence records, communicating our confidence in the quality of data. Variant analysis revealed wide circulation of Delta in the fall of 2021 and captured the sweep of Omicron and subsequent diversification of this lineage through the end of the sampling period. This project successfully achieved two important goals for the FDA's GenomeTrakr program: first, contributing timely genomic data for the SARS-CoV-2 pandemic response, and second, establishing both capacity and best practices for culture-independent, population-level environmental surveillance for other pathogens of interest to the FDA. IMPORTANCE: This paper serves two primary objectives. First, it summarizes the genomic and contextual data collected during a Covid-19 pandemic response project, which utilized the FDA's laboratory network, traditionally employed for sequencing foodborne pathogens, for sequencing SARS-CoV-2 from wastewater samples. Second, it outlines best practices for gathering and organizing population-level next generation sequencing (NGS) data collected for culture-free, surveillance of pathogens sourced from environmental samples.

Elucidating the evolutionary history and expression patterns of nucleoside phosphorylase paralogs (vegetative storage proteins) in Populus and the plant kingdom.

BACKGROUND: Nucleoside phosphorylases (NPs) have been extensively investigated in human and bacterial systems for their role in metabolic nucleotide salvaging and links to oncogenesis. In plants, NP-like proteins have not been comprehensively studied, likely because there is no evidence of a metabolic function in nucleoside salvage. However, in the forest trees genus Populus a family of NP-like proteins function as an important ecophysiological adaptation for inter- and intra-seasonal nitrogen storage and cycling. RESULTS: We conducted phylogenetic analyses to determine the distribution and evolution of NP-like proteins in plants. These analyses revealed two major clusters of NP-like proteins in plants. Group I proteins were encoded by genes across a wide range of plant taxa while proteins encoded by Group II genes were dominated by species belonging to the order Malpighiales and included the Populus Bark Storage Protein (BSP) and WIN4-like proteins. Additionally, we evaluated the NP-like genes in Populus by examining the transcript abundance of the 13 NP-like genes found in the Populus genome in various tissues of plants exposed to long-day (LD) and short-day (SD) photoperiods. We found that all 13 of the Populus NP-like genes belonging to either Group I or II are expressed in various tissues in both LD and SD conditions. Tests of natural selection and expression evolution analysis of the Populus genes suggests that divergence in gene expression may have occurred recently during the evolution of Populus, which supports the adaptive maintenance models. Lastly, in silico analysis of cis-regulatory elements in the promoters of the 13 NP-like genes in Populus revealed common regulatory elements known to be involved in light regulation, stress/pathogenesis and phytohormone responses. CONCLUSION: In Populus, the evolution of the NP-like protein and gene family has been shaped by duplication events and natural selection. Expression data suggest that previously uncharacterized NP-like proteins may function in nutrient sensing and/or signaling. These proteins are members of Group I NP-like proteins, which are widely distributed in many plant taxa. We conclude that NP-like proteins may function in plants, although this function is undefined.

Baseline survey of the anatomical microbial ecology of an important food plant: Solanum lycopersicum (tomato).

BACKGROUND: Research to understand and control microbiological risks associated with the consumption of fresh fruits and vegetables has examined many environments in the farm to fork continuum. An important data gap however, that remains poorly studied is the baseline description of microflora that may be associated with plant anatomy either endemically or in response to environmental pressures. Specific anatomical niches of plants may contribute to persistence of human pathogens in agricultural environments in ways we have yet to describe. Tomatoes have been implicated in outbreaks of Salmonella at least 17 times during the years spanning 1990 to 2010. Our research seeks to provide a baseline description of the tomato microbiome and possibly identify whether or not there is something distinctive about tomatoes or their growing ecology that contributes to persistence of Salmonella in this important food crop. RESULTS: DNA was recovered from washes of epiphytic surfaces of tomato anatomical organs; leaves, stems, roots, flowers and fruits of Solanum lycopersicum (BHN602), grown at a site in close proximity to commercial farms previously implicated in tomato-Salmonella outbreaks. DNA was amplified for targeted 16S and 18S rRNA genes and sheared for shotgun metagenomic sequencing. Amplicons and metagenomes were used to describe "native" bacterial microflora for diverse anatomical parts of Virginia-grown tomatoes. CONCLUSIONS: Distinct groupings of microbial communities were associated with different tomato plant organs and a gradient of compositional similarity could be correlated to the distance of a given plant part from the soil. Unique bacterial phylotypes (at 95% identity) were associated with fruits and flowers of tomato plants. These include Microvirga, Pseudomonas, Sphingomonas, Brachybacterium, Rhizobiales, Paracocccus, Chryseomonas and Microbacterium. The most frequently observed bacterial taxa across aerial plant regions were Pseudomonas and Xanthomonas. Dominant fungal taxa that could be identified to genus with 18S amplicons included Hypocrea, Aureobasidium and Cryptococcus. No definitive presence of Salmonella could be confirmed in any of the plant samples, although 16S sequences suggested that closely related genera were present on leaves, fruits and roots.