Power law analysis of the human microbiome

Temporal turnover in species' ranks can explain variation in Taylor's slope for ecological timeseries

The scaling exponent relating the mean and variance of the density of individual organisms in space (i.e., Taylor's slope: zspace) is well studied in ecology, but the analogous scaling exponent for temporal datasets (ztime) is underdeveloped. Previous theory suggests the narrow distribution of ztime (e.g., typically 1-2) could be due to interspecific competition. Here, using 1694 communities time series, we show that ztime can exceed 2, and reaffirm how this can affect our inference about the stabilizing effect of biodiversity. We also develop a new theory, based on temporal change in the ranks of species abundances, to help account for the observed ztime distribution. Specifically, we find that communities with minimal turnover in species' rank abundances are more likely to have higher ztime. Our analysis shows how species-level variability affects our inference about the stability of ecological communities.

Temporal dynamics and metagenomics of phosphorothioate epigenomes in the human gut microbiome

Background

Epigenetic regulation of gene expression and host defense is well established in microbial communities, with dozens of DNA modifications comprising the epigenomes of prokaryotes and bacteriophage. Phosphorothioation (PT) of DNA, in which a chemically-reactive sulfur atom replaces a non-bridging oxygen in the sugar-phosphate backbone, is catalyzed by dnd and ssp gene families widespread in bacteria and archaea. However, little is known about the role of PTs or other microbial epigenetic modifications in the human microbiome. Here we optimized and applied fecal DNA extraction, mass spectrometric, and metagenomics technologies to characterize the landscape and temporal dynamics of gut microbes possessing PT modifications.

Results

Exploiting the nuclease-resistance of PTs, mass spectrometric analysis of limit digests of PT-containing DNA reveals PT dinucleotides as part of genomic consensus sequences, with 16 possible dinucleotide combinations. Analysis of mouse fecal DNA revealed a highly uniform spectrum of 11 PT dinucleotides in all littermates, with PTs estimated to occur in 5-10% of gut microbes. Though at similar levels, PT dinucleotides in fecal DNA from 11 healthy humans possessed signature combinations and levels of individual PTs. Comparison with a widely distributed microbial epigenetic mark, m6dA, suggested temporal dynamics consistent with expectations for gut microbial communities based on Taylor's Power Law. Application of PT-seq for site-specific metagenomic analysis of PT-containing bacteria in one fecal donor revealed the larger consensus sequences for the PT dinucleotides in Bacteroidota, Firmicutes, Actinobacteria, and Proteobacteria, which differed from unbiased metagenomics and suggested that the abundance of PT-containing bacteria did not simply mirror the spectrum of gut bacteria. PT-seq further revealed low abundance PT sites not detected as dinucleotides by mass spectrometry, attesting to the complementarity of the technologies.

Conclusions

The results of our studies provide a benchmark for understanding the behavior of an abundant and chemically-reactive epigenetic mark in the human gut microbiome, with implications for inflammatory conditions of the gut.

Towards a unified medical microbiome ecology of the OMU for metagenomes and the OTU for microbes

BACKGROUND

Metagenomic sequencing technologies offered unprecedented opportunities and also challenges to microbiology and microbial ecology particularly. The technology has revolutionized the studies of microbes and enabled the high-profile human microbiome and earth microbiome projects. The terminology-change from microbes to microbiomes signals that our capability to count and classify microbes (microbiomes) has achieved the same or similar level as we can for the biomes (macrobiomes) of plants and animals (macrobes). While the traditional investigations of macrobiomes have usually been conducted through naturalists' (Linnaeus & Darwin) naked eyes, and aerial and satellite images (remote-sensing), the large-scale investigations of microbiomes have been made possible by DNA-sequencing-based metagenomic technologies. Two major types of metagenomic sequencing technologies-amplicon sequencing and whole-genome (shotgun sequencing)-respectively generate two contrastingly different categories of metagenomic reads (data)-OTU (operational taxonomic unit) tables representing microorganisms and OMU (operational metagenomic unit), a new term coined in this article to represent various cluster units of metagenomic genes.

RESULTS

The ecological science of microbiomes based on the OTU representing microbes has been unified with the classic ecology of macrobes (macrobiomes), but the unification based on OMU representing metagenomes has been rather limited. In a previous series of studies, we have demonstrated the applications of several classic ecological theories (diversity, composition, heterogeneity, and biogeography) to the studies of metagenomes. Here I push the envelope for the unification of OTU and OMU again by demonstrating the applications of metacommunity assembly and ecological networks to the metagenomes of human gut microbiomes. Specifically, the neutral theory of biodiversity (Sloan's near neutral model), Ning et al.stochasticity framework, core-periphery network, high-salience skeleton network, special trio-motif, and positive-to-negative ratio are applied to analyze the OMU tables from whole-genome sequencing technologies, and demonstrated with seven human gut metagenome datasets from the human microbiome project.

CONCLUSIONS

All of the ecological theories demonstrated previously and in this article, including diversity, composition, heterogeneity, stochasticity, and complex network analyses, are equally applicable to OMU metagenomic analyses, just as to OTU analyses. Consequently, I strongly advocate the unification of OTU/OMU (microbiomes) with classic ecology of plants and animals (macrobiomes) in the context of medical ecology.

A new hypothesis on BV etiology: dichotomous and crisscrossing categorization of complex versus simple on healthy versus BV vaginal microbiomes

IMPORTANCE

BV may influence as many as one-third of women, but its etiology remains unclear. A traditional view is that dominance by Lactobacillus is the hallmark of a healthy vaginal microbiome and lack of dominance may make women BV-prone. Recent studies show that the human VMs can be classified into five major types, four of which possess type-specific dominant species of Lactobacillus. The remaining one (type IV) is not dominated by Lactobacillus and contains a handful of strictly anaerobic bacteria. Nevertheless, exceptions to the first hypothesis have been noticed from the very beginning, and there is not a definite relationship, suggested yet, between the five VM types and BV status. Here, we propose and test a novel hypothesis that assumes the existence of four VM types from dichotomous crisscrossing of "complex versus simple (high diversity or low dominance versus low diversity or high dominance)" on "healthy versus BV." Consequently, there are simple BV versus complex BV.

Helminth Communities of Common Fish Species in the Coastal Zone off Crimea: Species Composition, Diversity, and Structure

In this paper, we analyzed the diversity and structure of helminth communities of 12 common fish species from the coastal zone of Crimea. A total of 53 helminth species were found. The total number of parasite species per host fish ranged from 3 to 18. Species richness at the infracommunity and component community levels were from 1.4-4.2 to 1.7-7, respectively. The Brillouin index for the infracommunites was 0.1-1, while the Shannon index for the component communities was 0.3-1.2. Component communities demonstrated a bi- or tri-modal distribution of the parasite prevalence and positive correlations between the prevalence and log-transformed abundance indices, thus following the "core-satellite" conception. Overall, the prevalence and abundance index of the dominant parasite in the component communities ranged from 18 to 80% and from 0.6 to 61.5 ind. per fish, respectively. The structure of the helminth component communities demonstrated good accordance with the nestedness mode where the rarest species occurred in the most diverse infracommunities, while the poorest infracommunities were composed of a few dominating species. More than two-thirds of the studied helminth species had an aggregated distribution indicating well-structured and developed communities. Our data provide a basis for further research and may be used for fish resource monitoring and management.

Interspecific and intraspecific Taylor's laws for frog skin microbes

Amphibians are known to have an abundance of microorganisms colonizing their skin, and these symbionts often protect the host from disease. There are now many comprehensive studies on amphibian skin microbes, but the interspecific and intraspecific abundance distributions (or abundance heterogeneity) of amphibian skin microbes remain unclear. Furthermore, we have a very limited understanding of how the abundance and heterogeneity of microbial communities relate to the body size (or more specifically, skin surface area) of amphibian hosts. In this study, we evaluated the interspecific and intraspecific abundance distribution patterns of amphibian skin microbes and evaluated whether the symbiotic skin microbes of different anuran species share a fundamental heterogeneity scaling parameter. If scaling invariance exists, we hypothesize that a fundamental heterogeneity scaling value also exists. A total of 358 specimens of 10 amphibian host species were collected, and we used Type-I and III Taylor's power law expansions (TPLE) to assess amphibian skin microbial heterogeneity at the community and mixed-species population levels, respectively. The obtained results showed that, at the community scale, a high aggregation of the microbial abundance distribution on the skin barely changed with host size. In a mixed-species population (i.e., a community context), the abundance distribution pattern of mixed microbial species populations also does not change with host size and always remains highly aggregated. These findings suggest that while amphibian skin microbiomes located in different hosts may have different environmental conditions, they share a fundamental heterogeneity scaling parameter, and thus, scale invariance exists. Finally, we found that microhabitat area provided by the host skin is vital to the stability of the symbiotic microbial community.

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).

Influences of Helicobacter pylori infection on diversity, heterogeneity, and composition of human gastric microbiomes across stages of gastric cancer development

BACKGROUND

About a half of the world's population is infected with Helicobacter pylori (H. pylori), but only 1%-3% of them develop gastric cancer. As a primary risk factor for gastric cancer, the relationship between H. pylori infection and gastric microbiome has been a focus in recent years.

MATERIALS AND METHODS

We reanalyze 11 human gastric microbiome datasets with or without H. pylori, covering the healthy control (HC) and four disease stages (chronic gastritis (CG), atrophic gastritis (AG), intestinal metaplasia (IM), and gastric cancer (GC)) of gastric cancer development to quantitatively compare the influences of the H. pylori infection and disease stages on the diversity, heterogeneity, and composition of gastric microbiome. Four medical ecology approaches including (i) diversity analysis with Hill numbers, (ii) heterogeneity analysis with Taylor's power law extensions (TPLE), (iii) diversity scaling analysis with diversity-area relationship (DAR) model, and (iv) shared species analysis were applied to fulfill the data reanalysis.

RESULTS

(i) The influences of H. pylori infection on the species diversity, spatial heterogeneity, and potential diversity of gastric microbiome seem to be more prevalent than the influences of disease stages during gastric cancer development. (ii) The influences of H. pyloriinfection on diversity, heterogeneity, and composition of gastric microbiomes in HC, CG, IM, and GC stages appear more prevalent than those in AG stage.

CONCLUSION

Our study confirmed the impact of H. pylori infection on human gastric microbiomes: The influences of H. pylori infection on the diversity, heterogeneity, and composition of gastric microbiomes appear to be disease-stage dependent.

Ecology of Human Medical Enterprises: From Disease Ecology of Zoonoses, Cancer Ecology Through to Medical Ecology of Human Microbiomes

In nature, the interaction between pathogens and their hosts is only one of a handful of interaction relationships between species, including parasitism, predation, competition, symbiosis, commensalism, and among others. From a non-anthropocentric view, parasitism has relatively fewer essential differences from the other relationships; but from an anthropocentric view, parasitism and predation against humans and their well-beings and belongings are frequently related to heinous diseases. Specifically, treating (managing) diseases of humans, crops and forests, pets, livestock, and wildlife constitute the so-termed medical enterprises (sciences and technologies) humans endeavor in biomedicine and clinical medicine, veterinary, plant protection, and wildlife conservation. In recent years, the significance of ecological science to medicines has received rising attentions, and the emergence and pandemic of COVID-19 appear accelerating the trend. The facts that diseases are simply one of the fundamental ecological relationships in nature, and the study of the relationships between species and their environment is a core mission of ecology highlight the critical importance of ecological science. Nevertheless, current studies on the ecology of medical enterprises are highly fragmented. Here, we (i) conceptually overview the fields of disease ecology of wildlife, cancer ecology and evolution, medical ecology of human microbiome-associated diseases and infectious diseases, and integrated pest management of crops and forests, across major medical enterprises. (ii) Explore the necessity and feasibility for a unified medical ecology that spans biomedicine, clinical medicine, veterinary, crop (forest and wildlife) protection, and biodiversity conservation. (iii) Suggest that a unified medical ecology of human diseases is both necessary and feasible, but laissez-faire terminologies in other human medical enterprises may be preferred. (iv) Suggest that the evo-eco paradigm for cancer research can play a similar role of evo-devo in evolutionary developmental biology. (v) Summarized 40 key ecological principles/theories in current disease-, cancer-, and medical-ecology literatures. (vi) Identified key cross-disciplinary discovery fields for medical/disease ecology in coming decade including bioinformatics and computational ecology, single cell ecology, theoretical ecology, complexity science, and the integrated studies of ecology and evolution. Finally, deep understanding of medical ecology is of obvious importance for the safety of human beings and perhaps for all living things on the planet.

Variance in Landscape Connectivity Shifts Microbial Population Scaling

Understanding mechanisms shaping distributions and interactions of soil microbes is essential for determining their impact on large scale ecosystem services, such as carbon sequestration, climate regulation, waste decomposition, and nutrient cycling. As the functional unit of soil ecosystems, we focus our attention on the spatial structure of soil macroaggregates. Emulating this complex physico-chemical environment as a patchy habitat landscape we investigate on-chip the effect of changing the connectivity features of this landscape as Escherichia coli forms a metapopulation. We analyze the distributions of E. coli occupancy using Taylor's law, an empirical law in ecology which asserts that the fluctuations in populations is a power law function of the mean. We provide experimental evidence that bacterial metapopulations in patchy habitat landscapes on microchips follow this law. Furthermore, we find that increased variance of patch-corridor connectivity leads to a qualitative transition in the fluctuation scaling. We discuss these results in the context of the spatial ecology of microbes in soil.

Bacterial low-abundant taxa are key determinants of a healthy airway metagenome in the early years of human life

The default removal of low-abundance (rare) taxa from microbial community analyses may lead to an incomplete picture of the taxonomic and functional microbial potential within the human habitat. Publicly available shotgun metagenomics data of healthy children and children with cystic fibrosis (CF) were reanalysed to study the development of the rare species biosphere, which was here defined by either the 15th, 25th or 35th species abundance percentile. We found that healthy children contained an age-independent network of abundant (core) and rare species with both entities being essential in maintaining the network structure. The protein sequence usage for more than 100 bacterial metabolic pathways differed between the core and rare species biosphere. In CF children, the background structure was underdeveloped and random forest bootstrapping based on all constituents of the early airway metagenome and host-associated factors indicated that rare taxa were the most important variables in deciding whether a child was healthy or suffered from the life-limiting CF disease. Attempts failed to make the age-independent CF network as robust as the healthy structure when an increasing number of bacterial taxa from the healthy network was incorporated into the CF structure by computer-based model simulations. However, the transfer of a key combination of taxa from the healthy to the CF network structure with high species diversity and low species dominance, correlated with a more robust CF network and a topological approximation of CF and healthy graph structures. Rothia mucilaginosa, Streptococci and rare species were essential in improving the underdeveloped CF network.

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.

An online time-to-event dashboard comparing the effective control of COVID-19 among continents using the inflection point on an ogive curve: Observational study

BACKGROUND

During the COVID-19 pandemic, one of the frequently asked questions is which countries (or continents) are severely hit. Aside from using the number of confirmed cases and the fatality to measure the impact caused by COVID-19, few adopted the inflection point (IP) to represent the control capability of COVID-19. How to determine the IP days related to the capability is still unclear. This study aims to (i) build a predictive model based on item response theory (IRT) to determine the IP for countries, and (ii) compare which countries (or continents) are hit most.

METHODS

We downloaded COVID-19 outbreak data of the number of confirmed cases in all countries as of October 19, 2020. The IRT-based predictive model was built to determine the pandemic IP for each country. A model building scheme was demonstrated to fit the number of cumulative infected cases. Model parameters were estimated using the Solver add-in tool in Microsoft Excel. The absolute advantage coefficient (AAC) was computed to track the IP at the minimum of incremental points on a given ogive curve. The time-to-event analysis (a.k.a. survival analysis) was performed to compare the difference in IPs among continents using the area under the curve (AUC) and the respective 95% confidence intervals (CIs). An online comparative dashboard was created on Google Maps to present the epidemic prediction for each country.

RESULTS

The top 3 countries that were hit severely by COVID-19 were France, Malaysia, and Nepal, with IP days at 263, 262, and 262, respectively. The top 3 continents that were hit most based on IP days were Europe, South America, and North America, with their AUCs and 95% CIs at 0.73 (0.61-0.86), 0.58 (0.31-0.84), and 0.54 (0.44-0.64), respectively. An online time-event result was demonstrated and shown on Google Maps, comparing the IP probabilities across continents.

CONCLUSION

An IRT modeling scheme fitting the epidemic data was used to predict the length of IP days. Europe, particularly France, was hit seriously by COVID-19 based on the IP days. The IRT model incorporated with AAC is recommended to determine the pandemic IP.

Heterogeneity-disease relationship in the human microbiome-associated diseases

Space is a critical and also challenging frontier in human microbiome research. It has been found that lack of consideration of scales beyond individual and ignoring of microbe dispersal are two crucial roadblocks in preventing deep understanding of the spatial heterogeneity of human microbiome. Assessing and interpreting the heterogeneity and dispersal in microbiomes explicitly are particularly challenging, but implicit approaches such as Taylor's power law (TPL) can be rather effective. Based on TPL, which achieved a rare status of ecological laws, we introduce a general methodology for characterizing the spatial heterogeneity of microbiome (i.e. characterization of microbial spatial distribution) and further apply it for investigating the heterogeneity-disease relationship (HDR) via analyzing a big dataset of 26 MAD (microbiome-associated disease) studies covering nearly all high-profile MADs including obesity, diabetes and gout. It was found that in majority of the MAD cases, the microbiome was sufficiently resilient to endure the disease disturbances. Specifically, in ∼10-16% cases, disease effects were significant-the healthy and diseased cohorts exhibited statistically significant differences in the TPL heterogeneity parameters. We further compared HDR with classic diversity-disease relationship (DDR) and explained their mechanistic differences. Both HDR and DDR cross-verified remarkable resilience of the human microbiomes against MADs.

Predicting the Outbreak Risks and Inflection Points of COVID-19 Pandemic with Classic Ecological Theories

Predicting the outbreak risks and/or the inflection (turning or tipping) points of COVID-19 can be rather challenging. Here, it is addressed by modeling and simulation approaches guided by classic ecological theories and by treating the COVID-19 pandemic as a metapopulation dynamics problem. Three classic ecological theories are harnessed, including TPL (Taylor's power-law) and Ma's population aggregation critical density (PACD) for spatiotemporal aggregation/stability scaling, approximating virus metapopulation dynamics with Hubbell's neutral theory, and Ma's diversity-time relationship adapted for the infection-time relationship. Fisher-Information for detecting critical transitions and tipping points are also attempted. It is discovered that: (i) TPL aggregation/stability scaling parameter (b > 2), being significantly higher than the b-values of most macrobial and microbial species including SARS, may interpret the chaotic pandemic of COVID-19. (ii) The infection aggregation critical threshold (M 0) adapted from PACD varies with time (outbreak-stage), space (region) and public-health interventions. Exceeding M 0, local contagions may become aggregated and connected regionally, leading to epidemic/pandemic. (iii) The ratio of fundamental dispersal to contagion numbers can gauge the relative importance between local contagions vs. regional migrations in spreading infections. (iv) The inflection (turning) points, pair of maximal infection number and corresponding time, are successfully predicted in more than 80% of Chinese provinces and 68 countries worldwide, with a precision >80% generally.

Does Ulcerative Colitis Influence the Inter-individual Heterogeneity of the Human Intestinal Mucosal Microbiome?

The dysbiosis of the gut microbiome associated with ulcerative colitis (UC) has been extensively studied in recent years. However, the question of whether UC influences the spatial heterogeneity of the human gut mucosal microbiome has not been addressed. Spatial heterogeneity (specifically, the inter-individual heterogeneity in microbial species abundances) is one of the most important characterizations at both population and community scales, and can be assessed and interpreted by Taylor’s power law (TPL) and its community-scale extensions (TPLEs). Due to the high mobility of microbes, it is difficult to investigate their spatial heterogeneity explicitly; however, TPLE offers an effective approach to implicitly analyze the microbial communities. Here, we investigated the influence of UC on the spatial heterogeneity of the gut microbiome with intestinal mucosal microbiome samples collected from 28 UC patients and healthy controls. Specifically, we applied Type-I TPLE for measuring community spatial heterogeneity and Type-III TPLE for measuring mixed-species population heterogeneity to evaluate the heterogeneity changes of the mucosal microbiome induced by UC at both the community and species scales. We further used permutation test to determine the possible differences between UC patients and healthy controls in heterogeneity scaling parameters. Results showed that UC did not significantly influence gut mucosal microbiome heterogeneity at either the community or mixed-species levels. These findings demonstrated significant resilience of the human gut microbiome and confirmed a prediction of TPLE: that the inter-subject heterogeneity scaling parameter of the gut microbiome is an intrinsic property to humans, invariant with UC disease.

Critical Network Structures and Medical Ecology Mechanisms Underlying Human Microbiome-Associated Diseases

A fundamental problem in studies on human microbiome-associated diseases (MADs) is to understand the relationships between microbiome structures and health status of hosts. For example, species diversity metrics have been routinely evaluated in virtually all studies on MADs, yet a recent meta-analysis revealed that, in only approximately one-third of the cases, diversity and diseases were related. In this study, we ask whether Hubbell's neutral theory (supplemented with the normalized stochasticity ratio [NSR]) or critical microbiome network structures may offer better alternatives. Whereas neutral theory and NSR focus on stochastic processes, we use core/periphery and high-salience skeleton networks to evaluate deterministic, asymmetrical niche effects, assuming that all species or their interactions were not “born” equal and focusing on non-neutral, critical network structures. We found that properties of critical network structures are more indicative of disease effects. Finally, seven findings (mechanisms, interpretations, and postulations) regarding medical ecology mechanisms underlying MADs were summarized.

•

Seven findings (mechanisms/interpretations/postulations) of medical ecology proposed

•

Critical network structures more indicative of disease effects than ecology metrics

•

One-third seems ceiling of diversity-disease relations, half to two-thirds of network structures

•

Super resilience (unexplained one-third to half gap) is likely attributed to host genome

Estimating the Optimum Coverage and Quality of Amplicon Sequencing With Taylor's Power Law Extensions

Theoretical analysis of DNA sequencing coverage problem has been investigated with complex mathematical models such as Lander-Waterman expectation theory and Stevens' theorem for randomly covering a domain. In the field of metagenomics sequencing, several approaches have been developed to estimate the coverage of whole-genome shotgun sequencing, but surprisingly few studies addressed the coverage problem for marker-gene amplicon sequencing, for which arguably the biggest challenge is the complexity or heterogeneity of microbial communities. Overall, much of the practice still relies variously on speculation, semi-empirical and ad hoc heuristic models. Conservatively raising coverage may ensure the success of sequencing project, but often with unduly cost. In this study, we borrow the principles and approaches of optimum sampling methodology originated in applied entomology, achieved equal success in plant pathology and parasitology, and plays a critical role in the decision-making for global crop and forest protection against economic pests since 1970s when the pesticide crisis and food safety concerns forced the reduction of pesticide usages, which in turn requires reliable sampling techniques for monitoring pest populations. We realized that sequencing coverage is essentially an optimum sampling problem. Perhaps the only essential difference between sampling insects and sampling microbiome is the "instrument" used. In traditional entomology, it is usually humans that visually count the numbers of insects, occasionally aided by binocular microscope. In the metagenomics research, it is the DNA sequencers that count the number of DNA reads. Furthermore, a key theoretical foundation for sampling insect pest populations, i.e., Taylor's power law, which achieved rare status of ecological law and captures the population aggregation, has been recently extended to the community level for describing community heterogeneity and stability, namely, Taylor's power law extensions (TPLEs). This theoretical advance enabled us to develop a novel approach to assessing the quality and determining optimum reads (coverage) of amplicon sequencing operations. Specifically, two applications were developed: one is, in hindsight, to assess the quality of amplicon sequencing operation in terms of the precision and confidence levels. Another is, prior to sequencing operation, to determine the minimum sequencing efforts for a sequencing project to achieve preset precision and confidence levels.

Assessing and Interpreting the Metagenome Heterogeneity With Power Law

There are two major sequencing technologies for investigating the microbiome: the amplicon sequencing that generates the OTU (operational taxonomic unit) tables of marker genes (e.g., bacterial 16S-rRNA), and the metagenomic shotgun sequencing that generates metagenomic gene abundance (MGA) tables. The OTU table is the counterpart of species abundance tables in macrobial ecology of plants and animals, and has been the target of numerous ecological and network analyses in recent gold rush for microbiome research and in great efforts for establishing an inclusive theoretical ecology. Nevertheless, MGA analyses have been largely limited to bioinformatics pipelines and ad hoc statistical methods, and systematic approaches to MGAs guided by classic ecological theories are still few. Here, we argue that, the difference between “gene kinds” and “gene species” are nominal, and the metagenome that a microbiota carries is essentially a ‘community’ of metagenomic genes (MGs). Each row of a MGA table represents a metagenome of a microbiota, and the whole MGA table represents a ‘meta-metagenome’ (or an assemblage of metagenomes) of N microbiotas (microbiome samples). Consequently, the same ecological/network analyses used in OTU analyses should be equally applicable to MGA tables. Here we choose to analyze the heterogeneity of metagenome by introducing classic Taylor’s power law (TPL) and its recent extensions in community ecology. Heterogeneity is a fundamental property of metagenome, particularly in the context of human microbiomes. Recent studies have shown that the heterogeneity of human metagenomes is far more significant than that of human genomes. Therefore, without deep understanding of the human metagenome heterogeneity, personalized medicine of the human microbiome-associated diseases is hardly feasible. The TPL extensions have been successfully applied to measure the heterogeneity of human microbiome based on amplicon-sequencing reads of marker genes (e.g., 16s-rRNA). In this article, we demonstrate the analysis of the metagenomic heterogeneity of human gut microbiome at whole metagenome scale (with type-I power law extension) and metagenomic gene scale (type-III), as well as the heterogeneity of gene clusters, respectively. We further examine the influences of obesity, IBD and diabetes on the heterogeneity, which is of important ramifications for the diagnosis and treatment of human microbiome-associated diseases.

Macroecological dynamics of gut microbiota

The gut microbiota is now widely recognized as a dynamic ecosystem that plays an important role in health and disease. Although current sequencing technologies make it possible to explore how relative abundances of host-associated bacteria change over time, the biological processes governing microbial dynamics remain poorly understood. Therefore, as in other ecological systems, it is important to identify quantitative relationships describing various aspects of gut microbiota dynamics. In the present study, we use multiple high-resolution time series data obtained from humans and mice to demonstrate that, despite their inherent complexity, gut microbiota dynamics can be characterized by several robust scaling relationships. Interestingly, the observed patterns are highly similar to those previously identified across diverse ecological communities and economic systems, including the temporal fluctuations of animal and plant populations and the performance of publicly traded companies. Specifically, we find power-law relationships describing short- and long-term changes in gut microbiota abundances, species residence and return times, and the correlation between the mean and the temporal variance of species abundances. The observed scaling laws are altered in mice receiving different diets and are affected by context-specific perturbations in humans. We use the macroecological relationships to reveal specific bacterial taxa, the dynamics of which are substantially perturbed by dietary and environmental changes. Overall, our results suggest that a quantitative macroecological framework will be important for characterizing and understanding the complex dynamics of diverse microbial communities.

How and Why Men and Women Differ in Their Microbiomes: Medical Ecology and Network Analyses of the Microgenderome

Microgenderome or sexual dimorphism in microbiome refers to the bidirectional interactions between microbiotas, sex hormones, and immune systems, and it is highly relevant to disease susceptibility. A critical step in exploring microgenderome is to dissect the sex differences in key community ecology properties, which has not been systematically analyzed. This study aims at filling the gap by reanalyzing the Human Microbiome Project datasets with two objectives: (i) dissecting the sex differences in community diversity and their intersubject scaling, species composition, core/periphery species, and high-salience skeletons (species interactions); (ii) offering mechanistic interpretations for (i). Conceptually, the Vellend-Hanson synthesis of community ecology that stipulates selection, drift, speciation, and dispersal as the four processes driving community dynamics is followed. Methodologically, seven approaches reflecting the state-of-the-art research in medical ecology of human microbiomes are harnessed to achieve the objectives. It is postulated that the revealed microgenderome characteristics (categorized as seven aspects of differences/similarities) exert far reaching influences on disease susceptibility, and are primarily due to the sex difference in selection effects (deterministic fitness differences in microbial species and/or species interactions with each other or with their hosts), which are, in turn, shaped/modulated by host physiology (immunity, hormones, gut-brain communications, etc.).

The Microbiota-Gut-Brain Axis

The importance of the gut-brain axis in maintaining homeostasis has long been appreciated. However, the past 15 yr have seen the emergence of the microbiota (the trillions of microorganisms within and on our bodies) as one of the key regulators of gut-brain function and has led to the appreciation of the importance of a distinct microbiota-gut-brain axis. This axis is gaining ever more traction in fields investigating the biological and physiological basis of psychiatric, neurodevelopmental, age-related, and neurodegenerative disorders. The microbiota and the brain communicate with each other via various routes including the immune system, tryptophan metabolism, the vagus nerve and the enteric nervous system, involving microbial metabolites such as short-chain fatty acids, branched chain amino acids, and peptidoglycans. Many factors can influence microbiota composition in early life, including infection, mode of birth delivery, use of antibiotic medications, the nature of nutritional provision, environmental stressors, and host genetics. At the other extreme of life, microbial diversity diminishes with aging. Stress, in particular, can significantly impact the microbiota-gut-brain axis at all stages of life. Much recent work has implicated the gut microbiota in many conditions including autism, anxiety, obesity, schizophrenia, Parkinson’s disease, and Alzheimer’s disease. Animal models have been paramount in linking the regulation of fundamental neural processes, such as neurogenesis and myelination, to microbiome activation of microglia. Moreover, translational human studies are ongoing and will greatly enhance the field. Future studies will focus on understanding the mechanisms underlying the microbiota-gut-brain axis and attempt to elucidate microbial-based intervention and therapeutic strategies for neuropsychiatric disorders.

Comparative power law analysis for the spatial heterogeneity scaling of the hot-spring microbiomes

Spatial heterogeneity is a fundamental property of any natural ecosystems, including hot spring and human microbiomes. Two important scales that spatial heterogeneity exhibits are population and community scales, and Taylor's power law (PL) and its extensions (PLEs) offer ideal quantitative models to assess population- and community-level heterogeneities. Here we analyse 165 hot spring microbiome samples at the global scale that cover a wide range of temperatures (7.5-99°C) and pH levels (3.3-9). We explore a question of fundamental importance for measuring the spatial heterogeneity of the hot-spring microbiome and further discuss their ecological implications: How do critical environmental factors such as temperature and pH influence the scaling of community spatial heterogeneity? We are particularly interested in the existence of a universal scaling model that is independent of environmental gradients. By applying PL and PLEs, we were able to obtain such scaling parameters of the hot spring at both community and population levels, which are temperature- and pH-invariant. These findings suggest that while the hot-spring microbiomes located at different regions may have different environmental conditions, they share a fundamental heterogeneity scaling parameter, analogically similar to the gravitational acceleration on Earth, which may vary slightly depending on altitude and latitude, but is invariant overall. In contrast, similar to the physics of the Moon and Earth, which have different gravitational accelerations, the hot spring and human microbiomes can have different scaling parameters as demonstrated in this study.

Optimal Microbiome Networks: Macroecology and Criticality

The human microbiome is an extremely complex ecosystem considering the number of bacterial species, their interactions, and its variability over space and time. Here, we untangle the complexity of the human microbiome for the Irritable Bowel Syndrome (IBS) that is the most prevalent functional gastrointestinal disorder in human populations. Based on a novel information theoretic network inference model, we detected potential species interaction networks that are functionally and structurally different for healthy and unhealthy individuals. Healthy networks are characterized by a neutral symmetrical pattern of species interactions and scale-free topology versus random unhealthy networks. We detected an inverse scaling relationship between species total outgoing information flow, meaningful of node interactivity, and relative species abundance (RSA). The top ten interacting species are also the least relatively abundant for the healthy microbiome and the most detrimental. These findings support the idea about the diminishing role of network hubs and how these should be defined considering the total outgoing information flow rather than the node degree. Macroecologically, the healthy microbiome is characterized by the highest Pareto total species diversity growth rate, the lowest species turnover, and the smallest variability of RSA for all species. This result challenges current views that posit a universal association between healthy states and the highest absolute species diversity in ecosystems. Additionally, we show how the transitory microbiome is unstable and microbiome criticality is not necessarily at the phase transition between healthy and unhealthy states. We stress the importance of considering portfolios of interacting pairs versus single node dynamics when characterizing the microbiome and of ranking these pairs in terms of their interactions (i.e., species collective behavior) that shape transition from healthy to unhealthy states. The macroecological characterization of the microbiome is useful for public health and disease diagnosis and etiognosis, while species-specific analyses can detect beneficial species leading to personalized design of pre- and probiotic treatments and microbiome engineering.

Sketching the Human Microbiome Biogeography with DAR (Diversity-Area Relationship) Profiles

SAR (species area relationship) is a classic ecological theory that has been extensively investigated and applied in the studies of global biogeography and biodiversity conservation in macro-ecology. It has also found important applications in microbial ecology in recent years thanks to the breakthroughs in metagenomic sequencing technology. Nevertheless, SAR has a serious limitation for practical applications-ignoring the species abundance and treating all species as equally abundant. This study aims to explore the biogeography discoveries of human microbiome over 18 sites of 5 major microbiome habitats, establish the baseline DAR (diversity-area scaling relationship) parameters, and perform comparisons with the classic SAR. The extension from SAR to DAR by adopting the Hill numbers as diversity measures not only overcomes the previously mentioned flaw of SAR but also allows for obtaining a series of important findings on the human microbiome biodiversity and biogeography. Specifically, two types of DAR models were built, the traditional power law (PL) and power law with exponential cutoff (PLEC), using comprehensive datasets from the HMP (human microbiome project). Furthermore, the biogeography "maps" for 18 human microbiome sites using their DAR profiles for assessing and predicting the diversity scaling across individuals, PDO profiles (pair-wise diversity overlap) for measuring diversity overlap (similarity), and MAD profile (for predicting the maximal accrual diversity in a population) were sketched out. The baseline biogeography maps for the healthy human microbiome diversity can offer guidelines for conserving human microbiome diversity and investigating the health implications of the human microbiome diversity and heterogeneity.

The gut microbiota heterogeneity and assembly changes associated with the IBD

Inflammatory bowel disease (IBD) is an immunologically mediated disease and may be caused by abnormal immunological response to gut microbes. Although several studies on the ecological changes associated with IBD, such as community diversities, were reported, no previous studies have investigated the changes in the spatial heterogeneity and the mechanism of community assembly of the gut microbiota associated with IBD. In the present study, we first applied the Taylor’s power law extensions to compare the community spatial heterogeneity between the gut microbial communities of the IBD patients and those of the healthy individuals. We found that the community spatial heterogeneity of gut microbiota in IBD patients is slightly lower than in the healthy individuals. This finding suggests that IBD may lower the spatial heterogeneity of gut microbiota, possibly via lowering the abundance of dominant species. We further applied the neutral theory of biodiversity to comparatively investigate the community assembly and diversity maintenance of the gut microbiota with and without IBD, and our application suggested that deterministic factors such as host immunity should be dominant forces shaping gut microbiota assembly, and diseases such as IBD may not be strong enough to change the trend set by the deterministic host factors.

DAR (diversity-area relationship): Extending classic SAR (species-area relationship) for biodiversity and biogeography analyses

I extend the classic SAR, which has achieved status of ecological law and plays a critical role in global biodiversity and biogeography analyses, to general DAR (diversity-area relationship). The extension was aimed to remedy a serious application limitation of the traditional SAR that only addressed one aspect of biodiversity scaling-species richness scaling over space, but ignoring species abundance information. The extension was further inspired by a recent consensus that Hill numbers offer the most appropriate measures for alpha-diversity and multiplicative beta-diversity. In particular, Hill numbers are essentially a series of Renyi's entropy values weighted differently along the rare-common-dominant spectrum of species abundance distribution and are in the units of effective number of species (or species equivalents such as OTUs). I therefore postulate that Hill numbers should follow the same or similar law of the traditional SAR. I test the postulation with the American gut microbiome project (AGP) dataset of 1,473 healthy North American individuals. I further propose three new concepts and develop their statistical estimation formulae based on the new DAR extension, including: (i) DAR profile-z-q relationship (DAR scaling parameter z at different diversity order q), (ii) PDO (pair-wise diversity overlap) profile-g-q relationship (PDO parameter g at order q, and (iii) MAD (maximal accrual diversity: D max) profile-D max-q. While the classic SAR is a special case of our new DAR profile, the PDO and MAD profiles offer novel tools for analyzing biodiversity (including alpha-diversity and beta-diversity) and biogeography over space.

Assessing and Interpreting the Within-Body Biogeography of Human Microbiome Diversity

A human body hosts a relatively independent microbiome including five major regional biomes (i.e., airway, oral, gut, skin, and urogenital). Each of them may possess different regional characteristics with important implications to our health and diseases (i.e., so-termed microbiome associated diseases). Nevertheless, these regional microbiomes are connected with each other through diffusions and migrations. Here, we investigate the within-body (intra-individual) distribution feature of microbiome diversity via diversity area relationship (DAR) modeling, which, to the best of our knowledge, has not been systematically studied previously. We utilized the Hill numbers for measuring alpha and beta-diversities and built 1,200 within-body DAR models with to date the most comprehensive human microbiome datasets of 18 sites from the human microbiome project (HMP) cohort. We established the intra-DAR profile (z-q pattern: the diversity scaling parameter z of the power law (PL) at diversity order q = 0-3), intra-PDO (pair-wise diversity overlap) profile (g-q), and intra-MAD (maximal accrual diversity) profile (D max-q) for the within-body biogeography of the human microbiome. These profiles constitute the "maps" of the within-body biogeography, and offer important insights on the within-body distribution of the human microbiome. Furthermore, we investigated the heterogeneity among individuals in their biogeography parameters and found that there is not an "average Joe" that can represent majority of individuals in a cohort or population. For example, we found that most individuals in the HMP cohort have relatively lower maximal accrual diversity (MAD) or in the "long tail" of the so-termed power law distribution. In the meantime, there are a small number of individuals in the cohort who possess disproportionally higher MAD values. These findings may have important implications for personalized medicine of the human microbiome associated diseases in practice, besides their theoretical significance in microbiome research such as establishing the baseline for the conservation of human microbiome.

Diversity time-period and diversity-time-area relationships exemplified by the human microbiome

We extend the ecological laws of species-time relationship (STR) and species-time-area relationship (STAR) to general diversity time-period relationship (DTR) and diversity-time-area relationship (DTAR), and test the extensions with the human vaginal microbiome datasets by building 1460 DTR/DTAR models. Our extensions were inspired by the observation that Hill numbers, well regarded as the most appropriate measure of alpha-diversity and also particularly suitable for multiplicative beta-diversity partitioning, are actually in the units of effective species, and therefore, should be able to substitute for species in the STR and STAR. We found that the traditional power law (PL) model is only applicable for DTR at diversity order zero (i.e., species richness); at higher diversity orders (q = 1–4), the power law with exponent cutoff (PLEC) and power law with inverse exponent cutoff (PLIEC) are more appropriate. In particular, PLEC has an advantage over PLIEC in predicting maximal accumulation diversity (MAD) over time. In fact, with the DTR extensions, we can construct DTR and MAD profiles. To the best of our knowledge, this is the first comprehensive investigation of the DTR/DTAR in human microbiome. Methodologically, our DTR/DTAR profiles can characterize general diversity scaling beyond species richness, covering both alpha- and beta-diversity regimes across different diversity orders.

The P/N (Positive-to-Negative Links) Ratio in Complex Networks-A Promising In Silico Biomarker for Detecting Changes Occurring in the Human Microbiome

Relatively little progress in the methodology for differentiating between the healthy and diseased microbiomes, beyond comparing microbial community diversities with traditional species richness or Shannon index, has been made. Network analysis has increasingly been called for the task, but most currently available microbiome datasets only allows for the construction of simple species correlation networks (SCNs). The main results from SCN analysis are a series of network properties such as network degree and modularity, but the metrics for these network properties often produce inconsistent evidence. We propose a simple new network property, the P/N ratio, defined as the ratio of positive links to the number of negative links in the microbial SCN. We postulate that the P/N ratio should reflect the balance between facilitative and inhibitive interactions among microbial species, possibly one of the most important changes occurring in diseased microbiome. We tested our hypothesis with five datasets representing five major human microbiome sites and discovered that the P/N ratio exhibits contrasting differences between healthy and diseased microbiomes and may be harnessed as an in silico biomarker for detecting disease-associated changes in the human microbiome, and may play an important role in personalized diagnosis of the human microbiome-associated diseases.

A large-scale survey of the postmortem human microbiome, and its potential to provide insight into the living health condition

The microbiome plays many roles in human health, often through the exclusive lens of clinical interest. The inevitable end point for all living hosts, death, has its own altered microbiome configurations. However, little is understood about the ecology and changes of microbial communities after death, or their potential utility for understanding the health condition of the recently living. Here we reveal distinct postmortem microbiomes of human hosts from a large-scale survey of death cases representing a predominantly urban population, and demonstrated these microbiomes reflected antemortem health conditions within 24–48 hours of death. Our results characterized microbial community structure and predicted function from 188 cases representing a cross-section of an industrial-urban population. We found strong niche differentiation of anatomic habitat and microbial community turnover based on topographical distribution. Microbial community stability was documented up to two days after death. Additionally, we observed a positive relationship between cell motility and time since host death. Interestingly, we discovered evidence that microbial biodiversity is a predictor of antemortem host health condition (e.g., heart disease). These findings improve the understanding of postmortem host microbiota dynamics, and provide a robust dataset to test the postmortem microbiome as a tool for assessing health conditions in living populations.

Temporal Stability of the Human Skin Microbiome

Biogeography and individuality shape the structural and functional composition of the human skin microbiome. To explore these factors' contribution to skin microbial community stability, we generated metagenomic sequence data from longitudinal samples collected over months and years. Analyzing these samples using a multi-kingdom, reference-based approach, we found that despite the skin's exposure to the external environment, its bacterial, fungal, and viral communities were largely stable over time. Site, individuality, and phylogeny were all determinants of stability. Foot sites exhibited the most variability; individuals differed in stability; and transience was a particular characteristic of eukaryotic viruses, which showed little site-specificity in colonization. Strain and single-nucleotide variant-level analysis showed that individuals maintain, rather than reacquire, prevalent microbes from the environment. Longitudinal stability of skin microbial communities generates hypotheses about colonization resistance and empowers clinical studies exploring alterations observed in disease states.