Wednesday, September 13, 2017

unweighted UniFrac distance metric

UniFrac

From Wikipedia, the free encyclopedia
UniFrac is a distance metric used for comparing biological communities. It differs from dissimilarity measures such as Bray-Curtis dissimilarity in that it incorporates information on the relative relatedness of community members by incorporating phylogenetic distances between observed organisms in the computation. Both weighted (quantitative) and unweighted (qualitative) variants of UniFrac are widely used in microbial ecology, where the former accounts for abundance of observed organisms, while the latter only considers their presence or absence. The method was devised by Catherine Lozupone, when she was working under Rob Knight of the University of Colorado at Boulder in 2005.

Research methods

The distance is calculated between pairs of samples (each sample represents an organismal community). All taxa found in one or both samples are placed on a phylogenetic tree. A branch leading to taxa from both samples is marked as "shared" and branches leading to taxa which appears only in one sample are marked as "unshared". The distance between the two samples is then calculated as (the sum of "unshared" branch lengths)/(the sum of all tree branch lengths (= shared+unshared)), i.e. the fraction of total branch length which is unshared. This definition satisfies the requirements of a distance metric, being non-negative, zero only when entities are identical, transitive, and conformant to the triangle inequality.
If there are several different samples, a distance matrix can be created, by making a tree for each pair of samples and calculating their UniFrac measure. Later, standard multivariate statistical, methods such as data clustering and principal co-ordinates analysis can be used.
One can determine the statistical significance of the Unifrac distance between two samples using Monte Carlo simulations. By randomizing the sample classification of each taxa on the tree (leaving the branch structure unchanged) and creating a distribution of UniFrac distance values, one can obtain a distribution of UniFrac values. From this, a p-value can be given to the actual distance between the samples.
Additionally, there is a weighted version of the UniFrac metric which accounts for the relative abundance of each of the taxa within the communities. This is commonly used in metagenomic studies, where the number of metagenomic reads can be in the tens of thousands, and it is appropriate to 'bin' these reads into operational taxonomic units, or OTUs, which can then be dealt with as taxa within the UniFrac framework.
In 2012, a generalized UniFrac version, which unifies the weighted and unweighted UniFrac distance in a single framework, was proposed. The weighted and unweighted UniFrac distance place too much weight on either abundant lineages or rare lineages. Their power to detect environmental influence is limited under some setting, where the moderately abundant lineages are mostly affected. The generalized UniFrac distance corrects the limitation of the weighted/weighted UniFrac distance by down-weighting their emphasis on either abundant or rare lineages.

No comments: