 Methodology
 Open Access
 Published:
METAMVGL: a multiview graphbased metagenomic contig binning algorithm by integrating assembly and pairedend graphs
BMC Bioinformatics volume 22, Article number: 378 (2021)
Abstract
Background
Due to the complexity of microbial communities, de novo assembly on next generation sequencing data is commonly unable to produce complete microbial genomes. Metagenome assembly binning becomes an essential step that could group the fragmented contigs into clusters to represent microbial genomes based on contigs’ nucleotide compositions and read depths. These features work well on the long contigs, but are not stable for the short ones. Contigs can be linked by sequence overlap (assembly graph) or by the pairedend reads aligned to them (PE graph), where the linked contigs have high chance to be derived from the same clusters.
Results
We developed METAMVGL, a multiview graphbased metagenomic contig binning algorithm by integrating both assembly and PE graphs. It could strikingly rescue the short contigs and correct the binning errors from dead ends. METAMVGL learns the two graphs’ weights automatically and predicts the contig labels in a uniform multiview label propagation framework. In experiments, we observed METAMVGL made use of significantly more highconfidence edges from the combined graph and linked dead ends to the main graph. It also outperformed many stateoftheart contig binning algorithms, including MaxBin2, MetaBAT2, MyCC, CONCOCT, SolidBin and GraphBin on the metagenomic sequencing data from simulation, two mock communities and Sharon infant fecal samples.
Conclusions
Our findings demonstrate METAMVGL outstandingly improves the short contig binning and outperforms the other existing contig binning tools on the metagenomic sequencing data from simulation, mock communities and infant fecal samples.
Background
During longterm genetic evolution, animals, including humans, have formed complex ecosystems of symbiotic relationships with diverse microbes. The gut microbiome is a community with the highest microbial density in the human body, including thousands of microbial species mixed in varying proportions and constituting a dynamic system. Most gut microbes are difficult to be isolated and cultured in vitro. Metagenomic sequencing is designed to directly sequence a mixture of microbes and explore microbial compositions and abundances by data postprocessing.
Due to the paucity of highquality microbial reference genomes, current pipelines commonly target single genes or species using speciesspecific markers [1, 2]. But novel microbes would be lost by the alignmentbased approaches. Metagenome assembly is a promising strategy to explore the novel species by concatenating the shortreads into contigs. But these contigs could be fragmented and can be only regarded as pieces of the target genomes. Contig binning algorithms provide a supplement to genome assembly that group the contigs into clusters to represent the complete microbial genomes. This strategy has been widely adopted to explore the novel microbes from the human gut metagenomic sequencing data [3,4,5,6,7,8,9,10].
Many stateoftheart contig binning algorithms have been developed by considering contig nucleotide compositions (tetranucleotide frequencies (TNF), kmer frequencies) and read depths. MaxBin2 [11] uses Expectation–Maximization algorithm to maximize the probability of a contig belonging to the local cluster centers using TNF and read depth. These two types of information are also used in MetaBAT2 [12] to calculate the contig similarities. MetaBAT2 constructs a graph using contigs as vertices and their similarities as the edges’ weights, which is further partitioned into subgraphs by applying a modified label propagation algorithm. CONCOCT [13] applies Gaussian mixture models for contig clustering based on kmer frequencies and read depths across multiple samples. Besides considering TNF and read depths, MyCC [14] aggregates the contigs with complementary marker genes by affinity propagation; SolidBin [15] develops a spectral clustering algorithm using taxonomy alignments as mustlinks between contigs; BMC3C [16] applies codon usage in an ensemble clustering algorithm. All these methods are helpless in labeling short contigs (< 1 kb), because the limited number of nucleotides might lead to unstable TNF distributions and read depths. We observed a majority of the contigs (89.55%, Additional file 1: metaSPAdes assembly of Sharon dataset) in the assembly graph were shorter than 1 kb, which would be dropped by most of the existing binning algorithms.
To rescue those short contigs, Mallawaarachchi et al. developed GraphBin [17] to label the short contigs and correct the potential binning errors by employing label propagation on the assembly graph. In principle, the assembly graph should include s disconnected subgraphs, each representing one species. In practice, the subgraphs could be linked by the repeat sequences and some contigs are isolated from the main graph (the largest graph component) due to sequencing errors, imbalanced reads coverage, named dead ends. The performance of label propagation heavily relies on the number of edges and label density in the graph. The labels of short contigs would be significantly affected by the dead ends in two ways: (1) contigs are failed to be labeled if the dead end contains no label before propagation (Fig. 1 dead end 1); (2) erroneous labeling happens if only a small proportion of nodes are labeled in the dead end (Fig. 1 dead end 2).
Here we present METAMVGL (Fig. 2), a multiview graphbased metagenomic contig binning algorithm to address the abovementioned issues. METAMVGL not only considers the contig sequence overlaps from the assembly graph but also involves the pairedend graph (PE graph), representing the shared pairedend reads between two contigs. The two graphs are integrated together by autoweighting, where the weights together with the predicted contig labels are updated in a uniform framework [18] (Methods). Figure 1 gives a proofofconcept example on the simulated data, where the pairedend reads connect the two dead ends (dead end 1 and 2) to the main graph. Our experiments indicate METAMVGL substantially improves the binning performance of the stateoftheart algorithms, including MaxBin2, MetaBAT2, MyCC, CONCOCT, SolidBin and GraphBin in the simulated, mock and Sharon datasets (Figs. 3, 4, Additional files 3, 4, 5, 6: Figures). Comparing with assembly graph, the combined graph adds up to 8942.37% vertices and 15,114.06% edges to the main graph (Additional file 2: The assembly graph from MEGAHIT for Sharon dataset).
Methods
Figure 2 demonstrates the workflow of METAMVGL, which consists of two steps. In step 1, METAMVGL constructs the assembly graph and PE graph with contig labels generated by the existing binning tools. In step 2, we remove the ambiguous labels of contigs if their neighbors are labeled as belonging to the other binning groups. The two graphs are integrated based on the automatic weights and the unlabeled contigs will be further predicted by label propagation on this combined graph. Finally, METAMVGL removes the ambiguous labels and generates the binning results.
Step 1: Preprocessing
Construct assembly graph
We define the assembly graph as \(\mathcal {G}_1(\mathcal {V}, \mathcal {E}_1)\), where the vertex \(v_i \in \mathcal {V}\) represents the contig i, and an edge \(e_{i, j} \in \mathcal {E}_1\) exists if \(v_i\) and \(v_j\) are connected in the assembly graph and with \(k1\) mer (continuous nucleotide of length \(k1\)) overlap. In principle, the assembly graph should include s unconnected subgraphs, each representing one species and we can easily recognize contig binning groups. In practice, the subgraphs could be linked due to the interspecies repeat sequences and complicated by the sequencing errors and imbalanced genomic coverage. Commonly the assembly graph includes a main graph and several dead ends. Figure 2 illustrates an assembly graph with two dead ends (vertices 11 and 12). METAMVGL uses the assembly graph from metaSPAdes [19] or MEGAHIT [20]. The original assembly graph of metaSPAdes is a unitigbased graph, where each vertex represents a unitig. The contigs are sets of unitigs after resolving short repeats. Hence we convert the unitigbased graph to contigbased graph by adding the edge \(e_{i,j}\), if at least two unitigs connect to each other and belong to \(v_i\) and \(v_j\), respectively. MEGAHIT would not provide the assembly graph directly, so METAMVGL uses contig2fastg module in megahit_toolkit to generate the graph in fastg format.
Construct PE graph
In order to deal with the dead ends in assembly graph, METAMVGL constructs the PE graph by aligning pairedend reads to the contigs by BWAMEM[21]. For every two contigs \(v_i, v_j\), we maintain a read name set \(\mathcal {RS}_{i,j}\) to keep the names of read pairs, where the forward and reverse reads are aligned to the two contigs, respectively. The library insert size \(\mathcal {IS}\) is calculated based on the uniquely aligned pairedend reads in the same contigs. To alleviate the influence of chimeric reads, METAMVGL links \(v_i\) and \(v_j\) if at least half of the reads in \(\mathcal {RS}_{i,j}\) come from the two stretches with length \(\mathcal {IS}\) in \(v_i\) and \(v_j\), respectively [22].
We denote the PE graph as \(\mathcal {G}_2(\mathcal {V}, \mathcal {E}_2)\), where \(\mathcal {V}\) represents contigs, and \(\mathcal {E}_2\) the edges linked by the pairedend reads (PE links). According to our observation, PE graph is complementary to the assembly graph to some extent, because the edges in assembly graph (\(\mathcal {E}_1\)) only capture the overlaps between contigs, while the PE graph edges (\(\mathcal {E}_2\)) link the contigs with gaps. Figure 2 illustrates how dead ends of assembly graph can be linked to the main graph using PE links.
Initial binning
The contigs’ initial labels are generated by any existing binning tools. In experiments, we evaluated the performance of METAMVGL with the initial labels from MaxBin2, MetaBAT2, MyCC, CONCOCT and SolidBin in SolidBinSFS mode. We used the default parameters for these algorithms except the MetaBAT2, where the minimum contig length was set to 1.5kb to label more contigs.
Step 2: Autoweighted multiview binning
METAMVGL applies a multiview label propagation algorithm [18] to learn the weights of assembly and PE graphs automatically and predict the unlabeled contigs in a uniformed framework. We remove the ambiguous labels for two times before and after label propagation.
Remove ambiguous labels
The initial contig labels could be incorrect especially for the ones from the interspecies repeat sequences and their influence would be amplified in label propagation. METAMVGL computes the distance between two vertices as the length of shortest path between them. Let \(\mathcal {CLV}(v)\) be the set of labels from vertex v’s closest labeled neighbors in graph \(\mathcal {G}\) and v’s label is ambiguous if \(\mathcal {CLV}(v)\) contains a label that is different from v [17]. Let \(\mathcal {VA}(\mathcal {G})\) denote the set including all the vertices with ambiguous labels in graph \(\mathcal {G}\), and we remove the labels in \(\mathcal {VA}(\mathcal {G})\). In Fig. 2, the closest labeled vertices of \(v_6\) are \(\{v_4, v_7\}\). Because \(v_6\) and \(v_7\) have different labels, the \(v_6\)’s label is marked ambiguous. Algorithm 1 shows the procedure to remove ambiguous labels. As shown in Fig. 2, we applied Algorithm 1 on the assembly and combined graphs before and after label propagation, respectively. We only use assembly graph to mark ambiguous labels in the preprocessing step for keeping more labels before propagation.
Autoweighted multiview binning algorithm
Assume l contigs are initially labeled with s groups, denoted as \(Y_l=[y_1, y_2,\ldots , y_l]^T \in \mathbb {R}^{l \times s}\), where \(y_{ij} \in \{0, 1\}\), and \(y_{ij}=1\) indicates the contig \(v_i\) is labeled from group j. We define a indicator matrix \(F=[F_l;F_u] \in \mathbb {R}^{n \times s}\), where \(F_l=Y_l\) and \(F_u=[f_{l+1}, f_{l+2},\ldots , f_n]^T\) are labels to be inferred. Let \(D_i, W_i \in \mathbb {R}^{n \times n}\) denote the degree and adjacent matrices of \(\mathcal {G}_i\) (\(i\in \{1,2\}\)), respectively. The normalized Laplacian matrix of \(\mathcal {G}_i\) has the formulation \(\mathcal {L}_i = D_i^{1/2}(D_iW_i)D_i^{1/2}\). According to [18], the inference of \(F_u\) by label propagation can be modeled as the following optimization problem:
where \(\mathcal {TR}(\cdot )\) computes the trace of a matrix. The optimization problem is converted to
where \(\mathcal {L}=\sum _{i=1}^{2}\alpha _i\mathcal {L}_i\). \(\alpha _i\) is the weight of \(\mathcal {G}_i\), with initial values of 1/2. We partition \(\mathcal {L}\) from \((l+1)\)th row and column into four blocks as \([\mathcal {L}_{ll}, \mathcal {L}_{lu};\mathcal {L}_{ul},\mathcal {L}_{uu}]\). \(F_u\) and \(\alpha _i\) can be updated alternatively until convergence by the following equations [18]:
Equation 3 can be considered as performing label propagation in the combined graph with iteratively updated weight \(\alpha _i\), hence \(\alpha _i\) implies the confidence of each graph. After obtaining \(F_u\), we infer the labels of all the contigs by \(l_i = \arg \mathrm {max}_{j}\ F_{ij},i\in \{1, 2, \ldots , n\}, j\in \{1,2,\ldots ,s\}\). Algorithm 2 shows the procedure of autoweighted multiview binning, and Fig. 2 is an illustration of this algorithm.
Datasets
Simulated datasets
We simulated metagenomic sequencing data for a mixture of three strains with low, medium and high abundances. The components are:

Acinetobacter baumannii: 0.90%,

Streptococcus agalactiae: 9.01%,

Streptococcus mutans: 90.09%.
We downloaded the complete reference genomes of the three strains from the NCBI Nucleotide Database (Taxonomy ID: 400667, 208435, 210007). CAMISIM [23] generated shortreads for the three strains with corresponding abundances. Five simulated datasets were generated with read depths as 30x (SIM_30x), 50x (SIM_50x), 70x (SIM_70x), 90x (SIM_90x) and 110x (SIM_110x).
Mock datasets
We evaluated the performance of METAMVGL on the metagenomic sequencing from two mock communities:

BMock12 refers to the metagenomic sequencing for a mock community with 12 bacterial strains sequenced by Illumina HiSeq 2500 [24] (NCBI acc. no. SRX4901583). It contains 426.8 million 150 bp shortreads with a total size of 64.4G bases.

SYNTH64 is a metagenomic sequencing dataset for a synthetic community with 64 diverse bacterial and archaea species [25] (NCBI acc. no. SRX200676), sequenced by Illumina HiSeq 2000 with read length 101bp and total size 11.1G bases.
Real dataset
Sharon dataset [26] (NCBI acc. no. SRX144807) contains the metagenomic sequencing data of infant fecal samples from 18 time points, sequenced by Illumina HiSeq 2000 with a total of 274.4 million 100 bp shortreads. We combined all the 18 datasets for coassembly and referred them as Sharon.
Evaluation criteria
We annotated the potential species the contigs came from as ground truth to compare METAMVGL with the other tools. For the simulated and mock datasets, we aligned the contigs to the available reference genomes and selected the ones with unique alignments. For the Sharon dataset, we used Kraken2 [27] to annotate the contigs according to kmer similarities with the species from the buildin database.
Assume there are s ground truth species, and the binning result have k groups. To evaluate the binning result, we define the assessment matrix \([n_{i,j}]^{(k+1) \times (s+1)}\), where \(n_{i,j}\) represents the number of contigs in ith binning group that are annotated jth ground truth species. The \((k+1)\)th row denotes unbinned contigs. The \((s+1)\)th column indicates contigs without ground truth annotations. We applied (1) Precision, (2) Recall, (3) F1Score and (4) Adjusted Rand Index (ARI) to evaluate the performance of binning algorithms. Let \(N = \sum \nolimits _{i=1}^{k} \sum \nolimits _{j=1}^{s} n_{i,j}\) be the number of contigs; the four metrics are calculated as follows:

1.
Precision \(=\frac{1}{N}\sum \limits _{i=1}^{k} \underset{j \le s}{max}(n_{i,j})\),

2.
Recall \(=\frac{1}{N + \sum \limits _{j=1}^{s} n_{k+1,j}}\sum \limits _{j=1}^{s} \underset{i \le k}{max}(n_{i,j})\),

3.
F1Score \(=\frac{2\times Precision \times Recall}{Precision + Recall}\),

4.
ARI \(=\frac{\sum \nolimits _{i=1}^{k} \sum \nolimits _{j=1}^{s} \left( {\begin{array}{c}n_{i,j}\\ 2\end{array}}\right)  t}{\frac{1}{2}\bigg (\sum \nolimits _{i=1}^{k} \left( {\begin{array}{c}\sum \nolimits _{j=1}^{s} n_{i,j}\\ 2\end{array}}\right) + \sum \nolimits _{j=1}^{s} \left( {\begin{array}{c}\sum \nolimits _{i=1}^{k} n_{i,j}\\ 2\end{array}}\right) \bigg )  t}\),
where \(t = \frac{1}{\left( {\begin{array}{c}N\\ 2\end{array}}\right) }\sum \nolimits _{i=1}^{k} \left( {\begin{array}{c}\sum \nolimits _{j=1}^{s} n_{i,j}\\ 2\end{array}}\right) \sum \nolimits _{j=1}^{s} \left( {\begin{array}{c}\sum \nolimits _{i=1}^{k} n_{i,j}\\ 2\end{array}}\right)\).
Results
METAMVGL was compared to six binning tools, MaxBin2, MetaBAT2, MyCC, CONCOCT, SolidBin in SolidBinSFS mode and GraphBin. We analyzed their binning results on the five simulated datasets with various read depths, two mock communities (BMock12 and SYNTH64) and a real metagenomic sequencing dataset (Sharon dataset).
Evaluation on the simulated datasets
Figure 3 shows the binning results of the simulated datasets. The contigs and assembly graph were generated by MEGAHIT (Fig. 3a–d) and metaSPAdes (Fig. 3e–h). MaxBin2 was applied as the initial binning tool for GraphBin and METAMVGL.
All the three binning algorithms (MaxBin2, GraphBin and METAMVGL) yielded extremely high precision and ARI (Fig. 3a, d, e, h), due to the low complexity of the simulated datasets. Because of considering assembly and PE graphs jointly, METAMVGL labeled more contigs than GraphBin and MaxBin2 across various sequence depths, as shown in Fig. 3b, f. We also found both Recall and F1Score were improved as read depth became higher until SIM_70x (Fig. 3b, c, f, g). This observation was analog to the results from CAMISIM [23], suggesting too high read depth would introduce assembly noise even it might help in detecting the lowabundance microbes.
Evaluation on the mock communities
We illustrate the binning results for two mock communities with initial binning tool of MaxBin2 in Fig. 4a, b, d, e. In general, the graphbased methods (METAMVGL and GraphBin) were better than MaxBin2, but their performance would be influenced by the assembly graph. We observed the recalls could be significantly improved using the assembly graph from metaSPAdes (Fig. 4d, e), but the elevation became unobvious by the one from MEGAHIT (Fig. 4a, b). This was probably because metaSPAdes could generate more accurate and complete assembly graph than MEGAHIT. In the mock communities, METAMVGL was just slightly better than GraphBin, suggesting the PE graph was largely overlapped with the assembly graph (Additional file 2: BMock12 and SYNTH64 datasets). This observation only occurred if a perfect assembly graph was generated due to a low microbial complexity in the community. The results for the other initial binning tools (MetaBAT2, MyCC, CONCOCT and SolidBin) could be found in Additional files 3, 4, 5, 6: Figures, which were akin to the observations from MaxBin2.
Evaluation on Sharon datasets
Figure 4c, f describe the binning results of MaxBin2, GraphBin and METAMVGL in Sharon dataset. METAMVGL substantially improved the recalls comparing with GraphBin (2.12 and 2.46 times on the assembly graph from MEGAHIT and metaSPAdes, respectively) and MaxBin2 (2.29 and 3.24 times on the assembly graph from MEGAHIT and metaSPAdes, respectively). METAMVGL showed the highest precision on the assembly by metaSPAdes (Fig. 4f). All the three binning algorithms were comparable in ARI.
The outstanding recalls from METAMVGL validate the capability of PE graph to connect the dead ends to the main graph when the assembly graph is incomplete in the complex microbial community. We observed that MEGAHIT produced very fragmented assembly graph in the Sharon dataset, in which the main graph only had 59 vertices with 64 edges, while a total of 15,660 vertices existed in the whole graph (Additional file 2: MEGAHIT assembly of Sharon dataset). The fragmented assembly graph was also mentioned as a limitation of GraphBin [17]. With PE links, METAMVGL yielded 5335 vertices and 9737 edges in the main graph (Additional file 2: MEGAHIT assembly of Sharon dataset), rescuing a large number of unlabeled contigs from the dead ends (Fig. 4c). Although the assembly graph was more complete (23.69% vertices in the main graph) from metaSPAdes, the PE graph still added 28.97% edges to the main graph (Additional file 2: metaSPAdes assembly of Sharon dataset) and improved the recall substantially.
Discussion
De novo assembly together with contig binning offer a practical way to explore the novel microbes from metagenomic sequencing. But the current binning algorithms work stably merely on long contigs; the shorter ones are commonly neglected in the subsequent analysis. We observed a large proportion of contigs were shorter than 1 kb, which resulted in low completeness of the binning groups. A recent study [17] suggests the short contigs could be rescued from the assembly graph by considering their connections with the labeled ones. Assembly graph is accurate, but its connectivity relies heavily on the complexity of microbial community. Extremely low or high read depth, sequencing errors and imbalanced coverage could generate considerable dead ends, which would introduce both missing labels and labeling errors (Fig. 1).
In experiments, we observed a slightly lower ARI of METAMVGL comparing with MaxBin2 and GraphBin when read depth was low (Fig. 3d, h). It might because METAMVGL could only retrieve very few and lowconfidence PE links. First, the label propagation would perform poorly on the contigs with low read depth. Comparing with GraphBin, METAMVGL included more edges in the graph, but these edges were sparse and cannot guarantee the good performance of label propagation, e.g. the unlabeled contigs might only have one neighbor. Second, the label refinement after label propagation could remove a majority of erroneous labels generated by METAMVGL based on our experience. Due to the paucity of edges, this step also performed inefficiently. Third, the quality of the contigs assembled from the sequencing data with low read depth was poor, making difficulties in aligning pairedend reads correctly.
In this paper, we developed METAMVGL, a multiview graphbased contig binning algorithm to integrate both assembly and PE graphs to label short contigs and correct initial labeling errors. PE graph could link the dead ends to the main graph and increase the graph connectivity. METAMVGL automatically weights the two graphs and performs label propagation to label the short contigs. In experiments, we observed METAMVGL could substantially improve the recalls without loosing any precision comparing with the existing contig binning tools, especially for the metagenomic data from the complex microbial community (Fig. 4c, f). We also evaluated METAMVGL: 1. on the assembly graphs from metaSPAdes and MEGAHIT; and 2. using the initial binning labels from different tools. All these results support METAMVGL outperform GraphBin in different experimental configurations. On average, METAMVGL could finish the contig binning in 3.38 min and requires 2.81 Gb RAM to store the two graphs and perform label propagation. It requires a little bit more computational resources than GraphBin due to the analysis of more complex and complete graph (Additional file 7). Sometimes, we found the combined graph was still incomplete after incorporating both assembly and PE graphs and there still required to consider other information to reveal contig longrange connectivity from various longfragment sequencing (PacBio and Oxford Nanopore sequencing) or linkedread sequencing (Tellseq and stLFR sequencing) technologies.
Conclusion
Metagenomic sequencing has been proved as an efficient technology to explore and recognize the novel microbes in the environmental and human fecal samples. Due to the scarcity of the reference genomes, the genomes of novel species could be obtained by de novo assembly. Because only fragmented contigs can be assembled from the mainstream shortread sequencing technologies, the interests rise quickly in developing efficient contig binning algorithms. But most of the available algorithms can only handle long contigs based on their sequence contexts and read depths. In this study, we developed METAMVGL, a multiview graphbased contig binning algorithm to integrate both assembly and PE graphs to label short contigs and correct initial labeling errors. METAMVGL could weight the two graphs automatically and connect the dead ends to the main graph efficiently. Our experiments proved it can significantly improve the recalls without loosing any precision comparing with the existing contig binning tools on the metagenomic sequencing data from simulation, mock communities and infant fecal samples. We believe METAMVGL would attract more interests of the fastgrowing metagenomic research field and pave the way to future understanding the microbial genome dark matter.
Availability of data and materials
The source code of METAMVGL is publicly available at https://github.com/ZhangZhenmiao/METAMVGL. The Illumina shortreads of BMock12, SYNTH64 and Sharon datasets are available in NCBI Sequence Read Archive (SRA); the accession numbers are SRX4901583, SRX200676 and SRX144807, respectively.
About this supplement
This article has been published as part of BMC Bioinformatics Volume 22 Supplement 10 2021: Selected articles from the 19th Asia Pacific Bioinformatics Conference (APBC 2021): bioinformatics. The full contents of the supplement are available at https://bmcbioinformatics.biomedcentral.com/articles/supplements/volume22supplement10.
Abbreviations
 METAMVGL:

Metagenomic contig binning using multiview graph learning
 PE graph:

Pairedend graph
 PE link:

Pairedend link
 TNF:

Tetranucleotide frequencies
 ARI:

Adjusted rand index
 RAM:

Randomaccess memory
References
Li J, Jia H, Cai X, Zhong H, Feng Q, Sunagawa S, Arumugam M, Kultima JR, Prifti E, Nielsen T, et al. An integrated catalog of reference genes in the human gut microbiome. Nat Biotechnol. 2014;32(8):834–41.
Truong DT, Franzosa EA, Tickle TL, Scholz M, Weingart G, Pasolli E, Tett A, Huttenhower C, Segata N. Metaphlan2 for enhanced metagenomic taxonomic profiling. Nat Methods. 2015;12(10):902–3.
Almeida A, Mitchell AL, Boland M, Forster SC, Gloor GB, Tarkowska A, Lawley TD, Finn RD. A new genomic blueprint of the human gut microbiota. Nature. 2019;568(7753):499–504.
Almeida A, Nayfach S, Boland M, Strozzi F, Beracochea M, Shi ZJ, Pollard KS, Sakharova E, Parks DH, Hugenholtz P, et al. A unified catalog of 204,938 reference genomes from the human gut microbiome. Nat Biotechnol. 2020; 1–10.
Poyet M, Groussin M, Gibbons S, AvilaPacheco J, Jiang X, Kearney S, Perrotta A, Berdy B, Zhao S, Lieberman T, et al. A library of human gut bacterial isolates paired with longitudinal multiomics data enables mechanistic microbiome research. Nat Med. 2019;25(9):1442–52.
Nayfach S, Shi ZJ, Seshadri R, Pollard KS, Kyrpides NC. New insights from uncultivated genomes of the global human gut microbiome. Nature. 2019;568(7753):505–10.
Zou Y, Xue W, Luo G, Deng Z, Qin P, Guo R, Sun H, Xia Y, Liang S, Dai Y, et al. 1520 reference genomes from cultivated human gut bacteria enable functional microbiome analyses. Nat Biotechnol. 2019;37(2):179–85.
Forster SC, Kumar N, Anonye BO, Almeida A, Viciani E, Stares MD, Dunn M, Mkandawire TT, Zhu A, Shao Y, et al. A human gut bacterial genome and culture collection for improved metagenomic analyses. Nat Biotechnol. 2019;37(2):186–92.
Consortium HMJRS, et al. A catalog of reference genomes from the human microbiome. Science. 2010;328(5981):994–9.
Pasolli E, Asnicar F, Manara S, Zolfo M, Karcher N, Armanini F, Beghini F, Manghi P, Tett A, Ghensi P, et al. Extensive unexplored human microbiome diversity revealed by over 150,000 genomes from metagenomes spanning age, geography, and lifestyle. Cell. 2019;176(3):649–62.
Wu YW, Simmons BA, Singer SW. Maxbin 2.0: an automated binning algorithm to recover genomes from multiple metagenomic datasets. Bioinformatics. 2016;32(4):605–7.
Kang DD, Li F, Kirton E, Thomas A, Egan R, An H, Wang Z. Metabat 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ. 2019;7:7359.
Alneberg J, Bjarnason BS, De Bruijn I, Schirmer M, Quick J, Ijaz UZ, Lahti L, Loman NJ, Andersson AF, Quince C. Binning metagenomic contigs by coverage and composition. Nat Methods. 2014;11(11):1144–6.
Lin HH, Liao YC. Accurate binning of metagenomic contigs via automated clustering sequences using information of genomic signatures and marker genes. Sci Rep. 2016;6:24175.
Wang Z, Wang Z, Lu YY, Sun F, Zhu S. Solidbin: improving metagenome binning with semisupervised normalized cut. Bioinformatics. 2019;35(21):4229–38.
Yu G, Jiang Y, Wang J, Zhang H, Luo H. Bmc3c: binning metagenomic contigs using codon usage, sequence composition and read coverage. Bioinformatics. 2018;34(24):4172–9.
Mallawaarachchi V, Wickramarachchi A, Lin Y. Graphbin: refined binning of metagenomic contigs using assembly graphs. Bioinformatics. 2020;36(11):3307–13.
Nie F, Li J, Li X, et al. Parameterfree autoweighted multiple graph learning: a framework for multiview clustering and semisupervised classification. In: IJCAI, 2016; p. 1881–7.
Nurk S, Meleshko D, Korobeynikov A, Pevzner PA. metaspades: a new versatile metagenomic assembler. Genome Res. 2017;27(5):824–34.
Li D, Liu CM, Luo R, Sadakane K, Lam TW. Megahit: an ultrafast singlenode solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics. 2015;31(10):1674–6.
Li H. Aligning sequence reads, clone sequences and assembly contigs with bwamem. arXiv preprint arXiv:1303.3997 (2013).
Bishara A, Moss EL, Kolmogorov M, Parada AE, Weng Z, Sidow A, Dekas AE, Batzoglou S, Bhatt AS. Highquality genome sequences of uncultured microbes by assembly of read clouds. Nat Biotechnol. 2018;36(11):1067–75.
Fritz A, Hofmann P, Majda S, Dahms E, Dröge J, Fiedler J, Lesker TR, Belmann P, DeMaere MZ, Darling AE, et al. Camisim: simulating metagenomes and microbial communities. Microbiome. 2019;7(1):1–12.
Sevim V, Lee J, Egan R, Clum A, Hundley H, Lee J, Everroad RC, Detweiler AM, Bebout BM, PettRidge J, et al. Shotgun metagenome data of a defined mock community using oxford nanopore, pacbio and illumina technologies. Sci data. 2019;6(1):1–9.
Shakya M, Quince C, Campbell JH, Yang ZK, Schadt CW, Podar M. Comparative metagenomic and rrna microbial diversity characterization using archaeal and bacterial synthetic communities. Environ Microbiol. 2013;15(6):1882–99.
Sharon I, Morowitz MJ, Thomas BC, Costello EK, Relman DA, Banfield JF. Time series community genomics analysis reveals rapid shifts in bacterial species, strains, and phage during infant gut colonization. Genome Res. 2013;23(1):111–20.
Wood DE, Lu J, Langmead B. Improved metagenomic analysis with kraken 2. Genome Biol. 2019;20(1):257.
Funding
Publication costs are funded by Research Grant Council Early Career Scheme (HKBU 22201419). This work is supported by Research Grant Council Early Career Scheme (HKBU 22201419), IRCMS HKBU (No. IRCMS/1920/D02), Guangdong Basic and Applied Basic Research Foundation (No. 2019A1515011046). The funders did not play any role in the design of the study, the collection, analysis, and interpretation of data, or in writing of the manuscript.
Author information
Authors and Affiliations
Contributions
LZ conceived the study. ZMZ implemented METAMVGL. ZMZ and LZ analyzed the results. ZMZ and LZ wrote the paper. Both authors have read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors have no conflicts of interest.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Additional file 1.
The statistics of assemblies from MEGAHIT andmetaSPAdes for the simulated, BMock12, SYNTH64 and Sharon datasets.
Additional file 2.
The statistics of the largest component (main graph) inthe assembly graph, PE graph, and combined graph for the simulated,BMock12, SYNTH64 and Sharon datasets.
Additional file 3.
The performance of MetaBAT2, GraphBin andMETAMVGL on the BMock12, SYNTH64 and Sharon datasets: (a) and(d) for BMock12 dataset; (b) and (e) for SYNTH64 dataset; (c) and (f)for Sharon dataset. MEGAHIT and metaSPAdes are used to generate theassembly graphs. The initial binning tool is MetaBAT2.
Additional file 4.
The performance of MyCC, GraphBin and METAMVGLon the BMock12, SYNTH64 and Sharon datasets: (a) and (d) forBMock12 dataset; (b) and (e) for SYNTH64 dataset; (c) and (f) forSharon dataset. MEGAHIT and metaSPAdes are used to generate theassembly graphs. The initial binning tool is MyCC.
Additional file 5.
The performance of CONCOCT, GraphBin andMETAMVGL on the BMock12, SYNTH64 and Sharon datasets: (a) and(d) for BMock12 dataset; (b) and (e) for SYNTH64 dataset; (c) and (f)for Sharon dataset. MEGAHIT and metaSPAdes are used to generate theassembly graphs. The initial binning tool is CONCOCT.
Additional file 6.
The performance of SolidBin, GraphBin and METAMVGLon the BMock12, SYNTH64 and Sharon datasets: (a) and (d) forBMock12 dataset; (b) and (e) for SYNTH64 dataset; (c) and (f) forSharon dataset. MEGAHIT and metaSPAdes are used to generate theassembly graphs. The initial binning tool is SolidBin.
Additional file 7.
Running time and memory usage of GraphBin andMETAMVGL.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Zhang, Z., Zhang, L. METAMVGL: a multiview graphbased metagenomic contig binning algorithm by integrating assembly and pairedend graphs. BMC Bioinformatics 22, 378 (2021). https://doi.org/10.1186/s12859021042844
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12859021042844
Keywords
 Contig binning
 Assembly graph
 Pairedend graph
 Dead end
 Multiview label propagation