Stochastic models and numerical algorithms for a class of regulatory gene networks.

Regulatory gene networks contain generic modules, like those involving feedback loops, which are essential for the regulation of many biological functions (Guido et al. in Nature 439:856-860, 2006). We consider a class of self-regulated genes which are the building blocks of many regulatory gene networks, and study the steady-state distribution of the associated Gillespie algorithm by providing efficient numerical algorithms. We also study a regulatory gene network of interest in gene therapy, using mean-field models with time delays. Convergence of the related time-nonhomogeneous Markov chain is established for a class of linear catalytic networks with feedback loops.

Constraining ACT-R models of decision strategies: An experimental paradigm

It has been repeatedly debated which strategies people rely on in inference. These debates have been difficult to resolve, partially because hypotheses about the decision processes assumed by these strategies have typically been formulated qualitatively, making it hard to test precise quantitative predictions about response times and other behavioral data. One way to increase the precision of strategies is to implement them in cognitive architectures such as ACT-R. Often, however, a given strategy can be implemented in several ways, with each implementation yielding different behavioral predictions. We present and report a study with an experimental paradigm that can help to identify the correct implementations of classic compensatory and non-compensatory strategies such as the take-the-best and tallying heuristics, and the weighted-linear model.

Modelling hydrodynamics in the Rio Parana, Argentina : an evaluation and inter-comparison of reduced-complexity and physics based models applied to a large sand-bed river

Depth-averaged velocities and unit discharges within a 30 km reach of one of the world's largest rivers, the Rio Parana, Argentina, were simulated using three hydrodynamic models with different process representations: a reduced complexity (RC) model that neglects most of the physics governing fluid flow, a two-dimensional model based on the shallow water equations, and a three-dimensional model based on the Reynolds-averaged Navier-Stokes equations. Row characteristics simulated using all three models were compared with data obtained by acoustic Doppler current profiler surveys at four cross sections within the study reach. This analysis demonstrates that, surprisingly, the performance of the RC model is generally equal to, and in some instances better than, that of the physics based models in terms of the statistical agreement between simulated and measured flow properties. In addition, in contrast to previous applications of RC models, the present study demonstrates that the RC model can successfully predict measured flow velocities. The strong performance of the RC model reflects, in part, the simplicity of the depth-averaged mean flow patterns within the study reach and the dominant role of channel-scale topographic features in controlling the flow dynamics. Moreover, the very low water surface slopes that typify large sand-bed rivers enable flow depths to be estimated reliably in the RC model using a simple fixed-lid planar water surface approximation. This approach overcomes a major problem encountered in the application of RC models in environments characterised by shallow flows and steep bed gradients. The RC model is four orders of magnitude faster than the physics based models when performing steady-state hydrodynamic calculations. However, the iterative nature of the RC model calculations implies a reduction in computational efficiency relative to some other RC models. A further implication of this is that, if used to simulate channel morphodynamics, the present RC model may offer only a marginal advantage in terms of computational efficiency over approaches based on the shallow water equations. These observations illustrate the trade off between model realism and efficiency that is a key consideration in RC modelling. Moreover, this outcome highlights a need to rethink the use of RC morphodynamic models in fluvial geomorphology and to move away from existing grid-based approaches, such as the popular cellular automata (CA) models, that remain essentially reductionist in nature. In the case of the world's largest sand-bed rivers, this might be achieved by implementing the RC model outlined here as one element within a hierarchical modelling framework that would enable computationally efficient simulation of the morphodynamics of large rivers over millennial time scales. (C) 2012 Elsevier B.V. All rights reserved.

Integrating species distribution models (SDMs) and phylogeography for two species of Alpine Primula.

The major intention of the present study was to investigate whether an approach combining the use of niche-based palaeodistribution modeling and phylo-geography would support or modify hypotheses about the Quaternary distributional history derived from phylogeographic methods alone. Our study system comprised two closely related species of Alpine Primula. We used species distribution models based on the extant distribution of the species and last glacial maximum (LGM) climate models to predict the distribution of the two species during the LGM. Phylogeographic data were generated using amplified fragment length polymorphisms (AFLPs). In Primula hirsuta, models of past distribution and phylogeographic data are partly congruent and support the hypothesis of widespread nunatak survival in the Central Alps. Species distribution models (SDMs) allowed us to differentiate between alpine regions that harbor potential nunatak areas and regions that have been colonized from other areas. SDMs revealed that diversity is a good indicator for nunataks, while rarity is a good indicator for peripheral relict populations that were not source for the recolonization of the inner Alps. In P. daonensis, palaeo-distribution models and phylogeographic data are incongruent. Besides the uncertainty inherent to this type of modeling approach (e.g., relatively coarse 1-km grain size), disagreement of models and data may partly be caused by shifts of ecological niche in both species. Nevertheless, we demonstrate that the combination of palaeo-distribution modeling with phylogeographical approaches provides a more differentiated picture of the distributional history of species and partly supports (P. hirsuta) and partly modifies (P. daonensis and P. hirsuta) hypotheses of Quaternary distributional history. Some of the refugial area indicated by palaeodistribution models could not have been identified with phylogeographic data.

Phenotypic and molecular characterization of the onset of neuropathy in rodent models of diabetes mellitus

Abstract : The principal focus of this work was to study the molecular changes leading to the development of diabetic peripheral neuropathy (DPN). DPN is the most common complication associated with both type I and II diabetes mellitus (DM). This pathology is the leading cause of non-traumatic amputations. Even though the pathological and morphological changes underlying DPN are relatively well described, the implicated molecular mechanisms remain poorly understood. The following two approaches were developed to study the development of DPN in a rodent model of DM type I. As a first approach, we studied the implication of lipid metabolism in DPN phenotype, concentrating on Sterol Response Element Binding Protein (SREBP)-lc which is the key regulator of storage lipid metabolism. We showed that SREBP-1c was expressed in peripheral nerves and that its expression profile followed the expression of genes involved in storage lipid metabolism. In addition, the expression of SREBP-1c in the endoneurium of peripheral nerves was dependant upon nutritional status and this expression was also perturbed in type I diabetes. In line with this, we showed that insulin elevated the expression of SREBP-1c in primary cultured Schwann cells by activating the SREBP-1c promoter. Taken together, these findings reveal that SREBP-1c expression in Schwann cells responds to metabolic stimuli including insulin and that this response is affected in type I diabetes mellitus. This suggests that disturbed SREBP-1c regulated lipid metabolism may contribute to the pathophysiology of DPN. As a second approach, we performed a comprehensive analysis of the molecular changes associated with DPN in the Akital~1~+ mouse which is a model of spontaneous early-onset type I diabetes mellitus. This mouse expresses a mutated non-functional isoform of insulin, leading to hypoinsulinemia and hyperglycaemia. To determine the onset of DPN, weight, blood glucose and motor nerve conduction velocity (MNCV) were measured in Akital+/+ mice during the first three months of life. A decrease in MNCV was evident akeady one week after the onset of hyperglycemia. To explore the molecular changes associated with the development of DPN in these mice, we performed gene expression profiling using sciatic nerve endoneurium and dorsal root ganglia (DRG) isolated from early diabetic male Akita+/+ mice and sex-matched littermate controls. No major transcriptional changes were detected either in the DRG or in the sciatic nerve endoneurium. This experiment indicates that the phenotypic changes observed during the development of DPN are not correlated with major transcriptional alterations, but mainly with alterations at the protein level. Résumé Lors ce travail, nous nous sommes intéressés aux changements moléculaires aboutissant aux neuropathies périphériques dues au diabète (NPD). Les NPD sont la complication la plus commune du diabète de type I et de type II. Cette pathologie est une cause majeure d'amputations. Même si les changements pathologiques et morphologiques associés aux NPD sont relativement bien décrits, les mécanismes moléculaires provoquant cette pathologie sont mal connus. Deux approches ont principalement été utilisées pour étudier le développement des NPD dans des modèles murins du diabète de type I. Nous avons d'abord étudié l'impact du métabolisme des lipides sur le développement des NPD en nous concentrant sur Sterol Response Element Binding Protein (SREBP)-1 c qui est un régulateur clé des lipides de stockage. Nous avons montré que SREBP-1 c est exprimé dans les nerfs périphériques et que son profil d'expression suit celui de gènes impliqués dans le métabolisme des lipides de stockage. De plus, l'expression de SREBP-1c dans l'endoneurium des nerfs périphériques est dépendante du statut nutritionnel et est dérégulée lors de diabète de type I. Nous avons également pu montrer que l'insuline augmente l'expression de SREBP-1c dans des cultures primaires de cellules de Schwann en activant le promoteur de SREBP-1c. Ses résultats démontrent que l'expression de SREBP-1c dans les cellules de Schwann est contrôlée par des stimuli métaboliques comme l'insuline et que cette réponse est affectée dans le cas d'un diabète de type I. Ces données suggèrent que la dérégulation de l'expression de SREBP-1c lors du diabète pourrait affecter le métabolisme des lipides et ainsi contribuer à la pathophysiologie des NPD. Comme seconde approche, nous avons réalisé une analyse globale des changements moléculaires associés au développement des NPD chez les souris Akita+/+, un modèle de diabète de type I. Cette souris exprime une forme mutée et non fonctionnelle de l'insuline provoquant une hypoinsulinémie et une hyperglycémie. Afin de déterminer le début du développement de la NPD, le poids, le niveau de glucose sanguin et la vitesse de conduction nerveuse (VCN) ont été mesurés durant les 3 premiers mois de vie. Une diminution de la VCN a été détectée une semaine seulement après le développement de l'hyperglycémie. Pour explorer les changements moléculaires associés avec le développement des NPD, nous avons réalisé un profil d'expression de l'endoneurium du nerf sciatique et des ganglions spinaux isolés à partir de souris Akital+/+ et de souris contrôles Akita+/+. Aucune altération transcriptionnelle majeure n'a été détectée dans nos échantillons. Cette expérience suggère que les changements phénotypiques observés durant le développement des NPD ne sont pas corrélés avec des changements importants au niveau transcriptionnel, mais plutôt avec des altérations au niveau protéique. Résumé : Lors ce travail, nous nous sommes intéressés aux changements moléculaires aboutissant aux neuropathies périphériques dues au diabète (NPD). Les NPD sont la complication la plus commune du diabète de type I et de type II. Cette pathologie est une cause majeure d'amputations. Même si les changements pathologiques et morphologiques associés aux NPD sont relativement bien décrits, les mécanismes moléculaires provoquant cette pathologie sont mal connus. Deux approches ont principalement été utilisées pour étudier le développement des NPD dans des modèles murins du diabète de type I. Nous avons d'abord étudié l'impact du métabolisme des lipides sur le développement des NPD en nous concentrant sur Sterol Response Element Binding Protein (SREBP)-1c qui est un régulateur clé des lipides de stockage. Nous avons montré que SREBP-1 c est exprimé dans les nerfs périphériques et que son profil d'expression suit celui de gènes impliqués dans le métabolisme des lipides de stockage. De plus, l'expression de SREBP-1c dans l'endoneurium des nerfs périphériques est dépendante du statut nutritionnel et est dérégulée lors de diabète de type I. Nous avons également pu montrer que l'insuline augmente l'expression de SREBP-1c dans des cultures primaires de cellules de Schwann en activant le promoteur de SREBP-1c. Ses résultats démontrent que l'expression de SREBP-1c dans les cellules de Schwann est contrôlée par des stimuli métaboliques comme l'insuline et que cette réponse est affectée dans le cas d'un diabète de type I. Ces données suggèrent que la dérégulation de l'expression de SREBP-1c lors du diabète pourrait affecter le métabolisme des lipides et ainsi contribuer à la pathophysiologie des NPD. Comme seconde approche, nous avons réalisé une analyse globale des changements moléculaires associés au développement des NPD chez les souris Akita~~Z~+, un modèle de diabète de type I. Cette souris exprime une forme mutée et non fonctionnelle de l'insuline provoquant une hypoinsulinémie et une hyperglycémie. Afin de déterminer le début du développement de la NPD, le poids, le niveau de glucose sanguin et la vitesse de conduction nerveuse (VCN) ont été mesurés durant les 3 premiers mois de vie. Une diminution de la VCN a été détectée une semaine seulement après le développement de l'hyperglycémie. Pour explorer les changements moléculaires associés avec le développement des NPD, nous avons réalisé un profil d'expression de l'endoneurium du nerf sciatique et des ganglions spinaux isolés à partir de souris Akital+/+ et de souris contrôles Akita+/+. Aucune altération transcriptionnelle majeure n'a été détectée dans nos échantillons. Cette expérience suggère que les changements phénotypiques observés durant le développement des NPD ne sont pas corrélés avec des changements importants au niveau transcriptionnel, mais plutôt avec des altérations au niveau protéique.

Effects of alternative sets of climatic predictors on species distribution models and associated estimates of extinction risk: a test with plants in an arid environment

Aim To evaluate the effects of using distinct alternative sets of climatic predictor variables on the performance, spatial predictions and future projections of species distribution models (SDMs) for rare plants in an arid environment. . Location Atacama and Peruvian Deserts, South America (18º30'S - 31º30'S, 0 - 3 000 m) Methods We modelled the present and future potential distributions of 13 species of Heliotropium sect. Cochranea, a plant group with a centre of diversity in the Atacama Desert. We developed and applied a sequential procedure, starting from climate monthly variables, to derive six alternative sets of climatic predictor variables. We used them to fit models with eight modelling techniques within an ensemble forecasting framework, and derived climate change projections for each of them. We evaluated the effects of using these alternative sets of predictor variables on performance, spatial predictions and projections of SDMs using Generalised Linear Mixed Models (GLMM). Results The use of distinct sets of climatic predictor variables did not have a significant effect on overall metrics of model performance, but had significant effects on present and future spatial predictions. Main conclusion Using different sets of climatic predictors can yield the same model fits but different spatial predictions of current and future species distributions. This represents a new form of uncertainty in model-based estimates of extinction risk that may need to be better acknowledged and quantified in future SDM studies.

Integrating nuclear receptor mobility in models of gene regulation.

The mode of action of nuclear receptors in living cells is an actively investigated field but much remains hypothetical due to the lack, until recently, of methods allowing the assessment of molecular mechanisms in vivo. However, these last years, the development of fluorescence microscopy methods has allowed initiating the dissection of the molecular mechanisms underlying gene regulation by nuclear receptors directly in living cells or organisms. Following our analyses on peroxisome proliferator activated receptors (PPARs) in living cells, we discuss here the different models arising from the use of these tools, that attempt to link mobility, DNA binding or chromatin interaction, and transcriptional activity.

Evolutionarily stable learning schedules and cumulative culture in discrete generation models.

Individual learning (e.g., trial-and-error) and social learning (e.g., imitation) are alternative ways of acquiring and expressing the appropriate phenotype in an environment. The optimal choice between using individual learning and/or social learning may be dictated by the life-stage or age of an organism. Of special interest is a learning schedule in which social learning precedes individual learning, because such a schedule is apparently a necessary condition for cumulative culture. Assuming two obligatory learning stages per discrete generation, we obtain the evolutionarily stable learning schedules for the three situations where the environment is constant, fluctuates between generations, or fluctuates within generations. During each learning stage, we assume that an organism may target the optimal phenotype in the current environment by individual learning, and/or the mature phenotype of the previous generation by oblique social learning. In the absence of exogenous costs to learning, the evolutionarily stable learning schedules are predicted to be either pure social learning followed by pure individual learning ("bang-bang" control) or pure individual learning at both stages ("flat" control). Moreover, we find for each situation that the evolutionarily stable learning schedule is also the one that optimizes the learned phenotype at equilibrium.

Global habitat suitability models of terrestrial mammals.

Detailed large-scale information on mammal distribution has often been lacking, hindering conservation efforts. We used the information from the 2009 IUCN Red List of Threatened Species as a baseline for developing habitat suitability models for 5027 out of 5330 known terrestrial mammal species, based on their habitat relationships. We focused on the following environmental variables: land cover, elevation and hydrological features. Models were developed at 300 m resolution and limited to within species' known geographical ranges. A subset of the models was validated using points of known species occurrence. We conducted a global, fine-scale analysis of patterns of species richness. The richness of mammal species estimated by the overlap of their suitable habitat is on average one-third less than that estimated by the overlap of their geographical ranges. The highest absolute difference is found in tropical and subtropical regions in South America, Africa and Southeast Asia that are not covered by dense forest. The proportion of suitable habitat within mammal geographical ranges correlates with the IUCN Red List category to which they have been assigned, decreasing monotonically from Least Concern to Endangered. These results demonstrate the importance of fine-resolution distribution data for the development of global conservation strategies for mammals.

Using niche-based models to improve the sampling of rare species

Because data on rare species usually are sparse, it is important to have efficient ways to sample additional data. Traditional sampling approaches are of limited value for rare species because a very large proportion of randomly chosen sampling sites are unlikely to shelter the species. For these species, spatial predictions from niche-based distribution models can be used to stratify the sampling and increase sampling efficiency. New data sampled are then used to improve the initial model. Applying this approach repeatedly is an adaptive process that may allow increasing the number of new occurrences found. We illustrate the approach with a case study of a rare and endangered plant species in Switzerland and a simulation experiment. Our field survey confirmed that the method helps in the discovery of new populations of the target species in remote areas where the predicted habitat suitability is high. In our simulations the model-based approach provided a significant improvement (by a factor of 1.8 to 4 times, depending on the measure) over simple random sampling. In terms of cost this approach may save up to 70% of the time spent in the field.