Propagation of genetic variation in gene regulatory networks.

A future quantitative genetics theory should link genetic variation to phenotypic variation in a causally cohesive way based on how genes actually work and interact. We provide a theoretical framework for predicting and understanding the manifestation of genetic variation in haploid and diploid regulatory networks with arbitrary feedback structures and intra-locus and inter-locus functional dependencies. Using results from network and graph theory, we define propagation functions describing how genetic variation in a locus is propagated through the network, and show how their derivatives are related to the network's feedback structure. Similarly, feedback functions describe the effect of genotypic variation of a locus on itself, either directly or mediated by the network. A simple sign rule relates the sign of the derivative of the feedback function of any locus to the feedback loops involving that particular locus. We show that the sign of the phenotypically manifested interaction between alleles at a diploid locus is equal to the sign of the dominant feedback loop involving that particular locus, in accordance with recent results for a single locus system. Our results provide tools by which one can use observable equilibrium concentrations of gene products to disclose structural properties of the network architecture. Our work is a step towards a theory capable of explaining the pleiotropy and epistasis features of genetic variation in complex regulatory networks as functions of regulatory anatomy and functional location of the genetic variation.

Astrocytic mechanisms explaining neural-activity-induced shrinkage of extraneuronal space.

Neuronal stimulation causes approximately 30% shrinkage of the extracellular space (ECS) between neurons and surrounding astrocytes in grey and white matter under experimental conditions. Despite its possible implications for a proper understanding of basic aspects of potassium clearance and astrocyte function, the phenomenon remains unexplained. Here we present a dynamic model that accounts for current experimental data related to the shrinkage phenomenon in wild-type as well as in gene knockout individuals. We find that neuronal release of potassium and uptake of sodium during stimulation, astrocyte uptake of potassium, sodium, and chloride in passive channels, action of the Na/K/ATPase pump, and osmotically driven transport of water through the astrocyte membrane together seem sufficient for generating ECS shrinkage as such. However, when taking into account ECS and astrocyte ion concentrations observed in connection with neuronal stimulation, the actions of the Na(+)/K(+)/Cl(-) (NKCC1) and the Na(+)/HCO(3) (-) (NBC) cotransporters appear to be critical determinants for achieving observed quantitative levels of ECS shrinkage. Considering the current state of knowledge, the model framework appears sufficiently detailed and constrained to guide future key experiments and pave the way for more comprehensive astroglia-neuron interaction models for normal as well as pathophysiological situations.

Analysis of gene regulatory network models with graded and binary transcriptional responses.

The steep sigmoid framework developed by Plahte and Kjøglum [Plahte, E., Kjøglum, S., 2005. Analysis and generic properties of gene regulatory networks with graded response functions. Phys. D 201, 150-176, doi:10.1016/j.physd.2004.11.014] provides a uniform description of gene regulatory networks in which there may be both graded and binary transcriptional responses, as well as a method for analysing the models developed. Here we extend this framework. We show that there is a relation between the location of steady states and the feedback structure of a system, thus generalising existing results for Boolean type models. In addition, we justify underlying assumptions and generic features of the modelling framework in terms of biology and generalise the overall approach to take into account that each transcription factor only regulates one gene at a given threshold. By this assumption, the analysis of the models are greatly simplified.

Targeted reduction of complex models with time scale hierarchy--a case study.

With the increasing flow of biological data there is a growing demand for mathematical tools whereby essential aspects of complex causal dynamic models can be captured and detected by simpler mathematical models without sacrificing too much of the realism provided by the original ones. Given the presence of a time scale hierarchy, singular perturbation techniques represent an elegant method for making such minimised mathematical representations. Any reduction of a complex model by singular perturbation methods is a targeted reduction by the fact that one has to pick certain mechanisms, processes or aspects thought to be essential in a given explanatory context. Here we illustrate how such a targeted reduction of a complex model of melanogenesis in mammals recently developed by the authors provides a way to improve the understanding of how the melanogenic system may behave in a switch-like manner between production of the two major types of melanins. The reduced model is shown by numerical means to be in good quantitative agreement with the original model. Furthermore, it is shown how the reduced model discloses hidden robustness features of the full model, and how the making of a reduced model represents an efficient analytical sensitivity analysis. In addition to yielding new insights concerning the melanogenic system, the paper provides an illustration of a protocol that could be followed to make validated simplifications of complex biological models possessing time scale hierarchies.

Monotonicity is a key feature of genotype-phenotype maps.

It was recently shown that monotone gene action, i.e., order-preservation between allele content and corresponding genotypic values in the mapping from genotypes to phenotypes, is a prerequisite for achieving a predictable parent-offspring relationship across the whole allele frequency spectrum. Here we test the consequential prediction that the design principles underlying gene regulatory networks are likely to generate highly monotone genotype-phenotype maps. To this end we present two measures of the monotonicity of a genotype-phenotype map, one based on allele substitution effects, and the other based on isotonic regression. We apply these measures to genotype-phenotype maps emerging from simulations of 1881 different 3-gene regulatory networks. We confirm that in general, genotype-phenotype maps are indeed highly monotonic across network types. However, regulatory motifs involving incoherent feedforward or positive feedback, as well as pleiotropy in the mapping between genotypes and gene regulatory parameters, are clearly predisposed for generating non-monotonicity. We present analytical results confirming these deep connections between molecular regulatory architecture and monotonicity properties of the genotype-phenotype map. These connections seem to be beyond reach by the classical distinction between additive and non-additive gene action.

Allele interaction--single locus genetics meets regulatory biology.

BACKGROUND: Since the dawn of genetics, additive and dominant gene action in diploids have been defined by comparison of heterozygote and homozygote phenotypes. However, these definitions provide little insight into the underlying intralocus allelic functional dependency and thus cannot serve directly as a mediator between genetics theory and regulatory biology, a link that is sorely needed. METHODOLOGY/PRINCIPAL FINDINGS: We provide such a link by distinguishing between positive, negative and zero allele interaction at the genotype level. First, these distinctions disclose that a biallelic locus can display 18 qualitatively different allele interaction sign motifs (triplets of +, - and 0). Second, we show that for a single locus, Mendelian dominance is not related to heterozygote allele interaction alone, but is actually a function of the degrees of allele interaction in all the three genotypes. Third, we demonstrate how the allele interaction in each genotype is directly quantifiable in gene regulatory models, and that there is a unique, one-to-one correspondence between the sign of autoregulatory feedback loops and the sign of the allele interactions. CONCLUSION/SIGNIFICANCE: The concept of allele interaction refines single locus genetics substantially, and it provides a direct link between classical models of gene action and gene regulatory biology. Together with available empirical data, our results indicate that allele interaction can be exploited experimentally to identify and explain intricate intra- and inter-locus feedback relationships in eukaryotes.

The genotype-phenotype relationship in multicellular pattern-generating models--the neglected role of pattern descriptors.

BACKGROUND: A deep understanding of what causes the phenotypic variation arising from biological patterning processes, cannot be claimed before we are able to recreate this variation by mathematical models capable of generating genotype-phenotype maps in a causally cohesive way. However, the concept of pattern in a multicellular context implies that what matters is not the state of every single cell, but certain emergent qualities of the total cell aggregate. Thus, in order to set up a genotype-phenotype map in such a spatiotemporal pattern setting one is actually forced to establish new pattern descriptors and derive their relations to parameters of the original model. A pattern descriptor is a variable that describes and quantifies a certain qualitative feature of the pattern, for example the degree to which certain macroscopic structures are present. There is today no general procedure for how to relate a set of patterns and their characteristic features to the functional relationships, parameter values and initial values of an original pattern-generating model. Here we present a new, generic approach for explorative analysis of complex patterning models which focuses on the essential pattern features and their relations to the model parameters. The approach is illustrated on an existing model for Delta-Notch lateral inhibition over a two-dimensional lattice. RESULTS: By combining computer simulations according to a succession of statistical experimental designs, computer graphics, automatic image analysis, human sensory descriptive analysis and multivariate data modelling, we derive a pattern descriptor model of those macroscopic, emergent aspects of the patterns that we consider of interest. The pattern descriptor model relates the values of the new, dedicated pattern descriptors to the parameter values of the original model, for example by predicting the parameter values leading to particular patterns, and provides insights that would have been hard to obtain by traditional methods. CONCLUSION: The results suggest that our approach may qualify as a general procedure for how to discover and relate relevant features and characteristics of emergent patterns to the functional relationships, parameter values and initial values of an underlying pattern-generating mathematical model.

Threshold-dominated regulation hides genetic variation in gene expression networks.

BACKGROUND: In dynamical models with feedback and sigmoidal response functions, some or all variables have thresholds around which they regulate themselves or other variables. A mathematical analysis has shown that when the dose-response functions approach binary or on/off responses, any variable with an equilibrium value close to one of its thresholds is very robust to parameter perturbations of a homeostatic state. We denote this threshold robustness. To check the empirical relevance of this phenomenon with response function steepnesses ranging from a near on/off response down to Michaelis-Menten conditions, we have performed a simulation study to investigate the degree of threshold robustness in models for a three-gene system with one downstream gene, using several logical input gates, but excluding models with positive feedback to avoid multistationarity. Varying parameter values representing functional genetic variation, we have analysed the coefficient of variation (CV) of the gene product concentrations in the stable state for the regulating genes in absolute terms and compared to the CV for the unregulating downstream gene. The sigmoidal or binary dose-response functions in these models can be considered as phenomenological models of the aggregated effects on protein or mRNA expression rates of all cellular reactions involved in gene expression. RESULTS: For all the models, threshold robustness increases with increasing response steepness. The CVs of the regulating genes are significantly smaller than for the unregulating gene, in particular for steep responses. The effect becomes less prominent as steepnesses approach Michaelis-Menten conditions. If the parameter perturbation shifts the equilibrium value too far away from threshold, the gene product is no longer an effective regulator and robustness is lost. Threshold robustness arises when a variable is an active regulator around its threshold, and this function is maintained by the feedback loop that the regulator necessarily takes part in and also is regulated by. In the present study all feedback loops are negative, and our results suggest that threshold robustness is maintained by negative feedback which necessarily exists in the homeostatic state. CONCLUSION: Threshold robustness of a variable can be seen as its ability to maintain an active regulation around its threshold in a homeostatic state despite external perturbations. The feedback loop that the system necessarily possesses in this state, ensures that the robust variable is itself regulated and kept close to its threshold. Our results suggest that threshold regulation is a generic phenomenon in feedback-regulated networks with sigmoidal response functions, at least when there is no positive feedback. Threshold robustness in gene regulatory networks illustrates that hidden genetic variation can be explained by systemic properties of the genotype-phenotype map.

Nonlinear regulation enhances the phenotypic expression of trans-acting genetic polymorphisms.

BACKGROUND: Genetic variation explains a considerable part of observed phenotypic variation in gene expression networks. This variation has been shown to be located both locally (cis) and distally (trans) to the genes being measured. Here we explore to which degree the phenotypic manifestation of local and distant polymorphisms is a dynamic feature of regulatory design. RESULTS: By combining mathematical models of gene expression networks with genetic maps and linkage analysis we find that very different network structures and regulatory motifs give similar cis/trans linkage patterns. However, when the shape of the cis-regulatory input functions is more nonlinear or threshold-like, we observe for all networks a dramatic increase in the phenotypic expression of distant compared to local polymorphisms under otherwise equal conditions. CONCLUSION: Our findings indicate that genetic variation affecting the form of cis-regulatory input functions may reshape the genotype-phenotype map by changing the relative importance of cis and trans variation. Our approach combining nonlinear dynamic models with statistical genetics opens up for a systematic investigation of how functional genetic variation is translated into phenotypic variation under various systemic conditions.

The regulatory basis of melanogenic switching.

Melanin produced in follicular melanocytes is the major basis for pigmentation of hair and wool in mammals. Two major types of melanin may be synthesized, the black/brown eumelanin and the reddish/yellow pheomelanin. Based on available cell biological evidence and reasonable assumptions, a mathematical model is developed to improve our understanding of melanogenic switching, i.e. the switching between eumelanin and pheomelanin production depending on the extracellular signalling context. In 1993, Ito proposed that melanogenic switching is due to the covalent binding of the intermediate DOPAquinone to the enzyme glutathione reductase. We were only able to obtain a good fit to available experimental data on the relation between pheomelanin levels and the activity of the key enzyme tyrosinase by taking Ito's hypothesis into account. Thus, our results support Ito's hypothesis, and suggest that melanogenic switching may be due to a jump between two stable production pattern states when the tyrosinase activity varies between two bifurcation levels. This implies that small changes in the levels of external regulatory factors may cause an accentuated change in the proportion of the produced colour pigments and may explain the fact that mammalian coat patterns often exhibit sharply delimited patches of either black or reddish colour.