UniFrac

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:

$\left({\frac {sum~of~unshared~branch~lengths}{sum~of~all~tree~branch~lengths}}\right){=}fraction~of~total~unshared~branch~lengths$

This definition satisfies the requirements of a distance metric, being non-negative, zero only when entities are identical, transitive, and conforming to the triangle inequality.

Three examples of the triangle inequality for triangles with sides of lengths x, y, z. The top example shows a case where z is much less than the sum x + y of the other two sides, and the bottom example shows a case where the side z is only slightly less than x + y.

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. Subsequently, 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 taxon 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 settings, 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.

External links

• UniFrac online
• Knight Lab website
• Description of UniFrac, with worked examples