Variance adjusted weighted UniFrac: a powerful beta diversity measure for comparing communities based on phylogeny

Abstract

BackgroundBeta diversity, which involves the assessment of differences between communities, is an important problem in ecological studies. Many statistical methods have been developed to quantify beta diversity, and among them, UniFrac and weighted-UniFrac (W-UniFrac) are widely used. The W-UniFrac is a weighted sum of branch lengths in a phylogenetic tree of the sequences from the communities. However, W-UniFrac does not consider the variation of the weights under random sampling resulting in less power detecting the differences between communities.ResultsWe develop a new statistic termed variance adjusted weighted UniFrac (VAW-UniFrac) to compare two communities based on the phylogenetic relationships of the individuals. The VAW-UniFrac is used to test if the two communities are different. To test the power of VAW-UniFrac, we first ran a series of simulations which revealed that it always outperforms W-UniFrac, as well as UniFrac when the individuals are not uniformly distributed. Next, all three methods were applied to analyze three large 16S rRNA sequence collections, including human skin bacteria, mouse gut microbial communities, microbial communities from hypersaline soil and sediments, and a tropical forest census data. Both simulations and applications to real data show that VAW-UniFrac can satisfactorily measure differences between communities, considering not only the species composition but also abundance information.ConclusionsVAW-UniFrac can recover biological insights that cannot be revealed by other beta diversity measures, and it provides a novel alternative for comparing communities.

Associating microbiome composition with environmental covariates using generalized UniFrac distances

“…A recently proposed variance adjusted UniFrac distance ( d VAW ) (Chang et al ., 2011) is also compared in this setting. d VAW was developed to moderate the branch proportion difference by its variance and was shown to increase the power of detecting the difference between two microbial communities.…”

Characterization of oral bacterial diversity of irradiated patients by high-throughput sequencing

“…By comparing the compositions of oral communities from different time periods, we can learn how specific factors affect community development and how species or individuals associate with each other. 28 Cluster analysis using weighted UniFrac metrics (Figure 4) found that samples before radiotherapy (PT) and the early stages of radiotherapy (10 Gy, 20 Gy and 30 Gy) formed a cluster distinct from the time points at later stages of radiotherapy (40 Gy, 50 Gy and 60 Gy), indicating that there existed temporal variations in microbial communities through the course of treatment. Comparing the cluster analysis with the phylogenetic tree generated by the R program showed similar results for the distances based on genus composition (Figure 2).…”

Insect exclusion limits variation in bacterial microbiomes of tomato flowers and fruit

“…Beta diversity metrics weighted by phylogenetic distance were used to generate distance matrices for further analysis (Lozupone and Knight ; Chang et al . ). NMDS plots were created in R using ggplot2 and vegan (Oksanen et al .…”