1
|
Ulrich JU, Renard BY. Fast and space-efficient taxonomic classification of long reads with hierarchical interleaved XOR filters. Genome Res 2024; 34:914-924. [PMID: 38886068 PMCID: PMC11293544 DOI: 10.1101/gr.278623.123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2023] [Accepted: 05/23/2024] [Indexed: 06/20/2024]
Abstract
Metagenomic long-read sequencing is gaining popularity for various applications, including pathogen detection and microbiome studies. To analyze the large data created in those studies, software tools need to taxonomically classify the sequenced molecules and estimate the relative abundances of organisms in the sequenced sample. Because of the exponential growth of reference genome databases, the current taxonomic classification methods have large computational requirements. This issue motivated us to develop a new data structure for fast and memory-efficient querying of long reads. Here, we present Taxor as a new tool for long-read metagenomic classification using a hierarchical interleaved XOR filter data structure for indexing and querying large reference genome sets. Taxor implements several k-mer-based approaches, such as syncmers, for pseudoalignment to classify reads and an expectation-maximization algorithm for metagenomic profiling. Our results show that Taxor outperforms state-of-the-art tools regarding precision while having a similar recall for long-read taxonomic classification. Most notably, Taxor reduces the memory requirements and index size by >50% and is among the fastest tools regarding query times. This enables real-time metagenomics analysis with large reference databases on a small laptop in the field.
Collapse
Affiliation(s)
- Jens-Uwe Ulrich
- Data Analytics and Computational Statistics, Hasso Plattner Institute, Digital Engineering Faculty, University of Potsdam, 14482 Potsdam, Germany;
- Phylogenomics Unit, Center for Artificial Intelligence in Public Health Research, Robert Koch Institute, 15745 Wildau, Germany
- Department of Mathematics and Computer Science, Free University of Berlin, 14195 Berlin, Germany
| | - Bernhard Y Renard
- Data Analytics and Computational Statistics, Hasso Plattner Institute, Digital Engineering Faculty, University of Potsdam, 14482 Potsdam, Germany;
| |
Collapse
|
2
|
Shi A, Nguyen H, Kuo CB, Beringer PM. Drug-induced liver injury associated with elexacaftor/tezacaftor/ivacaftor: A pharmacovigilance analysis of the FDA adverse event reporting system (FAERS). J Cyst Fibros 2024; 23:566-572. [PMID: 38233246 DOI: 10.1016/j.jcf.2024.01.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2023] [Revised: 11/28/2023] [Accepted: 01/04/2024] [Indexed: 01/19/2024]
Abstract
INTRODUCTION The efficacy and safety of elexacaftor/tezacaftor/ivacaftor (ETI) have been established in prospective clinical trials. Liver function test elevations were observed in a greater proportion of patients receiving ETI compared with placebo; however, the relatively small number of patients and short duration of study preclude detection of rare but clinically significant associations with drug-induced liver injury (DILI). To address this gap, we assessed the real-world risk of DILI associated with ETI through data mining of the FDA adverse event reporting system (FAERS). METHODS Disproportionality analyses were conducted on FAERS data from the fourth quarter of 2019 through the third quarter of 2022. Comparative patient demographics, onset time and outcomes for ETI-DILI were also obtained. RESULTS 452 reports of DILI associated with ETI were found, representing 2.1 % of all adverse event reports for ETI. All disproportionality measures were significant for ETI-DILI at p < 0.05; the reporting odds ratio (ROR) (2.82) was comparable to that of drugs classified by FDA as "Most-DILI concern". The most notable demographic finding was a male majority (5:4 male to female ratio) for ETI-DILI compared to a female majority (4:5 male to female ratio) for non ETI-DILI. Median ETI-DILI onset time was 50.5 days, and hospitalization was the second most common complication. CONCLUSION Using FAERS data, ETI was found to be disproportionately associated with DILI. Future research is needed to investigate the hepatotoxic mechanisms and assess potential mitigation strategies for ETI-induced hepatotoxicity.
Collapse
Affiliation(s)
- Alan Shi
- Department of Clinical Pharmacy, Mann School of Pharmacy and Pharmaceutical Sciences, University of Southern California, 1985 Zonal Ave, Los Angeles, CA 90033, USA
| | - Harold Nguyen
- Department of Clinical Pharmacy, Mann School of Pharmacy and Pharmaceutical Sciences, University of Southern California, 1985 Zonal Ave, Los Angeles, CA 90033, USA
| | - C Benson Kuo
- Department of Regulatory and Quality Sciences, Mann School of Pharmacy and Pharmaceutical Sciences, University of Southern California, 1985 Zonal Ave, Los Angeles, CA 90033, USA
| | - Paul M Beringer
- Department of Clinical Pharmacy, Mann School of Pharmacy and Pharmaceutical Sciences, University of Southern California, 1985 Zonal Ave, Los Angeles, CA 90033, USA; USC Anton Yelchin CF Clinic, 1510 San Pablo St, Los Angeles, CA 90033, USA.
| |
Collapse
|
3
|
Cadenas-Castrejón E, Verleyen J, Boukadida C, Díaz-González L, Taboada B. Evaluation of tools for taxonomic classification of viruses. Brief Funct Genomics 2023; 22:31-41. [PMID: 36335985 DOI: 10.1093/bfgp/elac036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2022] [Revised: 09/21/2022] [Accepted: 09/28/2022] [Indexed: 11/09/2022] Open
Abstract
Viruses are the most abundant infectious agents on earth, and they infect living organisms such as bacteria, plants and animals, among others. They play an important role in the balance of different ecosystems by modulating microbial populations. In humans, they are responsible for some common diseases and may cause severe illnesses. Viral metagenomic studies have become essential and offer the possibility to understand and extend the knowledge of virus diversity and functionality. For these approaches, an essential step is the classification of viral sequences. In this work, 11 taxonomic classification tools were compared by analysing their performances, in terms of sensitivity and precision, to classify reads at the species and family levels using the same (viral and nonviral) datasets and evaluation metrics, as well as their processing times and memory requirements. The results showed that factors such as richness (numbers of viral species in samples), taxonomic level in the classification and read length influence tool performance. High values of viral richness in samples decreased the performances of most tools. Additionally, the classifications were better at higher taxonomic levels, such as families, compared to lower taxonomic levels, such as species, and were more evident in short reads. The results also indicated that BLAST and Kraken2 were the best tools for classifying all types of reads, while FastViromeExplorer and VirusFinder were only good when used for long reads and Centrifuge, DIAMOND, and One Codex when used for short reads. Regarding nonviral datasets (human and bacterial), all tools correctly classified them as nonviral.
Collapse
|
4
|
Shen W, Xiang H, Huang T, Tang H, Peng M, Cai D, Hu P, Ren H. KMCP: accurate metagenomic profiling of both prokaryotic and viral populations by pseudo-mapping. Bioinformatics 2022; 39:6965021. [PMID: 36579886 PMCID: PMC9828150 DOI: 10.1093/bioinformatics/btac845] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2022] [Revised: 12/17/2022] [Accepted: 12/28/2022] [Indexed: 12/30/2022] Open
Abstract
MOTIVATION The growing number of microbial reference genomes enables the improvement of metagenomic profiling accuracy but also imposes greater requirements on the indexing efficiency, database size and runtime of taxonomic profilers. Additionally, most profilers focus mainly on bacterial, archaeal and fungal populations, while less attention is paid to viral communities. RESULTS We present KMCP (K-mer-based Metagenomic Classification and Profiling), a novel k-mer-based metagenomic profiling tool that utilizes genome coverage information by splitting the reference genomes into chunks and stores k-mers in a modified and optimized Compact Bit-Sliced Signature Index for fast alignment-free sequence searching. KMCP combines k-mer similarity and genome coverage information to reduce the false positive rate of k-mer-based taxonomic classification and profiling methods. Benchmarking results based on simulated and real data demonstrate that KMCP, despite a longer running time than all other methods, not only allows the accurate taxonomic profiling of prokaryotic and viral populations but also provides more confident pathogen detection in clinical samples of low depth. AVAILABILITY AND IMPLEMENTATION The software is open-source under the MIT license and available at https://github.com/shenwei356/kmcp. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Wei Shen
- To whom correspondence should be addressed. or or or
| | - Hongyan Xiang
- Key Laboratory of Molecular Biology for Infectious Diseases (Ministry of Education), Department of Infectious Diseases, Institute for Viral Hepatitis, The Second Affiliated Hospital, Chongqing Medical University, Chongqing 400010, China
| | - Tianquan Huang
- Key Laboratory of Molecular Biology for Infectious Diseases (Ministry of Education), Department of Infectious Diseases, Institute for Viral Hepatitis, The Second Affiliated Hospital, Chongqing Medical University, Chongqing 400010, China
| | - Hui Tang
- Key Laboratory of Molecular Biology for Infectious Diseases (Ministry of Education), Department of Infectious Diseases, Institute for Viral Hepatitis, The Second Affiliated Hospital, Chongqing Medical University, Chongqing 400010, China
| | - Mingli Peng
- Key Laboratory of Molecular Biology for Infectious Diseases (Ministry of Education), Department of Infectious Diseases, Institute for Viral Hepatitis, The Second Affiliated Hospital, Chongqing Medical University, Chongqing 400010, China
| | - Dachuan Cai
- Key Laboratory of Molecular Biology for Infectious Diseases (Ministry of Education), Department of Infectious Diseases, Institute for Viral Hepatitis, The Second Affiliated Hospital, Chongqing Medical University, Chongqing 400010, China
| | - Peng Hu
- To whom correspondence should be addressed. or or or
| | - Hong Ren
- To whom correspondence should be addressed. or or or
| |
Collapse
|
5
|
Abstract
Assigning taxonomy remains a challenging topic in microbiome studies, due largely to ambiguity of reads which overlap multiple reference genomes. With the Web of Life (WoL) reference database hosting 10,575 reference genomes and growing, the percentage of ambiguous reads will only increase. The resulting artifacts create both the illusion of co-occurrence and a long tail end of extraneous reference hits that confound interpretation. We introduce genome cover, the fraction of reference genome overlapped by reads, to distinguish these artifacts. We show how to dynamically predict genome cover by read count and examine our model in Staphylococcus aureus monoculture. Our modeling cleanly separates both S. aureus and true contaminants from the false artifacts of reference overlap. We next introduce saturated genome cover, the true fraction of a reference genome overlapped by sample contents. Genome cover may not saturate for low abundance or low prevalence bacteria. We assuage this worry with examination of a large human fecal data set. By compositing the metric across like samples, genome cover saturates even for rare species. We note that it is a threshold on saturated genome cover, not genome cover itself, which indicates a spurious reference hit or distant relative. We present Zebra, a method to compute and threshold the genome cover metric across like samples, a recurrence to estimate genome cover and confirm saturation, and provide guidance for choosing cover thresholds in real world scenarios. Standalone genome cover and integration into Woltka are available: https://github.com/biocore/zebra_filter, https://github.com/qiyunzhu/woltka. IMPORTANCE Taxonomic assignment, assigning sequences to specific taxonomic units, is a crucial processing step in microbiome analyses. Issues in taxonomic assignment affect interpretation of what microbes are present in each sample and may be associated with specific environmental or clinical conditions. Assigning importance to a particular taxon relies strongly on independence of assigned counts. The false inclusion of thousands of correlated taxa makes interpretation ambiguous, leading to underconstrained results which cannot be reproduced. The importance sometimes attached to implausible artifacts such as anthrax or bubonic plague is especially problematic. We show that the Zebra filter retrieves only the nearest relatives of sample contents enabling more reproducible and biologically plausible interpretation of metagenomic data.
Collapse
|
6
|
PathoLive—Real-Time Pathogen Identification from Metagenomic Illumina Datasets. Life (Basel) 2022; 12:life12091345. [PMID: 36143382 PMCID: PMC9505849 DOI: 10.3390/life12091345] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2022] [Revised: 08/24/2022] [Accepted: 08/24/2022] [Indexed: 11/18/2022] Open
Abstract
Over the past years, NGS has become a crucial workhorse for open-view pathogen diagnostics. Yet, long turnaround times result from using massively parallel high-throughput technologies as the analysis can only be performed after sequencing has finished. The interpretation of results can further be challenged by contaminations, clinically irrelevant sequences, and the sheer amount and complexity of the data. We implemented PathoLive, a real-time diagnostics pipeline for the detection of pathogens from clinical samples hours before sequencing has finished. Based on real-time alignment with HiLive2, mappings are scored with respect to common contaminations, low-entropy areas, and sequences of widespread, non-pathogenic organisms. The results are visualized using an interactive taxonomic tree that provides an easily interpretable overview of the relevance of hits. For a human plasma sample that was spiked in vitro with six pathogenic viruses, all agents were clearly detected after only 40 of 200 sequencing cycles. For a real-world sample from Sudan, the results correctly indicated the presence of Crimean-Congo hemorrhagic fever virus. In a second real-world dataset from the 2019 SARS-CoV-2 outbreak in Wuhan, we found the presence of a SARS coronavirus as the most relevant hit without the novel virus reference genome being included in the database. For all samples, clinically irrelevant hits were correctly de-emphasized. Our approach is valuable to obtain fast and accurate NGS-based pathogen identifications and correctly prioritize and visualize them based on their clinical significance: PathoLive is open source and available on GitLab and BioConda.
Collapse
|
7
|
Gut Microbial Signatures of Distinct Trimethylamine N-Oxide Response to Raspberry Consumption. Nutrients 2022; 14:nu14081656. [PMID: 35458219 PMCID: PMC9027468 DOI: 10.3390/nu14081656] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2022] [Revised: 04/11/2022] [Accepted: 04/13/2022] [Indexed: 02/07/2023] Open
Abstract
The aim of this exploratory study was to evaluate the gut microbial signatures of distinct trimethylamine N-oxide (TMAO) responses following raspberry consumption. Investigations were carried out in 24 subjects at risk of developing metabolic syndrome who received 280 g/day of frozen raspberries for 8 weeks. Blood and stool samples were collected at weeks 0 and 8. Inter-individual variability in plasma TMAO levels was analyzed, 7 subjects were excluded due to noninformative signals and 17 subjects were kept for analysis and further stratified according to their TMAO response. Whole-metagenome shotgun sequencing analysis was used to determine the impact of raspberry consumption on gut microbial composition. Before the intervention, the relative abundance of Actinobacteriota was significantly higher in participants whose TMAO levels increased after the intervention (p = 0.03). The delta TMAO (absolute differences of baseline and week 8 levels) was positively associated with the abundance of gut bacteria such as Bilophila wadsworthia (p = 0.02; r2 = 0.37), from the genus Granulicatella (p = 0.03; r2 = 0.48) or the Erysipelotrichia class (p = 0.03; r2 = 0.45). Changes in the gut microbial ecology induced by raspberry consumption over an 8-week period presumably impacted quaternary amines-utilizing activity and thus plasma TMAO levels.
Collapse
|
8
|
Franck M, de Toro-Martín J, Varin TV, Garneau V, Pilon G, Roy D, Couture P, Couillard C, Marette A, Vohl MC. Raspberry consumption: identification of distinct immune-metabolic response profiles by whole blood transcriptome profiling. J Nutr Biochem 2022; 101:108946. [PMID: 35016998 DOI: 10.1016/j.jnutbio.2022.108946] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2021] [Revised: 11/02/2021] [Accepted: 12/22/2021] [Indexed: 01/03/2023]
Abstract
Numerous studies have reported that diets rich in phenolic compounds are beneficial to immune-metabolic health, yet these effects are heterogeneous and the underlying mechanisms are poorly understood. To investigate the inter-individual variability of the immune-metabolic response to raspberry consumption, whole-blood RNAseq data from 24 participants receiving 280g/day of raspberries for 8 weeks were used for the identification of responsiveness subgroups by using partial least squares-discriminant analysis (PLSDA) and hierarchical clustering. Transcriptomic-based clustering regrouped participants into two distinct subgroups of 13 and 11 participants, so-called responders and non-responders, respectively. Following raspberry consumption, a significant decrease in triglycerides, cholesterol and C-reactive protein levels were found in responders, as compared to non-responders. Two major gene expression components of 100 and 220 genes were identified by sparse PLSDA as those better discriminating responders from non-responders, and functional analysis identified pathways related to cytokine production, leukocyte activation and immune response as significantly enriched with most discriminant genes. As compared to non-responders, the plasma lipidomic profile of responders was characterized by a significant decrease in triglycerides and an increase in phosphatidylcholines following raspberry consumption. Prior to the intervention, a distinct metagenomic profile was identified by PLSDA between responsiveness subgroups, and the Firmicutes-to-Bacteroidota ratio was found significantly lower in responders, as compared to non-responders. Findings point to this transcriptomic-based clustering approach as a suitable tool to identify distinct responsiveness subgroups to raspberry consumption. This approach represents a promising framework to tackle the issue of inter-individual variability in the understanding of the impact of foods on immune-metabolic health.
Collapse
Affiliation(s)
- Maximilien Franck
- Institute of Nutrition and Functional Foods (INAF), Université Laval, Québec, QC G1V 0A6, Canada; Centre Nutrition, santé et société (NUTRISS), Université Laval, Québec, QC G1V 0A6, Canada; School of Nutrition, Université Laval, Québec, QC G1V 0A6, Canada
| | - Juan de Toro-Martín
- Institute of Nutrition and Functional Foods (INAF), Université Laval, Québec, QC G1V 0A6, Canada; Centre Nutrition, santé et société (NUTRISS), Université Laval, Québec, QC G1V 0A6, Canada; School of Nutrition, Université Laval, Québec, QC G1V 0A6, Canada
| | - Thibault V Varin
- Institute of Nutrition and Functional Foods (INAF), Université Laval, Québec, QC G1V 0A6, Canada; Centre Nutrition, santé et société (NUTRISS), Université Laval, Québec, QC G1V 0A6, Canada; Quebec Heart and Lung Institute (IUCPQ) Research Center, Québec, QC G1V 4G5, Canada
| | - Véronique Garneau
- Institute of Nutrition and Functional Foods (INAF), Université Laval, Québec, QC G1V 0A6, Canada; Centre Nutrition, santé et société (NUTRISS), Université Laval, Québec, QC G1V 0A6, Canada; School of Nutrition, Université Laval, Québec, QC G1V 0A6, Canada
| | - Geneviève Pilon
- Institute of Nutrition and Functional Foods (INAF), Université Laval, Québec, QC G1V 0A6, Canada; Centre Nutrition, santé et société (NUTRISS), Université Laval, Québec, QC G1V 0A6, Canada; Quebec Heart and Lung Institute (IUCPQ) Research Center, Québec, QC G1V 4G5, Canada
| | - Denis Roy
- Institute of Nutrition and Functional Foods (INAF), Université Laval, Québec, QC G1V 0A6, Canada
| | - Patrick Couture
- Institute of Nutrition and Functional Foods (INAF), Université Laval, Québec, QC G1V 0A6, Canada; Centre Nutrition, santé et société (NUTRISS), Université Laval, Québec, QC G1V 0A6, Canada; Endocrinology and Nephrology Unit, CHU de Quebec Research Center, Québec, QC G1V 4G2, Canada
| | - Charles Couillard
- Institute of Nutrition and Functional Foods (INAF), Université Laval, Québec, QC G1V 0A6, Canada; Centre Nutrition, santé et société (NUTRISS), Université Laval, Québec, QC G1V 0A6, Canada; School of Nutrition, Université Laval, Québec, QC G1V 0A6, Canada
| | - André Marette
- Institute of Nutrition and Functional Foods (INAF), Université Laval, Québec, QC G1V 0A6, Canada; Centre Nutrition, santé et société (NUTRISS), Université Laval, Québec, QC G1V 0A6, Canada; Quebec Heart and Lung Institute (IUCPQ) Research Center, Québec, QC G1V 4G5, Canada
| | - Marie-Claude Vohl
- Institute of Nutrition and Functional Foods (INAF), Université Laval, Québec, QC G1V 0A6, Canada; Centre Nutrition, santé et société (NUTRISS), Université Laval, Québec, QC G1V 0A6, Canada; School of Nutrition, Université Laval, Québec, QC G1V 0A6, Canada.
| |
Collapse
|
9
|
Zheng H, Xie T, Li S, Qiao X, Lu Y, Feng Y. Analysis of oral microbial dysbiosis associated with early childhood caries. BMC Oral Health 2021; 21:181. [PMID: 33827534 PMCID: PMC8028703 DOI: 10.1186/s12903-021-01543-x] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2020] [Accepted: 03/30/2021] [Indexed: 11/22/2022] Open
Abstract
Background "Core microbes" play a key role in the development of caries and lead to microbial disorders. Our goal was to detect the core microbes associated with the microbiota imbalance in early childhood caries (ECC). Methods Fifteen caries-free children and fifteen high-caries (DMFT ≥ 10) children aged 4–6 years old were recruited according to the diagnostic criteria of caries suggested by the WHO. The 16S rRNA genes from samples of plaque in saliva were amplified by PCR, and the PCR products were sequenced by the Illumina Miseq platform. The sequencing results were analyzed by professional software to determine the composition and structure of the saliva microorganisms. Results There were statistically significant differences between the groups regarding the relative abundance of Streptococcus mutans (Wilcoxon rank-sum test, P < 0.05). No significant difference was found between the groups regarding other species or functional genes. Conclusion S. mutans, together with other pathogens, may play a prominent role and act as "core microbes" in the occurrence and development of early childhood caries.
Collapse
Affiliation(s)
- Hui Zheng
- Stomatological Key Laboratory of Fujian College and University and Department of Orthodontics, School and Hospital of Stomatology, Fujian Medical University, Fuzhou, China
| | - Tengfei Xie
- Fujian Provincial Engineering Research Center of Oral Biomaterial and Department of Preventive Dentistry, School and Hospital of Stomatology, Fujian Medical University, Fuzhou, China
| | - Shaokai Li
- School and Hospital of Stomatology, Fujian Medical University, Fuzhou, China
| | - Xiaotong Qiao
- School and Hospital of Stomatology, Fujian Medical University, Fuzhou, China
| | - Youguang Lu
- Institute of Stomatology and Department of Preventive Dentistry, School and Hospital of Stomatology, Fujian Medical University, 246 Middle Yangqiao Road, Fuzhou, 350000, China.
| | - Yan Feng
- Fujian Key Laboratory of Oral Diseases and Department of Preventive Dentistry, School and Hospital of Stomatology, Fujian Medical University, 246 Middle Yangqiao Road, Fuzhou, 350000, China.
| |
Collapse
|
10
|
Irinyi L, Hu Y, Hoang MTV, Pasic L, Halliday C, Jayawardena M, Basu I, McKinney W, Morris AJ, Rathjen J, Stone E, Chen S, Sorrell TC, Schwessinger B, Meyer W. Long-read sequencing based clinical metagenomics for the detection and confirmation of Pneumocystis jirovecii directly from clinical specimens: A paradigm shift in mycological diagnostics. Med Mycol 2021; 58:650-660. [PMID: 31758176 DOI: 10.1093/mmy/myz109] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2019] [Revised: 08/26/2019] [Accepted: 11/19/2019] [Indexed: 11/14/2022] Open
Abstract
The advent of next generation sequencing technologies has enabled the characterization of the genetic content of entire communities of organisms, including those in clinical specimens, without prior culturing. The MinION from Oxford Nanopore Technologies offers real-time, direct sequencing of long DNA fragments directly from clinical samples. The aim of this study was to assess the ability of unbiased, genome-wide, long-read, shotgun sequencing using MinION to identify Pneumocystis jirovecii directly from respiratory tract specimens and to characterize the associated mycobiome. Pneumocystis pneumonia (PCP) is a life-threatening fungal disease caused by P. jirovecii. Currently, the diagnosis of PCP relies on direct microscopic or real-time quantitative polymerase chain reaction (PCR) examination of respiratory tract specimens, as P. jirovecii cannot be cultured readily in vitro. P. jirovecii DNA was detected in bronchoalveolar lavage (BAL) and induced sputum (IS) samples from three patients with confirmed PCP. Other fungi present in the associated mycobiome included known human pathogens (Aspergillus, Cryptococcus, Pichia) as well as commensal species (Candida, Malassezia, Bipolaris). We have established optimized sample preparation conditions for the generation of high-quality data, curated databases, and data analysis tools, which are key to the application of long-read MinION sequencing leading to a fundamental new approach in fungal diagnostics.
Collapse
Affiliation(s)
- Laszlo Irinyi
- Molecular Mycology Research Laboratory, Centre for Infectious Diseases and Microbiology, Faculty of Medicine and Health, Sydney Medical School, Westmead Clinical School, The University of Sydney, Sydney, NSW, Australia.,Marie Bashir Institute for Infectious Diseases and Biosecurity, The University of Sydney, Sydney, NSW, Australia.,Westmead Institute for Medical Research, Westmead, NSW Australia
| | - Yiheng Hu
- Research School of Biology, Australian National University, Canberra, ACT, Australia
| | - Minh Thuy Vi Hoang
- Molecular Mycology Research Laboratory, Centre for Infectious Diseases and Microbiology, Faculty of Medicine and Health, Sydney Medical School, Westmead Clinical School, The University of Sydney, Sydney, NSW, Australia.,Marie Bashir Institute for Infectious Diseases and Biosecurity, The University of Sydney, Sydney, NSW, Australia.,Westmead Institute for Medical Research, Westmead, NSW Australia
| | - Lana Pasic
- Molecular Mycology Research Laboratory, Centre for Infectious Diseases and Microbiology, Faculty of Medicine and Health, Sydney Medical School, Westmead Clinical School, The University of Sydney, Sydney, NSW, Australia.,Marie Bashir Institute for Infectious Diseases and Biosecurity, The University of Sydney, Sydney, NSW, Australia.,Westmead Institute for Medical Research, Westmead, NSW Australia
| | - Catriona Halliday
- Centre for Infectious Diseases and Microbiology-Laboratory Services, Institute for Clinical Pathology and Medical Research, NSW Health Pathology, Westmead, NSW, Australia
| | - Menuk Jayawardena
- Centre for Infectious Diseases and Microbiology-Laboratory Services, Institute for Clinical Pathology and Medical Research, NSW Health Pathology, Westmead, NSW, Australia
| | - Indira Basu
- Microbiology Department, LabPLUS, Auckland City Hospital, Auckland, New Zealand
| | - Wendy McKinney
- Microbiology Department, LabPLUS, Auckland City Hospital, Auckland, New Zealand
| | - Arthur J Morris
- Microbiology Department, LabPLUS, Auckland City Hospital, Auckland, New Zealand
| | - John Rathjen
- Research School of Biology, Australian National University, Canberra, ACT, Australia
| | - Eric Stone
- Research School of Biology, Australian National University, Canberra, ACT, Australia.,ANU-CSIRO Centre for Genomics, Metabolomics and Bioinformatics, Canberra, ACT, Australia
| | - Sharon Chen
- Molecular Mycology Research Laboratory, Centre for Infectious Diseases and Microbiology, Faculty of Medicine and Health, Sydney Medical School, Westmead Clinical School, The University of Sydney, Sydney, NSW, Australia.,Marie Bashir Institute for Infectious Diseases and Biosecurity, The University of Sydney, Sydney, NSW, Australia.,Centre for Infectious Diseases and Microbiology-Laboratory Services, Institute for Clinical Pathology and Medical Research, NSW Health Pathology, Westmead, NSW, Australia
| | - Tania C Sorrell
- Molecular Mycology Research Laboratory, Centre for Infectious Diseases and Microbiology, Faculty of Medicine and Health, Sydney Medical School, Westmead Clinical School, The University of Sydney, Sydney, NSW, Australia.,Marie Bashir Institute for Infectious Diseases and Biosecurity, The University of Sydney, Sydney, NSW, Australia
| | - Benjamin Schwessinger
- Research School of Biology, Australian National University, Canberra, ACT, Australia
| | - Wieland Meyer
- Molecular Mycology Research Laboratory, Centre for Infectious Diseases and Microbiology, Faculty of Medicine and Health, Sydney Medical School, Westmead Clinical School, The University of Sydney, Sydney, NSW, Australia.,Marie Bashir Institute for Infectious Diseases and Biosecurity, The University of Sydney, Sydney, NSW, Australia.,Westmead Institute for Medical Research, Westmead, NSW Australia.,Westmead Hospital (Research and Education Network), Westmead, NSW, Australia
| |
Collapse
|
11
|
Morissette A, Kropp C, Songpadith JP, Junges Moreira R, Costa J, Mariné-Casadó R, Pilon G, Varin TV, Dudonné S, Boutekrabt L, St-Pierre P, Levy E, Roy D, Desjardins Y, Raymond F, Houde VP, Marette A. Blueberry proanthocyanidins and anthocyanins improve metabolic health through a gut microbiota-dependent mechanism in diet-induced obese mice. Am J Physiol Endocrinol Metab 2020; 318:E965-E980. [PMID: 32228321 DOI: 10.1152/ajpendo.00560.2019] [Citation(s) in RCA: 59] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Blueberry consumption can prevent obesity-linked metabolic diseases, and it has been proposed that the polyphenol content of blueberries may contribute to these effects. Polyphenols have been shown to favorably impact metabolic health, but the role of specific polyphenol classes and whether the gut microbiota is linked to these effects remain unclear. We aimed to evaluate the impact of whole blueberry powder and blueberry polyphenols on the development of obesity and insulin resistance and to determine the potential role of gut microbes in these effects by using fecal microbiota transplantation (FMT). Sixty-eight C57BL/6 male mice were assigned to one of the following diets for 12 wk: balanced diet (Chow); high-fat, high-sucrose diet (HFHS); or HFHS supplemented with whole blueberry powder (BB), anthocyanidin (ANT)-rich extract, or proanthocyanidin (PAC)-rich extract. After 8 wk, mice were housed in metabolic cages, and an oral glucose tolerance test (OGTT) was performed. Sixty germ-free mice fed HFHS diet received FMT from one of the above groups biweekly for 8 wk, followed by an OGTT. PAC-treated mice were leaner than HFHS controls although they had the same energy intake and were more physically active. This observation was reproduced in germ-free mice receiving FMT from PAC-treated mice. PAC- and ANT-treated mice showed improved insulin responses during OGTT, and this finding was also reproduced in germ-free mice following FMT. These results show that blueberry PAC and ANT polyphenols can reduce diet-induced body weight and improve insulin sensitivity and that at least part of these beneficial effects are explained by modulation of the gut microbiota.
Collapse
Affiliation(s)
- Arianne Morissette
- Department of Medicine, Faculty of Medicine, Cardiology Axis of the Québec Heart and Lung Institute, Laval University, Hôpital Laval, Quebec City, Quebec, Canada
- Institute of Nutrition and Functional Foods, Laval University, Quebec City, Quebec, Canada
| | - Camille Kropp
- Department of Medicine, Faculty of Medicine, Cardiology Axis of the Québec Heart and Lung Institute, Laval University, Hôpital Laval, Quebec City, Quebec, Canada
| | - Jean-Philippe Songpadith
- Department of Medicine, Faculty of Medicine, Cardiology Axis of the Québec Heart and Lung Institute, Laval University, Hôpital Laval, Quebec City, Quebec, Canada
- Institute of Nutrition and Functional Foods, Laval University, Quebec City, Quebec, Canada
| | - Rafael Junges Moreira
- Department of Medicine, Faculty of Medicine, Cardiology Axis of the Québec Heart and Lung Institute, Laval University, Hôpital Laval, Quebec City, Quebec, Canada
| | - Janice Costa
- Department of Medicine, Faculty of Medicine, Cardiology Axis of the Québec Heart and Lung Institute, Laval University, Hôpital Laval, Quebec City, Quebec, Canada
| | - Roger Mariné-Casadó
- Department of Medicine, Faculty of Medicine, Cardiology Axis of the Québec Heart and Lung Institute, Laval University, Hôpital Laval, Quebec City, Quebec, Canada
| | - Geneviève Pilon
- Department of Medicine, Faculty of Medicine, Cardiology Axis of the Québec Heart and Lung Institute, Laval University, Hôpital Laval, Quebec City, Quebec, Canada
| | - Thibault V Varin
- Institute of Nutrition and Functional Foods, Laval University, Quebec City, Quebec, Canada
| | - Stéphanie Dudonné
- Institute of Nutrition and Functional Foods, Laval University, Quebec City, Quebec, Canada
| | - Lemia Boutekrabt
- Institute of Nutrition and Functional Foods, Laval University, Quebec City, Quebec, Canada
| | - Philippe St-Pierre
- Department of Medicine, Faculty of Medicine, Cardiology Axis of the Québec Heart and Lung Institute, Laval University, Hôpital Laval, Quebec City, Quebec, Canada
| | - Emile Levy
- Institute of Nutrition and Functional Foods, Laval University, Quebec City, Quebec, Canada
- Research Centre, CHU-Sainte-Justine, Montreal, Quebec, Canada
- Department of Nutrition, Université de Montréal, Montreal, Quebec, Canada
| | - Denis Roy
- Institute of Nutrition and Functional Foods, Laval University, Quebec City, Quebec, Canada
| | - Yves Desjardins
- Institute of Nutrition and Functional Foods, Laval University, Quebec City, Quebec, Canada
| | - Frédéric Raymond
- Institute of Nutrition and Functional Foods, Laval University, Quebec City, Quebec, Canada
| | - Vanessa P Houde
- Department of Medicine, Faculty of Medicine, Cardiology Axis of the Québec Heart and Lung Institute, Laval University, Hôpital Laval, Quebec City, Quebec, Canada
| | - André Marette
- Department of Medicine, Faculty of Medicine, Cardiology Axis of the Québec Heart and Lung Institute, Laval University, Hôpital Laval, Quebec City, Quebec, Canada
- Institute of Nutrition and Functional Foods, Laval University, Quebec City, Quebec, Canada
| |
Collapse
|
12
|
Kobus R, Abuín JM, Müller A, Hellmann SL, Pichel JC, Pena TF, Hildebrandt A, Hankeln T, Schmidt B. A big data approach to metagenomics for all-food-sequencing. BMC Bioinformatics 2020; 21:102. [PMID: 32164527 PMCID: PMC7069206 DOI: 10.1186/s12859-020-3429-6] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2019] [Accepted: 02/24/2020] [Indexed: 01/31/2023] Open
Abstract
BACKGROUND All-Food-Sequencing (AFS) is an untargeted metagenomic sequencing method that allows for the detection and quantification of food ingredients including animals, plants, and microbiota. While this approach avoids some of the shortcomings of targeted PCR-based methods, it requires the comparison of sequence reads to large collections of reference genomes. The steadily increasing amount of available reference genomes establishes the need for efficient big data approaches. RESULTS We introduce an alignment-free k-mer based method for detection and quantification of species composition in food and other complex biological matters. It is orders-of-magnitude faster than our previous alignment-based AFS pipeline. In comparison to the established tools CLARK, Kraken2, and Kraken2+Bracken it is superior in terms of false-positive rate and quantification accuracy. Furthermore, the usage of an efficient database partitioning scheme allows for the processing of massive collections of reference genomes with reduced memory requirements on a workstation (AFS-MetaCache) or on a Spark-based compute cluster (MetaCacheSpark). CONCLUSIONS We present a fast yet accurate screening method for whole genome shotgun sequencing-based biosurveillance applications such as food testing. By relying on a big data approach it can scale efficiently towards large-scale collections of complex eukaryotic and bacterial reference genomes. AFS-MetaCache and MetaCacheSpark are suitable tools for broad-scale metagenomic screening applications. They are available at https://muellan.github.io/metacache/afs.html (C++ version for a workstation) and https://github.com/jmabuin/MetaCacheSpark (Spark version for big data clusters).
Collapse
Affiliation(s)
- Robin Kobus
- Department of Computer Science, Johannes Gutenberg University, Mainz, 55099 Germany
| | - José M. Abuín
- IPCA, Polytechnic Institute of Cávado and Ave, Barcelos, 4750-810 Portugal
- CiTIUS, Universidade de Santiago de Compostela, Santiago de Compostela, 15782 Spain
| | - André Müller
- Department of Computer Science, Johannes Gutenberg University, Mainz, 55099 Germany
| | - Sören Lukas Hellmann
- Molecular Genetics and Genome Analysis, Institute of Organismal and Molecular Evolution, Johannes Gutenberg University, Mainz, 55099 Germany
| | - Juan C. Pichel
- CiTIUS, Universidade de Santiago de Compostela, Santiago de Compostela, 15782 Spain
| | - Tomás F. Pena
- CiTIUS, Universidade de Santiago de Compostela, Santiago de Compostela, 15782 Spain
| | - Andreas Hildebrandt
- Department of Computer Science, Johannes Gutenberg University, Mainz, 55099 Germany
| | - Thomas Hankeln
- Molecular Genetics and Genome Analysis, Institute of Organismal and Molecular Evolution, Johannes Gutenberg University, Mainz, 55099 Germany
| | - Bertil Schmidt
- Department of Computer Science, Johannes Gutenberg University, Mainz, 55099 Germany
| |
Collapse
|
13
|
Liu Y, Bible PW, Zou B, Liang Q, Dong C, Wen X, Li Y, Ge X, Li X, Deng X, Ma R, Guo S, Liang J, Chen T, Pan W, Liu L, Chen W, Wang X, Wei L. CSMD: a computational subtraction-based microbiome discovery pipeline for species-level characterization of clinical metagenomic samples. Bioinformatics 2020; 36:1577-1583. [PMID: 31626280 DOI: 10.1093/bioinformatics/btz790] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2019] [Revised: 09/22/2019] [Accepted: 10/16/2019] [Indexed: 11/12/2022] Open
Abstract
MOTIVATION Microbiome analyses of clinical samples with low microbial biomass are challenging because of the very small quantities of microbial DNA relative to the human host, ubiquitous contaminating DNA in sequencing experiments and the large and rapidly growing microbial reference databases. RESULTS We present computational subtraction-based microbiome discovery (CSMD), a bioinformatics pipeline specifically developed to generate accurate species-level microbiome profiles for clinical samples with low microbial loads. CSMD applies strategies for the maximal elimination of host sequences with minimal loss of microbial signal and effectively detects microorganisms present in the sample with minimal false positives using a stepwise convergent solution. CSMD was benchmarked in a comparative evaluation with other classic tools on previously published well-characterized datasets. It showed higher sensitivity and specificity in host sequence removal and higher specificity in microbial identification, which led to more accurate abundance estimation. All these features are integrated into a free and easy-to-use tool. Additionally, CSMD applied to cell-free plasma DNA showed that microbial diversity within these samples is substantially broader than previously believed. AVAILABILITY AND IMPLEMENTATION CSMD is freely available at https://github.com/liuyu8721/csmd. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Yu Liu
- Department of Statistical Science, School of Mathematics, Sun Yat-Sen University, Guangzhou 510275, China.,State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou 510060, China
| | - Paul W Bible
- State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou 510060, China.,College of Arts and Sciences, Marian University, Indianapolis, IN 46222, USA
| | - Bin Zou
- State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou 510060, China
| | - Qiaoxing Liang
- State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou 510060, China
| | - Cong Dong
- College of Chemistry, Sun Yat-Sen University, Guangzhou 510275, China
| | - Xiaofeng Wen
- State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou 510060, China
| | - Yan Li
- State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou 510060, China
| | - Xiaofei Ge
- State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou 510060, China
| | - Xifang Li
- State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou 510060, China
| | - Xiuli Deng
- State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou 510060, China
| | - Rong Ma
- State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou 510060, China
| | - Shixin Guo
- State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou 510060, China
| | - Juanran Liang
- State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou 510060, China
| | - Tingting Chen
- State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou 510060, China
| | - Wenliang Pan
- Department of Statistical Science, School of Mathematics, Sun Yat-Sen University, Guangzhou 510275, China
| | - Lixin Liu
- College of Chemistry, Sun Yat-Sen University, Guangzhou 510275, China
| | - Wei Chen
- Department of Biostatistics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, PA 15261, USA.,Division of Pulmonary Medicine, Department of Pediatrics, Children's Hospital of Pittsburgh of UPMC, University of Pittsburgh, Pittsburgh, PA 15224, USA
| | - Xueqin Wang
- Department of Statistical Science, School of Mathematics, Sun Yat-Sen University, Guangzhou 510275, China.,Southern China Research Center of Statistical Science, Sun Yat-Sen University, Guangzhou 510275, China
| | - Lai Wei
- State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou 510060, China
| |
Collapse
|
14
|
Abstract
The bacterial communities that live within the human gut have been linked to health and disease. However, we are still just beginning to understand how those bacteria interact and what potential interventions to our gut microbiome can make us healthier. Here, we present a mathematical modeling framework (named MICOM) that can recapitulate the growth rates of diverse bacterial species in the gut and can simulate metabolic interactions within microbial communities. We show that MICOM can unravel the ecological rules that shape the microbial landscape in our gut and that a given dietary or probiotic intervention can have widely different effects in different people. Compositional changes in the gut microbiota have been associated with a variety of medical conditions such as obesity, Crohn’s disease, and diabetes. However, connecting microbial community composition to ecosystem function remains a challenge. Here, we introduce MICOM, a customizable metabolic model of the human gut microbiome. By using a heuristic optimization approach based on L2 regularization, we were able to obtain a unique set of realistic growth rates that corresponded well with observed replication rates. We integrated adjustable dietary and taxon abundance constraints to generate personalized metabolic models for individual metagenomic samples. We applied MICOM to a balanced cohort of metagenomes from 186 people, including a metabolically healthy population and individuals with type 1 and type 2 diabetes. Model results showed that individual bacterial genera maintained conserved niche structures across humans, while the community-level production of short-chain fatty acids (SCFAs) was heterogeneous and highly individual specific. Model output revealed complex cross-feeding interactions that would be difficult to measure in vivo. Metabolic interaction networks differed somewhat consistently between healthy and diabetic subjects. In particular, MICOM predicted reduced butyrate and propionate production in a diabetic cohort, with restoration of SCFA production profiles found in healthy subjects following metformin treatment. Overall, we found that changes in diet or taxon abundances have highly personalized effects. We believe MICOM can serve as a useful tool for generating mechanistic hypotheses for how diet and microbiome composition influence community function. All methods are implemented in an open-source Python package, which is available at https://github.com/micom-dev/micom. IMPORTANCE The bacterial communities that live within the human gut have been linked to health and disease. However, we are still just beginning to understand how those bacteria interact and what potential interventions to our gut microbiome can make us healthier. Here, we present a mathematical modeling framework (named MICOM) that can recapitulate the growth rates of diverse bacterial species in the gut and can simulate metabolic interactions within microbial communities. We show that MICOM can unravel the ecological rules that shape the microbial landscape in our gut and that a given dietary or probiotic intervention can have widely different effects in different people.
Collapse
|
15
|
Kauter A, Epping L, Semmler T, Antao EM, Kannapin D, Stoeckle SD, Gehlen H, Lübke-Becker A, Günther S, Wieler LH, Walther B. The gut microbiome of horses: current research on equine enteral microbiota and future perspectives. Anim Microbiome 2019; 1:14. [PMID: 33499951 PMCID: PMC7807895 DOI: 10.1186/s42523-019-0013-3] [Citation(s) in RCA: 57] [Impact Index Per Article: 11.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2019] [Accepted: 09/09/2019] [Indexed: 12/26/2022] Open
Abstract
Understanding the complex interactions of microbial communities including bacteria, archaea, parasites, viruses and fungi of the gastrointestinal tract (GIT) associated with states of either health or disease is still an expanding research field in both, human and veterinary medicine. GIT disorders and their consequences are among the most important diseases of domesticated Equidae, but current gaps of knowledge hinder adequate progress with respect to disease prevention and microbiome-based interventions. Current literature on enteral microbiomes mirrors a vast data and knowledge imbalance, with only few studies tackling archaea, viruses and eukaryotes compared with those addressing the bacterial components.Until recently, culture-dependent methods were used for the identification and description of compositional changes of enteral microorganisms, limiting the outcome to cultivatable bacteria only. Today, next generation sequencing technologies provide access to the entirety of genes (microbiome) associated with the microorganisms of the equine GIT including the mass of uncultured microbiota, or "microbial dark matter".This review illustrates methods commonly used for enteral microbiome analysis in horses and summarizes key findings reached for bacteria, viruses and fungi so far. Moreover, reasonable possibilities to combine different explorative techniques are described. As a future perspective, knowledge expansion concerning beneficial compositions of microorganisms within the equine GIT creates novel possibilities for early disorder diagnostics as well as innovative therapeutic approaches. In addition, analysis of shotgun metagenomic data enables tracking of certain microorganisms beyond species barriers: transmission events of bacteria including pathogens and opportunists harboring antibiotic resistance factors between different horses but also between humans and horses will reach new levels of depth concerning strain-level distinctions.
Collapse
Affiliation(s)
- Anne Kauter
- Advanced Light and Electron Microscopy (ZBS-4), Robert Koch Institute, Seestraße 10, 13353, Berlin, Germany
| | - Lennard Epping
- Microbial Genomics (NG1), Robert Koch Institute, Berlin, Germany
| | - Torsten Semmler
- Microbial Genomics (NG1), Robert Koch Institute, Berlin, Germany
| | | | - Dania Kannapin
- Equine Clinic, Surgery and Radiology, Freie Universität Berlin, Berlin, Germany
| | - Sabita D Stoeckle
- Equine Clinic, Surgery and Radiology, Freie Universität Berlin, Berlin, Germany
| | - Heidrun Gehlen
- Equine Clinic, Surgery and Radiology, Freie Universität Berlin, Berlin, Germany
| | - Antina Lübke-Becker
- Institute of Microbiology and Epizootics, Centre for Infection Medicine, Freie Universität Berlin, Berlin, Germany
| | - Sebastian Günther
- Pharmaceutical Biology Institute of Pharmacy, Universität Greifswald, Greifswald, Germany
| | | | - Birgit Walther
- Advanced Light and Electron Microscopy (ZBS-4), Robert Koch Institute, Seestraße 10, 13353, Berlin, Germany.
| |
Collapse
|
16
|
Dadi TH, Siragusa E, Piro VC, Andrusch A, Seiler E, Renard BY, Reinert K. DREAM-Yara: an exact read mapper for very large databases with short update time. Bioinformatics 2019; 34:i766-i772. [PMID: 30423080 DOI: 10.1093/bioinformatics/bty567] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
Motivation Mapping-based approaches have become limited in their application to very large sets of references since computing an FM-index for very large databases (e.g. >10 GB) has become a bottleneck. This affects many analyses that need such index as an essential step for approximate matching of the NGS reads to reference databases. For instance, in typical metagenomics analysis, the size of the reference sequences has become prohibitive to compute a single full-text index on standard machines. Even on large memory machines, computing such index takes about 1 day of computing time. As a result, updates of indices are rarely performed. Hence, it is desirable to create an alternative way of indexing while preserving fast search times. Results To solve the index construction and update problem we propose the DREAM (Dynamic seaRchablE pArallel coMpressed index) framework and provide an implementation. The main contributions are the introduction of an approximate search distributor via a novel use of Bloom filters. We combine several Bloom filters to form an interleaved Bloom filter and use this new data structure to quickly exclude reads for parts of the databases where they cannot match. This allows us to keep the databases in several indices which can be easily rebuilt if parts are updated while maintaining a fast search time. The second main contribution is an implementation of DREAM-Yara a distributed version of a fully sensitive read mapper under the DREAM framework. Availability and implementation https://gitlab.com/pirovc/dream_yara/.
Collapse
Affiliation(s)
| | - Enrico Siragusa
- Computational Genomics, IBM Thomas J Watson Research Center, Yorktown Heights, NY, USA
| | - Vitor C Piro
- Bioinformatics Unit (MF1), Robert Koch Institute, Berlin, Germany.,CAPES Foundation, Ministry of Education of Brazil, Brasília DF, Brazil
| | - Andreas Andrusch
- Centre for Biological Threats and Special Pathogens (ZBS1), Robert Koch Institute, Berlin, Germany
| | - Enrico Seiler
- Algorithmic Bioinformatics, Institute for Bioinformatics, FU Berlin, Berlin, Germany
| | | | - Knut Reinert
- Algorithmic Bioinformatics, Institute for Bioinformatics, FU Berlin, Berlin, Germany
| |
Collapse
|
17
|
Velsko IM, Perez MS, Richards VP. Resolving Phylogenetic Relationships for Streptococcus mitis and Streptococcus oralis through Core- and Pan-Genome Analyses. Genome Biol Evol 2019; 11:1077-1087. [PMID: 30847473 PMCID: PMC6461889 DOI: 10.1093/gbe/evz049] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/05/2019] [Indexed: 02/06/2023] Open
Abstract
Taxonomic and phylogenetic relationships of Streptococcus mitis and Streptococcus oralis have been difficult to establish biochemically and genetically. We used core-genome analyses of S. mitis and S. oralis, as well as the closely related species Streptococcus pneumoniae and Streptococcus parasanguinis, to clarify the phylogenetic relationships between S. mitis and S. oralis, as well as within subclades of S. oralis. All S. mitis (n = 67), S. oralis (n = 89), S. parasanguinis (n = 27), and 27 S. pneumoniae genome assemblies were downloaded from NCBI and reannotated. All genes were delineated into homologous clusters and maximum-likelihood phylogenies built from putatively nonrecombinant core gene sets. Population structure was determined using Bayesian genome clustering, and patristic distance was calculated between populations. Population-specific gene content was assessed using a phylogenetic-based genome-wide association approach. Streptococcus mitis and S. oralis formed distinct clades, but species mixing suggests taxonomic misassignment. Patristic distance between populations suggests that S. oralis subsp. dentisani is a distinct species, whereas S. oralis subsp. tigurinus and subsp. oralis are supported as subspecies, and that S. mitis comprises two subspecies. None of the genes within the pan-genomes of S. mitis and S. oralis could be statistically correlated with either, and the dispensable genomes showed extensive variation among isolates. These are likely important factors contributing to established overlap in biochemical characteristics for these taxa. Based on core-genome analysis, the substructure of S. oralis and S. mitis should be redefined, and species assignments within S. oralis and S. mitis should be made based on whole-genome analysis to be robust to misassignment.
Collapse
Affiliation(s)
| | - Megan S Perez
- Department of Biological Sciences, Clemson University
- Department of Arts and Sciences, LeTourneau University
| | | |
Collapse
|
18
|
O'Connell CM, Brochu H, Girardi J, Harrell E, Jones A, Darville T, Seña AC, Peng X. Simultaneous profiling of sexually transmitted bacterial pathogens, microbiome, and concordant host response in cervical samples using whole transcriptome sequencing analysis. MICROBIAL CELL 2019; 6:177-183. [PMID: 30854394 PMCID: PMC6402362 DOI: 10.15698/mic2019.03.672] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Pelvic inflammatory disease (PID) is a female upper genital tract inflammatory disorder that arises after sexually transmitted bacterial infections (STI). Factors modulating risk for reproductive sequelae include co-infection, microbiota, host genetics and physiology. In a pilot study of cervical samples obtained from women at high risk for STIs, we examined the potential for unbiased characterization of host, pathogen and microbiome interactions using whole transcriptome sequencing analysis of ribosomal RNA-depleted total RNAs (Total RNA-Seq). Only samples from women with STI infection contained pathogen-specific sequences (3 to 38% transcriptome coverage). Simultaneously, we identified and quantified their active microbial communities. After integration with host-derived reads from the same data, we detected clustering of host transcriptional profiles that reflected microbiome differences and STI infection. Together, our study suggests that total RNA profiling will advance understanding of the interplay of pathogen, host and microbiota during natural infection and may reveal novel, outcome-relevant biomarkers.
Collapse
Affiliation(s)
- Catherine M O'Connell
- Department of Pediatrics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| | - Hayden Brochu
- Department of Molecular Biomedical Sciences, North Carolina State University, Raleigh, North Carolina, USA
| | - Jenna Girardi
- Department of Pediatrics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| | - Erin Harrell
- Department of Molecular Biomedical Sciences, North Carolina State University, Raleigh, North Carolina, USA
| | - Aiden Jones
- Department of Molecular Biomedical Sciences, North Carolina State University, Raleigh, North Carolina, USA
| | - Toni Darville
- Department of Pediatrics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| | - Arlene C Seña
- Department of Medicine, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| | - Xinxia Peng
- Department of Molecular Biomedical Sciences, North Carolina State University, Raleigh, North Carolina, USA.,Bioinformatics Research Center, North Carolina State University, Raleigh, North Carolina, USA
| |
Collapse
|
19
|
Breitwieser FP, Baker DN, Salzberg SL. KrakenUniq: confident and fast metagenomics classification using unique k-mer counts. Genome Biol 2018; 19:198. [PMID: 30445993 PMCID: PMC6238331 DOI: 10.1186/s13059-018-1568-0] [Citation(s) in RCA: 204] [Impact Index Per Article: 34.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2018] [Accepted: 10/18/2018] [Indexed: 12/30/2022] Open
Abstract
False-positive identifications are a significant problem in metagenomics classification. We present KrakenUniq, a novel metagenomics classifier that combines the fast k-mer-based classification of Kraken with an efficient algorithm for assessing the coverage of unique k-mers found in each species in a dataset. On various test datasets, KrakenUniq gives better recall and precision than other methods and effectively classifies and distinguishes pathogens with low abundance from false positives in infectious disease samples. By using the probabilistic cardinality estimator HyperLogLog, KrakenUniq runs as fast as Kraken and requires little additional memory. KrakenUniq is freely available at https://github.com/fbreitwieser/krakenuniq .
Collapse
Affiliation(s)
- F P Breitwieser
- Center for Computational Biology, McKusick-Nathans Institute of Genetic Medicine, Johns Hopkins School of Medicine, Baltimore, MD, USA.
| | - D N Baker
- Center for Computational Biology, McKusick-Nathans Institute of Genetic Medicine, Johns Hopkins School of Medicine, Baltimore, MD, USA
- Department of Computer Science, Johns Hopkins University, Baltimore, MD, USA
| | - S L Salzberg
- Center for Computational Biology, McKusick-Nathans Institute of Genetic Medicine, Johns Hopkins School of Medicine, Baltimore, MD, USA.
- Department of Computer Science, Johns Hopkins University, Baltimore, MD, USA.
- Departments of Biomedical Engineering and Biostatistics, Johns Hopkins University, Baltimore, MD, USA.
| |
Collapse
|
20
|
Walsh AM, Crispie F, O'Sullivan O, Finnegan L, Claesson MJ, Cotter PD. Species classifier choice is a key consideration when analysing low-complexity food microbiome data. MICROBIOME 2018; 6:50. [PMID: 29554948 PMCID: PMC5859664 DOI: 10.1186/s40168-018-0437-0] [Citation(s) in RCA: 46] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/01/2017] [Accepted: 03/05/2018] [Indexed: 05/03/2023]
Abstract
BACKGROUND The use of shotgun metagenomics to analyse low-complexity microbial communities in foods has the potential to be of considerable fundamental and applied value. However, there is currently no consensus with respect to choice of species classification tool, platform, or sequencing depth. Here, we benchmarked the performances of three high-throughput short-read sequencing platforms, the Illumina MiSeq, NextSeq 500, and Ion Proton, for shotgun metagenomics of food microbiota. Briefly, we sequenced six kefir DNA samples and a mock community DNA sample, the latter constructed by evenly mixing genomic DNA from 13 food-related bacterial species. A variety of bioinformatic tools were used to analyse the data generated, and the effects of sequencing depth on these analyses were tested by randomly subsampling reads. RESULTS Compositional analysis results were consistent between the platforms at divergent sequencing depths. However, we observed pronounced differences in the predictions from species classification tools. Indeed, PERMANOVA indicated that there was no significant differences between the compositional results generated by the different sequencers (p = 0.693, R2 = 0.011), but there was a significant difference between the results predicted by the species classifiers (p = 0.01, R2 = 0.127). The relative abundances predicted by the classifiers, apart from MetaPhlAn2, were apparently biased by reference genome sizes. Additionally, we observed varying false-positive rates among the classifiers. MetaPhlAn2 had the lowest false-positive rate, whereas SLIMM had the greatest false-positive rate. Strain-level analysis results were also similar across platforms. Each platform correctly identified the strains present in the mock community, but accuracy was improved slightly with greater sequencing depth. Notably, PanPhlAn detected the dominant strains in each kefir sample above 500,000 reads per sample. Again, the outputs from functional profiling analysis using SUPER-FOCUS were generally accordant between the platforms at different sequencing depths. Finally, and expectedly, metagenome assembly completeness was significantly lower on the MiSeq than either on the NextSeq (p = 0.03) or the Proton (p = 0.011), and it improved with increased sequencing depth. CONCLUSIONS Our results demonstrate a remarkable similarity in the results generated by the three sequencing platforms at different sequencing depths, and, in fact, the choice of bioinformatics methodology had a more evident impact on results than the choice of sequencer did.
Collapse
Affiliation(s)
- Aaron M Walsh
- Teagasc Food Research Centre, Moorepark, Fermoy, Co. Cork, Ireland
- APC Microbiome Institute, University College Cork, Co. Cork, Ireland
- Microbiology Department, University College Cork, Co. Cork, Ireland
| | - Fiona Crispie
- Teagasc Food Research Centre, Moorepark, Fermoy, Co. Cork, Ireland
- APC Microbiome Institute, University College Cork, Co. Cork, Ireland
| | - Orla O'Sullivan
- Teagasc Food Research Centre, Moorepark, Fermoy, Co. Cork, Ireland
- APC Microbiome Institute, University College Cork, Co. Cork, Ireland
| | - Laura Finnegan
- Teagasc Food Research Centre, Moorepark, Fermoy, Co. Cork, Ireland
- APC Microbiome Institute, University College Cork, Co. Cork, Ireland
| | - Marcus J Claesson
- APC Microbiome Institute, University College Cork, Co. Cork, Ireland
- Microbiology Department, University College Cork, Co. Cork, Ireland
| | - Paul D Cotter
- Teagasc Food Research Centre, Moorepark, Fermoy, Co. Cork, Ireland.
- APC Microbiome Institute, University College Cork, Co. Cork, Ireland.
| |
Collapse
|
21
|
Identification of fungi in shotgun metagenomics datasets. PLoS One 2018; 13:e0192898. [PMID: 29444186 PMCID: PMC5812651 DOI: 10.1371/journal.pone.0192898] [Citation(s) in RCA: 47] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2017] [Accepted: 01/31/2018] [Indexed: 01/22/2023] Open
Abstract
Metagenomics uses nucleic acid sequencing to characterize species diversity in different niches such as environmental biomes or the human microbiome. Most studies have used 16S rRNA amplicon sequencing to identify bacteria. However, the decreasing cost of sequencing has resulted in a gradual shift away from amplicon analyses and towards shotgun metagenomic sequencing. Shotgun metagenomic data can be used to identify a wide range of species, but have rarely been applied to fungal identification. Here, we develop a sequence classification pipeline, FindFungi, and use it to identify fungal sequences in public metagenome datasets. We focus primarily on animal metagenomes, especially those from pig and mouse microbiomes. We identified fungi in 39 of 70 datasets comprising 71 fungal species. At least 11 pathogenic species with zoonotic potential were identified, including Candida tropicalis. We identified Pseudogymnoascus species from 13 Antarctic soil samples initially analyzed for the presence of bacteria capable of degrading diesel oil. We also show that Candida tropicalis and Candida loboi are likely the same species. In addition, we identify several examples where contaminating DNA was erroneously included in fungal genome assemblies.
Collapse
|
22
|
Grześkowiak Ł, Martínez-Vallespín B, Dadi TH, Radloff J, Amasheh S, Heinsen FA, Franke A, Reinert K, Vahjen W, Zentek J, Pieper R. Formula Feeding Predisposes Neonatal Piglets to Clostridium difficile Gut Infection. J Infect Dis 2017; 217:1442-1452. [DOI: 10.1093/infdis/jix567] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2017] [Accepted: 10/24/2017] [Indexed: 12/20/2022] Open
Affiliation(s)
| | | | - Temesgen H Dadi
- Department of Mathematics and Computer Science, Institute of Computer Science, Kiel, Germany
- Max Planck Institute for Molecular Genetics, Berlin, Kiel, Germany
| | - Judith Radloff
- Institute of Veterinary Physiology, Freie Universität Berlin, Kiel, Germany
| | - Salah Amasheh
- Institute of Veterinary Physiology, Freie Universität Berlin, Kiel, Germany
| | - Femke-Anouska Heinsen
- Institute of Clinical Molecular Biology, Christian-Albrechts-University of Kiel, University Hospital Schleswig Holstein, Kiel, Germany
| | - Andre Franke
- Institute of Clinical Molecular Biology, Christian-Albrechts-University of Kiel, University Hospital Schleswig Holstein, Kiel, Germany
| | - Knut Reinert
- Department of Mathematics and Computer Science, Institute of Computer Science, Kiel, Germany
| | | | | | | |
Collapse
|