Next Article in Journal
A Similarity Regression Fusion Model for Integrating Multi-Omics Data to Identify Cancer Subtypes
Previous Article in Journal
The Case of X and Y Localization of Nucleolus Organizer Regions (NORs) in Tragulus javanicus (Cetartiodactyla, Mammalia)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

CoreProbe: A Novel Algorithm for Estimating Relative Abundance Based on Metagenomic Reads

1
School of Mathematics and Physics, University of Science and Technology Beijing, Beijing 100083, China
2
Sinotech Genomics, Shanghai 200120, China
3
Department of Medicine, Stanford University School of Medicine, 269 Campus Dr., Stanford, CA 94305, USA
*
Authors to whom correspondence should be addressed.
Submission received: 11 April 2018 / Revised: 12 June 2018 / Accepted: 13 June 2018 / Published: 20 June 2018
(This article belongs to the Section Technologies and Resources for Genetics)

Abstract

:
With the rapid development of high-throughput sequencing technology, the analysis of metagenomic sequencing data and the accurate and efficient estimation of relative microbial abundance have become important ways to explore the microbial composition and function of microbes. In addition, the accuracy and efficiency of the relative microbial abundance estimation are closely related to the algorithm and the selection of the reference sequence for sequence alignment. We introduced the microbial core genome as the reference sequence for potential microbes in a metagenomic sample, and we constructed a finite mixture and latent Dirichlet models and used the Gibbs sampling algorithm to estimate the relative abundance of microorganisms. The simulation results showed that our approach can improve the efficiency while maintaining high accuracy and is more suitable for high-throughput metagenomic data. The new approach was implemented in our CoreProbe package which provides a pipeline for an accurate and efficient estimation of the relative abundance of microbes in a community. This tool is available free of charge from the CoreProbe’s website: Access the Docker image with the following instruction: sudo docker pull panhongfei/coreprobe:1.0.

Graphical Abstract

1. Introduction

Microbial organisms are ubiquitous in virtually all the natural environments of the earth’s biosphere. They play integral and unique roles in ecosystems [1], they are involved in the biogeochemical cycling of the earth [2], and they have a great impact on human health. There is a dynamic equilibrium between the intestinal microflora, the host, and the environment. Once the structure, composition, and function of the microbiota cause alterations in the metabolites of the intestinal microbiota, there can be host diseases such as obesity [3], malnutrition [4] and diabetes [5]; intestinal flora disorders and irritable bowel syndrome [6], ulcerative colitis and Crohn’s disease [7,8] and other chronic bowel diseases, colon cancer [9,10] and gastric cancer [11]. Thus, profiling the taxonomic composition using the microbial abundances of related communities is critical for understanding the microbial ecology of the environment and for human health.
Recent innovations in metagenomic shotgun sequencing have made it possible to characterize microbial contents in uncultured samples by yielding billions of short reads from metagenomes. Many algorithms have been merged to estimate a community’s taxonomic composition by analyzing the metagenomics sequencing data. These algorithms can be classified into two categories, alignment-based or composition-based, according to their different resolutions of taxonomic binning.
The composition-based algorithms classify metagenomic reads mostly according to their k-mer frequencies. A number of unsupervised methods of this kind have been employed for clustering the reads generated from similar taxonomies. For example, TETRA [12] clusters reads based on Pearson’s correlation coefficients of z-scores, which evaluate the divergence between the observed and expected reads of tetra nucleotide frequencies. CompostBin [13] and SCIMM [14] step up by using a weighted PCA algorithm and an interpolated Markov model, respectively, to bin the reads. MetaCluster [15] involves two phases of clustering to guarantee the accuracy of the read binning. Some supervised methods have also been proposed to consider the available genomic information and to assign taxonomic labels. Various machine learning methods such as the support vector machine (SVM) classifier, naive Bayes classifier, and Gaussian kernel function are employed in Phylopythia [16], NBC [17] and TACOA, respectively [18]. In addition, RAIphy [19] assigns reads taxonomically according to each read’s k-mer index value sum based on its Relative Abundance Index model. Other composition-based methods directly explore the k-mer sequences of the reads, such as LMAT [20], Kraken [21] and CLARK [22], which assigns taxonomic labels using the lowest common ancestor (LCA) approach based on the matches of the reads’ k-mer sequences with their differently self-built k-mer databases. PhymmBL [23] and RITA [24] use hybrid information from BLAST results to improve their accuracy. MetaTopics [25] can efficiently extract the latent microbial communities by using the topic model, and it reflects the intrinsic relations or interactions among several major microbes.
Old composition-based methods such as TETRA and TACOA show sharply reduced accuracy during the abundance estimation when the lengths of the reads are less than 800 bp or the community complexity of the samples is high. Although it is claimed that newly merged composition-based methods can be applied quickly and accurately to metagenomic datasets from next-generation sequencing with read lengths of less than 300 bp, these methods still do not appear to be widely used in real metagenomic analyses because there is currently an open question about how k-mer compositions are associated with microbial phylogeny. More evidence must be discovered until the compositional characteristics of DNA sequences become solid phylogenetic signals.
In the alignment-based approaches, alignment and mapping tools, such as BLAST [26], BWA [27] and Bowtie2 [28], are used to find similarity hits in metagenomic reads to reference sequences and then assign the mapped reads to the most plausible microbial lineages. Early alignment-based methods simply estimate the microbial abundances according to the mapped counts. Most of these methods focus primarily on precisely filtering mapped reads and employing the lowest common ancestor algorithm. Methods such as MEGAN [29], CARMA [30], PaPaRa [31] and MTR [32] all have different strict mapping and filtering procedures or preprocessing procedures to achieve an accurate estimation. Then, the reference length normalization and probabilistic model are introduced into the metagenomic reads analysis. GAAS [33] also improves the microbial abundance estimates by introducing similarity weighting based on alignment e-values. GRAMMy [34] accurately estimates the microbial species abundances based on the mixture model theory and EM algorithm, while Pathoscope [35] presents a complete framework of the metagenomic composition analysis. GASic [36] and MetaMix [37] add additional probabilistic procedures such as similarity correction by Least Absolute Shrinkage and Selection Operator (LASSO) or model comparison by Markov-Chain Monte Carlo (MCMC) to improve the estimation.
Another way to perform microbial abundance estimations other than introducing the probabilistic model to alignment-based approaches is to substitute traditional complete genome sequences with more condensed and informative marker sequences as references. MetaPhlAn [38] and mOTU [39], which employ clade-specific marker genes as reference sequences, have greatly improved the speed and accuracy of analyzing extremely large metagenomic datasets. Clade-specific marker genes can be identified in almost every clade at different microbial taxonomic levels, and thus they are much more representative than 16S rRNA in metagenomic analysis. The use of marker genes as references only accounts for approximately 4% of the sequenced microbial genes, which leads to a significant conservation in alignment time and storage. However, this approach will discard large numbers of reads in the meantime, because most reads will hardly be mapped to marker genes, which is a regretful loss of information from the metagenomic datasets. Moreover, some microbial species may have a rather small percentage of marker genes because of frequent mutations.
Here, we recommend using core genomes as the reference sequences of potential species and estimating the microbial abundances in metagenomic analysis by using the probabilistic model. According to the pan-genome concepts motivated by Tettelin [40], a clade’s core-genome, which contains genes shared by all the strains within the clade, typically includes the genes responsible for the major phenotypic traits, which account for nearly 8% of the genetic repertoire. In involving clade-specific marker genes as a subset, the core-genome generally contains some genes that may be shared by different clades. The reasons for using core-genomes as references are as follows: First, the coexistence of different subpopulations of a microbial species may be a general feature in highly mixed habitats. Evidence has been discovered of both selected cultured isolates and wild uncultured populations, from marine environments [41] to human body sites [42]. The existence of multiple subpopulations or strains, the abundances of which are very likely to vary greatly, often brings extensive genomic diversity and community complexity and leads to inaccurate results when using only several complete genomes or a small percentage of marker genes as references for a clade. However, using core-genomes that consist of the genes shared by all the microbial strains that were studied as references provides us with an opportunity to estimate the relative abundance of a species with a complex substrain composition quickly and accurately. Second, variational subpopulations of a microbial species that have never been sequenced before, the strains of one species that may have high similarity with another species [43], and some mechanisms such as horizontal gene transfer and lysogeny, which currently tend to be considered to occur more frequently in natural environments than previously thought [44], are all obstacles to alignment accuracy when using the strains’ complete reference genomes. Otherwise, using core-genomes as references should include the general new strains of the species because most currently identified core genes should also be shared by those strains. Different microbial species with high similarity or exotic sequences would also be addressed because they would differ significantly in their core-genomes. Third, compared to using marker genes as clade references, the shared or redundant genes between different clades in the core-genomes form the basis for introducing probabilistic models, which create a balance between mapping quickly and precisely to the references and utilizing the information hidden in ambiguous reads.
In this paper, we introduce CoreProbe, a relative abundance estimation framework for microbes that employs microbial core-genomes as references for metagenomic analysis. In addition, CoreProbe takes advantage of the mixture model theory and describes the sequencing procedure for metagenomic reads as a generative model to accurately estimate microbial abundances with the Gibbs sampling algorithm [45,46,47]. We first tested CoreProbe both in our own simulated metagenomic read sets using MetaSim [48] and in third-party synthetic communities [49]. From these experiments, we can observe that using core-genomes as references in CoreProbe sharply outperforms the use of ordinary complete microbial genomes, whether the specific strains of the reference genomes are “in” the metagenomic datasets or not. Compared to other methods including GRAMMy [34], Pathoscope [35], MetaPhlAn [38] and Kraken [21], CoreProbe also shows improved accuracy in abundance estimations. We then analyzed 25 real metagenomic read sets from Human Microbiome Project (HMP) (https://portal.hmpdacc.org/), yielding new insights into microbiomes from different human body sites. Finally, we implemented CoreProbe in C++, and we accessed the Docker image with the following instruction: sudo docker pull panhongfei/coreprobe:1.0.

2. Methods

2.1. A Finite Mixture and Latent Dirichlet Model

To estimate the relative abundance of reference species accurately according to a metagenomic dataset, we describe the sampling and sequencing procedure of metagenomic reads as a generative model: First, in using M microbial organisms as the reference species, the metagenome M ¯ in a metagenomic dataset can be denoted as
M ¯ = i = 1 M θ i g i
where { g 1 , g 2 , , g M } G represents the reference sequences of M known species. These reference sequences can be contigs, complete genomes, pan-genomes, core-genomes, etc. ( θ 1 , θ 2 , , θ M ) θ denotes the mixture parameters of those reference sequences (or reference species). It should be noted that
i = 1 M θ i = 1
and each θ i is proportional to its relative species abundance a i and the corresponding reference sequence’s base length l i , i.e., θ i a i l i according to Xia et al. [33]. Here, those organisms are subject to a particular discrete distribution ( θ 1 , θ 2 , , θ M ) in the metagenomic generative model. Second, from metagenome M , we randomly chose reference species g i , which is subject to the multinomial probability θ i , because each read must be generated from the biological sequence of a particular species. Third, given the chosen genome g i , we randomly generated read r k . The generation of reads from the reference species g i is subject to g i ’s read-composition distribution
( φ r 1 , g i , φ r 2 , g i , , φ r K , g i ) φ i
where φ r k , g i denotes the probability of generating a particular read r k from genome g i , i.e., p { r e a d _ k | G = g i } . Here, we assume K , the total number of different reads generated from metagenome M ¯ . We denote the set of the read-composition distributions for all the reference species as { φ g 1 , φ g 2 , , φ g M } Φ .
Based on the assumptions in the metagenomic generative model, we can easily induce the formula of the relative abundances for reference species when their reference sequence lengths { l 1 , l 2 , , l M } and mixture parameters ( θ 1 , θ 2 , , θ M ) are known. Noting that i = 1 M θ i = 1 and θ i a i l i , the relative abundance formula of the reference species under the metagenomic generative model is as follows:
a i = θ i l i s = 1 M θ s l s
The procedure in the metagenomic generative model will be repeated N times to obtain N metagenomic reads. We can then use the obtained reads to infer the mixture parameters and then the relative abundances of the reference species based on the generative model. As we can observe, both sampling procedures for the second and third steps are subject to multinomial distributions, i.e., g i M u l t ( θ ) and r k M u l t ( φ g i ) , respectively. For the sake of the calculation, we followed the suggestion of Pritchard [45] and used the Dirichlet distribution as the prior distribution of species mixture parameters in the metagenome because of the Dirichlet-multinomial conjugacy. However, we do not assume prior distributions for the read-composition distributions Φ of those reference species, because we can approximate their full conditional probabilities during parameter inference via read-to-reference-sequence alignment results. The above generative model for metagenomic reads is shown as the pseudo-code in Algorithm 1.
Algorithm 1 (A finite mixture and latent Dirichlet model for metagenomics).
Require the hyperparameter α , the total number of reads N , the species set G , the read-composition distributions Φ
Ensure
the read dataset R
sample the species mixture parameters θ D i r ( α ) for a metagenome
repeat
1. sample species g i M u l t ( θ )
2. sample read r k M u l t ( φ g i )
untilthe total number of metagenomic reads N is reached
returnthe read dataset { r 1 , r 2 , , r N } R

2.2. Mixture Parameter Inference and Gibbs Sampling

Ideally, we would need exact knowledge about the occurrence counts of the read-origin species for each read in a metagenomic dataset to estimate the mixture parameters θ accurately. However, this task is basically impossible to perform through sequence alignment due to the inherent ambiguity of relatively short NGS reads and the complex microbial communities, which generally consist of species with similar reference genomes. Hence, we employ a parameter inference procedure involving Gibbs sampling to infer the mixture parameter.
In this section, we introduce hidden variables ( z 1 , z 2 , z N ) z for the first time, where z i is an index denoting the reference species that generates the corresponding read r i . We develop an approximate inference algorithm of Gibbs sampling to emulate the probability distribution of z i given the observations of metagenomic reads R conditioned on prior probability α and the read-composition distributions Φ of those reference species; i.e., p { z | R ; α , Φ } . We can then statistically estimate the counts of the reference species using the samples of p { z | R ; α , Φ } after the burn-in period of Gibbs sampling, and herewith obtain the estimation of parameters θ and then the relative species abundances α . Specifically, Gibbs sampling generates an instance of each dimension z i of z in turn, subject to their full conditional p { z | z ¬ i , R ; α , Φ } , where z ¬ i denotes all other dimensions of z except z i . It can be shown [50] that under this condition, the sequence of samples { z 1 , z 2 , } constitutes a Markov chain whose stationary distribution is p { z | R ; α , Φ } . During the real application of metagenomic data, we can extract a certain number of samples of z after the burn-in period, and we can calculate the average counts of reference species to infer the parameters.
Next, we derive the full conditional p { z | R ; α , Φ } for Gibbs sampling. First, we calculate the probability p { z | α } . Starting with the probability z as conditioned to the species mixture probability θ , and noting that the hidden species indices are generated as multinomial trials, we have
p { z | θ } = i = 1 M θ i n i
where n i refers to the number of reads whose corresponding species index is i . Noting that we assume p { θ | α } is subject to Dirichlet distribution, we have
p { θ | α } d θ = 1 ( α ) i = 1 M θ i α i 1 d θ = 1
where
( α ) = k = 1 M Γ ( α k ) Γ ( k = 1 M α k )
Then, by using the above two formulas, and integrating out θ in the following conditional probability formula, we obtain
p { z | α } = p { z | θ } p { θ | α } d θ = 1 ( α ) i = 1 M θ i n i + α i 1 d θ = ( n + α ) ( α )
Second, we obtain the probability of reads conditioned on the hidden variables z with the knowledge of read-composition distribution set Φ as
p { r | z ; Φ } = i = 1 N p { r i | z i , Φ } = i = 1 N ϕ r i , g i = k = 1 K m = 1 M ( ϕ r k , g m ) n [ k , m ]
where we assume that the generation of each read is independent of other reads, K is the number of different reads, M is the number of reference species, and n [ k , m ] is the count of reads that have the same sequence as r k and whose corresponding species is g m . In real applications, we can estimate ϕ r i , g i by finding the ratio of high-quality hits for r i to all the high-quality read hits on the target reference sequence g i from the alignment result; that is,
ϕ r i , g i #   o f   r i   t h a t   h i t   g i   w i t h   h i g h   q u a l i t y #   o f   r e a d s   t h a t   h i t   g i   w i t h   h i g h   q u a l i t y
Thus, we have
p { z i = t | z ¬ i , R ; α , Φ } = p { z , r } p { z ¬ i , r } = p { r | z } p { z } p { r ¬ i | z ¬ i } p { z ¬ i } p { r i } p { r | z } p { r ¬ i | z ¬ i } p { z } p { z ¬ i } = k = 1 K m = 1 M ( ϕ r k , g m ) n [ k , m ] k = 1 K m = 1 M ( ϕ r k , g m ) n [ k , m ] , ¬ i ( α + n ) ( α + n ¬ i ) = ϕ r i , g t ( n t , ¬ i + α i ) j = 1 M ( n j , ¬ i + α j ) ϕ r i , g t ( n t , ¬ i + α i )
Here, we leave the priors out to simplify the notations. The relation Γ ( a + 1 ) = a Γ ( a ) is used in the above formula, and ¬ i indicates that the number is counted exclusive of the read r i . n [ k , m ] is the count of reads that have the same sequence as r k and whose corresponding species is g m . n [ k , m ] and n [ k , m ] , ¬ i are different only when k = i , m = t in the above formula, and n [ i , t ] = n [ i , t ] , ¬ i + 1 . n indicates the occurrence counts of read-origin species in the metagenomic reads, whose element
n j = n [ · , j ] = s = 1 M n [ s , j ]
Similarly, n and n ¬ i are the same except for their t th element in the above formula, and n t = n t , ¬ i + 1 . The above formula shows that read r i is more likely to be assigned to reference species g t if the probability that g t generates r i is large and if there are many reads in the metagenomic datasets that have been assigned to g t .
Finally, we can infer the genome mixture parameters θ using the counts vector n . For the Dirichlet-multinomial conjugacy, we have
p { θ | z , α } = D i r ( θ | n + α )
and we can estimate θ from the expectation of its distribution as follows:
θ k = n k + α k i = 1 M ( n i + α i )
n can be statistically estimated by finding the mean of a certain number of samples after the burn-in period of Gibbs sampling described above to overcome the ambiguity of the reads. The abundances can then be estimated. The pseudo-code of Gibbs sampling is shown in Algorithm 2.
Algorithm 2 (Gibbs Sampling Algorithm for the Metagenomic Model).
Require:references species M , metagenomic reads R , hyperparameter α
Global data:count statistics n , read-composition distributions Φ , memory for full conditionals p { z i | z ¬ i , R ; α , Φ }
Ensure:mixture parameters
θ //initialization: obtain read-composition distributions Φ according to alignment results zero all count statistics n
for i = 1 to N do sample the species index z i = m M u l t ( M ) increment sampled species count n m = n m + 1
end for//Gibbs sampling
whilenot finished do
for i = 1 to N dodecrement target species count n m = n m 1
sample a new species index z i = m ˜ p { z i | z ¬ i , R ; α , Φ }
increment sampled species’ count n m ˜ = n m ˜ + 1
end for
ifconverged and a given number of samples generated then
returnmixture parameter θ according to the equation
end if
end while

3. Results

3.1. The CoreProbe Framework

The primary contents of this article can be summarized as shown in Figure 1.
We then created a brief introduction to the specific CoreProbe process, which can be found in Figure 2.

3.2. Simulation Result

To evaluate the performance of the CoreProbe framework and compare it with the existing methods, we generated 90 simulated metagenomic datasets. We chose 40 microbial species as embed microorganisms in simulated metagenomic datasets. We aimed to estimate the relative abundances of 10 species, leaving the others as “unknown” species, and we compared the precision of the estimations among the existing methods. We downloaded the current completely sequenced microbial genomes of those species in which sequences may coexist, representing different strains of one species. We then employed MetaSim [48] to generate simulated reads of the selected genomes with preset relative abundances. Specifically, we built an Empirical Error Model based on the study of the Illumina’s sequencing technology [51,52,53] to generate metagenomic datasets that contained 1000, 2000, 5000, 10,000, 20,000, 50,000, 100,000, 200,000, and 500,000 single-ended reads with an average length of approximately 100 bp, and each dataset had 10 replicates.
When applying CoreProbe to those simulated datasets, we first downloaded all the comprehensive non-redundant reference core gene catalogs, which contained classification information and their corresponding sequences, and they also included over 2800 microbial species. Second, we selected certain microbial species as reference species, and we extracted all the core gene sequences for each reference species in the downloaded core gene catalogue. Third, we concatenated the core gene sequences of a reference species together into one FASTA file as a reference sequence of the species. We then used BWA MEM to map the simulated metagenomic reads against the reference core-genome sequences. All parameters were set as the defaults, and the output consisted of all the alignments in a SAM file.
The output SAM files were then passed to our program to estimate the relative abundances. In addition, we applied the three widely used algorithms for an abundance estimation, including GRAMMy [34], MetaPhlAn [38] and Kraken [21], and then used the microbial whole genome sequences as references for the metagenomic analysis to process the above 90 groups of metagenome simulation data and to estimate the relative abundance of 10 species of microorganisms. In addition, we compared the results using the CoreProbe algorithm to compare the accuracy and speed, to support the effectiveness of this algorithm.
It should be noted that when estimating the relative abundance of microorganisms, two types of whole genome sequenced strains were selected as the reference genome of the above algorithm, namely, (1) the selected strain used as a reference genome is contained in a sequence of strains used to generate simulated metagenomic datasets as g 1 ; and (2) the selected strain used as a reference genome is not included in the strain that is used to generate the simulated metagenomic datasets, and it is recorded as g 2 . This classification is because there are often new strains of certain microbial species in the natural environment, and the sequence has not been effectively sequenced. The use of as a reference genome is used to simulate these cases.

3.2.1. Comparison of Algorithm Accuracy

We compare the accuracy by using the relative mean square error of the relative abundance of each algorithm, and the results are shown in Figure 3a,b.
There are some common characteristics of the two graphs. With the increased capacity of the metagenome reads (i.e., the improvement of the depth of the sequencing), the mean variance of the relative abundance estimation of each algorithm decreases, indicating that the increase in the sequencing depth can improve the accuracy of each algorithm. MetaPhlAn [38] is rather special, in that the mean variance in the relative abundance values obtained here will begin to rise and then quickly decrease with the increase in the sequencing depth. When the metagenome reading capacity is less than 10,000, the relative abundance of the estimated variance is even more than 100%, indicating that with the increased depth of the sequencing, the accuracy will decrease; and when the reading capacity continues to increase, its mean square error will decline rapidly. Especially when the capacity of reads reaches up to 50,000, MetaPhlAn [38] is more accurate than Kraken [21]. Thus, we can observe that the MetaPhlAn algorithm [38] is not stable from the point of view of the sequencing depth, though it achieves high accuracy when the sequencing depth is sufficient. Kraken [21] is rarely affected by sequencing, with a mean square deviation remaining between 60% and 70%, indicating that the accuracy of the Kraken [21] is not high; there is a certain gap between the calculation results and the actual value. GRAMMy [34] has the highest accuracy; its variance is below 40%, and decreases occur with the increasing sequencing depth. The CoreProbe algorithm mentioned in this paper also uses the whole genome for the reference sequence. Its accuracy is better than that of MetaPhlAn [38] and Kraken, but there is still a gap with the GRAMMy algorithm [34]. However, when CoreProbe uses the core-genome sequence as the reference sequence, the accuracy is significantly improved. When the metagenome reads capacity is increased to 50,000, the accuracy is very close to that of GRAMMy [34], and the mean square error in the preset value is less than 10%.
However, by comparing Figure 3a and Figure 3b, the accuracy of each algorithm was found to decrease when the reference genome changes from g 1 to g 2 . Although the difference is not large, the selected reference genome sequence corresponds to the strains present in the sample environment or not, and they will have some impact on the algorithm. In this paper, CoreProbe uses the core gene sequence as the reference sequence, and thus, the variations in the strain in the real environment and the unknown strain that affect the algorithm are relatively small. Therefore, the relative abundance estimation is made using the microbial core genome as the reference. When the actual environment has not yet sequenced new strains, the algorithm estimates of the relative abundance for the accuracy can be more credible.

3.2.2. Comparison of Algorithm Speeds

The various algorithms mentioned in this paper are used in different algorithms, and during the process of dealing with data, there are different sequence mapping algorithms; for example, Kraken [21], MetaPhlAn [38], and GRAMMy [34] correspond to BLAST [26], Bowtie2 [28], and BWA [27]. Thus, the time spent on the calculation process (the sum of time used by the sequence alignment and the statistical algorithm) will be different, and we used the line graph to show the efficiency of each algorithm, as shown in Figure 4a,b.
The overall situation of the two figures is similar. As shown in Figure 4a,b, with the increased capacity of the metagenome (i.e., the improvement of the depth of the sequencing), the relative abundance estimation time of each algorithm is increased. However, the Kraken [21] is special, because the time it consumes does not increase as the depth of the sequencing increases, and the time was almost maintained at approximately ten seconds. GRAMMy [34] is the most accurate, but it is time-consuming, and especially when the sequencing depth increases, the growth is faster than it is in other algorithms. MetaPhlAn [38] is significantly more time-consuming than the other algorithms, and the computational accuracy is lower than that of other algorithms. When the sequencing depth increases, its required time starts to increase significantly, and its accuracy also starts to improve significantly. In addition, our method, CoreProbe, is the least time-consuming, and it only needs more time than Kraken [21] when the number of reads goes up to 200,000. At the same time, by comparing Figure 4a and Figure 4b when the reference genome changes from to , the speed of each algorithm is almost unaffected. This finding shows that each algorithm is relatively stable, and when the amount of data is equivalent, the calculation process will not change the necessary time. Because the microbial core genome is the gene that is common to all the strains in the clade, the number is much smaller than that of the whole genome, and thus the computational time can be dramatically reduced and the relative abundance can be estimated more efficiently. As shown in the figure, when the core genome was used as a reference sequence, the time consumption was less than that of the metagenomes g 1 and g 2 .

3.3. Real Metagenomic Datasets Analysis

As an immediate application, we applied this CoreProbe pipeline to a set of actual data. These data come from the HMP, and its sub-database HMIWGS contains 764 groups of high-throughput metagenomic sample data from 16 different human tissues via the Illumina platforms. In this paper, five different human tissues were selected, including the anterior nares, buccal mucosa, posterior fornix, right retro-auricular crease, stool, and five different genomic samples that were randomly selected from each organization. The metagenome data in this paper have passed the quality control (QC) test of the HMP. We intended to use these data to estimate the relative abundance of microorganisms and to classify the metagenomic data. Thus, no additional preprocessing of the above metagenome data was performed.
In this section, 150 common microbes were selected, and the gene sequences belonging to their core genomes were selected from the MetaRef database. Their core gene sequences were grouped together to form a FASTA formatted file as the reference genome sequence of the corresponding microorganism. In this paper, we used the BWA MEM program to compare the sequence of the actual genome reading with the reference genome. This parameter was used to locate the default value and output all the matching results. Finally, the CoreProbe algorithm implemented in this paper was used to estimate the relative abundance of microbes in each sample.
On this basis, the relative abundance of the obtained microbes was analyzed by heat map. In Figure 5, the buccal mucosa, vaginal dome and stool samples were basically clustered into one class, indicating that the relative abundance of microbes showed similarities within their respective tissues and had significant differences from other human tissues, which is similar to the results of other studies [54]. In addition, the anterior nares and the right retro-auricular crease corresponding to the metagenomic samples have been doped, and they cannot be clustered well, indicating that the relative microbial abundance between the two human tissues is of great similarity, which may result from close contact between the anterior nares and the right retro-auricular crease and the lesser secretion of body fluid. In addition, there are more white cells relative to the posterior fornix, indicating that the microbial species of this tissue is less than that of others; a small number of cells presented a violet color, and the rest were generally shallow, indicating that the tissue often consists of a small number of microorganisms as the dominant population, and the remaining microbial content is scant. Studies have shown that this finding is due to the important position of the female vagina in the human reproductive system, and its need to inhibit the growth of microbial species. In addition, according to its dominant microbial species, they can be divided into three types of microbial community types [55]. In the classification tree on the left side of the thermogram, it is also clear that the posterior fornix sample is divided into three sub-trees.
Through the above analysis, we can see that the CoreProbe algorithm can accurately and quickly estimate the relative abundance of microbes and analyze them by heat map, which can distinguish among the different genomic samples from different tissues, and it can be used to find the differences from and links in microbial distributions between different tissues. It is of great significance to understand the relationship between human health and microorganisms.

4. Discussion

We developed the CoreProbe framework to estimate the relative abundance of genomes. This approach has three unique advantages. First, a complete probability model is established for the relative abundance of microbes. The distribution and structure of microbes are simulated on the basis of this finding; thus, the relative abundance of various microorganisms can be estimated more accurately.
Second, CoreProbe uses the stochastic algorithm of the Gibbs sampling algorithm as an alternative to the deterministic algorithm of statistical reasoning (such as the expectation maximization algorithm) to achieve a better fuzzy allocation of the readings, which makes the method particularly suitable for short read data sets. The distribution of ambiguous sources may be sequencing errors, genetic variation, horizontal gene transfer or closely related genomes. The Gibbs sampling algorithm is a good solution to this series of problems.
Third, due to the changes in microorganisms in the actual environment to adapt to the changes in the surrounding ecology, an abnormal mutation mechanism is evolved, and the genomic polymorphism at the strain level is not uncommon. An environment often contains different strains of a certain microorganism, and the relative abundance of these strains is inconsistent. For the strains that have never been found and sequenced, these factors have greatly affected the accuracy of the relative abundance statistics algorithm for the genome sequence of the few sequenced strains as the reference genome. The core genome is a collection of all the genes of all the strains, which maintain high stability during microbial variation [56].
In using the microbial core genome as a reference genome, the effect of genomic polymorphism and strain diversity on the sequence alignment process can be eliminated to a greater extent, and the accuracy of the relative abundance estimation can be improved. In addition, the core genome can be used as a common gene set of all the strains, accounting for only approximately 8% of all the gene sets of the microbial species. Compared with the algorithm that uses the whole genome of the sequenced strain as a reference gene, the efficiency of the comparison is greatly improved.
In summary, the CoreProbe method we provide is likely to provide a more accurate and efficient estimate of the relative abundance estimates of microbes and to uncover a new direction for estimating abundance.

Author Contributions

D.A., R.H. and L.C.X. conceived and designed the experiments; H.P. and R.H. performed the experiments and analyzed the data; and D.A., H.P. and R.H. wrote the paper.

Funding

This work was supported by grants from the National Natural Science Foundation of China (61370131).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Barberán, A.; Bates, S.T.; Casamayor, E.O.; Fierer, N. Using network analysis to explore co-occurrence patterns in soil microbial communities. ISME J. 2012, 6, 343–351. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Allison, S.D.; Lu, Y.; Weihe, C.; Goulden, M.L.; Martiny, A.C.; Treseder, K.K.; Martiny, J.B. Microbial abundance and composition influence litter decomposition response to environmental change. Ecology 2013, 94, 714–725. [Google Scholar] [CrossRef] [PubMed]
  3. Li, X.; Yan, Q.; Xie, S.; Hu, W.; Yu, Y.; Hu, Z. Gut microbiota contributes to the growth of fast-growing transgenic common carp (Cyprinus carpio L.). PLoS ONE 2013, 8, e64577. [Google Scholar] [CrossRef] [PubMed]
  4. Smith, M.I.; Yatsunenko, T.; Manary, M.J.; Trehan, I.; Mkakosya, R.; Cheng, J.; Kau, A.L.; Rich, S.S.; Concannon, P.; Mychaleckyj, J.C. Gut microbiomes of Malawian twin pairs discordant for kwashiorkor. Science 2013, 339, 548–554. [Google Scholar] [CrossRef] [PubMed]
  5. Vaarala, O. Human intestinal microbiota and type 1 diabetes. Curr. Diabetes Rep. 2013, 13, 601–607. [Google Scholar] [CrossRef] [PubMed]
  6. Ghoshal, U.C.; Shukla, R.; Ghoshal, U.; Gwee, K.-A.; Ng, S.C.; Quigley, E.M. The gut microbiota and irritable bowel syndrome: Friend or foe? Int. J. Inflamm. 2012, 2012. [Google Scholar] [CrossRef] [PubMed]
  7. Manichanh, C.; Borruel, N.; Casellas, F.; Guarner, F. The gut microbiota in IBD. Nat. Rev. Gastroenterol. Hepatol. 2012, 9, 599–608. [Google Scholar] [CrossRef] [PubMed]
  8. Li, Q.; Wang, C.; Tang, C.; Li, N.; Li, J. Molecular-phylogenetic characterization of the microbiota in ulcerated and non-ulcerated regions in the patients with Crohn’s disease. PLoS ONE 2012, 7, e34939. [Google Scholar] [CrossRef] [PubMed]
  9. Kostic, A.D.; Gevers, D.; Pedamallu, C.S.; Michaud, M.; Duke, F.; Earl, A.M.; Ojesina, A.I.; Jung, J.; Bass, A.J.; Tabernero, J. Genomic analysis identifies association of Fusobacterium with colorectal carcinoma. Genome Res. 2012, 22, 292–298. [Google Scholar] [CrossRef] [PubMed]
  10. Kostic, A.D.; Chun, E.; Robertson, L.; Glickman, J.N.; Gallini, C.A.; Michaud, M.; Clancy, T.E.; Chung, D.C.; Lochhead, P.; Hold, G.L. Fusobacterium nucleatum potentiates intestinal tumorigenesis and modulates the tumor-immune microenvironment. Cell Host Microbe 2013, 14, 207–215. [Google Scholar] [CrossRef] [PubMed]
  11. Marshall, B.J.; Windsor, H.M. The relation of Helicobacter pylori to gastric adenocarcinoma and lymphoma: Pathophysiology, epidemiology, screening, clinical presentation, treatment, and prevention. Med. Clin. 2005, 89, 313–344. [Google Scholar] [CrossRef] [PubMed]
  12. Teeling, H.; Waldmann, J.; Lombardot, T.; Bauer, M.; Glöckner, F.O. TETRA: A web-service and a stand-alone program for the analysis and comparison of tetranucleotide usage patterns in DNA sequences. BMC Bioinform. 2004, 5, 163. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Chatterji, S.; Yamazaki, I.; Bai, Z.; Eisen, J. CompostBin: A DNA composition-based algorithm for binning environmental shotgun reads. In Research in Computational Molecular Biology; Springer: Berlin, Germany, 2008; pp. 17–28. [Google Scholar]
  14. Kelley, D.R.; Salzberg, S.L. Clustering metagenomic sequences with interpolated Markov models. BMC Bioinform. 2010, 11, 544. [Google Scholar] [CrossRef] [PubMed]
  15. Leung, H.C.; Yiu, S.-M.; Yang, B.; Peng, Y.; Wang, Y.; Liu, Z.; Chen, J.; Qin, J.; Li, R.; Chin, F.Y. A robust and accurate binning algorithm for metagenomic sequences with arbitrary species abundance ratio. Bioinformatics 2011, 27, 1489–1495. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. McHardy, A.C.; Martin, H.G.; Tsirigos, A.; Hugenholtz, P.; Rigoutsos, I. Accurate phylogenetic classification of variable-length DNA fragments. Nat. Methods 2007, 4, 63–72. [Google Scholar] [CrossRef] [PubMed]
  17. Rosen, G.; Garbarine, E.; Caseiro, D.; Polikar, R.; Sokhansanj, B. Metagenome Fragment Classification Using N-Mer Frequency Profiles. Adv. Bioinf. 2008. [Google Scholar] [CrossRef] [PubMed]
  18. Diaz, N.N.; Krause, L.; Goesmann, A.; Niehaus, K.; Nattkemper, T.W. TACOA—Taxonomic classification of environmental genomic fragments using a kernelized nearest neighbor approach. BMC Bioinform. 2009, 10, 56. [Google Scholar] [CrossRef] [PubMed]
  19. Nalbantoglu, O.U.; Way, S.F.; Hinrichs, S.H.; Sayood, K. RAIphy: Phylogenetic classification of metagenomics samples using iterative refinement of relative abundance index profiles. BMC Bioinform. 2011, 12, 41. [Google Scholar] [CrossRef] [PubMed]
  20. Ames, S.K.; Hysom, D.A.; Gardner, S.N.; Lloyd, G.S.; Gokhale, M.B.; Allen, J.E. Scalable metagenomic taxonomy classification using a reference genome database. Bioinformatics 2013, 29, 2253–2260. [Google Scholar] [CrossRef] [PubMed]
  21. Wood, D.E.; Salzberg, S.L. Kraken: Ultrafast metagenomic sequence classification using exact alignments. Genome Biol. 2014, 15, R46. [Google Scholar] [CrossRef] [PubMed]
  22. Ounit, R.; Wanamaker, S.; Close, T.J.; Lonardi, S. CLARK: Fast and accurate classification of metagenomic and genomic sequences using discriminative k-mers. BMC Genom. 2015, 16, 236. [Google Scholar] [CrossRef] [PubMed]
  23. Brady, A.; Salzberg, S.L. Phymm and PhymmBL: Metagenomic phylogenetic classification with interpolated Markov models. Nat. Methods 2009, 6, 673–676. [Google Scholar] [CrossRef] [PubMed]
  24. MacDonald, N.J.; Parks, D.H.; Beiko, R.G. Rapid identification of high-confidence taxonomic assignments for metagenomic data. Nucleic Acids Res. 2012, 40, e111. [Google Scholar] [CrossRef] [PubMed]
  25. Yan, J.; Chuai, G.; Qi, T.; Shao, F.; Zhou, C.; Zhu, C.; Yang, J.; Yu, Y.; Shi, C.; Kang, N. MetaTopics: An integration tool to analyze microbial community profile by topic model. BMC Genom. 2017, 18, 962. [Google Scholar] [CrossRef] [PubMed]
  26. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  27. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows—Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [PubMed]
  28. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Huson, D.H.; Auch, A.F.; Qi, J.; Schuster, S.C. MEGAN analysis of metagenomic data. Genome Res. 2007, 17, 377–386. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Krause, L.; Diaz, N.N.; Goesmann, A.; Kelley, S.; Nattkemper, T.W.; Rohwer, F.; Edwards, R.A.; Stoye, J. Phylogenetic classification of short environmental DNA fragments. Nucleic Acids Res. 2008, 36, 2230–2239. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Berger, S.A.; Stamatakis, A. Aligning short reads to reference alignments and trees. Bioinformatics 2011, 27, 2068–2075. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Gori, F.; Folino, G.; Jetten, M.S.; Marchiori, E. MTR: Taxonomic annotation of short metagenomic reads using clustering at multiple taxonomic ranks. Bioinformatics 2010, 27, 196–203. [Google Scholar] [CrossRef] [PubMed]
  33. Angly, F.E.; Willner, D.; Prieto-Davó, A.; Edwards, R.A.; Schmieder, R.; Vega-Thurber, R.; Antonopoulos, D.A.; Barott, K.; Cottrell, M.T.; Desnues, C. The GAAS metagenomic tool and its estimations of viral and microbial average genome size in four major biomes. PLoS Comput. Biol. 2009, 5, e1000593. [Google Scholar] [CrossRef] [PubMed]
  34. Xia, L.C.; Cram, J.A.; Chen, T.; Fuhrman, J.A.; Sun, F. Accurate genome relative abundance estimation based on shotgun metagenomic reads. PLoS ONE 2011, 6, e27992. [Google Scholar] [CrossRef] [PubMed]
  35. Hong, C.; Manimaran, S.; Shen, Y.; Perez-Rogers, J.F.; Byrd, A.L.; Castro-Nallar, E.; Crandall, K.A.; Johnson, W.E. PathoScope 2.0: A complete computational framework for strain identification in environmental or clinical sequencing samples. Microbiome 2014, 2, 33. [Google Scholar] [CrossRef] [PubMed]
  36. Lindner, M.S.; Renard, B.Y. Metagenomic abundance estimation and diagnostic testing on species level. Nucleic Acids Res. 2012, 41, e10. [Google Scholar] [CrossRef] [PubMed]
  37. Morfopoulou, S.; Plagnol, V. Bayesian mixture analysis for metagenomic community profiling. Bioinformatics 2015, 31, 2930–2938. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Segata, N.; Waldron, L.; Ballarini, A.; Narasimhan, V.; Jousson, O.; Huttenhower, C. Metagenomic microbial community profiling using unique clade-specific marker genes. Nat. Methods 2012, 9, 811–814. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Sunagawa, S.; Mende, D.R.; Zeller, G.; Izquierdo-Carrasco, F.; Berger, S.A.; Kultima, J.R.; Coelho, L.P.; Arumugam, M.; Tap, J.; Nielsen, H.B. Metagenomic species profiling using universal phylogenetic marker genes. Nat. Methods 2013, 10, 1196–1199. [Google Scholar] [CrossRef] [PubMed]
  40. Tettelin, H.; Riley, D.; Cattuto, C.; Medini, D. Comparative genomics: The bacterial pan-genome. Curr. Opin. Microbiol. 2008, 11, 472–477. [Google Scholar] [CrossRef] [PubMed]
  41. Kashtan, N.; Roggensack, S.E.; Rodrigue, S.; Thompson, J.W.; Biller, S.J.; Coe, A.; Ding, H.; Marttinen, P.; Malmstrom, R.R.; Stocker, R. Single-cell genomics reveals hundreds of coexisting subpopulations in wild Prochlorococcus. Science 2014, 344, 416–420. [Google Scholar] [CrossRef] [PubMed]
  42. Ley, R.E.; Peterson, D.A.; Gordon, J.I. Ecological and evolutionary forces shaping microbial diversity in the human intestine. Cell 2006, 124, 837–848. [Google Scholar] [CrossRef] [PubMed]
  43. Vernikos, G.; Medini, D.; Riley, D.R.; Tettelin, H. Ten years of pan-genome analyses. Curr. Opin. Microbiol. 2015, 23, 148–154. [Google Scholar] [CrossRef] [PubMed]
  44. Stecher, B.; Denzler, R.; Maier, L.; Bernet, F.; Sanders, M.J.; Pickard, D.J.; Barthel, M.; Westendorf, A.M.; Krogfelt, K.A.; Walker, A.W. Gut inflammation can boost horizontal gene transfer between pathogenic and commensal Enterobacteriaceae. Proc. Natl. Acad. Sci. USA 2012, 109, 1269–1274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 2000, 155, 945–959. [Google Scholar] [PubMed]
  46. Blei, D.M.; Ng, A.Y.; Jordan, M.I. Latent dirichlet allocation. J. Mach. Learn. Res. 2003, 3, 993–1022. [Google Scholar]
  47. Heinrich, G. Parameter Estimation for Text Analysis; Technical Report; University of Leipzig and Vsonix GmbH: Darmstadt, Germany, 2008. [Google Scholar]
  48. Richter DC, Ott F, Auch AF, Schmid R, Huson DH: MetaSim—A sequencing simulator for genomics and metagenomics. PLoS ONE 2008, 3, e3373.
  49. Shakya, M.; Quince, C.; Campbell, J.H.; Yang, Z.K.; Schadt, C.W.; Podar, M. Comparative metagenomic and rRNA microbial diversity characterization using archaeal and bacterial synthetic communities. Environ. Microbiol. 2013, 15, 1882–1899. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Gelman, A.; Carlin, J.B.; Stern, H.S.; Dunson, D.B.; Vehtari, A.; Rubin, D.B. Bayesian Data Analysis, 3rd ed.; CRC Press: Boca Raton, FL, USA, 2013. [Google Scholar]
  51. Schirmer, M.; Ijaz, U.Z.; D’Amore, R.; Hall, N.; Sloan, W.T.; Quince, C. Insight into biases and sequencing errors for amplicon sequencing with the Illumina MiSeq platform. Nucleic Acids Res. 2015, 43, e37. [Google Scholar] [CrossRef] [PubMed]
  52. Nakamura, K.; Oshima, T.; Morimoto, T.; Ikeda, S.; Yoshikawa, H.; Shiwa, Y.; Ishikawa, S.; Linak, M.C.; Hirai, A.; Takahashi, H. Sequence-specific error profile of Illumina sequencers. Nucleic Acids Res. 2011, 39, e90. [Google Scholar] [CrossRef] [PubMed]
  53. Minoche, A.E.; Dohm, J.C.; Himmelbauer, H. Evaluation of genomic high-throughput sequencing data generated on Illumina HiSeq and genome analyzer systems. Genome Biol. 2011, 12, R112. [Google Scholar] [CrossRef] [PubMed]
  54. Consortium, H.M.P. Structure, function and diversity of the healthy human microbiome. Nature 2012, 486, 207–214. [Google Scholar] [Green Version]
  55. Arumugam, M.; Raes, J.; Pelletier, E.; Le Paslier, D.; Yamada, T.; Mende, D.R.; Fernandes, G.R.; Tap, J.; Bruls, T.; Batto, J.-M. Addendum: Enterotypes of the human gut microbiome. Nature 2014, 506, 516. [Google Scholar] [CrossRef]
  56. Hu, B.; Xie, G.; Lo, C.-C.; Starkenburg, S.R.; Chain, P.S. Pathogen comparative genomics in the next-generation sequencing era: Genome alignments, pangenomics and metagenomics. Brief. Funct. Genom. 2011, 10, 322–333. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Overview of the overall content of the article. We first downloaded the reference sequence from the National Center for Biotechnology Information (NCBI) database (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/), the MetaRef database, and simulated the metagenomic data with MetaSim software [48]. We then used our CoreProbe to analyze the simulated data to estimate the genome relative abundance, to compare them with the existing methods to assess the accuracy and efficiency of our method. Last, our method was applied to the actual data obtained from Human Microbiome Project (HMP), and a preliminary analysis was made to support the practicability of our method.
Figure 1. Overview of the overall content of the article. We first downloaded the reference sequence from the National Center for Biotechnology Information (NCBI) database (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/), the MetaRef database, and simulated the metagenomic data with MetaSim software [48]. We then used our CoreProbe to analyze the simulated data to estimate the genome relative abundance, to compare them with the existing methods to assess the accuracy and efficiency of our method. Last, our method was applied to the actual data obtained from Human Microbiome Project (HMP), and a preliminary analysis was made to support the practicability of our method.
Genes 09 00313 g001
Figure 2. The CoreProbe framework. First, the metagenomic reads obtained by shotgun sequencing method was compared with the core genome reference sequence using BWA [27], and the results of the alignment were recorded. The results were then combined with our reads set and the reference sequences applied to the Gibbs sampling algorithm to obtain the relative abundance of the metagenome. Finally, if the taxonomic information for the input reference genomes was available, we could calculate a high level of taxonomic abundance combined with the strain level estimates.
Figure 2. The CoreProbe framework. First, the metagenomic reads obtained by shotgun sequencing method was compared with the core genome reference sequence using BWA [27], and the results of the alignment were recorded. The results were then combined with our reads set and the reference sequences applied to the Gibbs sampling algorithm to obtain the relative abundance of the metagenome. Finally, if the taxonomic information for the input reference genomes was available, we could calculate a high level of taxonomic abundance combined with the strain level estimates.
Genes 09 00313 g002
Figure 3. The relative mean square error of the relative abundance of each algorithm. (a) The selected strain as a reference genome is contained in a sequence of strains used to generate simulated metagenomic datasets as g 1 ; and (b) the selected strain as a reference genome is contained in a sequence of strains used to generate simulated metagenomic datasets as g 2 . The abscissa represents the different metagenome simulation data set of the reading capacity, and the ordinate is the mean square error of the algorithm relative abundance between the estimated value and the preset value.
Figure 3. The relative mean square error of the relative abundance of each algorithm. (a) The selected strain as a reference genome is contained in a sequence of strains used to generate simulated metagenomic datasets as g 1 ; and (b) the selected strain as a reference genome is contained in a sequence of strains used to generate simulated metagenomic datasets as g 2 . The abscissa represents the different metagenome simulation data set of the reading capacity, and the ordinate is the mean square error of the algorithm relative abundance between the estimated value and the preset value.
Genes 09 00313 g003
Figure 4. The time used for each algorithm. (a) The time used for each algorithm in different reads capacities with g 1 as the reference genome. (b) The time used for each algorithm in different reads capacities with g 2 as the reference genome. The abscissa represents the different metagenome simulation datasets of the reads capacity, and the vertical coordinate is the time spent on the calculation process of each method.
Figure 4. The time used for each algorithm. (a) The time used for each algorithm in different reads capacities with g 1 as the reference genome. (b) The time used for each algorithm in different reads capacities with g 2 as the reference genome. The abscissa represents the different metagenome simulation datasets of the reads capacity, and the vertical coordinate is the time spent on the calculation process of each method.
Genes 09 00313 g004
Figure 5. Thermal graph of microbial species relative abundance data in five different human tissues. The horizontal axis corresponds to the different microbial species, and the vertical axis corresponds to the metagenome samples of different human tissues. The white cells indicate that the microorganism corresponding to its vertical axis does not appear in the corresponding sample on the horizontal axis. The change in color from light blue to deep purple indicates that the relative abundance of microorganisms in the corresponding sample varies from low to high. The distance used by the cluster analysis of different metagenomic samples on the left side of the thermogram is derived from the Spearman correlation coefficient between the relative microbial abundance vectors (using the difference between unit 1 and the correlation coefficient as a distance).
Figure 5. Thermal graph of microbial species relative abundance data in five different human tissues. The horizontal axis corresponds to the different microbial species, and the vertical axis corresponds to the metagenome samples of different human tissues. The white cells indicate that the microorganism corresponding to its vertical axis does not appear in the corresponding sample on the horizontal axis. The change in color from light blue to deep purple indicates that the relative abundance of microorganisms in the corresponding sample varies from low to high. The distance used by the cluster analysis of different metagenomic samples on the left side of the thermogram is derived from the Spearman correlation coefficient between the relative microbial abundance vectors (using the difference between unit 1 and the correlation coefficient as a distance).
Genes 09 00313 g005

Share and Cite

MDPI and ACS Style

Ai, D.; Pan, H.; Huang, R.; Xia, L.C. CoreProbe: A Novel Algorithm for Estimating Relative Abundance Based on Metagenomic Reads. Genes 2018, 9, 313. https://0-doi-org.brum.beds.ac.uk/10.3390/genes9060313

AMA Style

Ai D, Pan H, Huang R, Xia LC. CoreProbe: A Novel Algorithm for Estimating Relative Abundance Based on Metagenomic Reads. Genes. 2018; 9(6):313. https://0-doi-org.brum.beds.ac.uk/10.3390/genes9060313

Chicago/Turabian Style

Ai, Dongmei, Hongfei Pan, Ruocheng Huang, and Li C. Xia. 2018. "CoreProbe: A Novel Algorithm for Estimating Relative Abundance Based on Metagenomic Reads" Genes 9, no. 6: 313. https://0-doi-org.brum.beds.ac.uk/10.3390/genes9060313

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop