Understanding dynamics using sensitivity analysis: caveat and solution

<p>Abstract</p> <p>Background</p> <p>Parametric sensitivity analysis (PSA) has become one of the most commonly used tools in computational systems biology, in which the sensitivity coefficients are used to study the parametric dependence of biological models. As many of these models describe dynamical behaviour of biological systems, the PSA has subsequently been used to elucidate important cellular processes that regulate this dynamics. However, in this paper, we show that the PSA coefficients are not suitable in inferring the mechanisms by which dynamical behaviour arises and in fact it can even lead to incorrect conclusions.</p> <p>Results</p> <p>A careful interpretation of parametric perturbations used in the PSA is presented here to explain the issue of using this analysis in inferring dynamics. In short, the PSA coefficients quantify the integrated change in the system behaviour due to persistent parametric perturbations, and thus the dynamical information of when a parameter perturbation matters is lost. To get around this issue, we present a new sensitivity analysis based on impulse perturbations on system parameters, which is named impulse parametric sensitivity analysis (iPSA). The inability of PSA and the efficacy of iPSA in revealing mechanistic information of a dynamical system are illustrated using two examples involving switch activation.</p> <p>Conclusions</p> <p>The interpretation of the PSA coefficients of dynamical systems should take into account the persistent nature of parametric perturbations involved in the derivation of this analysis. The application of PSA to identify the controlling mechanism of dynamical behaviour can be misleading. By using impulse perturbations, introduced at different times, the iPSA provides the necessary information to understand how dynamics is achieved, i.e. which parameters are essential and when they become important.</p

A role for the cell-wall protein silacidin in cell size of the diatom Thalassiosira pseudonana

Diatoms contribute 20% of global primary production and form the basis of many marine food webs. Although their species diversity correlates with broad diversity in cell size, there is also an intraspecific cell-size plasticity due to sexual reproduction and varying environmental conditions. However, despite the ecological significance of the diatom cell size for food-web structure and global biogeochemical cycles, our knowledge about genes underpinning the size of diatom cells remains elusive. Here, a combination of reverse genetics, experimental evolution and comparative RNA8 sequencing analyses enabled us to identify a previously unknown genetic control of cell size in the diatom Thalassiosira pseudonana. In particular, the targeted deregulation of the expression of the cell-wall protein silacidin caused a significant increase in valve diameter. Remarkably, the natural downregulation of the silacidin gene transcript due to experimental evolution under low temperature also correlated with cell-size increase. Our data give first evidence for a genetically controlled regulation of cell size in Thalassiosira pseudonana and possibly other centric diatoms as they also encode the silacidin gene in their genomes

Spontaneous Reaction Silencing in Metabolic Optimization

Metabolic reactions of single-cell organisms are routinely observed to become dispensable or even incapable of carrying activity under certain circumstances. Yet, the mechanisms as well as the range of conditions and phenotypes associated with this behavior remain very poorly understood. Here we predict computationally and analytically that any organism evolving to maximize growth rate, ATP production, or any other linear function of metabolic fluxes tends to significantly reduce the number of active metabolic reactions compared to typical non-optimal states. The reduced number appears to be constant across the microbial species studied and just slightly larger than the minimum number required for the organism to grow at all. We show that this massive spontaneous reaction silencing is triggered by the irreversibility of a large fraction of the metabolic reactions and propagates through the network as a cascade of inactivity. Our results help explain existing experimental data on intracellular flux measurements and the usage of latent pathways, shedding new light on microbial evolution, robustness, and versatility for the execution of specific biochemical tasks. In particular, the identification of optimal reaction activity provides rigorous ground for an intriguing knockout-based method recently proposed for the synthetic recovery of metabolic function.Comment: 34 pages, 6 figure

Metagenomics of the Deep Mediterranean, a Warm Bathypelagic Habitat

BACKGROUND: Metagenomics is emerging as a powerful method to study the function and physiology of the unexplored microbial biosphere, and is causing us to re-evaluate basic precepts of microbial ecology and evolution. Most marine metagenomic analyses have been nearly exclusively devoted to photic waters. METHODOLOGY/PRINCIPAL FINDINGS: We constructed a metagenomic fosmid library from 3,000 m-deep Mediterranean plankton, which is much warmer (approximately 14 degrees C) than waters of similar depth in open oceans (approximately 2 degrees C). We analyzed the library both by phylogenetic screening based on 16S rRNA gene amplification from clone pools and by sequencing both insert extremities of ca. 5,000 fosmids. Genome recruitment strategies showed that the majority of high scoring pairs corresponded to genomes from Rhizobiales within the Alphaproteobacteria, Cenarchaeum symbiosum, Planctomycetes, Acidobacteria, Chloroflexi and Gammaproteobacteria. We have found a community structure similar to that found in the aphotic zone of the Pacific. However, the similarities were significantly higher to the mesopelagic (500-700 m deep) in the Pacific than to the single 4000 m deep sample studied at this location. Metabolic genes were mostly related to catabolism, transport and degradation of complex organic molecules, in agreement with a prevalent heterotrophic lifestyle for deep-sea microbes. However, we observed a high percentage of genes encoding dehydrogenases and, among them, cox genes, suggesting that aerobic carbon monoxide oxidation may be important in the deep ocean as an additional energy source. CONCLUSIONS/SIGNIFICANCE: The comparison of metagenomic libraries from the deep Mediterranean and the Pacific ALOHA water column showed that bathypelagic Mediterranean communities resemble more mesopelagic communities in the Pacific, and suggests that, in the absence of light, temperature is a major stratifying factor in the oceanic water column, overriding pressure at least over 4000 m deep. Several chemolithotrophic metabolic pathways could supplement organic matter degradation in this most depleted habitat

The changing carbon cycle of the coastal ocean

The carbon cycle of the coastal ocean is a dynamic component of the global carbon budget. But the diverse sources and sinks of carbon and their complex interactions in these waters remain poorly understood. Here we discuss the sources, exchanges and fates of carbon in the coastal ocean and how anthropogenic activities have altered the carbon cycle. Recent evidence suggests that the coastal ocean may have become a net sink for atmospheric carbon dioxide during post-industrial times. Continued human pressures in coastal zones will probably have an important impact on the future evolution of the coastal ocean's carbon budget

Preservation of organic matter in mound-forming coral skeletons

This study demonstrates that intracrystalline organic matter in coral skeletons is well preserved over century timescales. The extent of preservation of organic matter in coral skeletons was investigated by measuring total organic carbon (TOC), total hydrolyzable amino acid (THAA), chloropigment, and lipid concentrations in 0-300 year old annual growth bands from Montastraea annularis (Florida Keys) and Porites lutea (Red Sea). Organic matter intrinsic to the calcium carbonate mineral (intracrystalline) was analyzed separately from total skeletal organic matter. The Red Sea coral had less TOC (0.02-0.04 wt%) than the Florida Keys coral (0.04-0.11 wt%), but a higher percent of intracrystalline organic matter in all annual bands measured. Carbon in the form of THAA, most likely from mineral-precipitating proteins, contributed 30-45% of the TOC in both corals. Carbon in lipids represented about 3% of the TOC in the coral skeletons. Chlorophyll-a and b were present in annual bands where endolithic algae were present, but these compounds were minor contributors to TOC. The distribution of specific organic compounds showed that organic matter was well preserved throughout the time period sampled in both the total and intracrystalline pools. Variations in THAA were not correlated with TOC over time, suggesting that organic matter that is involved in biomineralization, like amino acids, may be deposited in response to different environmental factors than are other components of skeletal organic matter. Differences in the quantity and composition of organic matter between the two corals investigated here were assessed using principal components analysis and suggest that location, species and skeletal structure may all influence organic matter content and possibly the degree of physical protection of organic matter by the coral skeleton. Further, our study suggests that intracrystalline organic matter may be better protected from diagenesis than non-intracrystalline organic matter and may therefore be a more reliable source of organic matter for paleoceanographic studies than total skeletal organic matter. © 2003 Elsevier Science Ltd

Preservation of organic matter in mound-forming coral skeletons

This study demonstrates that intracrystalline organic matter in coral skeletons is well preserved over century timescales. The extent of preservation of organic matter in coral skeletons was investigated by measuring total organic carbon (TOC), total hydrolyzable amino acid (THAA), chloropigment, and lipid concentrations in 0-300 year old annual growth bands from Montastraea annularis (Florida Keys) and Porites lutea (Red Sea). Organic matter intrinsic to the calcium carbonate mineral (intracrystalline) was analyzed separately from total skeletal organic matter. The Red Sea coral had less TOC (0.02-0.04 wt%) than the Florida Keys coral (0.04-0.11 wt%), but a higher percent of intracrystalline organic matter in all annual bands measured. Carbon in the form of THAA, most likely from mineral-precipitating proteins, contributed 30-45% of the TOC in both corals. Carbon in lipids represented about 3% of the TOC in the coral skeletons. Chlorophyll-a and b were present in annual bands where endolithic algae were present, but these compounds were minor contributors to TOC. The distribution of specific organic compounds showed that organic matter was well preserved throughout the time period sampled in both the total and intracrystalline pools. Variations in THAA were not correlated with TOC over time, suggesting that organic matter that is involved in biomineralization, like amino acids, may be deposited in response to different environmental factors than are other components of skeletal organic matter. Differences in the quantity and composition of organic matter between the two corals investigated here were assessed using principal components analysis and suggest that location, species and skeletal structure may all influence organic matter content and possibly the degree of physical protection of organic matter by the coral skeleton. Further, our study suggests that intracrystalline organic matter may be better protected from diagenesis than non-intracrystalline organic matter and may therefore be a more reliable source of organic matter for paleoceanographic studies than total skeletal organic matter. © 2003 Elsevier Science Ltd