{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,3,30]],"date-time":"2026-03-30T11:09:13Z","timestamp":1774868953526,"version":"3.50.1"},"reference-count":37,"publisher":"Oxford University Press (OUP)","issue":"5","license":[{"start":{"date-parts":[[2019,10,18]],"date-time":"2019-10-18T00:00:00Z","timestamp":1571356800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/academic.oup.com\/journals\/pages\/open_access\/funder_policies\/chorus\/standard_publication_model"}],"funder":[{"DOI":"10.13039\/501100012166","name":"National Basic Research Program of China","doi-asserted-by":"publisher","award":["2015CB964601"],"award-info":[{"award-number":["2015CB964601"]}],"id":[{"id":"10.13039\/501100012166","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100001809","name":"National Natural Science Foundation of China","doi-asserted-by":"publisher","award":["81570828"],"award-info":[{"award-number":["81570828"]}],"id":[{"id":"10.13039\/501100001809","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100001809","name":"National Natural Science Foundation of China","doi-asserted-by":"publisher","award":["11771462"],"award-info":[{"award-number":["11771462"]}],"id":[{"id":"10.13039\/501100001809","id-type":"DOI","asserted-by":"publisher"}]},{"name":"International Science and Technology Cooperation Program of Guangdong","award":["2016B050502007"],"award-info":[{"award-number":["2016B050502007"]}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2020,3,1]]},"abstract":"<jats:title>Abstract<\/jats:title><jats:sec><jats:title>Motivation<\/jats:title><jats:p>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.<\/jats:p><\/jats:sec><jats:sec><jats:title>Results<\/jats:title><jats:p>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.<\/jats:p><\/jats:sec><jats:sec><jats:title>Availability and implementation<\/jats:title><jats:p>CSMD is freely available at https:\/\/github.com\/liuyu8721\/csmd.<\/jats:p><\/jats:sec><jats:sec><jats:title>Supplementary information<\/jats:title><jats:p>Supplementary data are available at Bioinformatics online.<\/jats:p><\/jats:sec>","DOI":"10.1093\/bioinformatics\/btz790","type":"journal-article","created":{"date-parts":[[2019,10,16]],"date-time":"2019-10-16T22:20:14Z","timestamp":1571264414000},"page":"1577-1583","source":"Crossref","is-referenced-by-count":3,"title":["CSMD: a computational subtraction-based microbiome discovery pipeline for species-level characterization of clinical metagenomic samples"],"prefix":"10.1093","volume":"36","author":[{"given":"Yu","family":"Liu","sequence":"first","affiliation":[{"name":"Department of Statistical Science , School of Mathematics, Sun Yat-Sen University, Guangzhou 510275, China"},{"name":"State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center , Sun Yat-sen University, Guangzhou 510060, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Paul W","family":"Bible","sequence":"additional","affiliation":[{"name":"State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center , Sun Yat-sen University, Guangzhou 510060, China"},{"name":"College of Arts and Sciences, Marian University , Indianapolis, IN 46222, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-9243-9923","authenticated-orcid":false,"given":"Bin","family":"Zou","sequence":"additional","affiliation":[{"name":"State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center , Sun Yat-sen University, Guangzhou 510060, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Qiaoxing","family":"Liang","sequence":"additional","affiliation":[{"name":"State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center , Sun Yat-sen University, Guangzhou 510060, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Cong","family":"Dong","sequence":"additional","affiliation":[{"name":"College of Chemistry, Sun Yat-Sen University , Guangzhou 510275, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Xiaofeng","family":"Wen","sequence":"additional","affiliation":[{"name":"State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center , Sun Yat-sen University, Guangzhou 510060, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Yan","family":"Li","sequence":"additional","affiliation":[{"name":"State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center , Sun Yat-sen University, Guangzhou 510060, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Xiaofei","family":"Ge","sequence":"additional","affiliation":[{"name":"State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center , Sun Yat-sen University, Guangzhou 510060, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Xifang","family":"Li","sequence":"additional","affiliation":[{"name":"State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center , Sun Yat-sen University, Guangzhou 510060, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Xiuli","family":"Deng","sequence":"additional","affiliation":[{"name":"State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center , Sun Yat-sen University, Guangzhou 510060, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Rong","family":"Ma","sequence":"additional","affiliation":[{"name":"State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center , Sun Yat-sen University, Guangzhou 510060, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Shixin","family":"Guo","sequence":"additional","affiliation":[{"name":"State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center , Sun Yat-sen University, Guangzhou 510060, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Juanran","family":"Liang","sequence":"additional","affiliation":[{"name":"State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center , Sun Yat-sen University, Guangzhou 510060, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Tingting","family":"Chen","sequence":"additional","affiliation":[{"name":"State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center , Sun Yat-sen University, Guangzhou 510060, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Wenliang","family":"Pan","sequence":"additional","affiliation":[{"name":"Department of Statistical Science , School of Mathematics, Sun Yat-Sen University, Guangzhou 510275, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Lixin","family":"Liu","sequence":"additional","affiliation":[{"name":"College of Chemistry, Sun Yat-Sen University , Guangzhou 510275, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-7196-8703","authenticated-orcid":false,"given":"Wei","family":"Chen","sequence":"additional","affiliation":[{"name":"Department of Biostatistics, Graduate School of Public Health, University of Pittsburgh , Pittsburgh, PA 15261, USA"},{"name":"Division of Pulmonary Medicine, Department of Pediatrics, Children's Hospital of Pittsburgh of UPMC, University of Pittsburgh , Pittsburgh, PA 15224, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Xueqin","family":"Wang","sequence":"additional","affiliation":[{"name":"Department of Statistical Science , School of Mathematics, Sun Yat-Sen University, Guangzhou 510275, China"},{"name":"Southern China Research Center of Statistical Science, Sun Yat-Sen University , Guangzhou 510275, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Lai","family":"Wei","sequence":"additional","affiliation":[{"name":"State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center , Sun Yat-sen University, Guangzhou 510060, China"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2019,10,18]]},"reference":[{"key":"2023060910380799200_btz790-B1","doi-asserted-by":"crossref","first-page":"403","DOI":"10.1016\/S0022-2836(05)80360-2","article-title":"Basic local alignment search tool","volume":"215","author":"Altschul","year":"1990","journal-title":"J. Mol. Biol"},{"key":"2023060910380799200_btz790-B2","author":"Andrews","year":"2010"},{"key":"2023060910380799200_btz790-B3","doi-asserted-by":"crossref","first-page":"2114","DOI":"10.1093\/bioinformatics\/btu170","article-title":"Trimmomatic: a flexible trimmer for Illumina sequence data","volume":"30","author":"Bolger","year":"2014","journal-title":"Bioinformatics"},{"key":"2023060910380799200_btz790-B4","doi-asserted-by":"crossref","first-page":"59","DOI":"10.1038\/nmeth.3176","article-title":"Fast and sensitive protein alignment using DIAMOND","volume":"12","author":"Buchfink","year":"2015","journal-title":"Nat. Methods"},{"key":"2023060910380799200_btz790-B5","article-title":"The Human Oral Microbiome Database: a web accessible resource for investigating oral microbe taxonomic and genomic information","volume":"2010","author":"Chen","year":"2010","journal-title":"Database Oxford"},{"key":"2023060910380799200_btz790-B6","doi-asserted-by":"crossref","first-page":"e3138.","DOI":"10.7717\/peerj.3138","article-title":"SLIMM: species level identification of microorganisms from metagenomes","volume":"5","author":"Dadi","year":"2017","journal-title":"PeerJ"},{"key":"2023060910380799200_btz790-B7","doi-asserted-by":"crossref","first-page":"S17","DOI":"10.1186\/1471-2164-13-S8-S17","article-title":"Species-level classification of the vaginal microbiome","volume":"13 Suppl 8","author":"Fettweis","year":"2012","journal-title":"BMC Genomics"},{"key":"2023060910380799200_btz790-B8","doi-asserted-by":"crossref","first-page":"e19051.","DOI":"10.1371\/journal.pone.0019051","article-title":"CORE: a phylogenetically-curated 16S rDNA database of the core oral microbiome","volume":"6","author":"Griffen","year":"2011","journal-title":"PLoS One"},{"key":"2023060910380799200_btz790-B9","doi-asserted-by":"crossref","first-page":"1552","DOI":"10.1101\/gr.120618.111","article-title":"Integrative analysis of environmental sequences using MEGAN4","volume":"21","author":"Huson","year":"2011","journal-title":"Genome Res"},{"key":"2023060910380799200_btz790-B10","doi-asserted-by":"crossref","first-page":"357","DOI":"10.1038\/nmeth.3317","article-title":"HISAT: a fast spliced aligner with low memory requirements","volume":"12","author":"Kim","year":"2015","journal-title":"Nat. Methods"},{"key":"2023060910380799200_btz790-B151","doi-asserted-by":"crossref","first-page":"393","DOI":"10.1038\/nbt.1868","article-title":"PathSeq: software to identify or discover microbes by deep sequencing of human tissue","volume":"29","author":"Kostic","year":"2011","journal-title":"Nat. Biotechnol."},{"key":"2023060910380799200_btz790-B11","doi-asserted-by":"crossref","first-page":"9623","DOI":"10.1073\/pnas.1707009114","article-title":"Numerous uncharacterized and highly divergent microbes which colonize humans are revealed by circulating cell-free DNA","volume":"114","author":"Kowarsky","year":"2017","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"2023060910380799200_btz790-B12","doi-asserted-by":"crossref","first-page":"R25.","DOI":"10.1186\/gb-2009-10-3-r25","article-title":"Ultrafast and memory-efficient alignment of short DNA sequences to the human genome","volume":"10","author":"Langmead","year":"2009","journal-title":"Genome Biol"},{"key":"2023060910380799200_btz790-B13","doi-asserted-by":"crossref","first-page":"1754","DOI":"10.1093\/bioinformatics\/btp324","article-title":"Fast and accurate short read alignment with Burrows-Wheeler transform","volume":"25","author":"Li","year":"2009","journal-title":"Bioinformatics"},{"key":"2023060910380799200_btz790-B14","doi-asserted-by":"crossref","first-page":"1260","DOI":"10.1093\/bioinformatics\/btt147","article-title":"Analyzing genome coverage profiles with applications to quality control in metagenomics","volume":"29","author":"Lindner","year":"2013","journal-title":"Bioinformatics"},{"key":"2023060910380799200_btz790-B15","doi-asserted-by":"crossref","first-page":"e0117711.","DOI":"10.1371\/journal.pone.0117711","article-title":"Metagenomic profiling of known and unknown microbes with microbeGPS","volume":"10","author":"Lindner","year":"2015","journal-title":"PLoS One"},{"key":"2023060910380799200_btz790-B16","doi-asserted-by":"crossref","first-page":"182.","DOI":"10.1186\/s13059-017-1299-7","article-title":"Comprehensive benchmarking and ensemble approaches for metagenomic classifiers","volume":"18","author":"McIntyre","year":"2017","journal-title":"Genome Biol"},{"key":"2023060910380799200_btz790-B17","doi-asserted-by":"crossref","DOI":"10.1038\/ncomms11257","article-title":"Fast and sensitive taxonomic classification for metagenomics with Kaiju","volume":"7","author":"Menzel","year":"2016","journal-title":"Nat. Commun"},{"key":"2023060910380799200_btz790-B18","doi-asserted-by":"crossref","first-page":"51","DOI":"10.1186\/s13059-019-1646-y","article-title":"Assessing taxonomic metagenome profilers with OPAL","volume":"20","author":"Meyer","year":"2019","journal-title":"Genome Biol"},{"key":"2023060910380799200_btz790-B19","doi-asserted-by":"crossref","first-page":"1014","DOI":"10.1038\/s41467-019-08844-4","article-title":"Microbial abundance, activity and population genomic profiling with mOTUs2","volume":"10","author":"Milanese","year":"2019","journal-title":"Nat. Commun"},{"key":"2023060910380799200_btz790-B20","doi-asserted-by":"crossref","first-page":"e00218-17","DOI":"10.1128\/mSystems.00218-17","article-title":"KatharoSeq enables high-throughput microbiome analysis from low-biomass samples","volume":"3","author":"Minich","year":"2018","journal-title":"mSystems"},{"key":"2023060910380799200_btz790-B21","doi-asserted-by":"crossref","first-page":"3823","DOI":"10.1093\/bioinformatics\/btw542","article-title":"Higher classification sensitivity of short metagenomic reads with CLARK-S","volume":"32","author":"Ounit","year":"2016","journal-title":"Bioinformatics"},{"key":"2023060910380799200_btz790-B22","doi-asserted-by":"crossref","first-page":"236.","DOI":"10.1186\/s12864-015-1419-2","article-title":"CLARK: fast and accurate classification of metagenomic and genomic sequences using discriminative k-mers","volume":"16","author":"Ounit","year":"2015","journal-title":"BMC Genomics"},{"key":"2023060910380799200_btz790-B23","doi-asserted-by":"crossref","first-page":"649","DOI":"10.1016\/j.cell.2019.01.001","article-title":"Extensive unexplored human microbiome diversity revealed by over 150,000 genomes from metagenomes spanning age, geography, and lifestyle","volume":"176","author":"Pasolli","year":"2019","journal-title":"Cell"},{"key":"2023060910380799200_btz790-B24","doi-asserted-by":"crossref","first-page":"e1003107.","DOI":"10.1371\/journal.pcbi.1003107","article-title":"Bacteria-human somatic cell lateral gene transfer is enriched in cancer samples","volume":"9","author":"Riley","year":"2013","journal-title":"PLoS Comput. Biol"},{"key":"2023060910380799200_btz790-B25","doi-asserted-by":"crossref","first-page":"87","DOI":"10.1186\/s12915-014-0087-z","article-title":"Reagent and laboratory contamination can critically impact sequence-based microbiome analyses","volume":"12","author":"Salter","year":"2014","journal-title":"BMC Biol"},{"key":"2023060910380799200_btz790-B26","doi-asserted-by":"crossref","first-page":"269","DOI":"10.1002\/gepi.20575","article-title":"Efficient study design for next generation sequencing","volume":"35","author":"Sampson","year":"2011","journal-title":"Genet. Epidemiol"},{"key":"2023060910380799200_btz790-B27","doi-asserted-by":"crossref","first-page":"1063","DOI":"10.1038\/nmeth.4458","article-title":"Critical assessment of metagenome interpretation\u2014a benchmark of metagenomics software","volume":"14","author":"Sczyrba","year":"2017","journal-title":"Nat. Methods"},{"key":"2023060910380799200_btz790-B28","doi-asserted-by":"crossref","first-page":"811","DOI":"10.1038\/nmeth.2066","article-title":"Metagenomic microbial community profiling using unique clade-specific marker genes","volume":"9","author":"Segata","year":"2012","journal-title":"Nat. Methods"},{"key":"2023060910380799200_btz790-B29","doi-asserted-by":"crossref","first-page":"e00198.","DOI":"10.1128\/mBio.00198-16","article-title":"Changes in the eye microbiota associated with contact lens wearing","volume":"7","author":"Shin","year":"2016","journal-title":"mBio"},{"key":"2023060910380799200_btz790-B30","doi-asserted-by":"crossref","first-page":"6589","DOI":"10.1038\/s41598-017-06665-3","article-title":"Large-scale differences in microbial biodiversity discovery between 16S amplicon and shotgun sequencing","volume":"7","author":"Tessler","year":"2017","journal-title":"Sci. Rep.-UK"},{"key":"2023060910380799200_btz790-B31","doi-asserted-by":"crossref","first-page":"902","DOI":"10.1038\/nmeth.3589","article-title":"MetaPhlAn2 for enhanced metagenomic taxonomic profiling","volume":"12","author":"Truong","year":"2015","journal-title":"Nat. Methods"},{"key":"2023060910380799200_btz790-B32","doi-asserted-by":"crossref","first-page":"4287","DOI":"10.1093\/bioinformatics\/bty501","article-title":"GATK PathSeq: a customizable computational tool for the discovery and identification of microbial sequences in libraries from eukaryotic hosts","volume":"34","author":"Walker","year":"2018","journal-title":"Bioinformatics"},{"key":"2023060910380799200_btz790-B33","doi-asserted-by":"crossref","first-page":"6030","DOI":"10.1167\/iovs.17-22957","article-title":"The influence of age and sex on ocular surface microbiota in healthy adults","volume":"58","author":"Wen","year":"2017","journal-title":"Invest. Ophth. Vis. Sci"},{"key":"2023060910380799200_btz790-B34","doi-asserted-by":"crossref","first-page":"R46.","DOI":"10.1186\/gb-2014-15-3-r46","article-title":"Kraken: ultrafast metagenomic sequence classification using exact alignments","volume":"15","author":"Wood","year":"2014","journal-title":"Genome Biol"},{"key":"2023060910380799200_btz790-B35","doi-asserted-by":"crossref","first-page":"265.","DOI":"10.1186\/s13059-015-0821-z","article-title":"Identification of low abundance microbiome in clinical samples using whole genome sequencing","volume":"16","author":"Zhang","year":"2015","journal-title":"Genome Biol"},{"key":"2023060910380799200_btz790-B36","doi-asserted-by":"crossref","first-page":"128","DOI":"10.1167\/iovs.16-20231","article-title":"Conjunctival microbiome changes associated with soft contact lens and orthokeratology lens wearing","volume":"58","author":"Zhang","year":"2017","journal-title":"Invest. Ophth. Vis. Sci"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"http:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btz790\/30458583\/btz790.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/36\/5\/1577\/50552950\/bioinformatics_36_5_1577.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/36\/5\/1577\/50552950\/bioinformatics_36_5_1577.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2024,7,25]],"date-time":"2024-07-25T03:35:05Z","timestamp":1721878505000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/36\/5\/1577\/5593651"}},"subtitle":[],"editor":[{"given":"Jonathan","family":"Wren","sequence":"additional","affiliation":[],"role":[{"role":"editor","vocabulary":"crossref"}]}],"short-title":[],"issued":{"date-parts":[[2019,10,18]]},"references-count":37,"journal-issue":{"issue":"5","published-print":{"date-parts":[[2020,3,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btz790","relation":{},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"value":"1367-4803","type":"print"},{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2020,3]]},"published":{"date-parts":[[2019,10,18]]}}}