Inferring species interaction networks from species abundance data: A comparative evaluation of various statistical and machine learning methods

On the structure of species-function participation in multilayer ecological networks

Understanding how biotic interactions shape ecosystems and impact their functioning, resilience and biodiversity has been a sustained research priority in ecology. Yet, traditional assessments of ecological complexity typically focus on species-species interactions that mediate a particular function (e.g., pollination), overlooking both the synergistic effect that multiple functions might develop as well as the resulting species-function participation patterns that emerge in ecosystems that harbor multiple ecological functions. Here we propose a mathematical framework that integrates various types of biotic interactions observed between different species. Its application to recently collected data of an islet ecosystem-reporting 1537 interactions between 691 plants, animals and fungi across six different functions (pollination, herbivory, seed dispersal, decomposition, nutrient uptake, and fungal pathogenicity)-unveils a non-random, nested structure in the way plant species participate across different functions. The framework further allows us to identify a ranking of species and functions, where woody shrubs and fungal decomposition emerge as keystone actors whose removal have a larger-than-random effect on secondary extinctions. The dual insight-from species and functional perspectives-offered by the framework opens the door to a richer quantification of ecosystem complexity and to better calibrate the influence of multifunctionality on ecosystem functioning and biodiversity.

Signals of Potential Species Associations Offer Clues about Community Organisation of Stream Fish across Seasons

Simple Summary

Species interactions are one of the main factors affecting community assembly, yet the role of such interactions remains mostly unknown. Here, we investigated roles of potential species associations in fish community assembly in the Qiupu River, China. Our results suggested that potential species associations might have been underestimated in stream fish community assembly. The contribution of potential species associations to fish community assembly can be reflected by interaction network structures. Omnivorous species play an important role in maintaining network structure as they may have more associations with other species. This study highlights the importance of capturing species associations in river ecosystems across different geographical and environmental settings.

Abstract

Environmental filtering, spatial factors and species interactions are fundamental ecological mechanisms for community organisation, yet the role of such interactions across different environmental and spatial settings remains mostly unknown. In this study, we investigated fish community organisation scenarios and seasonal species-to-species associations potentially reflecting biotic associations along the Qiupu River (China). Based on a latent variable approach and a tree-based method, we compared the relative contribution of the abiotic environment, spatial covariates and potential species associations for variation in the community structure, and assessed whether different assembly scenarios were modulated by concomitant changes in the interaction network structure of fish communities across seasons. We found that potential species associations might have been underestimated in community-based assessments of stream fish. Omnivore species, since they have more associations with other species, were found to be key components sustaining fish interaction networks across different stream orders. Hence, we suggest that species interactions, such as predation and competition, likely played a key role in community structure. For instance, indices accounting for network structure, such as connectance and nestedness, were strongly correlated with the unexplained residuals from our latent variable approach, thereby re-emphasising that biotic signals, potentially reflecting species interactions, may be of primary importance in determining stream fish communities across seasons. Overall, our findings indicate that interaction network structures are a powerful tool to reflect the contribution of potential species associations to community assembly. From an applied perspective, this study should encourage freshwater ecologists to empirically capture and manage biotic constraints in stream ecosystems across different geographical and environmental settings, especially in the context of the ever-increasing impacts of human-induced local extinction debts and species invasions.

Machine Learning Advances in Microbiology: A Review of Methods and Applications

Microorganisms play an important role in natural material and elemental cycles. Many common and general biology research techniques rely on microorganisms. Machine learning has been gradually integrated with multiple fields of study. Machine learning, including deep learning, aims to use mathematical insights to optimize variational functions to aid microbiology using various types of available data to help humans organize and apply collective knowledge of various research objects in a systematic and scaled manner. Classification and prediction have become the main achievements in the development of microbial community research in the direction of computational biology. This review summarizes the application and development of machine learning and deep learning in the field of microbiology and shows and compares the advantages and disadvantages of different algorithm tools in four fields: microbiome and taxonomy, microbial ecology, pathogen and epidemiology, and drug discovery.

Perspectives in machine learning for wildlife conservation

Inexpensive and accessible sensors are accelerating data acquisition in animal ecology. These technologies hold great potential for large-scale ecological understanding, but are limited by current processing approaches which inefficiently distill data into relevant information. We argue that animal ecologists can capitalize on large datasets generated by modern sensors by combining machine learning approaches with domain knowledge. Incorporating machine learning into ecological workflows could improve inputs for ecological models and lead to integrated hybrid modeling tools. This approach will require close interdisciplinary collaboration to ensure the quality of novel approaches and train a new generation of data scientists in ecology and conservation.

Inferring species interactions from ecological survey data: A mechanistic approach to predict quantitative food webs of seed feeding by carabid beetles

Ecological networks are valuable for ecosystem analysis but their use is often limited by a lack of data because many types of ecological interaction, for example, predation, are short-lived and difficult to observe or detect. While there are different methods for inferring the presence of interactions, they have rarely been used to predict the interaction strengths that are required to construct weighted, or quantitative, ecological networks.Here, we develop a trait-based approach suitable for inferring weighted networks, that is, with varying interaction strengths. We developed the method for seed-feeding carabid ground beetles (Coleoptera: Carabidae) although the principles can be applied to other species and types of interaction.Using existing literature data from experimental seed-feeding trials, we predicted a per-individual interaction cost index based on carabid and seed size. This was scaled up to the population level to create inferred weighted networks using the abundance of carabids and seeds from empirical samples and energetic intake rates of carabids from the literature. From these weighted networks, we also derived a novel measure of expected predation pressure per seed type per network.This method was applied to existing ecological survey data from 255 arable fields with carabid data from pitfall traps and plant seeds from seed rain traps. Analysis of these inferred networks led to testable hypotheses about how network structure and predation pressure varied among fields.Inferred networks are valuable because (a) they provide null models for the structuring of food webs to test against empirical species interaction data, for example, DNA analysis of carabid gut regurgitates and (b) they allow weighted networks to be constructed whenever we can estimate interactions between species and have ecological census data available. This permits ecological network analysis even at times and in places when interactions were not directly assessed.

Enhanced inference of ecological networks by parameterizing ensembles of population dynamics models constrained with prior knowledge

Background

Accurate network models of species interaction could be used to predict population dynamics and be applied to manage real world ecosystems. Most relevant models are nonlinear, however, and data available from real world ecosystems are too noisy and sparsely sampled for common inference approaches. Here we improved the inference of generalized Lotka–Volterra (gLV) ecological networks by using a new optimization algorithm to constrain parameter signs with prior knowledge and a perturbation-based ensemble method.

Results

We applied the new inference to long-term species abundance data from the freshwater fish community in the Illinois River, United States. We constructed an ensemble of 668 gLV models that explained 79% of the data on average. The models indicated (at a 70% level of confidence) a strong positive interaction from emerald shiner (Notropis atherinoides) to channel catfish (Ictalurus punctatus), which we could validate using data from a nearby observation site, and predicted that the relative abundances of most fish species will continue to fluctuate temporally and concordantly in the near future. The network shows that the invasive silver carp (Hypophthalmichthys molitrix) has much stronger impacts on native predators than on prey, supporting the notion that the invader perturbs the native food chain by replacing the diets of predators.

Conclusions

Ensemble approaches constrained by prior knowledge can improve inference and produce networks from noisy and sparsely sampled time series data to fill knowledge gaps on real world ecosystems. Such network models could aid efforts to conserve ecosystems such as the Illinois River, which is threatened by the invasion of the silver carp.

Positive associations among rare species and their persistence in ecological assemblages

According to the competitive exclusion principle, species with low competitive abilities should be excluded by more efficient competitors; yet, they generally remain as rare species. Here, we describe the positive and negative spatial association networks of 326 disparate assemblages, showing a general organization pattern that simultaneously supports the primacy of competition and the persistence of rare species. Abundant species monopolize negative associations in about 90% of the assemblages. On the other hand, rare species are mostly involved in positive associations, forming small network modules. Simulations suggest that positive interactions among rare species and microhabitat preferences are the most probable mechanisms underpinning this pattern and rare species persistence. The consistent results across taxa and geography suggest a general explanation for the maintenance of biodiversity in competitive environments.

Predicting ecosystem components in the Gulf of Mexico and their responses to climate variability with a dynamic Bayesian network model

The Gulf of Mexico is an ecologically and economically important marine ecosystem that is affected by a variety of natural and anthropogenic pressures. These complex and interacting pressures, together with the dynamic environment of the Gulf, present challenges for the effective management of its resources. The recent adoption of Bayesian networks to ecology allows for the discovery and quantification of complex interactions from data after making only a few assumptions about observations of the system. In this study, we apply Bayesian network models, with different levels of structural complexity and a varying number of hidden variables to account for uncertainty when modeling ecosystem dynamics. From these models, we predict focal ecosystem components within the Gulf of Mexico. The predictive ability of the models varied with their structure. The model that performed best was parameterized through data-driven learning techniques and accounted for multiple ecosystem components’ associations and their interactions with human and natural pressures over time. Then, we altered sea surface temperature in the best performing model to explore the response of different ecosystem components to increased temperature. The magnitude and even direction of predicted responses varied by ecosystem components due to heterogeneity in driving factors and their spatial overlap. Our findings suggest that due to varying components’ sensitivity to drivers, changes in temperature will potentially lead to trade-offs in terms of population productivity. We were able to discover meaningful interactions between ecosystem components and their environment and show how sensitive these relationships are to climate perturbations, which increases our understanding of the potential future response of the system to increasing temperature. Our findings demonstrate that accounting for additional sources of variation, by incorporating multiple interactions and pressures in the model layout, has the potential for gaining deeper insights into the structure and dynamics of ecosystems.

Ecological Network Inference From Long-Term Presence-Absence Data

Ecological communities are characterized by complex networks of trophic and nontrophic interactions, which shape the dy-namics of the community. Machine learning and correlational methods are increasingly popular for inferring networks from co-occurrence and time series data, particularly in microbial systems. In this study, we test the suitability of these methods for inferring ecological interactions by constructing networks using Dynamic Bayesian Networks, Lasso regression, and Pear-son’s correlation coefficient, then comparing the model networks to empirical trophic and nontrophic webs in two ecological systems. We find that although each model significantly replicates the structure of at least one empirical network, no model significantly predicts network structure in both systems, and no model is clearly superior to the others. We also find that networks inferred for the Tatoosh intertidal match the nontrophic network much more closely than the trophic one, possibly due to the challenges of identifying trophic interactions from presence-absence data. Our findings suggest that although these methods hold some promise for ecological network inference, presence-absence data does not provide enough signal for models to consistently identify interactions, and networks inferred from these data should be interpreted with caution.

Next-Generation Global Biomonitoring: Large-scale, Automated Reconstruction of Ecological Networks

We foresee a new global-scale, ecological approach to biomonitoring emerging within the next decade that can detect ecosystem change accurately, cheaply, and generically. Next-generation sequencing of DNA sampled from the Earth's environments would provide data for the relative abundance of operational taxonomic units or ecological functions. Machine-learning methods would then be used to reconstruct the ecological networks of interactions implicit in the raw NGS data. Ultimately, we envision the development of autonomous samplers that would sample nucleic acids and upload NGS sequence data to the cloud for network reconstruction. Large numbers of these samplers, in a global array, would allow sensitive automated biomonitoring of the Earth's major ecosystems at high spatial and temporal resolution, revolutionising our understanding of ecosystem change.

Cross-biome comparison of microbial association networks

Clinical and environmental meta-omics studies are accumulating an ever-growing amount of microbial abundance data over a wide range of ecosystems. With a sufficiently large sample number, these microbial communities can be explored by constructing and analyzing co-occurrence networks, which detect taxon associations from abundance data and can give insights into community structure. Here, we investigate how co-occurrence networks differ across biomes and which other factors influence their properties. For this, we inferred microbial association networks from 20 different 16S rDNA sequencing data sets and observed that soil microbial networks harbor proportionally fewer positive associations and are less densely interconnected than host-associated networks. After excluding sample number, sequencing depth and beta-diversity as possible drivers, we found a negative correlation between community evenness and positive edge percentage. This correlation likely results from a skewed distribution of negative interactions, which take place preferentially between less prevalent taxa. Overall, our results suggest an under-appreciated role of evenness in shaping microbial association networks.

Incorporating uncertainty in predictive species distribution modelling

Motivated by the need to solve ecological problems (climate change, habitat fragmentation and biological invasions), there has been increasing interest in species distribution models (SDMs). Predictions from these models inform conservation policy, invasive species management and disease-control measures. However, predictions are subject to uncertainty, the degree and source of which is often unrecognized. Here, we review the SDM literature in the context of uncertainty, focusing on three main classes of SDM: niche-based models, demographic models and process-based models. We identify sources of uncertainty for each class and discuss how uncertainty can be minimized or included in the modelling process to give realistic measures of confidence around predictions. Because this has typically not been performed, we conclude that uncertainty in SDMs has often been underestimated and a false precision assigned to predictions of geographical distribution. We identify areas where development of new statistical tools will improve predictions from distribution models, notably the development of hierarchical models that link different types of distribution model and their attendant uncertainties across spatial scales. Finally, we discuss the need to develop more defensible methods for assessing predictive performance, quantifying model goodness-of-fit and for assessing the significance of model covariates.

Construction of gene regulatory networks using biclustering and Bayesian networks

Background

Understanding gene interactions in complex living systems can be seen as the ultimate goal of the systems biology revolution. Hence, to elucidate disease ontology fully and to reduce the cost of drug development, gene regulatory networks (GRNs) have to be constructed. During the last decade, many GRN inference algorithms based on genome-wide data have been developed to unravel the complexity of gene regulation. Time series transcriptomic data measured by genome-wide DNA microarrays are traditionally used for GRN modelling. One of the major problems with microarrays is that a dataset consists of relatively few time points with respect to the large number of genes. Dimensionality is one of the interesting problems in GRN modelling.

Results

In this paper, we develop a biclustering function enrichment analysis toolbox (BicAT-plus) to study the effect of biclustering in reducing data dimensions. The network generated from our system was validated via available interaction databases and was compared with previous methods. The results revealed the performance of our proposed method.

Conclusions

Because of the sparse nature of GRNs, the results of biclustering techniques differ significantly from those of previous methods.