Pan-cancer analyses reveal cancer-type-specific fungal ecologies and bacteriome interactions.

Cancer-microbe associations have been explored for centuries, but cancer-associated fungi have rarely been examined. Here, we comprehensively characterize the cancer mycobiome within 17,401 patient tissue, blood, and plasma samples across 35 cancer types in four independent cohorts. We report fungal DNA and cells at low abundances across many major human cancers, with differences in community compositions that differ among cancer types, even when accounting for technical background. Fungal histological staining of tissue microarrays supported intratumoral presence and frequent spatial association with cancer cells and macrophages. Comparing intratumoral fungal communities with matched bacteriomes and immunomes revealed co-occurring bi-domain ecologies, often with permissive, rather than competitive, microenvironments and distinct immune responses. Clinically focused assessments suggested prognostic and diagnostic capacities of the tissue and plasma mycobiomes, even in stage I cancers, and synergistic predictive performance with bacteriomes.

Individualized network analysis reveals a link between the gut microbiome, diet intervention and Gestational Diabetes Mellitus.

Gestational Diabetes Mellitus (GDM), a serious complication during pregnancy which is defined by abnormal glucose regulation, is commonly treated by diabetic diet and lifestyle changes. While recent findings place the microbiome as a natural mediator between diet interventions and diverse disease states, its role in GDM is still unknown. Here, based on observation data from healthy pregnant control group and GDM patients, we developed a new network approach using patterns of co-abundance of microorganism to construct microbial networks that represent human-specific information about gut microbiota in different groups. By calculating network similarity in different groups, we analyze the gut microbiome from 27 GDM subjects collected before and after two weeks of diet therapy compared with 30 control subjects to identify the health condition of microbial community balance in GDM subjects. Although the microbial communities remain similar after the diet phase, we find that the structure of their inter-species co-abundance network is significantly altered, which is reflected in that the ecological balance of GDM patients was not "healthier" after the diet intervention. In addition, we devised a method for individualized network analysis of the microbiome, thereby a pattern is found that GDM individuals whose microbial networks are with large deviations from the GDM group are usually accompanied by their abnormal glucose regulation. This approach may help the development of individualized diagnosis strategies and microbiome-based therapies in the future.

Universality of human microbial dynamics.

Human-associated microbial communities have a crucial role in determining our health and well-being, and this has led to the continuing development of microbiome-based therapies such as faecal microbiota transplantation. These microbial communities are very complex, dynamic and highly personalized ecosystems, exhibiting a high degree of inter-individual variability in both species assemblages and abundance profiles. It is not known whether the underlying ecological dynamics of these communities, which can be parameterized by growth rates, and intra- and inter-species interactions in population dynamics models, are largely host-independent (that is, universal) or host-specific. If the inter-individual variability reflects host-specific dynamics due to differences in host lifestyle, physiology or genetics, then generic microbiome manipulations may have unintended consequences, rendering them ineffective or even detrimental. Alternatively, microbial ecosystems of different subjects may exhibit universal dynamics, with the inter-individual variability mainly originating from differences in the sets of colonizing species. Here we develop a new computational method to characterize human microbial dynamics. By applying this method to cross-sectional data from two large-scale metagenomic studies--the Human Microbiome Project and the Student Microbiome Project--we show that gut and mouth microbiomes display pronounced universal dynamics, whereas communities associated with certain skin sites are probably shaped by differences in the host environment. Notably, the universality of gut microbial dynamics is not observed in subjects with recurrent Clostridium difficile infection but is observed in the same set of subjects after faecal microbiota transplantation. These results fundamentally improve our understanding of the processes that shape human microbial ecosystems, and pave the way to designing general microbiome-based therapies.

Inferring human microbial dynamics from temporal metagenomics data: Pitfalls and lessons.

The human gut microbiota is a very complex and dynamic ecosystem that plays a crucial role in health and well-being. Inferring microbial community structure and dynamics directly from time-resolved metagenomics data is key to understanding the community ecology and predicting its temporal behavior. Many methods have been proposed to perform the inference. Yet, as we point out in this review, there are several pitfalls along the way. Indeed, the uninformative temporal measurements and the compositional nature of the relative abundance data raise serious challenges in inference. Moreover, the inference results can be largely distorted when only focusing on highly abundant species by ignoring or grouping low-abundance species. Finally, the implicit assumptions in various regularization methods may not reflect reality. Those issues have to be seriously considered in ecological modeling of human gut microbiota.

On the Origins and Control of Community Types in the Human Microbiome.

Microbiome-based stratification of healthy individuals into compositional categories, referred to as "enterotypes" or "community types", holds promise for drastically improving personalized medicine. Despite this potential, the existence of community types and the degree of their distinctness have been highly debated. Here we adopted a dynamic systems approach and found that heterogeneity in the interspecific interactions or the presence of strongly interacting species is sufficient to explain community types, independent of the topology of the underlying ecological network. By controlling the presence or absence of these strongly interacting species we can steer the microbial ecosystem to any desired community type. This open-loop control strategy still holds even when the community types are not distinct but appear as dense regions within a continuous gradient. This finding can be used to develop viable therapeutic strategies for shifting the microbial composition to a healthy configuration.

Network impact of a single-time-point microbial sample.

The human microbiome plays a crucial role in determining our well-being and can significantly influence human health. The individualized nature of the microbiome may reveal host-specific information about the health state of the subject. In particular, the microbiome is an ecosystem shaped by a tangled network of species-species and host-species interactions. Thus, analysis of the ecological balance of microbial communities can provide insights into these underlying interrelations. However, traditional methods for network analysis require many samples, while in practice only a single-time-point microbial sample is available in clinical screening. Recently, a method for the analysis of a single-time-point sample, which evaluates its 'network impact' with respect to a reference cohort, has been applied to analyze microbial samples from women with Gestational Diabetes Mellitus. Here, we introduce different variations of the network impact approach and systematically study their performance using simulated 'samples' fabricated via the Generalized Lotka-Volttera model of ecological dynamics. We show that the network impact of a single sample captures the effect of the interactions between the species, and thus can be applied to anomaly detection of shuffled samples, which are 'normal' in terms of species abundance but 'abnormal' in terms of species-species interrelations. In addition, we demonstrate the use of the network impact in binary and multiclass classifications, where the reference cohorts have similar abundance profiles but different species-species interactions. Individualized analysis of the human microbiome has the potential to improve diagnosis and personalized treatments.

Model-free prediction of microbiome compositions.

BACKGROUND: The recent recognition of the importance of the microbiome to the host's health and well-being has yielded efforts to develop therapies that aim to shift the microbiome from a disease-associated state to a healthier one. Direct manipulation techniques of the species' assemblage are currently available, e.g., using probiotics or narrow-spectrum antibiotics to introduce or eliminate specific taxa. However, predicting the species' abundances at the new state remains a challenge, mainly due to the difficulties of deciphering the delicate underlying network of ecological interactions or constructing a predictive model for such complex ecosystems. RESULTS: Here, we propose a model-free method to predict the species' abundances at the new steady state based on their presence/absence configuration by utilizing a multi-dimensional k-nearest-neighbors (kNN) regression algorithm. By analyzing data from numeric simulations of ecological dynamics, we show that our predictions, which consider the presence/absence of all species holistically, outperform both the null model that uses the statistics of each species independently and a predictive neural network model. We analyze real metagenomic data of human-associated microbial communities and find that by relying on a small number of "neighboring" samples, i.e., samples with similar species assemblage, the kNN predicts the species abundance better than the whole-cohort average. By studying both real metagenomic and simulated data, we show that the predictability of our method is tightly related to the dissimilarity-overlap relationship of the training data. CONCLUSIONS: Our results demonstrate how model-free methods can prove useful in predicting microbial communities and may facilitate the development of microbial-based therapies. Video Abstract.

Interactions-based classification of a single microbial sample.

To address the limitation of overlooking crucial ecological interactions due to relying on single time point samples, we developed a computational approach that analyzes individual samples based on the interspecific microbial relationships. We verify, using both numerical simulations as well as real and shuffled microbial profiles from the human oral cavity, that the method can classify single samples based on their interspecific interactions. By analyzing the gut microbiome of people with autistic spectrum disorder, we found that our interaction-based method can improve the classification of individual subjects based on a single microbial sample. These results demonstrate that the underlying ecological interactions can be practically utilized to facilitate microbiome-based diagnosis and precision medicine.

Top-down identification of keystone taxa in the microbiome.

Keystone taxa in ecological communities are native taxa that play an especially important role in the stability of their ecosystem. However, we still lack an effective framework for identifying these taxa from the available high-throughput sequencing without the notoriously difficult step of reconstructing the detailed network of inter-specific interactions. In addition, while most microbial interaction models assume pair-wise relationships, it is yet unclear whether pair-wise interactions dominate the system, or whether higher-order interactions are relevant. Here we propose a top-down identification framework, which detects keystones by their total influence on the rest of the taxa. Our method does not assume a priori knowledge of pairwise interactions or any specific underlying dynamics and is appropriate to both perturbation experiments and metagenomic cross-sectional surveys. When applied to real high-throughput sequencing of the human gastrointestinal microbiome, we detect a set of candidate keystones and find that they are often part of a keystone module - multiple candidate keystone species with correlated occurrence. The keystone analysis of single-time-point cross-sectional data is also later verified by the evaluation of two-time-points longitudinal sampling. Our framework represents a necessary advancement towards the reliable identification of these key players of complex, real-world microbial communities.

Cascading failures in interdependent lattice networks: the critical role of the length of dependency links.

We study the cascading failures in a system composed of two interdependent square lattice networks A and B placed on the same Cartesian plane, where each node in network A depends on a node in network B randomly chosen within a certain distance r from the corresponding node in network A and vice versa. Our results suggest that percolation for small r below r(max)≈8 (lattice units) is a second-order transition, and for larger r is a first-order transition. For r

Complexity-stability trade-off in empirical microbial ecosystems.

May's stability theory, which holds that large ecosystems can be stable up to a critical level of complexity, a product of the number of resident species and the intensity of their interactions, has been a central paradigm in theoretical ecology. So far, however, empirically demonstrating this theory in real ecological systems has been a long-standing challenge with inconsistent results. Especially, it is unknown whether this theory is pertinent in the rich and complex communities of natural microbiomes, mainly due to the challenge of reliably reconstructing such large ecological interaction networks. Here we introduce a computational framework for estimating an ecosystem's complexity without relying on a priori knowledge of its underlying interaction network. By applying this method to human-associated microbial communities from different body sites and sponge-associated microbial communities from different geographical locations, we found that in both cases the communities display a pronounced trade-off between the number of species and their effective connectance. These results suggest that natural microbiomes are shaped by stability constraints, which limit their complexity.

Global coordination level in single-cell transcriptomic data.

Genes are linked by underlying regulatory mechanisms and by jointly implementing biological functions, working in coordination to apply different tasks in the cells. Assessing the coordination level between genes from single-cell transcriptomic data, without a priori knowledge of the map of gene regulatory interactions, is a challenge. A 'top-down' approach has recently been developed to analyze single-cell transcriptomic data by evaluating the global coordination level between genes (called GCL). Here, we systematically analyze the performance of the GCL in typical scenarios of single-cell RNA sequencing (scRNA-seq) data. We show that an individual anomalous cell can have a disproportionate effect on the GCL calculated over a cohort of cells. In addition, we demonstrate how the GCL is affected by the presence of clusters, which are very common in scRNA-seq data. Finally, we analyze the effect of the sampling size of the Jackknife procedure on the GCL statistics. The manuscript is accompanied by a description of a custom-built Python package for calculating the GCL. These results provide practical guidelines for properly pre-processing and applying the GCL measure in transcriptional data.

A top-down measure of gene-to-gene coordination for analyzing cell-to-cell variability.

Recent technological advances, such as single-cell RNA sequencing (scRNA-seq), allow the measurement of gene expression profiles of individual cells. These expression profiles typically exhibit substantial variations even across seemingly homogeneous populations of cells. Two main different sources contribute to this measured variability: actual differences between the biological activity of the cells and technical measurement errors. Analysis of the biological variability may provide information about the underlying gene regulation of the cells, yet distinguishing it from the technical variability is a challenge. Here, we apply a recently developed computational method for measuring the global gene coordination level (GCL) to systematically study the cell-to-cell variability in numerical models of gene regulation. We simulate 'biological variability' by introducing heterogeneity in the underlying regulatory dynamic of different cells, while 'technical variability' is represented by stochastic measurement noise. We show that the GCL decreases for cohorts of cells with increased 'biological variability' only when it is originated from the interactions between the genes. Moreover, we find that the GCL can evaluate and compare-for cohorts with the same cell-to-cell variability-the ratio between the introduced biological and technical variability. Finally, we show that the GCL is robust against spurious correlations that originate from a small sample size or from the compositionality of the data. The presented methodology can be useful for future analysis of high-dimensional ecological and biochemical dynamics.

Dependency-based targeted attacks in interdependent networks.

Modern large engineered network systems normally work in cooperation and incorporate dependencies between their components for purposes of efficiency and regulation. Such dependencies may become a major risk since they can cause small-scale failures to propagate throughout the system. Thus, the dependent nodes could be a natural target for malicious attacks that aim to exploit these vulnerabilities. Here we consider a type of targeted attack that is based on the dependencies between the networks. We study strategies of attacks that range from dependency-first to dependency-last, where a fraction 1-p of the nodes with dependency links, or nodes without dependency links, respectively, are initially attacked. We systematically analyze, both analytically and numerically, the percolation transition of partially interdependent networks, where a fraction q of the nodes in each network are dependent on nodes in the other network. We find that for a broad range of dependency strength q, the "dependency-first" attack strategy is actually less effective, in terms of lower critical percolation threshold p_{c}, compared with random attacks of the same size. In contrast, the "dependency-last" attack strategy is more effective, i.e., higher p_{c}, compared with a random attack. This effect is explained by exploring the dynamics of the cascading failures initiated by dependency-based attacks. We show that while "dependency-first" strategy increases the short-term impact of the initial attack, in the long term the cascade slows down compared with the case of random attacks and vice versa for "dependency-last." Our results demonstrate that the effectiveness of attack strategies over a system of interdependent networks should be evaluated not only by the immediate impact but mainly by the accumulated damage during the process of cascading failures. This highlights the importance of understanding the dynamics of avalanches that may occur due to different scenarios of failures in order to design resilient critical infrastructures.

Age-related loss of gene-to-gene transcriptional coordination among single cells.

A long-standing model holds that stochastic aberrations of transcriptional regulation play a key role in the process of ageing. While transcriptional dysregulation is observed in many cell types in the form of increased cell-to-cell variability, its generality to all cell types remains doubted. Here, we propose a new approach for analysing transcriptional regulation in single-cell RNA sequencing data by focusing on the global coordination between the genes rather than the variability of individual genes or correlations between pairs of genes. Consistently, across very different organisms and cell types, we find a decrease in the gene-to-gene transcriptional coordination in ageing cells. In addition, we find that loss of gene-to-gene transcriptional coordination is associated with high mutational load of a specific, age-related signature and with radiation-induced DNA damage. These observations suggest a general, potentially universal, stochastic attribute of transcriptional dysregulation in ageing.

Articulation points in complex networks.

An articulation point in a network is a node whose removal disconnects the network. Those nodes play key roles in ensuring connectivity of many real-world networks, from infrastructure networks to protein interaction networks and terrorist communication networks. Despite their fundamental importance, a general framework of studying articulation points in complex networks is lacking. Here we develop analytical tools to study key issues pertinent to articulation points, such as the expected number of them and the network vulnerability against their removal, in an arbitrary complex network. We find that a greedy articulation point removal process provides us a different perspective on the organizational principles of complex networks. Moreover, this process results in a rich phase diagram with two fundamentally different types of percolation transitions. Our results shed light on the design of more resilient infrastructure networks and the effective destruction of terrorist communication networks.

Localized attacks on spatially embedded networks with dependencies.

Many real world complex systems such as critical infrastructure networks are embedded in space and their components may depend on one another to function. They are also susceptible to geographically localized damage caused by malicious attacks or natural disasters. Here, we study a general model of spatially embedded networks with dependencies under localized attacks. We develop a theoretical and numerical approach to describe and predict the effects of localized attacks on spatially embedded systems with dependencies. Surprisingly, we find that a localized attack can cause substantially more damage than an equivalent random attack. Furthermore, we find that for a broad range of parameters, systems which appear stable are in fact metastable. Though robust to random failures-even of finite fraction-if subjected to a localized attack larger than a critical size which is independent of the system size (i.e., a zero fraction), a cascading failure emerges which leads to complete system collapse. Our results demonstrate the potential high risk of localized attacks on spatially embedded network systems with dependencies and may be useful for designing more resilient systems.

Network Physiology: How Organ Systems Dynamically Interact.

We systematically study how diverse physiologic systems in the human organism dynamically interact and collectively behave to produce distinct physiologic states and functions. This is a fundamental question in the new interdisciplinary field of Network Physiology, and has not been previously explored. Introducing the novel concept of Time Delay Stability (TDS), we develop a computational approach to identify and quantify networks of physiologic interactions from long-term continuous, multi-channel physiological recordings. We also develop a physiologically-motivated visualization framework to map networks of dynamical organ interactions to graphical objects encoded with information about the coupling strength of network links quantified using the TDS measure. Applying a system-wide integrative approach, we identify distinct patterns in the network structure of organ interactions, as well as the frequency bands through which these interactions are mediated. We establish first maps representing physiologic organ network interactions and discover basic rules underlying the complex hierarchical reorganization in physiologic networks with transitions across physiologic states. Our findings demonstrate a direct association between network topology and physiologic function, and provide new insights into understanding how health and distinct physiologic states emerge from networked interactions among nonlinear multi-component complex systems. The presented here investigations are initial steps in building a first atlas of dynamic interactions among organ systems.

Simultaneous first- and second-order percolation transitions in interdependent networks.

In a system of interdependent networks, an initial failure of nodes invokes a cascade of iterative failures that may lead to a total collapse of the whole system in the form of an abrupt first-order transition. When the fraction of initial failed nodes 1-p reaches criticality p = p(c), the abrupt collapse occurs by spontaneous cascading failures. At this stage, the giant component decreases slowly in a plateau form and the number of iterations in the cascade τ diverges. The origin of this plateau and its increasing with the size of the system have been unclear. Here we find that, simultaneously with the abrupt first-order transition, a spontaneous second-order percolation occurs during the cascade of iterative failures. This sheds light on the origin of the plateau and how its length scales with the size of the system. Understanding the critical nature of the dynamical process of cascading failures may be useful for designing strategies for preventing and mitigating catastrophic collapses.

Network physiology reveals relations between network topology and physiological function.

The human organism is an integrated network where complex physiological systems, each with its own regulatory mechanisms, continuously interact, and where failure of one system can trigger a breakdown of the entire network. Identifying and quantifying dynamical networks of diverse systems with different types of interactions is a challenge. Here we develop a framework to probe interactions among diverse systems, and we identify a physiological network. We find that each physiological state is characterized by a specific network structure, demonstrating a robust interplay between network topology and function. Across physiological states, the network undergoes topological transitions associated with fast reorganization of physiological interactions on time scales of a few minutes, indicating high network flexibility in response to perturbations. The proposed system-wide integrative approach may facilitate the development of a new field, Network Physiology.