The biological carbon pump is the process by which CO 2 is transformed to organic carbon via photosynthesis, exported through sinking particles, and finally sequestered in the deep ocean. While the intensity of the pump correlates with plankton community composition, the underlying ecosystem structure driving the process remains largely uncharacterised. Here we use environmental and metagenomic data gathered during the Tara Oceans expedition to improve our understanding of carbon export in the oligotrophic ocean. We show that specific plankton communities, from the surface and deep chlorophyll maximum, correlate with carbon export at 150 m and highlight unexpected taxa such as Radiolaria, alveolate parasites, as well as Synechococcus and their phages, as lineages most strongly associated with carbon export in the subtropical, nutrient-depleted, oligotrophic ocean. Additionally, we show that the relative abundance of just a few bacterial and viral genes can predict most of the variability in carbon export in these regions. Guidi et al. Page 2 Nature. Author manuscript; available in PMC 2016 September 22. Europe PMC Funders Author Manuscripts Europe PMC Funders Author ManuscriptsMarine planktonic photosynthetic organisms are responsible for approximately fifty percent of Earth's primary production and fuel the global ocean biological carbon pump 1 . The intensity of the pump is correlated to plankton community composition 2,3 , and controlled by the relative rates of primary production and carbon remineralisation 4 . About 10% of this newly produced organic carbon in the surface ocean is exported through gravitational sinking of particles. Finally, after multiple transformations, a fraction of the exported material reaches the deep ocean where it is sequestered over thousand-year timescales 5 .Like most biological systems, marine ecosystems in the sunlit upper layer of the ocean (denoted the euphotic zone) are complex 6,7 , characterised by a wide range of biotic and abiotic interactions [8][9][10] and in constant balance between carbon production, transfer to higher trophic levels, remineralisation, and export to the deep layers 11 . The marine ecosystem structure and its taxonomic and functional composition likely evolved to comply with this loss of energy by modifying organism turnover times and by the establishment of complex feedbacks between them 6 and the substrates they can exploit for metabolism 12 .Decades of groundbreaking research have focused on identifying independently the key players involved in the biological carbon pump. Among autotrophs, diatoms are commonly attributed to being important in carbon flux because of their large size and fast sinking rates 13-15 while small autotrophic picoplankton may contribute directly through subduction of surface water 16 or indirectly by aggregating with larger settling particles or consumption by organisms at higher trophic levels 17 . Among heterotrophs, zooplankton such as crustaceans impact carbon flux via production of fast-sinking fecal pellets...

Natural genetic diversity provides a powerful tool to study the complex interrelationship between metabolism and growth. Profiling of metabolic traits combined with network-based and statistical analyses allow the comparison of conditions and identification of sets of traits that predict biomass. However, it often remains unclear why a particular set of metabolites is linked with biomass and to what extent the predictive model is applicable beyond a particular growth condition. A panel of 97 genetically diverse Arabidopsis (Arabidopsis thaliana) accessions was grown in near-optimal carbon and nitrogen supply, restricted carbon supply, and restricted nitrogen supply and analyzed for biomass and 54 metabolic traits. Correlation-based metabolic networks were generated from the genotype-dependent variation in each condition to reveal sets of metabolites that show coordinated changes across accessions. The networks were largely specific for a single growth condition. Partial least squares regression from metabolic traits allowed prediction of biomass within and, slightly more weakly, across conditions (cross-validated Pearson correlations in the range of 0.27-0.58 and 0.21-0.51 and P values in the range of ,0.001-,0.13 and ,0.001-,0.023, respectively). Metabolic traits that correlate with growth or have a high weighting in the partial least squares regression were mainly condition specific and often related to the resource that restricts growth under that condition. Linear mixed-model analysis using the combined metabolic traits from all growth conditions as an input indicated that inclusion of random effects for the conditions improves predictions of biomass. Thus, robust prediction of biomass across a range of conditions requires condition-specific measurement of metabolic traits to take account of environment-dependent changes of the underlying networks.

Increasing amounts of sequence data are becoming available for a wide range of non-model organisms. Investigating and modelling the metabolic behaviour of those organisms is highly relevant to understand their biology and ecology. As sequences are often incomplete and poorly annotated, draft networks of their metabolism largely suffer from incompleteness. Appropriate gap-filling methods to identify and add missing reactions are therefore required to address this issue. However, current tools rely on phenotypic or taxonomic information, or are very sensitive to the stoichiometric balance of metabolic reactions, especially concerning the co-factors. This type of information is often not available or at least prone to errors for newly-explored organisms. Here we introduce , a tool dedicated to the topological gap-filling of genome-scale draft metabolic networks. reformulates gap-filling as a qualitative combinatorial optimization problem, omitting constraints raised by the stoichiometry of a metabolic network considered in other methods, and solves this problem using Answer Set Programming. Run on several artificial test sets gathering 10,800 degraded Escherichia coli networks was able to efficiently identify essential reactions missing in networks at high degradation rates, outperforming the stoichiometry-based tools in scalability. To demonstrate the utility of we applied it to two case studies. Its application to recent metabolic networks reconstructed for the brown algal model Ectocarpus siliculosus and an associated bacterium Candidatus Phaeomarinobacter ectocarpi revealed several candidate metabolic pathways for algal-bacterial interactions. Then was used to reconstruct, from transcriptomic and metabolomic data, the first metabolic network for the microalga Euglena mutabilis. These two case studies show that is a versatile tool to complete draft genome-scale metabolic networks produced from heterogeneous data, and to suggest relevant reactions that explain the metabolic capacity of a biological system.

Understanding the control of large-scale metabolic networks is central to biology and medicine. However, existing approaches either require specifying a cellular objective or can only be used for small networks. We introduce new coupling types describing the relations between reaction activities, and develop an efficient computational framework, which does not require any cellular objective for systematic studies of large-scale metabolism. We identify the driver reactions facilitating control of 23 metabolic networks from all kingdoms of life. We find that unicellular organisms require a smaller degree of control than multicellular organisms. Driver reactions are under complex cellular regulation in Escherichia coli, indicating their preeminent role in facilitating cellular control. In human cancer cells, driver reactions play pivotal roles in malignancy and represent potential therapeutic targets. The developed framework helps us gain insights into regulatory principles of diseases and facilitates design of engineering strategies at the interface of gene regulation, signaling, and metabolism.[Supplemental material is available for this article.]Understanding how cellular systems are controlled on a genomescale is a central issue in biology and medicine. Metabolic networks are at the center of systems biology approaches to unraveling cellular control, because metabolism carries the life-sustaining cellular functions shaping the molecular phenotype (Sweetlove and Ratcliffe 2011). The steady-state principle and physico-chemical constraints (e.g., mass balance and thermodynamics) have been employed to reduce the number of considered network states, facilitating the prediction of genotype-phenotype relationships and intervention strategies for biotechnological or medical purposes (McCloskey et al. 2013). In particular, flux balance analysis and variations thereof have been successfully applied to the metabolic networks of unicellular organisms to predict their metabolic and cellular phenotypes (Varma and Palsson 1994). Yet, those approaches are biased (Lewis et al. 2012) because they restrict the flux space to an a priori specified reference state by assuming a cellular objective to be optimized by the organism (Schuetz et al. 2007). While optimization of biomass yield has proven useful for unicellular organisms, identification of a suitable objective for multicellular organisms remains a nontrivial endeavor (Sweetlove and Ratcliffe 2011). Other approaches, e.g., elementary flux modes (Schuster and Schuster 1993) and extreme pathways analyses (Schilling et al. 2000), do not assume a cellular objective and hence are unbiased. However, despite extensive studies and recent advances (Terzer and Stelling 2008), these unbiased approaches are limited to rather small networks due to their intrinsic computational complexity. We still lack an unbiased computational approach for systematically studying the control of large-scale metabolic networks.Here, we develop such an approach by employing the flux coupling between reactio...

BackgroundFlux coupling analysis (FCA) has become a useful tool in the constraint-based analysis of genome-scale metabolic networks. FCA allows detecting dependencies between reaction fluxes of metabolic networks at steady-state. On the one hand, this can help in the curation of reconstructed metabolic networks by verifying whether the coupling between reactions is in agreement with the experimental findings. On the other hand, FCA can aid in defining intervention strategies to knock out target reactions.ResultsWe present a new method F2C2 for FCA, which is orders of magnitude faster than previous approaches. As a consequence, FCA of genome-scale metabolic networks can now be performed in a routine manner.ConclusionsWe propose F2C2 as a fast tool for the computation of flux coupling in genome-scale metabolic networks. F2C2 is freely available for non-commercial use at https://sourceforge.net/projects/f2c2/files/.

No abstract

Interplay within microbial communities impacts ecosystems on several scales, and elucidation of the consequent effects is a difficult task in ecology. In particular, the integration of genome-scale data within quantitative models of microbial ecosystems remains elusive. This study advocates the use of constraint-based modeling to build predictive models from recent high-resolution -omics datasets. Following recent studies that have demonstrated the accuracy of constraint-based models (CBMs) for simulating single-strain metabolic networks, we sought to study microbial ecosystems as a combination of single-strain metabolic networks that exchange nutrients. This study presents two multi-objective extensions of CBMs for modeling communities: multi-objective flux balance analysis (MO-FBA) and multi-objective flux variability analysis (MO-FVA). Both methods were applied to a hot spring mat model ecosystem. As a result, multiple trade-offs between nutrients and growth rates, as well as thermodynamically favorable relative abundances at community level, were emphasized. We expect this approach to be used for integrating genomic information in microbial ecosystems. Following models will provide insights about behaviors (including diversity) that take place at the ecosystem scale.

Highlights d Seipin is enriched at ER-MAMs d Seipin interacts with MAM calcium regulators in a nutritionally regulated manner d Adipocyte seipin deficiency impairs mitochondrial calcium import and ATP production d Inducible seipin removal from adipose tissue leads to rapid mitochondrial dysfunction