Multiscale Metabolic Modeling: Dynamic Flux Balance Analysis on a Whole-Plant Scale

Plant metabolism is characterized by a unique complexity on the cellular, tissue, and organ levels. On a whole-plant scale, changing source and sink relations accompanying plant development add another level of complexity to metabolism. With the aim of achieving a spatiotemporal resolution of source-sink interactions in crop plant metabolism, a multiscale metabolic modeling (MMM) approach was applied that integrates static organ-specific models with a whole-plant dynamic model. Allowing for a dynamic flux balance analysis on a whole-plant scale, the MMM approach was used to decipher the metabolic behavior of source and sink organs during the generative phase of the barley (Hordeum vulgare) plant. It reveals a sink-to-source shift of the barley stem caused by the senescence-related decrease in leaf source capacity, which is not sufficient to meet the nutrient requirements of sink organs such as the growing seed. The MMM platform represents a novel approach for the in silico analysis of metabolism on a whole-plant level, allowing for a systemic, spatiotemporally resolved understanding of metabolic processes involved in carbon partitioning, thus providing a novel tool for studying yield stability and crop improvement.Plants are of vital significance as a source of food (Grusak and DellaPenna, 1999; Rogalski and Carrer, 2011), feed (Lu et al., 2011), energy (Tilman et al., 2006; Parmar et al., 2011), and feedstocks for the chemical industry (Metzger and Bornscheuer, 2006; Kinghorn et al., 2011). Given the close connection between plant metabolism and the usability of plant products, there is a growing interest in understanding and predicting the behavior and regulation of plant metabolic processes. In order to increase crop quality and yield, there is a need for methods guiding the rational redesign of the plant metabolic network (Schwender, 2009).Mathematical modeling of plant metabolism offers new approaches to understand, predict, and modify complex plant metabolic processes. In plant research, the issue of metabolic modeling is constantly gaining attention, and different modeling approaches applied to plant metabolism exist, ranging from highly detailed quantitative to less complex qualitative approaches (for review, see Giersch, 2000; Morgan and Rhodes, 2002; Poolman et al., 2004; Rios-Estepa and Lange, 2007).A widely used modeling approach is flux balance analysis (FBA), which allows the prediction of metabolic capabilities and steady-state fluxes under different environmental and genetic backgrounds using (non)linear optimization (Orth et al., 2010). Assuming steady-state conditions, FBA has the advantage of not requiring the knowledge of kinetic parameters and, therefore, can be applied to model detailed, large-scale systems. In recent years, the FBA approach has been applied to several different plant species, such as maize (Zea mays; Dal’Molin et al., 2010; Saha et al., 2011), barley (Hordeum vulgare; Grafahrend-Belau et al., 2009b; Melkus et al., 2011; Rolletschek et al., 2011), rice (Oryza sativa; Lakshmanan et al., 2013), Arabidopsis (Arabidopsis thaliana; Poolman et al., 2009; de Oliveira Dal’Molin et al., 2010; Radrich et al., 2010; Williams et al., 2010; Mintz-Oron et al., 2012; Cheung et al., 2013), and rapeseed (Brassica napus; Hay and Schwender, 2011a, 2011b; Pilalis et al., 2011), as well as algae (Boyle and Morgan, 2009; Cogne et al., 2011; Dal’Molin et al., 2011) and photoautotrophic bacteria (Knoop et al., 2010; Montagud et al., 2010; Boyle and Morgan, 2011). These models have been used to study different aspects of metabolism, including the prediction of optimal metabolic yields and energy efficiencies (Dal’Molin et al., 2010; Boyle and Morgan, 2011), changes in flux under different environmental and genetic backgrounds (Grafahrend-Belau et al., 2009b; Dal’Molin et al., 2010; Melkus et al., 2011), and nonintuitive metabolic pathways that merit subsequent experimental investigations (Poolman et al., 2009; Knoop et al., 2010; Rolletschek et al., 2011). Although FBA of plant metabolic models was shown to be capable of reproducing experimentally determined flux distributions (Williams et al., 2010; Hay and Schwender, 2011b) and generating new insights into metabolic behavior, capacities, and efficiencies (Sweetlove and Ratcliffe, 2011), challenges remain to advance the utility and predictive power of the models.Given that many plant metabolic functions are based on interactions between different subcellular compartments, cell types, tissues, and organs, the reconstruction of organ-specific models and the integration of these models into interacting multiorgan and/or whole-plant models is a prerequisite to get insight into complex plant metabolic processes organized on a whole-plant scale (e.g. source-sink interactions). Almost all FBA models of plant metabolism are restricted to one cell type (Boyle and Morgan, 2009; Knoop et al., 2010; Montagud et al., 2010; Cogne et al., 2011; Dal’Molin et al., 2011), one tissue or one organ (Grafahrend-Belau et al., 2009b; Hay and Schwender, 2011a, 2011b; Pilalis et al., 2011; Mintz-Oron et al., 2012), and only one model exists taking into account the interaction between two cell types by specifying the interaction between mesophyll and bundle sheath cells in C4 photosynthesis (Dal’Molin et al., 2010). So far, no model representing metabolism at the whole-plant scale exists.Considering whole-plant metabolism raises the problem of taking into account temporal and environmental changes in metabolism during plant development and growth. Although classical static FBA is unable to predict the dynamics of metabolic processes, as the network analysis is based on steady-state solutions, time-dependent processes can be taken into account by extending the classical static FBA to a dynamic flux balance analysis (dFBA), as proposed by Mahadevan et al. (2002). The static (SOA) and dynamic optimization approaches introduced in this work provide a framework for analyzing the transience of metabolism by integrating kinetic expressions to dynamically constrain exchange fluxes. Due to the requirement of knowing or estimating a large number of kinetic parameters, so far dFBA has only been applied to a plant metabolic model once, to study the photosynthetic metabolism in the chloroplasts of C3 plants by a simplified model of five biochemical reactions (Luo et al., 2009). Integrating a dynamic model into a static FBA model is an alternative approach to perform dFBA.In this study, a multiscale metabolic modeling (MMM) approach was applied with the aim of achieving a spatiotemporal resolution of cereal crop plant metabolism. To provide a framework for the in silico analysis of the metabolic dynamics of barley on a whole-plant scale, the MMM approach integrates a static multiorgan FBA model and a dynamic whole-plant multiscale functional plant model (FPM) to perform dFBA. The performance of the novel whole-plant MMM approach was tested by studying source-sink interactions during the seed developmental phase of barley plants.     

The plant Selaginella moellendorffii possesses enzymes for synthesis and hydrolysis of the compatible solutes mannosylglycerate and glucosylglycerate

A mannosylglycerate synthase (MgS) gene detected in the genome of Selaginella moellendorffii was expressed in E. coli and the recombinant enzyme was purified and characterized. A remarkable and unprecedented feature of this enzyme was the ability to efficiently synthesize mannosylglycerate (MG) and glucosylglycerate (GG) alike, with maximal activity at 50 °C, pH 8.0 and with Mg²⁺ as reaction enhancer. We have also identified a novel glycoside hydrolase gene in this plant’s genome, which was functionally confirmed to be highly specific for the hydrolysis of MG and GG and named MG hydrolase (MgH), due to its homology with bacterial MgHs. The recombinant enzyme was maximally active at 40 °C and at pH 6.0–6.5. The activity was independent of cations, but Mn²⁺ was a strong stimulator. Regardless of these efficient enzymatic resources we could not detect MG or GG in S. moellendorffii or in the extracts of five additional Selaginella species. Herein, we describe the properties of the first eukaryotic enzymes for the synthesis and hydrolysis of the compatible solutes, MG and GG.     

Attenuated and Replication-Competent Vaccinia Virus Strains M65 and M101 with Distinct Biology and Immunogenicity as Potential Vaccine Candidates against Pathogens

Replication-competent poxvirus vectors with an attenuation phenotype and with a high immunogenic capacity of the foreign expressed antigen are being pursued as novel vaccine vectors against different pathogens. In this investigation, we have examined the replication and immunogenic characteristics of two vaccinia virus (VACV) mutants, M65 and M101. These mutants were generated after 65 and 101 serial passages of persistently infected Friend erythroleukemia (FEL) cells. In cultured cells of different origins, the mutants are replication competent and have growth kinetics similar to or slightly reduced in comparison with those of the parental Western Reserve (WR) virus strain. In normal and immune-suppressed infected mice, the mutants showed different levels of attenuation and pathogenicity in comparison with WR and modified vaccinia Ankara (MVA) strains. Wide genome analysis after deep sequencing revealed selected genomic deletions and mutations in a number of viral open reading frames (ORFs). Mice immunized in a DNA prime/mutant boost regimen with viral vectors expressing the LACK (Leishmania homologue for receptors of activated C kinase) antigen of Leishmania infantum showed protection or a delay in the onset of cutaneous leishmaniasis. Protection was similar to that triggered by MVA-LACK. In immunized mice, both polyfunctional CD4⁺ and CD8⁺ T cells with an effector memory phenotype were activated by the two mutants, but the DNA-LACK/M65-LACK protocol preferentially induced CD4⁺ whereas DNA-LACK/M101-LACK preferentially induced CD8⁺ T cell responses. Altogether, our findings showed the adaptive changes of the WR genome during long-term virus-host cell interaction and how the replication competency of M65 and M101 mutants confers distinct biological properties and immunogenicity in mice compared to those of the MVA strain. These mutants could have applicability for understanding VACV biology and as potential vaccine vectors against pathogens and tumors.     

A comparison of the strength of biodiversity effects across multiple functions

In order to predict which ecosystem functions are most at risk from biodiversity loss, meta-analyses have generalised results from biodiversity experiments over different sites and ecosystem types. In contrast, comparing the strength of biodiversity effects across a large number of ecosystem processes measured in a single experiment permits more direct comparisons. Here, we present an analysis of 418 separate measures of 38 ecosystem processes. Overall, 45 % of processes were significantly affected by plant species richness, suggesting that, while diversity affects a large number of processes not all respond to biodiversity. We therefore compared the strength of plant diversity effects between different categories of ecosystem processes, grouping processes according to the year of measurement, their biogeochemical cycle, trophic level and compartment (above- or belowground) and according to whether they were measures of biodiversity or other ecosystem processes, biotic or abiotic and static or dynamic. Overall, and for several individual processes, we found that biodiversity effects became stronger over time. Measures of the carbon cycle were also affected more strongly by plant species richness than were the measures associated with the nitrogen cycle. Further, we found greater plant species richness effects on measures of biodiversity than on other processes. The differential effects of plant diversity on the various types of ecosystem processes indicate that future research and political effort should shift from a general debate about whether biodiversity loss impairs ecosystem functions to focussing on the specific functions of interest and ways to preserve them individually or in combination.     

Foramen magnum position in bipedal mammals

The anterior position of the human foramen magnum is often explained as an adaptation for maintaining balance of the head atop the cervical vertebral column during bipedalism and the assumption of orthograde trunk postures. Accordingly, the relative placement of the foramen magnum on the basicranium has been used to infer bipedal locomotion and hominin status for a number of Mio-Pliocene fossil taxa. Nonetheless, previous studies have struggled to validate the functional link between foramen magnum position and bipedal locomotion. Here, we test the hypothesis that an anteriorly positioned foramen magnum is related to bipedalism through a comparison of basicranial anatomy between bipeds and quadrupeds from three mammalian clades: marsupials, rodents and primates. Additionally, we examine whether strepsirrhine primates that habitually assume orthograde trunk postures exhibit more anteriorly positioned foramina magna compared with non-orthograde strepsirrhines. Our comparative data reveal that bipedal marsupials and rodents have foramina magna that are more anteriorly located than those of quadrupedal close relatives. The foramen magnum is also situated more anteriorly in orthograde strepsirrhines than in pronograde or antipronograde strepsirrhines. Among the primates sampled, humans exhibit the most anteriorly positioned foramina magna. The results of this analysis support the utility of foramen magnum position as an indicator of bipedal locomotion in fossil hominins.     

Enhanced enzymatic cellulose degradation by cellobiohydrolases via product removal

Product inhibition by cellobiose decreases the rate of enzymatic cellulose degradation. The optimal reaction conditions for two Emericella (Aspergillus) nidulans-derived cellobiohydrolases I and II produced in Pichia pastoris were identified as CBHI: 52 °C, pH 4.5–6.5, and CBHII: 46 °C, pH 4.8. The optimum in a mixture of the two was 50 °C, pH 4.9. An almost fourfold increase in enzymatic hydrolysis yield was achieved with intermittent product removal of cellobiose with membrane filtration (2 kDa cut-off): The conversion of cotton cellulose after 72 h was ~19 % by weight, whereas the conversion in the parallel batch reaction was only ~5 % by weight. Also, a synergistic effect, achieving ~27 % substrate conversion, was obtained by addition of endo-1,4-β-d-glucanase. The synergistic effect was only obtained with product removal. By using pure, monoactive enzymes, the work illustrates the profound gains achievable by intermittent product removal during cellulose hydrolysis.     

Intraracial harassment on campus: explaining between- and within-group differences

Using the National Longitudinal Survey of Freshmen (NLSF), we examine both between- and within-group differences in the odds of feeling intraracially harassed. Specifically, we investigate the effects of colleges’ and universities’ racial composition as well as the nature of students’ associations with non-group members, including involvement in racially homogeneous campus organizations, ethnoracial diversity of friendship networks, and interracial dating. Our findings suggest that although college racial composition appears to have little effect on experiencing intraracial harassment, the nature of students’ involvement with other-race students matters a great deal. For all groups, interracial dating increased odds of harassment. Among black and white students, more diverse friendship networks did as well. And among Asian and Latino students, involvement in any racially homogeneous campus organization was associated with increases in reports of intraracial harassment. Thus, we propose a baseline theoretical model of intraracial harassment that highlights the nature of students’ associations with outgroups.     

Rapid Sample Processing for Detection of Food-Borne Pathogens via Cross-Flow Microfiltration

This paper reports an approach to enable rapid concentration and recovery of bacterial cells from aqueous chicken homogenates as a preanalytical step of detection. This approach includes biochemical pretreatment and prefiltration of food samples and development of an automated cell concentration instrument based on cross-flow microfiltration. A polysulfone hollow-fiber membrane module having a nominal pore size of 0.2 μm constitutes the core of the cell concentration instrument. The aqueous chicken homogenate samples were circulated within the cross-flow system achieving 500- to 1,000-fold concentration of inoculated Salmonella enterica serovar Enteritidis and naturally occurring microbiota with 70% recovery of viable cells as determined by plate counting and quantitative PCR (qPCR) within 35 to 45 min. These steps enabled 10 CFU/ml microorganisms in chicken homogenates or 10² CFU/g chicken to be quantified. Cleaning and sterilizing the instrument and membrane module by stepwise hydraulic and chemical cleaning (sodium hydroxide and ethanol) enabled reuse of the membrane 15 times before replacement. This approach begins to address the critical need for the food industry for detecting food pathogens within 6 h or less.