MetaCluster 4.0: a novel binning algorithm for NGS reads and huge number of species.
Next-generation sequencing (NGS) technologies allow the sequencing of microbial communities directly from the environment without prior culturing. The output of environmental DNA sequencing consists of many reads from genomes of different unknown species, making the clustering together reads from the same (or similar) species (also known as binning) a crucial step. The difficulties of the binning problem are due to the following four factors: (1) the lack of reference genomes; (2) uneven abundance ratio of species; (3) short NGS reads; and (4) a large number of species (can be more than a hundred). None of the existing binning tools can handle all four factors. No tools, including both AbundanceBin and MetaCluster 3.0, have demonstrated reasonable performance on a sample with more than 20 species. In this article, we introduce MetaCluster 4.0, an unsupervised binning algorithm that can accurately (with about 80% precision and sensitivity in all cases and at least 90% in some cases) and efficiently bin short reads with varying abundance ratios and is able to handle datasets with 100 species. The novelty of MetaCluster 4.0 stems from solving a few important problems: how to divide reads into groups by a probabilistic approach, how to estimate the 4-mer distribution of each group, how to estimate the number of species, and how to modify MetaCluster 3.0 to handle a large number of species. We show that Meta Cluster 4.0 is effective for both simulated and real datasets. Supplementary Material is available at www.liebertonline.com/cmb.
Department of Computer Science, The University of Hong Kong, Hong Kong, China.
Yi
Y
Wang
Henry C M
HC
Leung
S M
SM
Yiu
Francis Y L
FY
Chin
D000465:Algorithms; D001419:Bacteria / Q000235:genetics; D001483:Base Sequence; D016000:Cluster Analysis; D003627:Data Interpretation, Statistical; D016680:Genome, Bacterial; D059014:High-Throughput Nucleotide Sequencing; D015233:Models, Statistical; D017422:Sequence Analysis, DNA / Q000379:methods*; D012984:Software