Plasmids in the human gut reveal neutral dispersal and recombination that is overpowered by inflammatory diseases.

Plasmids are pivotal in driving bacterial evolution through horizontal gene transfer. Here, we investigated 3467 human gut microbiome samples across continents and disease states, analyzing 11,086 plasmids. Our analyses reveal that plasmid dispersal is predominantly stochastic, indicating neutral processes as the primary driver of their wide distribution. We find that only 20-25% of plasmid DNA is being selected in various disease states, constraining its distribution across hosts. Selective pressures shape specific plasmid segments with distinct ecological functions, influenced by plasmid mobilization lifestyle, antibiotic usage, and inflammatory gut diseases. Notably, these elements are more commonly shared within groups of individuals with similar health conditions, such as Inflammatory Bowel Disease (IBD), regardless of geographic location across continents. These segments contain essential genes such as iron transport mechanisms- a distinctive gut signature of IBD that impacts the severity of inflammation. Our findings shed light on mechanisms driving plasmid dispersal and selection in the human gut, highlighting their role as carriers of vital gene pools impacting bacterial hosts and ecosystem dynamics.

Ecological network analysis reveals cancer-dependent chaperone-client interaction structure and robustness.

Cancer cells alter the expression levels of metabolic enzymes to fuel proliferation. The mitochondrion is a central hub of metabolic reprogramming, where chaperones service hundreds of clients, forming chaperone-client interaction networks. How network structure affects its robustness to chaperone targeting is key to developing cancer-specific drug therapy. However, few studies have assessed how structure and robustness vary across different cancer tissues. Here, using ecological network analysis, we reveal a non-random, hierarchical pattern whereby the cancer type modulates the chaperones' ability to realize their potential client interactions. Despite the low similarity between the chaperone-client interaction networks, we highly accurately predict links in one cancer type based on another. Moreover, we identify groups of chaperones that interact with similar clients. Simulations of network robustness show that this group structure affects cancer-specific response to chaperone removal. Our results open the door for new hypotheses regarding the ecology and evolution of chaperone-client interaction networks and can inform cancer-specific drug development strategies.

Invasive species modulate the structure and stability of a multilayer mutualistic network.

Species interactions are critical for maintaining community structure and dynamics, but the effects of invasive species on multitrophic networks remain poorly understood. We leveraged an ongoing invasion scenario in Patagonia, Argentina, to explore how non-native ungulates affect multitrophic networks. Ungulates disrupt a hummingbird-mistletoe-marsupial keystone interaction, which alters community composition. We sampled pollination and seed dispersal interactions in intact and invaded sites. We constructed pollination and seed dispersal networks for each site, which we connected via shared plants. We calculated pollination-seed dispersal connectivity, identified clusters of highly connected species, and quantified species' roles in connecting species clusters. To link structural variation to stability, we quantified network tolerance to single random species removal (disturbance propagation) and sequential species removal (robustness) using a stochastic coextinction model. Ungulates reduced the connectivity between pollination and seed dispersal and produced fewer clusters with a skewed size distribution. Moreover, species shifted their structural role, fragmenting the network by reducing the 'bridges' among species clusters. These structural changes altered the dynamics of cascading effects, increasing disturbance propagation and reducing network robustness. Our results highlight invasive species' role in altering community structure and subsequent stability in multitrophic communities.

Multilayer networks of plasmid genetic similarity reveal potential pathways of gene transmission.

Antimicrobial resistance (AMR) is a significant threat to public health. Plasmids are principal vectors of AMR genes, significantly contributing to their spread and mobility across hosts. Nevertheless, little is known about the dynamics of plasmid genetic exchange across animal hosts. Here, we use theory and methodology from network and disease ecology to investigate the potential of gene transmission between plasmids using a data set of 21 plasmidomes from a single dairy cow population. We constructed a multilayer network based on pairwise plasmid genetic similarity. Genetic similarity is a signature of past genetic exchange that can aid in identifying potential routes and mechanisms of gene transmission within and between cows. Links between cows dominated the transmission network, and plasmids containing mobility genes were more connected. Modularity analysis revealed a network cluster where all plasmids contained a mobM gene, and one where all plasmids contained a beta-lactamase gene. Cows that contain both clusters also share transmission pathways with many other cows, making them candidates for super-spreading. In support, we found signatures of gene super-spreading in which a few plasmids and cows are responsible for most gene exchange. An agent-based transmission model showed that a new gene invading the cow population will likely reach all cows. Finally, we showed that edge weights contain a non-random signature for the mechanisms of gene transmission, allowing us to differentiate between dispersal and genetic exchange. These results provide insights into how genes, including those providing AMR, spread across animal hosts.

Conceptualizing microbe-plasmid communities as complex adaptive systems.

Plasmids shape microbial communities' diversity, structure, and function. Nevertheless, we lack a mechanistic understanding of how community structure and dynamics emerge from local microbe-plasmid interactions and coevolution. Addressing this gap is challenging because multiple processes operate simultaneously at multiple levels of organization. For example, immunity operates between a plasmid and a cell, but incompatibility mechanisms regulate coexistence between plasmids. Conceptualizing microbe-plasmid communities as complex adaptive systems is a promising approach to overcoming these challenges. I illustrate how agent-based evolutionary modeling, extended by network analysis, can be used to quantify the relative importance of local processes governing community dynamics. These theoretical developments can advance our understanding of plasmid ecology and evolution, especially when combined with empirical data.

Age-specific patterns of DBLα var diversity can explain why residents of high malaria transmission areas remain susceptible to Plasmodium falciparum blood stage infection throughout life.

Immunity to Plasmodium falciparum is non-sterilising, thus individuals residing in malaria-endemic areas are at risk of infection throughout their lifetime. Here we seek to find a genomic epidemiological explanation for why residents of all ages harbour blood stage infections despite lifelong exposure to P. falciparum in areas of high transmission. We do this by exploring, for the first known time, the age-specific patterns of diversity of variant antigen encoding (var) genes in the reservoir of infection. Microscopic and submicroscopic P. falciparum infections were analysed at the end of the wet and dry seasons in 2012-2013 for a cohort of 1541 residents aged from 1 to 91 years in an area characterised by high seasonal malaria transmission in Ghana. By sequencing the near ubiquitous Duffy-binding-like alpha domain (DBLα) that encodes immunogenic domains, we defined var gene diversity in an estimated 1096 genomes detected in sequential wet and dry season sampling of this cohort. Unprecedented var (DBLα) diversity was observed in all ages with 42,399 unique var types detected. There was a high degree of maintenance of types between seasons (>40% seen more than once), with many of the same types, especially upsA, appearing multiple times in isolates from different individuals. Children and adolescents were found to be significant reservoirs of var DBLα diversity compared with adults. Var repertoires within individuals were highly variable, with children having more related var repertoires compared to adolescents and adults. Individuals of all ages harboured multiple genomes with var repertoires unrelated to those infecting other hosts. High turnover of parasites with diverse isolate var repertoires was also observed in all ages. These age-specific patterns are best explained by variant-specific immune selection. The observed level of var diversity for the population was then used to simulate the development of variant-specific immunity to the diverse var types under conservative assumptions. Simulations showed that the extent of observed var diversity with limited repertoire relatedness was sufficient to explain why adolescents and adults in this community remain susceptible to blood stage infection, even with multiple genomes.

Habitat aridity as a determinant of the trade-off between water conservation and evaporative heat loss in bats.

The maintenance of water balance in arid environments might represent a formidable challenge for Chiroptera, since they have high surface-to-volume ratios. In deserts, bats conserve water, for example, using daily torpor, but they also might experience episodic heat bouts, when they may need to increase total evaporative water loss (TEWL) to thermoregulate. We hypothesized that in bats, habitat aridity and its variability determine a trade-off between water conservation and thermoregulation via evaporative means. To test this hypothesis, we collated data from the literature of 22 species of bats on TEWL, body temperature and resting metabolic rate, in torpor and euthermy. We also collected data on ambient temperature (Ta) and precipitation of the locations where bats were captured, calculated an aridity index, and built an index of variability of the environment. After correcting for phylogeny, we found that, as aridity and variability of the environment increased, bats had lower values of TEWL, but the rate at which TEWL increases with Ta was higher, supporting our hypothesis. These results suggest that at high Ta there is a trade-off between water conservation and evaporative heat loss in bats. The evolution of physiological mechanisms that allow water conservation and tolerance to conditions of high Ta without access to free water might thus be crucial to explain the distribution of desert bats.

The interplay between vaccination and social distancing strategies affects COVID19 population-level outcomes.

Social distancing is an effective population-level mitigation strategy to prevent COVID19 propagation but it does not reduce the number of susceptible individuals and bears severe social consequences-a dire situation that can be overcome with the recently developed vaccines. Although a combination of these interventions should provide greater benefits than their isolated deployment, a mechanistic understanding of the interplay between them is missing. To tackle this challenge we developed an age-structured deterministic model in which vaccines are deployed during the pandemic to individuals who do not show symptoms. The model allows for flexible and dynamic prioritization strategies with shifts between target groups. We find a strong interaction between social distancing and vaccination in their effect on the proportion of hospitalizations. In particular, prioritizing vaccines to elderly (60+) before adults (20-59) is more effective when social distancing is applied to adults or uniformly. In addition, the temporal reproductive number Rt is only affected by vaccines when deployed at sufficiently high rates and in tandem with social distancing. Finally, the same reduction in hospitalization can be achieved via different combination of strategies, giving decision makers flexibility in choosing public health policies. Our study provides insights into the factors that affect vaccination success and provides methodology to test different intervention strategies in a way that will align with ethical guidelines.

Frequency-Dependent Competition Between Strains Imparts Persistence to Perturbations in a Model of Plasmodium falciparum Malaria Transmission.

In high-transmission endemic regions, local populations of Plasmodium falciparum exhibit vast diversity of the var genes encoding its major surface antigen, with each parasite comprising multiple copies from this diverse gene pool. This strategy to evade the immune system through large combinatorial antigenic diversity is common to other hyperdiverse pathogens. It underlies a series of fundamental epidemiological characteristics, including large reservoirs of transmission from high prevalence of asymptomatics and long-lasting infections. Previous theory has shown that negative frequency-dependent selection (NFDS) mediated by the acquisition of specific immunity by hosts structures the diversity of var gene repertoires, or strains, in a pattern of limiting similarity that is both non-random and non-neutral. A combination of stochastic agent-based models and network analyses has enabled the development and testing of theory in these complex adaptive systems, where assembly of local parasite diversity occurs under frequency-dependent selection and large pools of variation. We show here the application of these approaches to theory comparing the response of the malaria transmission system to intervention when strain diversity is assembled under (competition-based) selection vs. a form of neutrality, where immunity depends only on the number but not the genetic identity of previous infections. The transmission system is considerably more persistent under NFDS, exhibiting a lower extinction probability despite comparable prevalence during intervention. We explain this pattern on the basis of the structure of strain diversity, in particular the more pronounced fraction of highly dissimilar parasites. For simulations that survive intervention, prevalence under specific immunity is lower than under neutrality, because the recovery of diversity is considerably slower than that of prevalence and decreased var gene diversity reduces parasite transmission. A Principal Component Analysis of network features describing parasite similarity reveals that despite lower overall diversity, NFDS is quickly restored after intervention constraining strain structure and maintaining patterns of limiting similarity important to parasite persistence. Given the described enhanced persistence under perturbation, intervention efforts will likely require longer times than the usual practice to eliminate P. falciparum populations. We discuss implications of our findings and potential analogies for ecological communities with non-neutral assembly processes involving frequency-dependence.

The network structure and eco-evolutionary dynamics of CRISPR-induced immune diversification.

As a heritable sequence-specific adaptive immune system, CRISPR-Cas is a powerful force shaping strain diversity in host-virus systems. While the diversity of CRISPR alleles has been explored, the associated structure and dynamics of host-virus interactions have not. We explore the role of CRISPR in mediating the interplay between host-virus interaction structure and eco-evolutionary dynamics in a computational model and compare the results with three empirical datasets from natural systems. We show that the structure of the networks describing who infects whom and the degree to which strains are immune, are respectively modular (containing groups of hosts and viruses that interact strongly) and weighted-nested (specialist hosts are more susceptible to subsets of viruses that in turn also infect the more generalist hosts with many spacers matching many viruses). The dynamic interplay between these networks influences transitions between dynamical regimes of virus diversification and host control. The three empirical systems exhibit weighted-nested immunity networks, a pattern our theory shows is indicative of hosts able to suppress virus diversification. Previously missing from studies of microbial host-pathogen systems, the immunity network plays a key role in the coevolutionary dynamics.

Pan-Cancer Analysis of Mitochondria Chaperone-Client Co-Expression Reveals Chaperone Functional Partitioning.

Metabolic reprogramming is a hallmark of cancer. Such reprogramming entails the up-regulation of the expression of specific mitochondrial proteins, thus increasing the burden on the mitochondrial protein quality control. However, very little is known about the specificity of interactions between mitochondrial chaperones and their clients, or to what extent the mitochondrial chaperone-client co-expression is coordinated. We hypothesized that a physical interaction between a chaperone and its client in mitochondria ought to be manifested in the co-expression pattern of both transcripts. Using The Cancer Genome Atlas (TCGA) gene expression data from 13 tumor entities, we constructed the mitochondrial chaperone-client co-expression network. We determined that the network is comprised of three distinct modules, each populated with unique chaperone-clients co-expression pairs belonging to distinct functional groups. Surprisingly, chaperonins HSPD1 and HSPE1, which are known to comprise a functional complex, each occupied a different module: HSPD1 co-expressed with tricarboxylic acid cycle cycle enzymes, while HSPE1 co-expressed with proteins involved in oxidative phosphorylation. Importantly, we found that the genes in each module were enriched for discrete transcription factor binding sites, suggesting the mechanism for the coordinated co-expression. We propose that our mitochondrial chaperone-client interactome can facilitate the identification of chaperones supporting specific mitochondrial pathways and bring forth a fundamental principle in metabolic adaptation.

Competition for hosts modulates vast antigenic diversity to generate persistent strain structure in Plasmodium falciparum.

In their competition for hosts, parasites with antigens that are novel to the host immune system will be at a competitive advantage. The resulting frequency-dependent selection can structure parasite populations into strains of limited genetic overlap. For the causative agent of malaria, Plasmodium falciparum, the high recombination rates and associated vast diversity of its highly antigenic and multicopy var genes preclude such clear clustering in endemic regions. This undermines the definition of strains as specific, temporally persisting gene variant combinations. We use temporal multilayer networks to analyze the genetic similarity of parasites in both simulated data and in an extensively and longitudinally sampled population in Ghana. When viewed over time, populations are structured into modules (i.e., groups) of parasite genomes whose var gene combinations are more similar within than between the modules and whose persistence is much longer than that of the individual genomes that compose them. Comparison to neutral models that retain parasite population dynamics but lack competition reveals that the selection imposed by host immunity promotes the persistence of these modules. The modular structure is, in turn, associated with a slower acquisition of immunity by individual hosts. Modules thus represent dynamically generated niches in host immune space, which can be interpreted as strains. Negative frequency-dependent selection therefore shapes the organization of the var diversity into parasite genomes, leaving a persistence signature over ecological time scales. Multilayer networks extend the scope of phylodynamics analyses by allowing quantification of temporal genetic structure in organisms that generate variation via recombination or other non-bifurcating processes. A strain structure similar to the one described here should apply to other pathogens with large antigenic spaces that evolve via recombination. For malaria, the temporal modular structure should enable the formulation of tractable epidemiological models that account for parasite antigenic diversity and its influence on intervention outcomes.

Networks of genetic similarity reveal non-neutral processes shape strain structure in Plasmodium falciparum.

Pathogens compete for hosts through patterns of cross-protection conferred by immune responses to antigens. In Plasmodium falciparum malaria, the var multigene family encoding for the major blood-stage antigen PfEMP1 has evolved enormous genetic diversity through ectopic recombination and mutation. With 50-60 var genes per genome, it is unclear whether immune selection can act as a dominant force in structuring var repertoires of local populations. The combinatorial complexity of the var system remains beyond the reach of existing strain theory and previous evidence for non-random structure cannot demonstrate immune selection without comparison with neutral models. We develop two neutral models that encompass malaria epidemiology but exclude competitive interactions between parasites. These models, combined with networks of genetic similarity, reveal non-neutral strain structure in both simulated systems and an extensively sampled population in Ghana. The unique population structure we identify underlies the large transmission reservoir characteristic of highly endemic regions in Africa.

The multilayer nature of ecological networks.

Although networks provide a powerful approach to study a large variety of ecological systems, their formulation does not typically account for multiple interaction types, interactions that vary in space and time, and interconnected systems such as networks of networks. The emergent field of 'multilayer networks' provides a natural framework for extending analyses of ecological systems to include such multiple layers of complexity, as it specifically allows one to differentiate and model 'intralayer' and 'interlayer' connectivity. The framework provides a set of concepts and tools that can be adapted and applied to ecology, facilitating research on high-dimensional, heterogeneous systems in nature. Here, we formally define ecological multilayer networks based on a review of previous, related approaches; illustrate their application and potential with analyses of existing data; and discuss limitations, challenges, and future applications. The integration of multilayer network theory into ecology offers largely untapped potential to investigate ecological complexity and provide new theoretical and empirical insights into the architecture and dynamics of ecological systems.

Asymmetric disease dynamics in multihost interconnected networks.

Epidemic spread in single-host systems strongly depends on the population's transmission network. However, little is known regarding the spread of epidemics across networks representing populations of multiple hosts. We explored cross-species transmission in a multilayer network where layers represent populations of two distinct hosts, and disease can spread across intralayer (within-host) and interlayer (between-host) edges. We developed an analytic framework for the SIR epidemic model to examine the effect of (i) source of infection and (ii) between-host asymmetry in infection probabilities, on disease risk. We measured risk as outbreak probability and outbreak size in a focal host, represented by one network layer. Numeric simulations were used to validate the analytic formulations. We found that outbreak probability is determined by a complex interaction between source of infection and between-host infection probabilities, whereas outbreak size is mainly affected by the non-focal host to focal host infection probability. Hence, inter-specific asymmetry in infection probabilities shapes disease dynamics in multihost networks. These results highlight the importance of considering multiple measures of disease risk and advance our understanding of disease spread in multihost systems. The study provides a flexible way to model disease dynamics in multiple hosts while considering contact heterogeneity within and between species. We strongly encourage empirical studies that include information on both cross-species infection rates and network structure of multiple hosts. Such studies are necessary to corroborate our theoretical results and to improve our understanding of multihost epidemiology.

The effect of water contamination and host-related factors on ectoparasite load in an insectivorous bat.

We examined the effects of sex, age, and reproductive state of the insectivorous bat Pipistrellus kuhlii on the abundance and prevalence of arthropod ectoparasites (Macronyssidae and Cimicidae) in habitats with either sewage-polluted or natural bodies of water, in the Negev Desert, Israel. We chose water pollution as an environmental factor because of the importance of water availability in desert environments, particularly for P. kuhlii, which needs to drink on a daily basis. We predicted that parasite infestation rates would be affected by both environment and demographic cohort of the host. We found that female bats in the polluted site harbored significantly more mites than female bats in the natural site and that juveniles in the polluted site harbored significantly more cimicid individuals than juveniles in the natural site. We further found that age and sex (host-related factors) affected ectoparasite prevalence and intensity (i.e., the abundance of parasites) in the polluted site. Our results may suggest that the interaction between host-related and environment-related factors affected parasite infestations, with females and young bats being more susceptible to ectoparasites when foraging over polluted water. This effect may be particularly important for bats that must drink or forage above water for other wildlife that depend on drinking water for survival.

Correlation-Based Network Analysis of Metabolite and Enzyme Profiles Reveals a Role of Citrate Biosynthesis in Modulating N and C Metabolism in Zea mays.

To investigate the natural variability of leaf metabolism and enzymatic activity in a maize inbred population, statistical and network analyses were employed on metabolite and enzyme profiles. The test of coefficient of variation showed that sugars and amino acids displayed opposite trends in their variance within the population, consistently with their related enzymes. The overall higher CV values for metabolites as compared to the tested enzymes are indicative for their greater phenotypic plasticity. H(2) tests revealed galactinol (1) and asparagine (0.91) as the highest scorers among metabolites and nitrate reductase (0.73), NAD-glutamate dehydrogenase (0.52), and phosphoglucomutase (0.51) among enzymes. The overall low H(2) scores for metabolites and enzymes are suggestive for a great environmental impact or gene-environment interaction. Correlation-based network generation followed by community detection analysis, partitioned the network into three main communities and one dyad, (i) reflecting the different levels of phenotypic plasticity of the two molecular classes as observed for the CV values and (ii) highlighting the concerted changes between classes of chemically related metabolites. Community 1 is composed mainly of enzymes and specialized metabolites, community 2' is enriched in N-containing compounds and phosphorylated-intermediates. The third community contains mainly organic acids and sugars. Cross-community linkages are supported by aspartate, by the photorespiration amino acids glycine and serine, by the metabolically related GABA and putrescine, and by citrate. The latter displayed the strongest node-betweenness value (185.25) of all nodes highlighting its fundamental structural role in the connectivity of the network by linking between different communities and to the also strongly connected enzyme aldolase.

Potential parasite transmission in multi-host networks based on parasite sharing.

Epidemiological networks are commonly used to explore dynamics of parasite transmission among individuals in a population of a given host species. However, many parasites infect multiple host species, and thus multi-host networks may offer a better framework for investigating parasite dynamics. We investigated the factors that influence parasite sharing--and thus potential transmission pathways--among rodent hosts in Southeast Asia. We focused on differences between networks of a single host species and networks that involve multiple host species. In host-parasite networks, modularity (the extent to which the network is divided into subgroups of rodents that interact with similar parasites) was higher in the multi-species than in the single-species networks. This suggests that phylogeny affects patterns of parasite sharing, which was confirmed in analyses showing that it predicted affiliation of individuals to modules. We then constructed "potential transmission networks" based on the host-parasite networks, in which edges depict the similarity between a pair of individuals in the parasites they share. The centrality of individuals in these networks differed between multi- and single-species networks, with species identity and individual characteristics influencing their position in the networks. Simulations further revealed that parasite dynamics differed between multi- and single-species networks. We conclude that multi-host networks based on parasite sharing can provide new insights into the potential for transmission among hosts in an ecological community. In addition, the factors that determine the nature of parasite sharing (i.e. structure of the host-parasite network) may impact transmission patterns.

Habitat fragmentation alters the properties of a host-parasite network: rodents and their helminths in South-East Asia.

1. While the effects of deforestation and habitat fragmentation on parasite prevalence or richness are well investigated, host-parasite networks are still understudied despite their importance in understanding the mechanisms of these major disturbances. Because fragmentation may negatively impact species occupancy, abundance and co-occurrence, we predict a link between spatiotemporal changes in habitat and the architecture of host-parasite networks. 2. For this, we used an extensive data set on 16 rodent species and 29 helminth species from seven localities of South-East Asia. We analysed the effects of rapid deforestation on connectance and modularity of helminth-parasite networks. We estimated both the degree of fragmentation and the rate of deforestation through the development of land uses and their changes through the last 20 to 30 years in order to take into account the dynamics of habitat fragmentation in our statistical analyses. 3. We found that rapid fragmentation does not affect helminth species richness per se but impacts host-parasite interactions as the rodent-helminth network becomes less connected and more modular. 4. Our results suggest that parasite sharing among host species may become more difficult to maintain with the increase of habitat disturbance.

Host-parasite network structure is associated with community-level immunogenetic diversity.

Genes of the major histocompatibility complex (MHC) encode proteins that recognize foreign antigens and are thus crucial for immune response. In a population of a single host species, parasite-mediated selection drives MHC allelic diversity. However, in a community-wide context, species interactions may modulate selection regimes because the prevalence of a given parasite in a given host may depend on its prevalence in other hosts. By combining network analysis with immunogenetics, we show that host species infected by similar parasites harbour similar alleles with similar frequencies. We further show, using a Bayesian approach, that the probability of mutual occurrence of a functional allele and a parasite in a given host individual is nonrandom and depends on other host-parasite interactions, driving co-evolution within subgroups of parasite species and functional alleles. Therefore, indirect effects among hosts and parasites can shape host MHC diversity, scaling it from the population to the community level.