RESUMEN
For decades, the 16S rRNA gene has been used to taxonomically classify prokaryotic species and to taxonomically profile microbial communities. However, the 16S rRNA gene has been criticized for being too conserved to differentiate between distinct species. We argue that the inability to differentiate between species is not a unique feature of the 16S rRNA gene. Rather, we observe the gradual loss of species-level resolution for other nearly-universal prokaryotic marker genes as the number of gene sequences increases in reference databases. This trend was strongly correlated with how represented a taxonomic group was in the database and indicates that, at the gene-level, the boundaries between many species might be fuzzy. Through our study, we argue that any approach that relies on a single marker to distinguish bacterial taxa is fraught even if some markers appear to be discriminative in current databases.
Asunto(s)
Bacterias , Bases de Datos Genéticas , ARN Ribosómico 16S , ARN Ribosómico 16S/genética , Bacterias/genética , Bacterias/clasificación , Marcadores Genéticos/genética , Filogenia , Biología Computacional/métodosRESUMEN
BACKGROUND: Oxford Nanopore provides high throughput sequencing platforms able to reconstruct complete bacterial genomes with 99.95% accuracy. However, even small levels of error can obscure the phylogenetic relationships between closely related isolates. Polishing tools have been developed to correct these errors, but it is uncertain if they obtain the accuracy needed for the high-resolution source tracking of foodborne illness outbreaks. RESULTS: We tested 132 combinations of assembly and short- and long-read polishing tools to assess their accuracy for reconstructing the genome sequences of 15 highly similar Salmonella enterica serovar Newport isolates from a 2020 onion outbreak. While long-read polishing alone improved accuracy, near perfect accuracy (99.9999% accuracy or ~ 5 nucleotide errors across the 4.8 Mbp genome, excluding low confidence regions) was only obtained by pipelines that combined both long- and short-read polishing tools. Notably, medaka was a more accurate and efficient long-read polisher than Racon. Among short-read polishers, NextPolish showed the highest accuracy, but Pilon, Polypolish, and POLCA performed similarly. Among the 5 best performing pipelines, polishing with medaka followed by NextPolish was the most common combination. Importantly, the order of polishing tools mattered i.e., using less accurate tools after more accurate ones introduced errors. Indels in homopolymers and repetitive regions, where the short reads could not be uniquely mapped, remained the most challenging errors to correct. CONCLUSIONS: Short reads are still needed to correct errors in nanopore sequenced assemblies to obtain the accuracy required for source tracking investigations. Our granular assessment of the performance of the polishing pipelines allowed us to suggest best practices for tool users and areas for improvement for tool developers.
Asunto(s)
Benchmarking , Brotes de Enfermedades , Genoma Bacteriano , Nanoporos , Secuenciación de Nanoporos/métodos , Secuenciación de Nucleótidos de Alto Rendimiento/métodos , Salmonella enterica/genética , Salmonella enterica/aislamiento & purificación , Humanos , FilogeniaRESUMEN
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.
Asunto(s)
Salmonella enterica , Serogrupo , Cebollas/genética , Granjas , Filogenia , Plásmidos/genética , Brotes de EnfermedadesRESUMEN
MOTIVATION: Microbial gene catalogs are data structures that organize genes found in microbial communities, providing a reference for standardized analysis of the microbes across samples and studies. Although gene catalogs are commonly used, they have not been critically evaluated for their effectiveness as a basis for metagenomic analyses. RESULTS: As a case study, we investigate one such catalog, the Integrated Gene Catalog (IGC), however, our observations apply broadly to most gene catalogs constructed to date. We focus on both the approach used to construct this catalog and on its effectiveness when used as a reference for microbiome studies. Our results highlight important limitations of the approach used to construct the IGC and call into question the broad usefulness of gene catalogs more generally. We also recommend best practices for the construction and use of gene catalogs in microbiome studies and highlight opportunities for future research. AVAILABILITY AND IMPLEMENTATION: All supporting scripts for our analyses can be found on GitHub: https://github.com/SethCommichaux/IGC.git. The supporting data can be downloaded from: https://obj.umiacs.umd.edu/igc-analysis/IGC_analysis_data.tar.gz. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Asunto(s)
Metagenoma , Microbiota , Microbiota/genética , MetagenómicaRESUMEN
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.
Asunto(s)
Listeria monocytogenes , Nanoporos , Genómica , Secuenciación de Nucleótidos de Alto Rendimiento , Listeria monocytogenes/genética , Análisis de Secuencia de ADN , Secuenciación Completa del GenomaRESUMEN
Seventeen Salmonella enterica serovar Schwarzengrund isolates from chicken (n = 9) and clinical samples including stool (n = 6), urine (n = 1), and gallbladder (n = 1) were sequenced and found to carry an IncFIB-IncFIC (FII) fusion plasmid of approximately 145 Kb. This information provides reference genomic data for comparative studies of S. Schwarzengrund pathogenicity and plasmid genetics.
RESUMEN
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.
RESUMEN
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.
RESUMEN
A total of 55 food and clinical S. Schwarzengrund isolates were assayed for plasmid content, among which an IncFIB-IncFIC(FII) fusion plasmid, conferring streptomycin resistance, was detected in 17 isolates. Among the 17 isolates, 9 were food isolates primarily collected from poultry meat, and 8 clinical isolates collected from stool, urine, and gallbladder. SNP-based phylogenetic analyses showed that the isolates carrying the fusion plasmid formed a subclade indicating the plasmid was acquired and is now maintained by the lineage. Phylogenetic analysis of the plasmid suggested it is derived from avian pathogenic plasmids and might confer an adaptive advantage to the S. Schwarzengrund isolates within birds. IncFIB-IncFIC(FII) fusion plasmids from all food and three clinical isolates were self-conjugative and successfully transferred into E. coli J53 by conjugation. Food and clinical isolates had similar virulome profiles and were able to invade human Caco-2 cells. However, the IncFIB-IncFIC(FII) plasmid did not significantly add to their invasion and persistence potential in human Caco-2 cells.
RESUMEN
For decades, the 16S rRNA gene has been used to taxonomically classify prokaryotic species and to taxonomically profile microbial communities. The 16S rRNA gene has been criticized for being too conserved to differentiate between distinct species. We argue that the inability to differentiate between species is not a unique feature of the 16S rRNA gene. Rather, we observe the gradual loss of species-level resolution for other marker genes as the number of gene sequences increases in reference databases. We demonstrate this effect through the analysis of three commonly used databases of nearly-universal prokaryotic marker genes: the SILVA 16S rRNA gene database, the Genome Taxonomy Database (GTDB), and a set of 40 taxonomically-informative single-copy genes. Our results reflect a more fundamental property of the taxonomies themselves and have broad implications for bioinformatic analyses beyond taxonomic classification. Effective solutions for fine-level taxonomic classification require a more precise, and operationally-relevant, definition of the taxonomic labels being sought, and the use of combinations of genomic markers in the classification process. Importance: The use of reference databases for assigning taxonomic labels to genomic and metagenomic sequences is a fundamental bioinformatic task in the characterization of microbial communities. The increasing accessibility of high throughput sequencing has led to a rapid increase in the size and number of sequences in databases. This has been beneficial for improving our understanding of the global microbial genetic diversity. However, there is evidence that as the microbial diversity is more densely sampled, increasingly longer genomic segments are needed to differentiate between distinct species. The scientific community needs to be aware of this issue and needs to develop methods that better account for it when assigning taxonomic labels to metagenomic sequences from microbial communities.
RESUMEN
Sprouts have been implicated in numerous foodborne illness outbreaks. To better understand baseline microbial profiles of irrigation water and subsequent spent irrigation water of alfalfa sprouts, DNA from water was extracted, sequenced, and annotated with CosmosID and a custom pipeline to provide bacterial, fungal, protist, and antimicrobial resistance gene profiles.
RESUMEN
Water from the Hickey Run Tributary of the Anacostia River is being collected quarterly (beginning August 2018) and analyzed to create high-resolution baseline taxonomic profiles of microbiota associated with this important aquatic ecosystem, which has a long history of exposure to residential and commercial effluents from Washington, DC. These United States National Arboretum Microbial Observatory data are available under NCBI BioProject number PRJNA498951.
RESUMEN
Accurate predictions across multiple fields of microbiome research have far-reaching benefits to society, but there are few widely accepted quantitative tools to make accurate predictions about microbial communities and their functions. More discussion is needed about the current state of microbiome analysis and the tools required to overcome the hurdles preventing development and implementation of predictive analyses. We summarize the ideas generated by participants of the Mid-Atlantic Microbiome Meet-up in January 2019. While it was clear from the presentations that most fields have advanced beyond simple associative and descriptive analyses, most fields lack essential elements needed for the development and application of accurate microbiome predictions. Participants stressed the need for standardization, reproducibility, and accessibility of quantitative tools as key to advancing predictions in microbiome analysis. We highlight hurdles that participants identified and propose directions for future efforts that will advance the use of prediction in microbiome research.
RESUMEN
Aquaculture plays an increasingly important role in the growing demand for seafood. Hatchery production of oyster larvae is an integral component of oyster farming, providing single seed for off-bottom farming or larvae for setting on larger substrates for on-bottom farming. Larvae from certain tanks in an established aquaculture enterprise were dying from an unknown etiologic agent. A metagenomic approach was used to examine oyster larvae and water from larval tanks with high and low survival rates to evaluate the epidemiological efficacy of this approach.