Power law analysis of the human milk microbiome

The human breast milk microbiome (HMM) has far reached health implications for both mothers and infants, and understanding the structure and dynamics of milk microbial communities is therefore of critical biomedical importance. Community heterogeneity, which has certain commonalities with familiar diversity but also with certain fundamental differences, is an important aspect of community structure and dynamics. Taylor's (1961) power law (TPL) (Nature, 1961) was discovered to govern the mean-variance power function relationship of population abundances and can be used to characterize population spatial aggregation (heterogeneity) and/or temporal stability. TPL was further extended to the community level to measure community spatial heterogeneity and/or temporal stability (Ma 2015, Molecular Ecology). Here, we applied TPL extensions (TPLE) to analyze the heterogeneity of the human milk microbiome by reanalyzing 12 datasets (2115 samples) of the healthy human milk microbiome. Our analysis revealed that the TPLE heterogeneity parameter (b) is rather stable across the 12 datasets, and there were approximately no statistically significant differences among ¾ of the datasets, which is consistent with the hypothesis that the heterogeneity scaling (i.e., change across individuals) of the human microbiome, including HMM, is rather stable or even constant. For this, we built a TPLE model for the pooled 12 datasets (b = 1.906), which can therefore represent the scaling rate of community-level spatial heterogeneity of HMM across individuals. Similarly, we also analyzed mixed-species ("averaged virtual species") level heterogeneity of HMM, and it was found that the mixed-species level heterogeneity was smaller than the heterogeneity at the previously mentioned community level (1.620 vs. 1.906).