A 3-dimensional mathematical model of microbial proliferation that generates the characteristic cumulative relative abundance distributions in gut microbiomes

The gut microbiome is highly variable among individuals, largely due to differences in host lifestyle and physiology. However, little is known about the underlying processes or rules that shape the complex microbial community. In this paper, we show that the cumulative relative abundance distribution (CRAD) of microbial species can be approximated by a power law function, and found that the power exponent of CRADs generated from 16S rRNA gene and metagenomic data for normal gut microbiomes of humans and mice was similar consistently with ∼0.9. A similarly robust power exponent was observed in CRADs of gut microbiomes during dietary interventions and several diseases. However, the power exponent was found to be ∼0.6 in CRADs from gut microbiomes characterized by lower species richness, such as those of human infants and the small intestine of mice. In addition, the CRAD of gut microbiomes of mice treated with antibiotics differed slightly from those of infants and the small intestines of mice. Based on these observations, in addition to data on the spatial distribution of microbes in the digestive tract, we developed a 3-dimensional mathematical model of microbial proliferation that reproduced the experimentally observed CRAD patterns. Our model indicated that the CRAD may be determined by the ratio of emerging to pre-existing species during non-uniform spatially competitive proliferation, independent of species composition.