Ashley E. Beck
Manuel Kleiner2- 1Department of Biological and Environmental Sciences, Carroll College, Helena, MT, United States
- 2Department of Plant and Microbial Biology, North Carolina State University, Raleigh, NC, United States
With a growing world population and increasing frequency of climate disturbance events, we are in dire need of methods to improve plant productivity, resilience, and resistance to both abiotic and biotic stressors, both for agriculture and conservation efforts. Microorganisms play an essential role in supporting plant growth, environmental response, and susceptibility to disease. However, understanding the specific mechanisms by which microbes interact with each other and with plants to influence plant phenotypes is a major challenge due to the complexity of natural communities, simultaneous competition and cooperation effects, signalling interactions, and environmental impacts. Synthetic communities are a major asset in reducing the complexity of these systems by simplifying to dominant components and isolating specific variables for controlled experiments, yet there still remains a large gap in our understanding of plant microbiome interactions. This perspectives article presents a brief review discussing ways in which metabolic modelling can be used in combination with synthetic communities to continue progress toward understanding the complexity of plant-microbe-environment interactions. We highlight the utility of metabolic models as applied to a community setting, identify different applications for both flux balance and elementary flux mode simulation approaches, emphasize the importance of ecological theory in guiding data interpretation, and provide ideas for how the integration of metabolic modelling techniques with big data may bridge the gap between simplified synthetic communities and the complexity of natural plant-microbe systems.
Introduction
Plants have evolved intricate signalling networks sensing and responding to environmental stimuli. In recent years, the complex network of interactions existing between plants, microorganisms, and their environment has been recognized as an important factor impacting plant health and productivity (Timm et al., 2018; Kumar and Verma, 2019; Harman et al., 2021). Plants rely on microbes for biologically available forms of essential nutrients (Richardson et al., 2009); beneficial microbes also aid in protecting against pathogens and pests and are involved in plant response to environmental conditions such as heat and drought (Bhattacharyya et al., 2016; Liu et al., 2020).
The numerous multi-directional interactions that occur among plants, colonizing microorganisms, and environmental variables create difficulty in determining cause-and-effect relationships. For example, soil contains thousands of different microbial taxa (Lennon et al., 2012), a subset of which are recruited for colonization by plant secretion of root exudates (e.g., Lebeis et al., 2015); in addition to chemical signalling factors from roots, environmental factors such as soil chemical properties, pH, moisture content, and temperature also affect microbial colonization (Andrew et al., 2012; Islam et al., 2020b), influencing both microbe-microbe and plant-microbe interactions. Synthetic communities provide a reductionist approach to help constrain biological factors, examining a limited number of microbial species with key roles in the community to better understand driving forces and their overall effect on the productivity and resilience of plant-microbe systems (Vorholt et al., 2017). In the last decade, synthetic communities have become widely applied in a variety of contexts including agriculture (De Roy et al., 2014; De Souza et al., 2020; Sgobba and Wendisch, 2020), recognizing their value in experimental hypothesis testing.
Even with the tools provided by synthetic communities, an accurate understanding of multi-species, cross-domain interactions is still complicated by the complexity of each individual member’s metabolic network. Metabolic modelling provides mathematical predictions of metabolic routes used by an organism under different biological and environmental constraints and can thus help quantify the function of individual members and identify potential roles in community interaction on a metabolic level (Bosi et al., 2017). This technique complements experimental findings with a computational aspect that can aid in elucidating interactions that may elude solely experimental approaches. This perspectives article provides insight into how metabolic modelling can be integrated with experiments, using ecological theories to mine the complexities of plant-microbe-environment interactions and illuminate underpinning rules that lead to increased productivity and robustness in these systems.
Overview of Synthetic Communities in Plant Systems
Synthetic communities have been developed for a variety of plants to provide more tractable systems for controlled lab-scale experiments. Microbes have been isolated and cultured from several different plant species of interest, including the model plants Arabidopsis thaliana and Lotus japonicus, and agriculturally important species such as rice, maize, tomato, potato, clover, and sugarcane (Liu et al., 2019; Wippel et al., 2021), identifying key functional microbial taxa that may influence plant-microbe interactions. Synthetic microbial communities have been developed for both rhizosphere and phyllosphere communities in several species: most notably Arabidopsis, as well as maize, sorghum, and duckweed, among others, with varying levels of phylogenetic representation (Table 1). Varied approaches have been taken for designing synthetic communities with the ultimate goal of developing a community representative of the function and interactions of a natural, complex microbiome (e.g., Timm et al., 2016; Wang et al., 2021). Co-occurrence network analyses are commonly used to identify core functional microbiomes – a reduced subset of microbial species that captures the main metabolic processes performed by a community (Toju et al., 2020), but other methods recognize the inherent complexity of non-equilibrium microbe-environment interactions and incorporate time-series data with empirical methods (Deyle et al., 2016; Suzuki et al., 2017). Metabolic modelling provides a theoretical approach to complement the experimental process of first selecting and building a synthetic community, as well as predicting physiological properties of the community.
Table 1. Highlighted synthetic microbial communities and community metabolic models.
What Is Metabolic Modelling?
A metabolic model is a mathematical representation of the stoichiometries of metabolic reactions occurring within an organism, ranging in complexity from prokaryotes to eukaryotes. Models are generated from genomic data using databases such as KEGG, MetaCyc, JGI IMG/M, RAST, and BRENDA to identify genes and associated metabolic reactions, collect metabolite stoichiometries for each reaction, and compile the information into a metabolite-reaction stoichiometry matrix via software tools such as KBase, Python, or CellNetAnalyzer (von Kamp et al., 2017; Arkin et al., 2018; Medlock et al., 2020). The process of building a metabolic model is a type of experiment in itself, as model accuracy is influenced by genome completeness, quality of genomic data and annotations (which can result in missing genes, reactions, or pathways), and availability of experimental data (such as biomass composition and maintenance energy requirements). The outcome of metabolic model construction is a hypothesis that can be used to predict physiological response under different environmental circumstances.
Simulation methods vary in their mathematical approach to solving the system of linear equations produced by the metabolite-reaction stoichiometric matrix, but all assume a steady state environment (Sv = 0, where S is the stoichiometric matrix and v is the flux vector). This system of equations is typically underdetermined due to a greater number of participating metabolites than number of reactions (Klamt et al., 2002); therefore, infinitely many flux solutions may satisfy the system. Elementary flux mode analysis and flux balance analysis are two of the most commonly employed simulation approaches and differ in computational intensity and scope of the solutions obtained, providing different advantages for plant microbiome applications.
Elementary flux mode analysis is a computationally intensive approach that enumerates all possible minimal pathways that balance the stoichiometric matrix (i.e., if a single reaction were removed from a minimal flux vector, the route through the network would be incomplete). This approach has been likened to finding the most efficient path when navigating a subway network (Zanghellini et al., 2013). From the entire set of elementary flux modes, any possible metabolism of the organism can be described by a linear combination of the elementary flux modes, serving as a comprehensive and unbiased examination of the organism’s physiological potential. The number of modes calculated for a given model depends on the model structure (reaction connectivity and reversibility, exchange reactions, etc.); but computational expense typically limits applications to small networks, though recent algorithm developments have enabled computation of 12 billion elementary flux modes (Buchner and Zanghellini, 2021). For complex communities, feasibility may be improved by performing elementary flux mode analysis for each species model separately, individually analyzing results for biologically and ecologically relevant pathways (e.g., Beck et al., 2017). The curated results can then be transferred into a less expensive optimization method (like flux balance analysis described below), similar to the hybrid analysis method demonstrated in Hunt et al. (2018), to enhance simulation tractability for plant-microbe systems.
Alternatively, flux balance analysis is a less computationally intensive approach using a linear optimization strategy to solve the stoichiometric matrix. By applying constraints to cellular uptake and secretion reactions along with an objective criterion (typically maximization of biomass production), the metabolic network is optimized under a restricted space. While this method is more easily applied to large networks, it presents bias toward the researcher-determined objective criterion, and it often does not provide a unique solution as more than one flux vector may achieve the same optimized objective value due to network redundancy. Flux balance analysis provides great utility in specifying uptake and secretion rates of nutrients and metabolic byproducts and can be used to test hypotheses of metabolite exchange in silico, which is key to elucidating plant-microbe and microbe-microbe interactions. Several algorithmic variations have extended this approach to examine more complex cellular environments, allowing consideration of other constraints on metabolism that are also useful to plant systems, such as cellular volume and biosynthetic costs of enzyme assembly (Mori et al., 2016), intracellular heterogeneity (Damiani et al., 2017), balance of resource usage (Goelzer and Fromion, 2011), and dynamic growth and integration with consumer-resource models (Mahadevan et al., 2002; Gowda et al., 2022).
How Can Metabolic Modelling Be Applied to Understand Multi-Organism Interactions?
While metabolic modelling was initially employed for singular species, increased availability of genomic data and development of more efficient computational algorithms has enabled application to community models. Stolyar et al. (2007) were the first to extend flux balance analysis to a multispecies context, using a multi-compartment flux balance model to predict metabolite exchange between the sulfate reducer Desulfovibrio vulgaris and the methanogen Methanococcus maripaludis, thereby highlighting the important role of hydrogen in syntrophy between the two microbes. Since this pioneering study, metabolic models have been used to investigate potential interactions within several different microbial community contexts, with a few recent key studies highlighted in Table 1.
Plants present a unique challenge for the construction of high-quality metabolic models due to large eukaryotic genomes, multiple differentiated tissues, and redundancy associated with polyploidy in many plant genomes. As more genomic data is collected, current efforts to improve databases, annotations, and primary and secondary pathways (e.g., PlantSEED; Seaver et al., 2018) provide automated platforms as a starting point for plant-specific modelling. To date, metabolic models have been developed for specialized tissue types of several different plant species, including soybean seed (Moreira et al., 2019), rapeseed (Pilalis et al., 2011), rice (Chatterjee and Kundu, 2015), potato (Botero et al., 2018), and maize (Saha et al., 2011). Integration of multiple compartments within a metabolic model has also allowed reconstruction of multi-tissue models for more accurate representation of plant metabolism (Shaw and Cheung, 2020); for example, one of the most comprehensive plant models developed thus far has encompassed six different key tissues for Arabidopsis to more accurately predict whole-plant physiological responses (Gomes de Oliveira Dal’Molin et al., 2015). Additionally, multi-scale models (e.g., incorporation of gene regulation or phenomic data) can improve accuracy and provide experimental validation, aiding in understanding physiological effects of factors such as genome redundancy and circadian rhythm (Jez et al., 2021; Krantz et al., 2021).
Despite these advances in modelling applications for both plants and microorganisms, we found no studies investigating plant-microbe interactions with a community metabolic model approach (i.e., integrating both plant and microbial models into the same simulation). Recent algal-bacterial community models are moving in this direction to study more complex phototroph-heterotroph and eukaryote-prokaryote interactions (Dittami et al., 2014; Gomez et al., 2021). Furthermore, there still exists a large gap between modelling, lab-scale experimental results, and field-scale results, offering abundant opportunity for advances in this area. The sections below outline four key areas in which modelling could be used to enhance plant microbiome studies harnessing the utility of synthetic communities, followed by overall challenges and opportunities moving forward.
Disentangling Multi-Directional Microbe–Microbe Interactions
Interspecies interactions can take a variety of forms, e.g., mutualistic, commensal, antagonistic, parasitic (Holland and DeAngelis, 2009). Often, interactions are based on metabolite exchange which may be the basis for either an obligate or facultative relationship (Zelezniak et al., 2015; Beck et al., 2016). Even in simplified synthetic communities with few members, interactions are still complex and multifaceted. Interactions can sometimes be teased apart with carefully designed experiments; however, experimental measurements are not always able to distinguish the specific donor and recipient of metabolite exchanges. Modelling therefore allows for in silico testing of a multitude of possible unidirectional and bidirectional interactions within a community in a more rapid time frame. For example, a systematic approach for dissecting community interactions may involve first selecting specific pairs of organisms and predicting potential metabolite exchange-based interactions through a dual member model. Comparing the pairwise interaction predictions with simulations of a larger community will then aid in understanding how interactions could change with additional community member(s) (Figure 1).
Figure 1. Schematic illustrating proposed metabolic modelling workflow with synthetic communities as input. The natural microbiome is a complex environment and is simplified to a synthetic community with a tractable number of microorganisms via experimental isolation techniques and/or metagenomics methods. Genomes for this limited number of microorganisms are used to construct metabolic models for each species, which can then be used with pathway simulation tools (flux balance analysis, FBA; elementary flux mode analysis, EFMA) to predict both pairwise and community interactions to observe changes due to additional community members and determine patterns of interaction. In the example shown, the blue microbe secretes a compound that is beneficial to both the red and green microbes and stabilizes their original negative pairwise interaction. Complementarily, culturing experiments can be used to examine the responses of microbial communities (both in planta and in vitro) to environmental variables such as limited nutrients or water stress. Collecting molecular data under specific experimental conditions provides data to integrate with modelling predictions, which can be interpreted with the aid of relevant ecological theory, such as resource ratio theory, resource allocation theory, or the maximum power principle. This process leads to iterative refinement of both models and experimental design, ultimately contributing to broader outcomes relevant to agricultural practice, such as the design of microbial inoculants to promote plant growth and resilience, an improved understanding of plant-microbe response to the environment under changing climate conditions, and the implementation of field-scale trials to further test interaction principles in the natural environment.
Modelling Plant-Specific Tissues and Community Assembly Processes
Rhizosphere and phyllosphere communities are complex assemblages of many species of organisms with the host plant, including not only bacteria but also archaea, fungi, and protozoa. The mechanism of community assembly is a major research question, including the effect of environmental stimuli on the assembly process, e.g., soil quality, nutrient availability, temperature, moisture, and presence of other organisms like pests and pathogens. Recruitment of microorganisms to a tissue involves metabolite exudation by the plant, but the ways in which plant-microbe and microbe-microbe interactions impact this process are currently poorly understood. Integration of plant and microbial metabolic models in a dynamic format (e.g., use of dynamic flux balance algorithms to investigate time-resolved interaction effects) will aid in better understanding and predicting the assembly process, particularly in understanding what drives colonization differences among plant species as well as how pathogens may disrupt the colonization process.
Interaction Effects Resulting From Abiotic and Biotic Stresses
Beyond initial community assembly, environmental stresses are also known to have profound impacts on plant physiological responses. The literature is filled with examples demonstrating effects on growth and productivity by factors such as drought, temperature, and nutrient limitation (Cramer et al., 2011; Fahad et al., 2017). However, research has focused primarily on plant molecular and physiological responses, whereas many of these responses are mediated by plant-microbe and microbe-microbe interactions. Our understanding of how interactions change under different environmental stressors can benefit from modelling investigations, which can test a particular environment as part of the research question. For example, in flux balance analysis, substrate availability or uptake can be restricted according to the desired environmental scenario. Elementary flux mode analysis commonly examines pathway efficiency, and under nutrient limitations, cells are typically presumed to use the limited nutrient more efficiently relative to other substrates, guiding the simulation analysis. Careful design of simulation environments can be used in conjunction with greenhouse or laboratory experiments to gain a better understanding of how environmental pressures impact plant-microbe systems via community interaction effects. Biotic stressors can be investigated similarly (e.g., expanding a community model to include a pathogen).
Incorporating Community Ecology Theory in Metabolic Modelling
Metabolic modelling uses genomic data combined with biochemical reactions governed by steady state conservation of mass to predict metabolite flux through a cell. However, this methodology further provides a platform for incorporating ecological theory in interpreting simulation results. For example, resource competition and allocation theories (Harcombe et al., 2014; Bernhardt et al., 2020) can be applied to predict how an organism might optimize its use of a certain resource, followed by comparison with efficient resource utilization predicted for competing species or the entire community. Similarly, resource ratio theory (de Mazancourt and Schwartz, 2010) examines tradeoffs between two key nutrients affecting growth, which can be expanded to examine resource usage within a community, as well (e.g., influence of microbial community composition on plant tissue allocation as observed in Henning et al., 2019). The maximum power principle, an alternative theory focused on energy usage, states that communities are evolutionarily organized to maximize energy consumption from the environment (Lotka, 1922; Sciubba, 2011). This principle has been less commonly applied in microbial communities but has been used to examine competition strategies (Cai et al., 2006; DeLong, 2008) and can provide another framework with which to evaluate the efficiency of community organization. Integrating these different community ecology theories can aid in elucidating the balance between competitive and cooperative interactions.
Key Challenges in Plant Microbiome Modelling
Community models present a unique challenge as compared to individual species models when designing model simulations, specifically with regard to determining the optimization criterion that best represents the community objective in a natural system, as accurately modelling the overall community “goal” can be challenging. Typically for single species models, growth rate is used as the optimization criterion to determine a metabolic route that promotes maximum biomass formation, which may or may not be reflective of actual cellular targets, depending upon the environment. Many natural environments are nutrient-poor and may not promote growth at the highest capacity; desirable metabolic routes may instead optimize for efficient use of a limiting resource like nitrogen or for efficient energy production. Many current community modelling configurations optimize the overall community growth rate (Zomorrodi and Maranas, 2012). According to ecological principles, however, communities may not always be operating to achieve the maximum total biomass; thus, the objectives relevant to the specific scenario must be carefully decided upon when designing a model simulation. Objectives should be based on ecological principles and should account for the conditions or treatments being investigated.
Experimental validation of model predictions is essential due to the complexity of microbial communities and the underdetermined nature of individual metabolic networks which allows multiple possible solutions. With the broader availability of omics technologies and the ability to generate large data sets at lower costs, it is now possible to obtain multiple layers of data to complement modelling studies. Actualization of genomic potential can be examined at many levels – e.g., metatranscriptomics (Bashiardes et al., 2016), metaproteomics (Kleiner, 2019; Salvato et al., 2021), metabolomics, and plant phenomics. These different expression levels are important for understanding experimentally how the intersection of multiple metabolic networks affects the plant, as using one data type alone may yield very different predictions (Blazier and Papin, 2012; Payne, 2015). A key area of algorithm development is the incorporation of omics data within a metabolic network to influence optimization results to more accurately reflect and predict the outcome of a particular experimental scenario (e.g., Jungreuthmayer et al., 2013; Gerstl et al., 2015; Tian and Reed, 2018).
Conclusion and Outlook
The above narrative addresses the current state of metabolic modelling as applied specifically to plant-microbe-environment interactions, detailing some of the key challenges in expanding modelling applications to these complex systems. Merging computational and experimental approaches will improve our limited understanding of multi-faceted plant-microbe-environment relationships (Figure 1). Advances will enable agricultural improvements, such as development of microbial inoculants to promote plant growth or resilience in the face of increasing global climate events. With these efforts, maintaining quality standards for modelling across the scientific community (Carey et al., 2020; Lieven et al., 2020) will continue to be essential when extending models and algorithms to more complex plant-microbe systems.
Envisioned future applications (Figure 1) involve starting from a relatively small synthetic microbial community (e.g., less than 10 members), developing a genome-enabled model for each member, validating individual models with experimental data (as much as possible), and investigating in silico pairwise interactions between species. Successively larger community models can be constructed by including additional member(s) and comparing with the previous pairwise models to observe how predictions change as more members are added. Laboratory and greenhouse experiments testing different community compositions, environmental factors, or other variables of interest can be used to verify and refine the community model, with the ultimate goal of moving to more complex systems (e.g., field-scale validation in agricultural settings). With continued effort and progress, the general principles and pattern of interactions uncovered through this systematic process is anticipated to be transferable to other plant systems; the microbes involved in the interactions may be different but likely follow similar governing interaction rules and can be used to engineer agricultural solutions in many different crops.
Data Availability Statement
The original contributions presented in the study are included in the article/supplementary material; further inquiries can be directed to the corresponding author.
Author Contributions
AB conceived the study, performed literature review, and prepared the initial draft of the manuscript. MK and A-KG performed literature review and curated the draft manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work is supported by Agriculture and Food Research Initiative Agricultural Microbiomes grant no. 2021-67013-34537 from the USDA National Institute of Food and Agriculture.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Andrew, D. R., Fitak, R. R., Munguia-Vega, A., Racolta, A., Martinson, V. G., and Dontsova, K. (2012). Abiotic factors shape microbial diversity in Sonoran Desert soils. Appl. Environ. Microbiol. 78, 7527–7537. doi: 10.1128/AEM.01459-12
Arkin, A. P., Cottingham, R. W., Henry, C. S., Harris, N. L., Stevens, R. L., Maslov, S., et al. (2018). KBase: the United States Department of Energy systems biology knowledgebase. Nat. Biotechnol. 36, 566–569. doi: 10.1038/nbt.4163
Bashiardes, S., Zilberman-Schapira, G., and Elinav, E. (2016). Use of metatranscriptomics in microbiome research. Bioinform. Biol. Insights. 10:BBI.S34610. doi: 10.4137/BBI.S34610
Beck, A. E., Bernstein, H. C., and Carlson, R. P. (2017). Stoichiometric network analysis of cyanobacterial acclimation to photosynthesis-associated stresses identifies heterotrophic niches. PRO 5:32. doi: 10.3390/pr5020032
Beck, A. E., Hunt, K. A., Bernstein, H. C., and Carlson, R. P. (2016). “Interpreting and designing microbial communities for bioprocess applications, from components to interactions to emergent properties,” in Biotechnology for Biofuel Production and Optimization. eds. C. A. Eckert and C. T. Trinh (Amsterdam, Netherlands: Elsevier), 407–432.
Bernhardt, J. R., Kratina, P., Pereira, A. L., Tamminen, M., Thomas, M. K., and Narwani, A. (2020). The evolution of competitive ability for essential resources. Philos. Trans. R. Soc. 375:247. doi: 10.1098/rstb.2019.0247
Bhattacharyya, P. N., Goswami, M. P., and Bhattacharyya, L. H. (2016). Perspective of beneficial microbes in agriculture under changing climatic scenario: a review. J. Phytology 8, 26–41. doi: 10.19071/jp.2016.v8.3022
Blazier, A. S., and Papin, J. A. (2012). Integration of expression data in genome-scale metabolic network reconstructions. Front. Physiol. 3:299. doi: 10.3389/fphys.2012.00299
Bosi, E., Bacci, G., Mengoni, A., and Fondi, M. (2017). Perspectives and challenges in microbial communities metabolic modeling. Front. Genet. 8:88. doi: 10.3389/fgene.2017.00088
Botero, K., Restrepo, S., and Pinzón, A. (2018). A genome-scale metabolic model of potato late blight suggests a photosynthesis suppression mechanism. BMC Genomics 19, 31–44. doi: 10.1186/s12864-018-5192-x
Buchner, B. A., and Zanghellini, J. (2021). EFMlrs: a Python package for elementary flux mode enumeration via lexicographic reverse search. BMC Bioinform. 22, 1–21. doi: 10.1186/s12859-021-04417-9
Cai, T. T., Montague, C. L., and Davis, J. S. (2006). The maximum power principle: an empirical investigation. Ecol. Model. 190, 317–335. doi: 10.1016/j.ecolmodel.2005.04.022
Carey, M. A., Dräger, A., Beber, M. E., Papin, J. A., and Yurkovich, J. T. (2020). Community standards to facilitate development and address challenges in metabolic modelling. Mol. Syst. Biol. 16:e9235. doi: 10.15252/msb.20199235
Carlström, C. I., Field, C. M., Bortfeld-Miller, M., Müller, B., Sunagawa, S., and Vorholt, J. A. (2019). Synthetic microbiota reveal priority effects and keystone strains in the Arabidopsis phyllosphere. Nat. Ecol. Evol. 3, 1445–1454. doi: 10.1038/s41559-019-0994-z
Chai, Y. N., Ge, Y., Stoerger, V., and Schachtman, D. P. (2021). High-resolution phenotyping of sorghum genotypic and phenotypic responses to low nitrogen and synthetic microbial communities. Plant Cell Environ. 44, 1611–1626. doi: 10.1111/pce.14004
Chan, S. H., Friedman, E. S., Wu, G. D., and Maranas, C. D. (2019). Predicting the longitudinally and radially varying gut microbiota composition using multi-scale microbial metabolic modelling. PRO 7:394. doi: 10.3390/pr7070394
Chatterjee, A., and Kundu, S. (2015). Revisiting the chlorophyll biosynthesis pathway using genome scale metabolic model of Oryza sativa japonica. Sci. Rep. 5, 1–15. doi: 10.1038/srep14975
Cramer, G. R., Urano, K., Delrot, S., Pezzotti, M., and Shinozaki, K. (2011). Effects of abiotic stress on plants: a systems biology perspective. BMC Plant Biol. 11:163. doi: 10.1186/1471-2229-11-163
Damiani, C., Di Filippo, M., Pescini, D., Maspero, D., Colombo, R., and Mauri, G. (2017). popFBA: tackling intratumour heterogeneity with flux balance analysis. Bioinformatics 33, i311–i318. doi: 10.1093/bioinformatics/btx251
De Mazancourt, C., and Schwartz, M. W. (2010). A resource ratio theory of cooperation. Ecol. Lett. 13, 349–359. doi: 10.1111/j.1461-0248.2009.01431.x
De Roy, K., Marzorati, M., Van den Abbeele, P., Van de Wiele, T., and Boon, N. (2014). Synthetic microbial ecosystems: an exciting tool to understand and apply microbial communities. Environ. Microbiol. 16, 1472–1481. doi: 10.1111/1462-2920.12343
De Souza, R. S. C., Armanhi, J. S. L., and Arruda, P. (2020). From microbiome to traits: designing synthetic microbial communities for improved crop resiliency. Front. Plant Sci. 11:1179. doi: 10.3389/fpls.2020.01179
DeLong, J. P. (2008). The maximum power principle predicts the outcomes of two-species competition experiments. Oikos 117, 1329–1336. doi: 10.1111/j.0030-1299.2008.16832.x
Deyle, E. R., May, R. M., Munch, S. B., and Sugihara, G. (2016). Tracking and forecasting ecosystem interactions in real time. Proc. R. Soc. B Biol. Sci. 283:2258. doi: 10.1098/rspb.2015.2258
Dittami, S. M., Eveillard, D., and Tonon, T. (2014). A metabolic approach to study algal–bacterial interactions in changing environments. Mol. Ecol. 23, 1656–1660. doi: 10.1111/mec.12670
Fahad, S., Bajwa, A. A., Nazir, U., Anjum, S. A., Farooq, A., Zohaib, A., et al. (2017). Crop production under drought and heat stress: plant responses and management options. Front. Plant Sci. 8:1147. doi: 10.3389/fpls.2017.01147
Finkel, O. M., Salas-González, I., Castrillo, G., Spaepen, S., Law, T. F., Teixeira, P. J. P. L., et al. (2019). The effects of soil phosphorus content on plant microbiota are driven by the plant phosphate starvation response. PLoS Biol. 17:e3000534. doi: 10.1371/journal.pbio.3000534
Gerstl, M. P., Jungreuthmayer, C., and Zanghellini, J. (2015). tEFMA: computing thermodynamically feasible elementary flux modes in metabolic networks. Bioinformatics 31, 2232–2234. doi: 10.1093/bioinformatics/btv111
Goelzer, A., and Fromion, V. (2011). Bacterial growth rate reflects a bottleneck in resource allocation. Biochim. Biophys. Acta Gen. Subj. 1810, 978–988. doi: 10.1016/j.bbagen.2011.05.014
Gomes de Oliveira Dal’Molin, C., Quek, L. E., Saa, P. A., and Nielsen, L. K. (2015). A multi-tissue genome-scale metabolic modelling framework for the analysis of whole plant systems. Front. Plant Sci. 6:4. doi: 10.3389/fpls.2015.00004
Gomez, J. A., Höffner, K., and Barton, P. I. (2021). Production of biofuels from sunlight and lignocellulosic sugars using microbial consortia. Chem. Eng. Sci. 239:116615. doi: 10.1016/j.ces.2021.116615
Gowda, K., Ping, D., Mani, M., and Kuehn, S. (2022). Genomic structure predicts metabolite dynamics in microbial communities. Cell 185, 530.e25–546.e25. doi: 10.1016/j.cell.2021.12.036
Harcombe, W. R., Riehl, W. J., Dukovski, I., Granger, B. R., Betts, A., Lang, A. H., et al. (2014). Metabolic resource allocation in individual microbes determines ecosystem interactions and spatial dynamics. Cell Rep. 7, 1104–1115. doi: 10.1016/j.celrep.2014.03.070
Harman, G., Khadka, R., Doni, F., and Uphoff, N. (2021). Benefits to plant health and productivity from enhancing plant microbial symbionts. Front. Plant Sci. 11:610065. doi: 10.3389/fpls.2020.610065
Henning, J. A., Weston, D. J., Pelletier, D. A., Timm, C. M., Jawdy, S. S., and Classen, A. T. (2019). Relatively rare root endophytic bacteria drive plant resource allocation patterns and tissue nutrient concentration in unpredictable ways. Am. J. Bot. 106, 1423–1434. doi: 10.1002/ajb2.1373
Holland, N. J., and DeAngelis, D. L. (2009). Consumer-resource theory predicts dynamic transitions between outcomes of interspecific interactions. Ecol. Lett. 12, 1357–1366. doi: 10.1111/j.1461-0248.2009.01390.x
Hunt, K. A., Jennings, R. M., Inskeep, W. P., and Carlson, R. P. (2018). Multiscale analysis of autotroph-heterotroph interactions in a high-temperature microbial community. PLoS Comput. Biol. 14:e1006431. doi: 10.1371/journal.pcbi.1006431
Ishizawa, H., Tada, M., Kuroda, M., Inoue, D., Futamata, H., and Ike, M. (2020). Synthetic bacterial community of duckweed: a simple and stable system to study plant-microbe interactions. Microbes Environ. 35:ME20112. doi: 10.1264/jsme2.ME20112
Islam, M. M., Le, T., Daggumati, S. R., and Saha, R. (2020a). Investigation of microbial community interactions between Lake Washington methanotrophs using genome-scale metabolic modelling. PeerJ 8:e9464. doi: 10.7717/peerj.9464
Islam, W., Noman, A., Naveed, H., Huang, Z., and Chen, H. Y. (2020b). Role of environmental factors in shaping the soil microbiome. Environ. Sci. Pollut. Res. 27, 41225–41247. doi: 10.1007/s11356-020-10471-2
Jez, J. M., Topp, C. N., Matthews, M. L., and Marshall-Colón, A. (2021). Multiscale plant modeling: From genome to phenome and beyond. Emerg. Top Life Sci. 5, 231–237. doi: 10.1042/ETLS20200276
Jungreuthmayer, C., Ruckerbauer, D. E., and Zanghellini, J. (2013). regEfmtool: speeding up elementary flux mode calculation using transcriptional regulatory rules in the form of three-state logic. Biosystems 113, 37–39. doi: 10.1016/j.biosystems.2013.04.002
Klamt, S., Schuster, S., and Gilles, E. D. (2002). Calculability analysis in underdetermined metabolic networks illustrated by a model of the central metabolism in purple nonsulfur bacteria. Biotechnol. Bioeng. 77, 734–751. doi: 10.1002/bit.10153
Kleiner, M. (2019). Metaproteomics: much more than measuring gene expression in microbial communities. mSystems 4, e00115–e00119. doi: 10.1128/mSystems.00115-19
Krantz, M., Zimmer, D., Adler, S. O., Kitashova, A., Klipp, E., Mühlhaus, T., et al. (2021). Data management and modeling in plant biology. Front. Plant Sci. 12:717958. doi: 10.3389/fpls.2021.717958
Kumar, R. K., Singh, N. K., Balakrishnan, S., Parker, C. W., Raman, K., and Venkateswaran, K. (2021). Metabolic modelling of the international space station microbiome reveals key microbial interactions. BioRxiv [Preprint]. doi: 10.1101/2021.09.03.458819
Kumar, A., and Verma, J. P. (2019). “The role of microbes to improve crop productivity and soil health,” in Ecological Wisdom Inspired Restoration Engineering eds. V. Achal and A. Mukherjee (Singapore: Springer), 249–265.
Lebeis, S. L., Paredes, S. H., Lundberg, D. S., Breakfield, N., Gehring, J., McDonald, M., et al. (2015). Salicylic acid modulates colonization of the root microbiome by specific bacterial taxa. Science 349, 860–864. doi: 10.1126/science.aaa8764
Lennon, J. T., Aanderud, Z. T., Lehmkuhl, B. K., and Schoolmaster, D. R. Jr. (2012). Mapping the niche space of soil microorganisms using taxonomy and traits. Ecology. 93, 1867–1879. doi: 10.1890/11-1745.1
Lieven, C., Beber, M. E., Olivier, B. G., Bergmann, F. T., Ataman, M., Babaei, P., et al. (2020). MEMOTE for standardized genome-scale metabolic model testing. Nat. Biotechnol. 38, 272–276. doi: 10.1038/s41587-020-0446-y
Liu, H., Brettell, L. E., Qiu, Z., and Singh, B. K. (2020). Microbiome-mediated stress resistance in plants. Trends Plant Sci. 25, 733–743. doi: 10.1016/j.tplants.2020.03.014
Liu, Y. X., Qin, Y., and Bai, Y. (2019). Reductionist synthetic community approaches in root microbiome research. Curr. Opin. Microbiol. 49, 97–102. doi: 10.1016/j.mib.2019.10.010
Lotka, A. J. (1922). Contribution to the energetics of evolution. PNAS 8, 147–151. doi: 10.1073/pnas.8.6.147
Mahadevan, R., Edwards, J. S., and Doyle, F. J. III (2002). Dynamic flux balance analysis of diauxic growth in Escherichia coli. Biophys. J. 83, 1331–1340. doi: 10.1016/S0006-3495(02)73903-9
Medlock, G. L., Moutinho, T. J., and Papin, J. A. (2020). Medusa: software to build and analyze ensembles of genome-scale metabolic network reconstructions. PLoS Comput. Biol. 16:e1007847. doi: 10.1371/journal.pcbi.1007847
Moreira, T. B., Shaw, R., Luo, X., Ganguly, O., Kim, H. S., Coelho, L. G. F., et al. (2019). A genome-scale metabolic model of soybean (Glycine max) highlights metabolic fluxes in seedlings. Plant Physiol. 180, 1912–1929. doi: 10.1104/pp.19.00122
Mori, M., Hwa, T., Martin, O. C., De Martino, A., and Marinari, E. (2016). Constrained allocation flux balance analysis. PLoS Comput. Biol. 12:e1004913. doi: 10.1371/journal.pcbi.1004913
Niu, B., Paulson, J. N., Zheng, X., and Kolter, R. (2017). Simplified and representative bacterial community of maize roots. PNAS 114, E2450–E2459. doi: 10.1073/pnas.1616148114
Payne, S. H. (2015). The utility of protein and mRNA correlation. Trends Biochem. Sci. 40, 1–3. doi: 10.1016/j.tibs.2014.10.010
Pilalis, E., Chatziioannou, A., Thomasset, B., and Kolisis, F. (2011). An in silico compartmentalized metabolic model of Brassica napus enables the systemic study of regulatory aspects of plant central metabolism. Biotechnol. Bioeng. 108, 1673–1682. doi: 10.1002/bit.23107
Richardson, A. E., Barea, J. M., McNeill, A. M., and Prigent-Combaret, C. (2009). Acquisition of phosphorus and nitrogen in the rhizosphere and plant growth promotion by microorganisms. Plant Soil 321, 305–339. doi: 10.1007/s11104-009-9895-2
Saha, R., Suthers, P. F., and Maranas, C. D. (2011). Zea mays iRS1563: a comprehensive genome-scale metabolic reconstruction of maize metabolism. PLoS One 6:e21784. doi: 10.1371/journal.pone.0021784
Salvato, F., Hettich, R. L., and Kleiner, M. (2021). “Five key aspects of metaproteomics as a tool to understand functional interactions in host associated microbiomes” in Advances in Clinical Immunology, Medical Microbiology, COVID-19, and Big Data (Singapore: Jenny Stanford Publishing), 647–660.
Sciubba, E. (2011). What did Lotka really say? A critical reassessment of the “maximum power principle”. Ecol. Model. 222, 1347–1353. doi: 10.1016/j.ecolmodel.2011.02.002
Seaver, S. M., Lerma-Ortiz, C., Conrad, N., Mikaili, A., Sreedasyam, A., Hanson, A. D., et al. (2018). PlantSEED enables automated annotation and reconstruction of plant primary metabolism with improved compartmentalization and comparative consistency. Plant J. 95, 1102–1113. doi: 10.1111/tpj.14003
Sgobba, E., and Wendisch, V. F. (2020). Synthetic microbial consortia for small molecule production. Curr. Opin. Biotechnol. 62, 72–79. doi: 10.1016/j.copbio.2019.09.011
Shaw, R., and Cheung, C. Y. (2020). Multi-tissue to whole plant metabolic modelling. Cell. Mol. Life Sci. 77, 489–495. doi: 10.1007/s00018-019-03384-y
Stolyar, S., Van Dien, S., Hillesland, K. L., Pinel, N., Lie, T. J., Leigh, J. A., et al. (2007). Metabolic modelling of a mutualistic microbial community. Mol. Syst. Biol. 3, 92. doi: 10.1038/msb4100131
Suzuki, K., Yoshida, K., Nakanishi, Y., and Fukuda, S. (2017). An equation-free method reveals the ecological interaction networks within complex microbial ecosystems. Methods Ecol. Evol. 8, 1774–1785. doi: 10.1111/2041-210X.12814
Tian, M., and Reed, J. L. (2018). Integrating proteomic or transcriptomic data into metabolic models using linear bound flux balance analysis. Bioinformatics 34, 3882–3888. doi: 10.1093/bioinformatics/bty445
Timm, C. M., Carter, K. R., Carrell, A. A., Jun, S. R., Jawdy, S. S., Vélez, J. M., et al. (2018). Abiotic stresses shift belowground Populus-associated bacteria toward a core stress microbiome. mSystems 3, e00070–e00017. doi: 10.1128/mSystems.00070-17
Timm, C. M., Pelletier, D. A., Jawdy, S. S., Gunter, L. E., Henning, J. A., Engle, N., et al. (2016). Two poplar-associated bacterial isolates induce additive favorable responses in a constructed plant-microbiome system. Front. Plant Sci. 7:497. doi: 10.3389/fpls.2016.00497
Toju, H., Abe, M. S., Ishii, C., Hori, Y., Fujita, H., and Fukuda, S. (2020). Scoring species for synthetic community design: network analyses of functional core microbiomes. Front. Microbiol. 11:1361. doi: 10.3389/fmicb.2020.01361
von Kamp, A., Thiele, S., Hädicke, O., and Klamt, S. (2017). Use of CellNetAnalyzer in biotechnology and metabolic engineering. J. Biotechnol. 261, 221–228. doi: 10.1016/j.jbiotec.2017.05.001
Vorholt, J. A., Vogel, C., Carlström, C. I., and Müller, D. B. (2017). Establishing causality: opportunities of synthetic communities for plant microbiome research. Cell Host Microbe 22, 142–155. doi: 10.1016/j.chom.2017.07.004
Wang, J., Carper, D. L., Burdick, L. H., Shrestha, H. K., Appidi, M. R., Abraham, P. E., et al. (2021). Formation, characterization and modeling of emergent synthetic microbial communities. Comput. Struct. Biotechnol. J. 19, 1917–1927. doi: 10.1016/j.csbj.2021.03.034
Weiss, A. S., Burrichter, A. G., Durai Raj, A. C., von Strempel, A., Meng, C., Kleigrewe, K., et al. (2021). In vitro interaction network of a synthetic gut bacterial community. ISME J., 1–15. doi: 10.1038/s41396-021-01153-z
Wippel, K., Tao, K., Niu, Y., Zgadzaj, R., Kiel, N., Guan, R., et al. (2021). Host preference and invasiveness of commensal bacteria in the Lotus and Arabidopsis root microbiota. Nat. Microbiol. 6, 1150–1162. doi: 10.1038/s41564-021-00941-9
Zelezniak, A., Andrejev, S., Ponomarova, O., Mende, D. R., Bork, P., and Patil, K. R. (2015). Metabolic dependencies drive species co-occurrence in diverse microbial communities. PNAS 112, 6449–6454. doi: 10.1073/pnas.1421834112
Zanghellini, J., Ruckerbauer, D. E., Hanscho, M., and Jungreuthmayer, C. (2013). Elementary flux modes in a nutshell: properties, calculation and applications. Biotechnol. J. 8, 1009–1016. doi: 10.1002/biot.201200269
Keywords: synthetic communities, plant microbiome, plant microbial interactions, metabolic modelling, flux balance analysis, elementary flux mode analysis
Citation: Beck AE, Kleiner M and Garrell A-K (2022) Elucidating Plant-Microbe-Environment Interactions Through Omics-Enabled Metabolic Modelling Using Synthetic Communities. Front. Plant Sci. 13:910377. doi: 10.3389/fpls.2022.910377
Edited by:
Zoran Nikoloski, Max Planck Institute of Molecular Plant Physiology, GermanyReviewed by:
Collin M. Timm, Johns Hopkins University, United StatesCopyright © 2022 Beck, Kleiner and Garrell. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ashley E. Beck, abeck@carroll.edu