Spatial heterogeneity analysis of the human virome with Taylor’s power law

Spatial heterogeneity is a fundamental characteristic of organisms from viruses to humans. Measuring heterogeneity is challenging, especially for naked-eye invisible viruses, but of obvious importance. For example, spatial heterogeneity of virus distribution may strongly influence infection spreading and outbreaks in the case of pathogenic viruses; the spatial distribution (i.e., the inter-subject heterogeneity) of commensal viruses within/on our bodies can influence the competition, coexistence, and dispersal of viruses within or between our bodies. Taylor’s power law (TPL) was first discovered in the 1960s to describe the spatial distributions of plant and/or animal populations, and since then it has been verified by numerous experimental and theoretical studies. Recently, TPL has been extended from population to community level and applied to bacterial communities. Here we report the first comprehensive testing of the TPL fitted to human virome datasets. It was found that the human virome follows the TPL as bacterial communities do. Furthermore, the TPL heterogeneity scaling parameter of human virome is virtually the same as that of the human bacterial microbiome (1.916 vs. 1.926). We postulate that the extreme closeness of human viruses and bacteria in heterogeneity scaling coefficients could be attributed to the fact that most of the viruses that were annotated in this study actually belong to bacteriophages (86% viral OTUs) that “piggyback” on their bacterial hosts, and their distributions are likely host-dependent. The scaling parameter, which measures the inter-subject heterogeneity changes, should be an innate property of human microbiomes including both bacteria and viruses. It is similar to the acceleration coefficient of the gravity (g = 9.8) as specified by Newton’s law, which is invariant on the earth. Nevertheless, we caution that our postulation is contingent on an implicit assumption that the proportion of bacteriophages to total virome may not change significantly when more virus species can be identified in future.