Background: Understanding the interaction among different species within a community and their responses to environmental changes is a central goal in ecology. However, defining the network structure in a microbial community is very challenging due to their extremely high diversity and as-yet uncultivated status. Although recent advance of metagenomic technologies, such as high throughout sequencing and functional gene arrays, provide revolutionary tools for analyzing microbial community structure, it is still difficult to examine network interactions in a microbial community based on high-throughput metagenomics data.Results: Here, we describe a novel mathematical and bioinformatics framework to construct ecological association networks named molecular ecological networks (MENs) through Random Matrix Theory (RMT)-based methods. Compared to other network construction methods, this approach is remarkable in that the network is automatically defined and robust to noise, thus providing excellent solutions to several common issues associated with high-throughput metagenomics data. We applied it to determine the network structure of microbial communities subjected to long-term experimental warming based on pyrosequencing data of 16 S rRNA genes. We showed that the constructed MENs under both warming and unwarming conditions exhibited topological features of scale free, small world and modularity, which were consistent with previously described molecular ecological networks. Eigengene analysis indicated that the eigengenes represented the module profiles relatively well. In consistency with many other studies, several major environmental traits including temperature and soil pH were found to be important in determining network interactions in the microbial communities examined. To facilitate its application by the scientific community, all these methods and statistical tools have been integrated into a comprehensive Molecular Ecological Network Analysis Pipeline (MENAP), which is open-accessible now (http://ieg2.ou.edu/MENA).Conclusions: The RMT-based molecular ecological network analysis provides powerful tools to elucidate network interactions in microbial communities and their responses to environmental changes, which are fundamentally important for research in microbial ecology and environmental microbiology.

Understanding the community assembly mechanisms controlling biodiversity patterns is a central issue in ecology. Although it is generally accepted that both deterministic and stochastic processes play important roles in community assembly, quantifying their relative importance is challenging. Here we propose a general mathematical framework to quantify ecological stochasticity under different situations in which deterministic factors drive the communities more similar or dissimilar than null expectation. An index, normalized stochasticity ratio (NST), was developed with 50% as the boundary point between more deterministic (50%) assembly. NST was tested with simulated communities by considering abiotic filtering, competition, environmental noise, and spatial scales. All tested approaches showed limited performance at large spatial scales or under very high environmental noise. However, in all of the other simulated scenarios, NST showed high accuracy (0.90 to 1.00) and precision (0.91 to 0.99), with averages of 0.37 higher accuracy (0.1 to 0.7) and 0.33 higher precision (0.0 to 1.8) than previous approaches. NST was also applied to estimate stochasticity in the succession of a groundwater microbial community in response to organic carbon (vegetable oil) injection. Our results showed that community assembly was shifted from more deterministic (NST = 21%) to more stochastic (NST = 70%) right after organic carbon input. As the vegetable oil was consumed, the community gradually returned to be more deterministic (NST = 27%). In addition, our results demonstrated that null model algorithms and community similarity metrics had strong effects on quantifying ecological stochasticity. Significance An ecological community is a dynamic complex system with a myriad of interacting species, which are controlled by various scale-dependent deterministic and stochastic forces. With rapid advances in genomics technologies, categorizing biological diversity, particularly microbial diversity, becomes relatively easy, but the great challenge is to disentangle the mechanisms controlling biological diversity. The general null model-based framework developed in this study provides an effective and robust tool to ecologists for quantitatively assessing ecological stochasticity. By highlighting the caveats such as model selection, similarity metrics, and spatial scales, this study provides guidance for appropriate use of null model-based approaches for examining community assembly processes. Although this framework was tested with microbial data, it should also be applicable to plant and animal ecology.

Diving into Deep Water: The Deepwater Horizon oil spill in the Gulf of Mexico was one of the largest oil spills on record. Its setting at the bottom of the sea floor posed an unanticipated risk as substantial amounts of hydrocarbons leaked into the deepwater column. Three separate cruises identified and sampled deep underwater hydrocarbon plumes that existed in May and June, 2010—before the well head was ultimately sealed. Camilli et al. (p. 201 ; published online 19 August) used an automated underwater vehicle to assess the dimensions of a stabilized, diffuse underwater plume of oil that was 22 miles long and estimated the daily quantity of oil released from the well, based on the concentration and dimensions of the plume. Hazen et al. (p. 204 ; published online 26 August) also observed an underwater plume at the same depth and found that hydrocarbon-degrading bacteria were enriched in the plume and were breaking down some parts of the oil. Finally, Valentine et al. (p. 208 ; published online 16 September) found that natural gas, including propane and ethane, were also present in hydrocarbon plumes. These gases were broken down quickly by bacteria, but primed the system for biodegradation of larger hydrocarbons, including those comprising the leaking crude oil. Differences were observed in dissolved oxygen levels in the plumes (a proxy for bacterial respiration), which may reflect differences in the location of sampling or the aging of the plumes.

Biodiversity and its responses to environmental changes are central issues in ecology and for society. Almost all microbial biodiversity research focuses on “species” richness and abundance but not on their interactions. Although a network approach is powerful in describing ecological interactions among species, defining the network structure in a microbial community is a great challenge. Also, although the stimulating effects of elevated CO ₂ (eCO ₂ ) on plant growth and primary productivity are well established, its influences on belowground microbial communities, especially microbial interactions, are poorly understood. Here, a random matrix theory (RMT)-based conceptual framework for identifying functional molecular ecological networks was developed with the high-throughput functional gene array hybridization data of soil microbial communities in a long-term grassland FACE (free air, CO ₂ enrichment) experiment. Our results indicate that RMT is powerful in identifying functional molecular ecological networks in microbial communities. Both functional molecular ecological networks under eCO ₂ and ambient CO ₂ (aCO ₂ ) possessed the general characteristics of complex systems such as scale free, small world, modular, and hierarchical. However, the topological structures of the functional molecular ecological networks are distinctly different between eCO ₂ and aCO ₂ , at the levels of the entire communities, individual functional gene categories/groups, and functional genes/sequences, suggesting that eCO ₂ dramatically altered the network interactions among different microbial functional genes/populations. Such a shift in network structure is also significantly correlated with soil geochemical variables. In short, elucidating network interactions in microbial communities and their responses to environmental changes is fundamentally important for research in microbial ecology, systems microbiology, and global change. IMPORTANCE Microorganisms are the foundation of the Earth's biosphere and play integral and unique roles in various ecosystem processes and functions. In an ecosystem, various microorganisms interact with each other to form complicated networks. Elucidating network interactions and their responses to environmental changes is difficult due to the lack of appropriate experimental data and an appropriate theoretical framework. This study provides a conceptual framework to construct interaction networks in microbial communities based on high-throughput functional gene array hybridization data. It also first documents that elevated carbon dioxide in the atmosphere dramatically alters the network interactions in soil microbial communities, which could have important implications in assessing the responses of ecosystems to climate change. The conceptual framework developed allows microbiologists to address research questions unapproachable previously by focusing on network interactions beyond the listing of, e.g., the number and abundance of species. Thus, this study could represent transformative research and a paradigm shift in microbial ecology.

Understanding the interactions among different species and their responses to environmental changes, such as elevated atmospheric concentrations of CO ₂ , is a central goal in ecology but is poorly understood in microbial ecology. Here we describe a novel random matrix theory (RMT)-based conceptual framework to discern phylogenetic molecular ecological networks using metagenomic sequencing data of 16S rRNA genes from grassland soil microbial communities, which were sampled from a long-term free-air CO ₂ enrichment experimental facility at the Cedar Creek Ecosystem Science Reserve in Minnesota. Our experimental results demonstrated that an RMT-based network approach is very useful in delineating phylogenetic molecular ecological networks of microbial communities based on high-throughput metagenomic sequencing data. The structure of the identified networks under ambient and elevated CO ₂ levels was substantially different in terms of overall network topology, network composition, node overlap, module preservation, module-based higher-order organization, topological roles of individual nodes, and network hubs, suggesting that the network interactions among different phylogenetic groups/populations were markedly changed. Also, the changes in network structure were significantly correlated with soil carbon and nitrogen contents, indicating the potential importance of network interactions in ecosystem functioning. In addition, based on network topology, microbial populations potentially most important to community structure and ecosystem functioning can be discerned. The novel approach described in this study is important not only for research on biodiversity, microbial ecology, and systems microbiology but also for microbial community studies in human health, global change, and environmental management. IMPORTANCE The interactions among different microbial populations in a community play critical roles in determining ecosystem functioning, but very little is known about the network interactions in a microbial community, owing to the lack of appropriate experimental data and computational analytic tools. High-throughput metagenomic technologies can rapidly produce a massive amount of data, but one of the greatest difficulties is deciding how to extract, analyze, synthesize, and transform such a vast amount of information into biological knowledge. This study provides a novel conceptual framework to identify microbial interactions and key populations based on high-throughput metagenomic sequencing data. This study is among the first to document that the network interactions among different phylogenetic populations in soil microbial communities were substantially changed by a global change such as an elevated CO ₂ level. The framework developed will allow microbiologists to address research questions which could not be approached previously, and hence, it could represent a new direction in microbial ecology research.

Disruption of healthy microbial communities has been linked to numerous diseases, yet microbial interactions are little understood. This is due in part to the large number of bacteria, and the much larger number of interactions (easily in the millions), making experimental investigation very difficult at best and necessitating the nascent field of computational exploration through microbial correlation networks. We benchmark the performance of eight correlation techniques on simulated and real data in response to challenges specific to microbiome studies: fractional sampling of ribosomal RNA sequences, uneven sampling depths, rare microbes and a high proportion of zero counts. Also tested is the ability to distinguish signals from noise, and detect a range of ecological and time-series relationships. Finally, we provide specific recommendations for correlation technique usage. Although some methods perform better than others, there is still considerable need for improvement in current techniques.