Low occurrence of multi-antimicrobial and heavy metal resistance in Salmonella enterica from wild birds in the United States.

Wild birds are common reservoirs of Salmonella enterica. Wild birds carrying resistant S. enterica may pose a risk to public health as they can spread the resistant bacteria across large spatial scales within a short time. Here, we whole-genome sequenced 375 S. enterica strains from wild birds collected in 41 U.S. states during 1978-2019 to examine bacterial resistance to antibiotics and heavy metals. We found that Typhimurium was the dominant S. enterica serovar, accounting for 68.3% (256/375) of the bird isolates. Furthermore, the proportions of the isolates identified as multi-antimicrobial resistant (multi-AMR: resistant to at least three antimicrobial classes) or multi-heavy metal resistant (multi-HMR: resistant to at least three heavy metals) were both 1.87% (7/375). Interestingly, all the multi-resistant S. enterica (n = 12) were isolated from water birds or raptors; none of them was isolated from songbirds. Plasmid profiling demonstrated that 75% (9/12) of the multi-resistant strains carried resistance plasmids. Our study indicates that wild birds do not serve as important reservoirs of multi-resistant S. enterica strains. Nonetheless, continuous surveillance for bacterial resistance in wild birds is necessary because the multi-resistant isolates identified in this study also showed close genetic relatedness with those from humans and domestic animals.

Evaluating the accuracy of Listeria monocytogenes assemblies from quasimetagenomic samples using long and short reads.

BACKGROUND: Whole genome sequencing of cultured pathogens is the state of the art public health response for the bioinformatic source tracking of illness outbreaks. Quasimetagenomics can substantially reduce the amount of culturing needed before a high quality genome can be recovered. Highly accurate short read data is analyzed for single nucleotide polymorphisms and multi-locus sequence types to differentiate strains but cannot span many genomic repeats, resulting in highly fragmented assemblies. Long reads can span repeats, resulting in much more contiguous assemblies, but have lower accuracy than short reads. RESULTS: We evaluated the accuracy of Listeria monocytogenes assemblies from enrichments (quasimetagenomes) of naturally-contaminated ice cream using long read (Oxford Nanopore) and short read (Illumina) sequencing data. Accuracy of ten assembly approaches, over a range of sequencing depths, was evaluated by comparing sequence similarity of genes in assemblies to a complete reference genome. Long read assemblies reconstructed a circularized genome as well as a 71 kbp plasmid after 24 h of enrichment; however, high error rates prevented high fidelity gene assembly, even at 150X depth of coverage. Short read assemblies accurately reconstructed the core genes after 28 h of enrichment but produced highly fragmented genomes. Hybrid approaches demonstrated promising results but had biases based upon the initial assembly strategy. Short read assemblies scaffolded with long reads accurately assembled the core genes after just 24 h of enrichment, but were highly fragmented. Long read assemblies polished with short reads reconstructed a circularized genome and plasmid and assembled all the genes after 24 h enrichment but with less fidelity for the core genes than the short read assemblies. CONCLUSION: The integration of long and short read sequencing of quasimetagenomes expedited the reconstruction of a high quality pathogen genome compared to either platform alone. A new and more complete level of information about genome structure, gene order and mobile elements can be added to the public health response by incorporating long read analyses with the standard short read WGS outbreak response.

Taxonomic and Functional Shifts in the Sprout Spent Irrigation Water Microbiome in Response to Salmonella Contamination of Alfalfa Seeds.

Despite recent advances in Salmonella-sprout research, little is known about the relationship between Salmonella and the sprout microbiome during sprouting. Sprout spent irrigation water (SSIW) provides an informative representation of the total microbiome of this primarily aquaponic crop. This study was designed to characterize the function and taxonomy of the most actively transcribed genes in SSIW from Salmonella enterica serovar Cubana-contaminated alfalfa seeds throughout the sprouting process. Genomic DNA and total RNA from SSIW was collected at regular intervals and sequenced using Illumina MiSeq and NextSeq platforms. Nucleic acid data were annotated using four different pipelines. Both metagenomic and metatranscriptomic analyses revealed a diverse and highly dynamic SSIW microbiome. A "core" SSIW microbiome comprised Klebsiella, Enterobacter, Pantoea, and Cronobacter The impact, however, of Salmonella contamination on alfalfa seeds influenced SSIW microbial community dynamics not only structurally but also functionally. Changes in genes associated with metabolism, genetic information processing, environmental information processing, and cellular processes were abundant and time dependent. At time points of 24 h, 48 h, and 96 h, totals of 541, 723, and 424 S Cubana genes, respectively, were transcribed at either higher or lower levels than at 0 h in SSIW during sprouting. An array of S Cubana genes (107) were induced at all three time points, including genes involved in biofilm formation and modulation, stress responses, and virulence and tolerance to antimicrobials. Taken together, these findings expand our understanding of the effect of Salmonella seed contamination on the sprout crop microbiome and metabolome.IMPORTANCE Interactions of human enteric pathogens like Salmonella with plants and plant microbiomes remain to be elucidated. The rapid development of next-generation sequencing technologies provides powerful tools enabling investigation of such interactions from broader and deeper perspectives. Using metagenomic and metatranscriptomic approaches, this study identified not only changes in microbiome structure of SSIW associated with sprouting but also changes in the gene expression patterns related to the sprouting process in response to Salmonella contamination of alfalfa seeds. This study advances our knowledge on Salmonella-plant (i.e., sprout) interaction.

HPLC-UV, Metabarcoding and Genome Skims of Botanical Dietary Supplements: A Case Study in Echinacea.

The use of DNA-based methods to authenticate botanical dietary supplements has been vigorously debated for a variety of reasons. More comparisons of DNA-based and chemical methods are needed, and concordant evaluation of orthogonal approaches on the same products will provide data to better understand the strengths and weaknesses of both approaches. The overall application of DNA-based methods is already firmly integrated into a wide array of continually modernizing stand alone and complementary authentication protocols. Recently, the use of full-length chloroplast genome sequences provided enhanced discriminatory capacity for closely related species of Echinacea compared to traditional DNA barcoding approaches (matK and rbcL). Here, two next-generation sequencing approaches were used: (1) genome skimming and (2) PCR amplicon (metabarcoding). The two genetic approaches were then combined with HPLC-UV to evaluate 20 commercially available dietary supplements of Echinacea representing "finished" products. The trade-offs involved in different DNA approaches were discussed, with a focus on how DNA methods support existing, accepted chemical methods. In most of the products (19/20), HPLC-UV suggested the presence of Echinacea spp. While metabarcoding was not useful with this genus and instead only resolved 7 products to the family level, genome skimming was able to resolve to species (9) or genus (1) with the 10/20 products where it was successful. Additional ingredients that HPLC-UV was unable to identify were also found in four products along with the relative sequence proportion of the constituents. Additionally, genome skimming was able to identify one product that was a different Echinacea species entirely.

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.

Using full chloroplast genomes of 'red' and 'yellow' Bixa orellana (achiote) for kmer based identification and phylogenetic inference.

BACKGROUND: Full chloroplast genomes provide high resolution taxonomic discrimination between closely related plant species and are quickly replacing single and multi-locus barcoding regions as reference materials of choice for DNA based taxonomic annotation of plants. Bixa orellana, commonly known as "achiote" and "annatto" is a plant used for both human and animal foods and was thus identified for full chloroplast sequencing for the Center for Veterinary Medicine (CVM) Complete Chloroplast Animal Feed database. This work was conducted in collaboration with the Instituto de Medicina Tradicional (IMET) in Iquitos, Peru. There is a wide range of color variation in pods of Bixa orellana for which genetic loci that distinguish phenotypes have not yet been identified. Here we apply whole chloroplast genome sequencing of "red" and "yellow" individuals of Bixa orellana to provide high quality reference genomes to support kmer database development for use identifying this plant from complex mixtures using shotgun data. Additionally, we describe chloroplast gene content, synteny and phylogeny, and identify an indel and snp that may be associated with seed pod color. RESULTS: Fully assembled chloroplast genomes were produced for both red and yellow Bixa orellana accessions (158,918 and 158,823 bp respectively). Synteny and gene content was identical to the only other previously reported full chloroplast genome of Bixa orellana (NC_041550). We observed a 17 base pair deletion at position 58,399-58,415 in both accessions, relative to NC_041550 and a 6 bp deletion at position 75,531-75,526 and a snp at position 86,493 in red Bixa orellana. CONCLUSIONS: Our data provide high quality reference genomes of individuals of red and yellow Bixa orellana to support kmer based identity markers for use with shotgun sequencing approaches for rapid, precise identification of Bixa orellana from complex mixtures. Kmer based phylogeny of full chloroplast genomes supports monophylly of Bixaceae consistent with alignment based approaches. A potentially discriminatory indel and snp were identified that may be correlated with the red phenotype.

Quasimetagenomic source tracking of Listeria monocytogenes from naturally contaminated ice cream.

BACKGROUND: The more quickly bacterial pathogens responsible for foodborne illness outbreaks can be linked to a vehicle of transmission or a source, the more illnesses can be prevented. Whole genome sequencing (WGS) based approaches to source tracking have greatly increased the speed and resolution with which public health response can pinpoint the vehicle and source of outbreaks. Traditionally, WGS approaches have focused on the culture of an individual isolate before proceeding to DNA extraction and sequencing. For Listeria monocytogenes (Lm), generation of an individual isolate for sequencing typically takes about 6 days. Here we demonstrate that a hybrid, "quasimetagenomic" approach ie; direct sequencing of microbiological enrichments (first step in pathogen detection and recovery) can provide high resolution source tracking sequence data, 5 days earlier than response that focuses on culture and sequencing of an individual isolate. This expedited approach could save lives, prevent illnesses and potentially minimize unnecessary destruction of food. METHODS: Naturally contaminated ice cream (from a 2015 outbreak) was enriched to recover Listeria monocytogenes following protocols outlined in the Bacteriological Analytic Manual (BAM). DNA from enriching microbiota was extracted and sequenced at incremental time-points during the first 48 h of pre-enrichment using the Illumina MiSeq platform (2 by 250), to evaluate genomic coverage of target pathogen, Listeria monocytogenes. RESULTS: Quasimetagenomic sequence data acquired from hour 20 were sufficient to discern whether or not Lm strain/s were part of the ongoing outbreak or not. Genomic data from hours 24, 28, 32, 36, 40, 44, and 48 of pre-enrichments all provided identical phylogenetic source tracking utility to the WGS of individual isolates (which require an additional 5 days to culture). CONCLUSIONS: The speed of this approach (more than twice as fast as current methods) has the potential to reduce the number of illnesses associated with any given outbreak by as many as 75% percent of total cases and potentially with continued optimization of the entire chain of response, contribute to minimized food waste.

Global transcriptomic analyses of Salmonella enterica in Iron-depleted and Iron-rich growth conditions.

BACKGROUND: Salmonella enterica possess several iron acquisition systems, encoded on the chromosome and plasmids. Recently, we demonstrated that incompatibility group (Inc) FIB plasmid-encoded iron acquisition systems (Sit and aerobactin) likely play an important role in persistence of Salmonella in human intestinal epithelial cells (Caco-2). In this study, we sought to determine global transcriptome analyses of S. enterica in iron-rich (IR) and iron-depleted (ID) growth conditions. RESULTS: The number of differentially-expressed genes were substantially higher for recipient (SE819) (n = 966) and transconjugant (TC) (n = 945) compared to the wild type (WT) (SE163A) (n = 110) strain in ID as compared to IR growth conditions. Several virulence-associated factors including T3SS, flagellin, cold-shock protein (cspE), and regulatory genes were upregulated in TC in ID compared to IR conditions. Whereas, IS1 and acrR/tetR transposases located on the IncFIB plasmid, ferritin and several regulatory genes were downregulated in TC in ID conditions. Enterobactin transporter (entS), iron ABC transporter (fepCD), colicin transporter, IncFIB-encoded enolase, cyclic di-GMP regulator (cdgR) and other regulatory genes of the WT strain were upregulated in ID compared to IR conditions. Conversely, ferritin, ferrous iron transport protein A (feoA), IncFIB-encoded IS1 and acrR/tetR transposases and ArtA toxin of WT were downregulated in ID conditions. SDS-PAGE coupled with LC-MS/MS analyses revealed that siderophore receptor proteins such as chromosomally-encoded IroN and, IncFIB-encoded IutA were upregulated in WT and TC in ID growth conditions. Both chromosome and IncFIB plasmid-encoded SitA was overexpressed in WT, but not in TC or recipient in ID conditions. Increased expression of flagellin was detected in recipient and TC, but not in WT in ID conditions. CONCLUSION: Iron concentrations in growth media influenced differential gene expressions both at transcriptional and translational levels, including genes encoded on the IncFIB plasmid. Limited iron availability within the host may promote pathogenic Salmonella to differentially express subsets of genes encoded by chromosome and/or plasmids, facilitating establishment of successful infection.

Metagenome tracking biogeographic agroecology: Phytobiota of tomatoes from Virginia, Maryland, North Carolina and California.

Describing baseline microbiota associated with agricultural commodities in the field is an important step towards improving our understanding of a wide range of important objectives from plant pathology and horticultural sustainability, to food safety. Environmental pressures on plants (wind, dust, drought, water, temperature) vary by geography and characterizing the impact of these variable pressures on phyllosphere microbiota will contribute to improved stewardship of fresh produce for both plant and human health. A higher resolution understanding of the incidence of human pathogens on food plants and co-occurring phytobiota using metagenomic approaches (metagenome tracking) may contribute to improved source attribution and risk assessment in cases where human pathogens become introduced to agro-ecologies. Between 1990 and 2007, as many as 1990 culture-confirmed Salmonella illnesses were linked to tomatoes from as many as 12 multistate outbreaks (Bell et al., 2012; Bell et al., 2015; Bennett et al., 2014; CDC, 2004; CDC, 2007; Greene et al., 2005a; Gruszynski et al., 2014). When possible, source attribution for these incidents revealed a biogeographic trend, most events were associated with eastern growing regions. To improve our understanding of potential biogeographically linked trends in contamination of tomatoes by Salmonella, we profiled microbiota from the surfaces of tomatoes from Virginia, Maryland, North Carolina and California. Bacterial profiles from California tomatoes were completely different than those of Maryland, Virginia and North Carolina (which were highly similar to each other). A statistically significant enrichment of Firmicutes taxa was observed in California phytobiota compared to the three eastern states. Rhizobiaceae, Sphingobacteriaceae and Xanthobacteraceae were the most abundant bacterial families associated with tomatoes grown in eastern states. These baseline metagenomic profiles of phyllosphere microbiota may contribute to improved understanding of how certain ecologies provide supportive resources for human pathogens on plants and how components of certain agro-ecologies may play a role in the introduction of human pathogens to plants.

Metagenomics of pasteurized and unpasteurized gouda cheese using targeted 16S rDNA sequencing.

BACKGROUND: The microbiome of cheese is diverse, even within a variety. The metagenomics of cheese is dependent on a vast array of biotic and abiotic factors. Biotic factors include the population of microbiota and their resulting cellular metabolism. Abiotic factors, including the pH, water activity, fat, salt, and moisture content of the cheese matrix, as well as environmental conditions (temperature, humidity, and location of aging), influence the biotic factors. This study assessed the metagenomics of commercial Gouda cheese prepared using pasteurized or unpasteurized cow milk or pasteurized goat milk via 16S rDNA sequencing. RESULTS: Results were analyzed and compared based on milk pasteurization and source, spatial variability (core, outer, and under the rind), and length of aging (2-4 up to 12-18 months). The dominant organisms in the Gouda cheeses, based on percentage of sequence reads identified at the family or genus levels, were Bacillaceae, Lactococcus, Lactobacillus, Streptococcus, and Staphylococcus. More genus- or family-level (e.g. Bacillaceae) identifications were observed in the Gouda cheeses prepared with unpasteurized cow milk (120) compared with those prepared with pasteurized cow milk (92). When assessing influence of spatial variability on the metagenomics of the cheese, more pronounced differences in bacterial genera were observed in the samples taken under the rind; Brachybacterium, Pseudoalteromonas, Yersinia, Klebsiella, and Weissella were only detected in these samples. Lastly, the aging length of the cheese greatly influenced the number of organisms observed. Twenty-seven additional genus-level identifications were observed in Gouda cheese aged for 12-18 months compared with cheese only aged 2-4 months. CONCLUSIONS: Collectively, the results of this study are important in determining the typical microbiota associated with Gouda cheese and how the microbiome plays a role in safety and quality.

Shifts in spinach microbial communities after chlorine washing and storage at compliant and abusive temperatures.

Fresh produce, like spinach, harbors diverse bacterial populations, including spoilage and potentially pathogenic bacteria. This study examined the effects of produce washing in chlorinated water and subsequent storage on the microbiota of spinach. Baby spinach leaves from a commercial fresh-cut produce processor were assessed before and after washing in chlorinated water, and then after one week's storage at 4, 10, and 15 °C. Microbial communities on spinach were analyzed by non-selective plating, qPCR, and 16S rDNA amplicon sequencing. Bacterial populations on spinach, averaging 6.12 ±â€¯0.61 log CFU/g, were reduced by 1.33 ±â€¯0.57 log after washing. However, populations increased by 1.77-3.24 log after storage, with larger increases occurring at higher temperature (15 > 10 > 4 °C). The predominant phylum identified on unwashed spinach leaves was Proteobacteria; dominant genera were Pseudomonas and Sphingomonas. Bacterial communities shifted significantly after chlorine washing and storage. Several Proteobacteria species, such as Stenotrophomonas sp. and Erwinia sp., were relatively tolerant of chlorine treatment, while species of Flavobacterium and Pedobacter (phylum Bacteroidetes) grew rapidly during storage, especially at abusive temperatures. Cupriavidus sp. and Ralstonia sp. showed significant increases after washing. After storage, microbial communities on spinach appeared to shift back toward the pre-washing distributions.

Development of a Reference Standard Library of Chloroplast Genome Sequences, GenomeTrakrCP.

Precise, species-level identification of plants in foods and dietary supplements is difficult. While the use of DNA barcoding regions (short regions of DNA with diagnostic utility) has been effective for many inquiries, it is not always a robust approach for closely related species, especially in highly processed products. The use of fully sequenced chloroplast genomes, as an alternative to short diagnostic barcoding regions, has demonstrated utility for closely related species. The U. S. Food and Drug Administration (FDA) has also developed species-specific DNA-based assays targeting plant species of interest by utilizing chloroplast genome sequences. Here, we introduce a repository of complete chloroplast genome sequences called GenomeTrakrCP, which will be publicly available at the National Center for Biotechnology Information (NCBI). Target species for inclusion are plants found in foods and dietary supplements, toxin producers, common contaminants and adulterants, and their close relatives. Publicly available data will include annotated assemblies, raw sequencing data, and voucher information with each NCBI accession associated with an authenticated reference herbarium specimen. To date, 40 complete chloroplast genomes have been deposited in GenomeTrakrCP (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA325670/), and this will be expanded in the future.

Tracing Origins of the Salmonella Bareilly Strain Causing a Food-borne Outbreak in the United States.

BACKGROUND: Using a novel combination of whole-genome sequencing (WGS) analysis and geographic metadata, we traced the origins of Salmonella Bareilly isolates collected in 2012 during a widespread food-borne outbreak in the United States associated with scraped tuna imported from India. METHODS: Using next-generation sequencing, we sequenced the complete genome of 100 Salmonella Bareilly isolates obtained from patients who consumed contaminated product, from natural sources, and from unrelated historically and geographically disparate foods. Pathogen genomes were linked to geography by projecting the phylogeny on a virtual globe and produced a transmission network. RESULTS: Phylogenetic analysis of WGS data revealed a common origin for outbreak strains, indicating that patients in Maryland and New York were infected from sources originating at a facility in India. CONCLUSIONS: These data represent the first report fully integrating WGS analysis with geographic mapping and a novel use of transmission networks. Results showed that WGS vastly improves our ability to delimit the scope and source of bacterial food-borne contamination events. Furthermore, these findings reinforce the extraordinary utility that WGS brings to global outbreak investigation as a greatly enhanced approach to protecting the human food supply chain as well as public health in general.

Enrichment dynamics of Listeria monocytogenes and the associated microbiome from naturally contaminated ice cream linked to a listeriosis outbreak.

BACKGROUND: Microbiota that co-enrich during efforts to recover pathogens from foodborne outbreaks interfere with efficient detection and recovery. Here, dynamics of co-enriching microbiota during recovery of Listeria monocytogenes from naturally contaminated ice cream samples linked to an outbreak are described for three different initial enrichment formulations used by the Food and Drug Administration (FDA), the International Organization of Standardization (ISO), and the United States Department of Agriculture (USDA). Enrichment cultures were analyzed using DNA extraction and sequencing from samples taken every 4 h throughout 48 h of enrichment. Resphera Insight and CosmosID analysis tools were employed for high-resolution profiling of 16S rRNA amplicons and whole genome shotgun data, respectively. RESULTS: During enrichment, other bacterial taxa were identified, including Anoxybacillus, Geobacillus, Serratia, Pseudomonas, Erwinia, and Streptococcus spp. Surprisingly, incidence of L. monocytogenes was proportionally greater at hour 0 than when tested 4, 8, and 12 h later with all three enrichment schemes. The corresponding increase in Anoxybacillus and Geobacillus spp.indicated these taxa co-enriched in competition with L. monocytogenes during early enrichment hours. L. monocytogenes became dominant after 24 h in all three enrichments. DNA sequences obtained from shotgun metagenomic data of Listeria monocytogenes at 48 h were assembled to produce a consensus draft genome which appeared to have a similar tracking utility to pure culture isolates of L. monocytogenes. CONCLUSIONS: All three methods performed equally well for enrichment of Listeria monocytogenes. The observation of potential competitive exclusion of L. mono by Anoxybacillus and Geobacillus in early enrichment hours provided novel information that may be used to further optimize enrichment formulations. Application of Resphera Insight for high-resolution analysis of 16S amplicon sequences accurately identified L. monocytogenes. Both shotgun and 16S rRNA data supported the presence of three slightly variable genomes of L. monocytogenes. Moreover, the draft assembly of a consensus genome of L. monocytogenes from shotgun metagenomic data demonstrated the potential utility of this approach to expedite trace-back of outbreak-associated strains, although further validation will be needed to confirm this utility.

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.

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.

Microbiota of laboratory channel catfish skin mucosa and aquaria water exposed to chloramine-T trihydrate.

Here, we describe the skin mucosa microbiome of channel catfish (Ictalurus punctatus) before and after exposure to chloramine-T trihydrate. We also describe the aquaria water microbiome after the post-treatment period. These data provide a unique baseline description of skin mucosa and aquaria water microbiome from catfish reared in research aquaria.

Fecal microbiomes of laboratory beagles receiving antiparasitic formulations in an experimental setting.

Here, we describe the fecal microbiome of laboratory beagles in a non-invasive experiment designed to contrast in vivo versus in vitro bioequivalence in response to antiparasitic drug administration. The experiment provided a unique opportunity to evaluate metagenomic profiles of canine feces before and after anti-parasitic drug exposure.

Paired metagenomic and chemical evaluation of aflatoxin-contaminated dog kibble.

Introduction: Identification of chemical toxins from complex or highly processed foods can present 'needle in the haystack' challenges for chemists. Metagenomic data can be used to guide chemical toxicity evaluations by providing DNA-based description of the wholistic composition (eukaryotic, bacterial, protozoal, viral, and antimicrobial resistance) of foods suspected to harbor toxins, allergens, or pathogens. This type of information can focus chemistry-based diagnostics, improve hazard characterization and risk assessment, and address data gaps. Additionally, there is increasing recognition that simultaneously co-occurring mycotoxins, either from single or multiple species, can impact dietary toxicity exposure. Metagenomic data provides a way to address data gaps related to co-occurrence of multiple fungal species. Methods: Paired metagenomic and chemical data were used to evaluate aflatoxin-contaminated kibble with known levels of specific mycotoxins. Kibble was ground to a fine powder for both chemical and molecular analyses. Chemical analyses were performed with Liquid Chromatography Mass Spectrometry (LCMS) and according to the AOAC Official method 2005.08: Aflatoxins in Corn, Raw Peanuts, and Peanut Butter using Liquid Chromatography with Post-Column Photochemical Derivatization. Metagenomes were created from DNA extracted from ground kibble and sequenced on an Illumina NextSeq 2000 with an average sequence depth of 180 million reads per replicate. Results and discussion: Metagenomic data demonstrated that the abundance of DNA from putative aflatoxigenic Aspergillus spp. correlated with the levels of aflatoxin quantified by LCMS. Metagenomic data also identified an expansive range of co-occurring fungal taxa which may produce additional mycotoxins. DNA data paired with chemical data provides a novel modality to address current data gaps surrounding dietary mycotoxin exposure, toxigenic fungal taxonomy, and mycotoxins of emerging concern.

Metagenomic survey of antimicrobial resistance (AMR) in Maryland surface waters differentiated by high and low human impact.

Here, we examine surface waters as a modality to better understand baseline antimicrobial resistance (AMR) across the environment to supplement existing AMR monitoring in pathogens associated with humans, foods, and animals. Data from metagenomic and quasimetagenomic (shotgun sequenced enrichments) are used to describe AMR in Maryland surface waters from high and low human impact classifications.