Systemic properties of metabolic networks lead to an epistasis-based model for heterosis

The genetic and molecular approaches to heterosis usually do not rely on any model of the genotype–phenotype relationship. From the generalization of Kacser and Burns’ biochemical model for dominance and epistasis to networks with several variable enzymes, we hypothesized that metabolic heterosis could be observed because the response of the flux towards enzyme activities and/or concentrations follows a multi-dimensional hyperbolic-like relationship. To corroborate this, we used the values of systemic parameters accounting for the kinetic behaviour of four enzymes of the upstream part of glycolysis, and simulated genetic variability by varying in silico enzyme concentrations. Then we “crossed” virtual parents to get 1,000 hybrids, and showed that best-parent heterosis was frequently observed. The decomposition of the flux value into genetic effects, with the help of a novel multilocus epistasis index, revealed that antagonistic additive-by-additive epistasis effects play the major role in this framework of the genotype–phenotype relationship. This result is consistent with various observations in quantitative and evolutionary genetics, and provides a model unifying the genetic effects underlying heterosis

Convergent genetic and expression data implicate immunity in Alzheimer's disease

Background Late–onset Alzheimer's disease (AD) is heritable with 20 genes showing genome wide association in the International Genomics of Alzheimer's Project (IGAP). To identify the biology underlying the disease we extended these genetic data in a pathway analysis. Methods The ALIGATOR and GSEA algorithms were used in the IGAP data to identify associated functional pathways and correlated gene expression networks in human brain. Results ALIGATOR identified an excess of curated biological pathways showing enrichment of association. Enriched areas of biology included the immune response (p = 3.27×10-12 after multiple testing correction for pathways), regulation of endocytosis (p = 1.31×10-11), cholesterol transport (p = 2.96 × 10-9) and proteasome-ubiquitin activity (p = 1.34×10-6). Correlated gene expression analysis identified four significant network modules, all related to the immune response (corrected p 0.002 – 0.05). Conclusions The immune response, regulation of endocytosis, cholesterol transport and protein ubiquitination represent prime targets for AD therapeutics

The Pitfalls of Heterosis Coefficients

Heterosis (hybrid vigour) is a universal phenomenon of crucial agro-economic and evolutionary importance. We show that the most common heterosis coefficients do not properly measure deviation from additivity because they include both a component accounting for “real” heterosis and a term that is not related to heterosis, since it is derived solely from parental values. Therefore, these coefficients are inadequate whenever the aim of the study is to compare heterosis levels between different traits, environments, genetic backgrounds, or developmental stages, as these factors may affect not only the level of non-additivity, but also parental values. The only relevant coefficient for such comparisons is the so-called “potence ratio”. Because most heterosis studies consider several traits/stages/environmental conditions, our observations support the use of the potence ratio, at least in non-agronomic contexts, because it is the only non-ambiguous heterosis coefficient

Data_Sheet_1_Heterosis Is a Systemic Property Emerging From Non-linear Genotype-Phenotype Relationships: Evidence From in Vitro Genetics and Computer Simulations.zip

<p>Heterosis, the superiority of hybrids over their parents for quantitative traits, represents a crucial issue in plant and animal breeding as well as evolutionary biology. Heterosis has given rise to countless genetic, genomic and molecular studies, but has rarely been investigated from the point of view of systems biology. We hypothesized that heterosis is an emergent property of living systems resulting from frequent concave relationships between genotypic variables and phenotypes, or between different phenotypic levels. We chose the enzyme-flux relationship as a model of the concave genotype-phenotype (GP) relationship, and showed that heterosis can be easily created in the laboratory. First, we reconstituted in vitro the upper part of glycolysis. We simulated genetic variability of enzyme activity by varying enzyme concentrations in test tubes. Mixing the content of “parental” tubes resulted in “hybrids,” whose fluxes were compared to the parental fluxes. Frequent heterotic fluxes were observed, under conditions that were determined analytically and confirmed by computer simulation. Second, to test this model in a more realistic situation, we modeled the glycolysis/fermentation network in yeast by considering one input flux, glucose, and two output fluxes, glycerol and acetaldehyde. We simulated genetic variability by randomly drawing parental enzyme concentrations under various conditions, and computed the parental and hybrid fluxes using a system of differential equations. Again we found that a majority of hybrids exhibited positive heterosis for metabolic fluxes. Cases of negative heterosis were due to local convexity between certain enzyme concentrations and fluxes. In both approaches, heterosis was maximized when the parents were phenotypically close and when the distributions of parental enzyme concentrations were contrasted and constrained. These conclusions are not restricted to metabolic systems: they only depend on the concavity of the GP relationship, which is commonly observed at various levels of the phenotypic hierarchy, and could account for the pervasiveness of heterosis.</p

Contribution des Insertion/Délétions à la variations des caractères, l’adaptation et aux performances hybrides au moyen d’une puce de génotypage haut débit

National audienc

Uranium perturbs signaling and iron uptake response in Arabidopsis thaliana roots.

International audienceUranium is a natural element which is mainly redistributed in the environment due to human activity, including accidents and spillages. Plants may be useful in cleaning up after incidents, although little is yet known about the relationship between metal speciation and plant response. Here, J-Chess modeling was used to predict U speciation and exposure conditions affecting U bioavailability for plants. The model was confirmed by exposing Arabidopsis thaliana plants to U under hydroponic conditions. The early root response was characterized using complete Arabidopsis transcriptome microarrays (CATMA). Expression of 111 genes was modified at the three timepoints studied. The associated biological processes were further examined by real-time quantitative RT-PCR. Annotation revealed that oxidative stress, cell wall and hormone biosynthesis, and signaling pathways (including phosphate signaling) were affected by U exposure. The main actors in iron uptake and signaling (IRT1, FRO2, AHA2, AHA7 and FIT1) were strongly down-regulated upon exposure to uranyl. A network calculated using IRT1, FRO2 and FIT1 as bait revealed a set of genes whose expression levels change under U stress. Hypotheses are presented to explain how U perturbs the iron uptake and signaling response. These results give preliminary insights into the pathways affected by uranyl uptake, which will be of interest for engineering plants to help clean areas contaminated with U

Simplified modelling of metabolic pathways for flux prediction and optimization: lessons from an in vitro reconstruction of the upper part of glycolysis

Explicit modelling of metabolic networks relies on well-known mathematical tools and specialized computer programs. However, identifying and estimating the values of the very numerous enzyme parameters inherent to the models remain a tedious and difficult task, and the rate equations of the reactions are usually not known in sufficient detail. A way to circumvent this problem is to use ‘non-mechanistic’ models, which may account for the behaviour of the systems with a limited number of parameters. Working on the first part of glycolysis reconstituted in vitro, we showed how to derive, from titration experiments, values of effective enzyme activity parameters that do not include explicitly any of the classical kinetic constants. With a maximum of only two parameters per enzyme, this approach produced very good estimates for the flux values, and enabled us to determine the optimization conditions of the system, i.e. to calculate the set of enzyme concentrations that maximizes the flux. This fast and easy method should be valuable in the context of integrative biology or for metabolic engineering, where the challenge is to deal with the dramatic increase in the number of parameters when the systems become complex

Convergent genetic and expression data implicate immunity in Alzheimer's disease

International audienceBackground Late‐onset Alzheimer's disease (AD) is heritable with 20 genes showing genome‐wide association in the International Genomics of Alzheimer's Project (IGAP). To identify the biology underlying the disease, we extended these genetic data in a pathway analysis. Methods The ALIGATOR and GSEA algorithms were used in the IGAP data to identify associated functional pathways and correlated gene expression networks in human brain. Results ALIGATOR identified an excess of curated biological pathways showing enrichment of association. Enriched areas of biology included the immune response ( P = 3.27 × 10 −12 after multiple testing correction for pathways), regulation of endocytosis ( P = 1.31 × 10 −11 ), cholesterol transport ( P = 2.96 × 10 −9 ), and proteasome‐ubiquitin activity ( P = 1.34 × 10 −6 ). Correlated gene expression analysis identified four significant network modules, all related to the immune response (corrected P = .002–.05). Conclusions The immune response, regulation of endocytosis, cholesterol transport, and protein ubiquitination represent prime targets for AD therapeutics