A hierarchical approach to model parameter optimization for developmental systems.

In the context of Systems Biology, computer simulations of gene regulatory networks provide a powerful tool to validate hypotheses and to explore possible system behaviors. Nevertheless, modeling a system poses some challenges of its own: especially the step of model calibration is often difficult due to insufficient data. For example when considering developmental systems, mostly qualitative data describing the developmental trajectory is available while common calibration techniques rely on high-resolution quantitative data. Focusing on the calibration of differential equation models for developmental systems, this study investigates different approaches to utilize the available data to overcome these difficulties. More specifically, the fact that developmental processes are hierarchically organized is exploited to increase convergence rates of the calibration process as well as to save computation time. Using a gene regulatory network model for stem cell homeostasis in Arabidopsis thaliana the performance of the different investigated approaches is evaluated, documenting considerable gains provided by the proposed hierarchical approach.

The application of Fick's law to evaluate the period of stay of mega-olivines in alkaline lavas

The integration of the differential equation of the second law of Fick applied to the diffusion of chemical elements in a semi-infinite solid made it easier to estimate the time of stay of olivine mega-cristals in contact with the host lava The results of this research show the existence of two groups of olivine. The first remained in contact with the magmatic liquid during 19 to 22 days, while the second remained so during only 5 to 9 days. This distinction is correlative to that based on the qualitative observation.

Ruin Theory with Excess of Loss Reinsurance and Reinstatements

The present paper studies the probability of ruin of an insurer, if excess of loss reinsurance with reinstatements is applied. In the setting of the classical Cramer-Lundberg risk model, piecewise deterministic Markov processes are used to describe the free surplus process in this more general situation. It is shown that the finite-time ruin probability is both the solution of a partial integro-differential equation and the fixed point of a contractive integral operator. We exploit the latter representation to develop and implement a recursive algorithm for numerical approximation of the ruin probability that involves high-dimensional integration. Furthermore we study the behavior of the finite-time ruin probability under various levels of initial surplus and security loadings and compare the efficiency of the numerical algorithm with the computational alternative of stochastic simulation of the risk process. (C) 2011 Elsevier Inc. All rights reserved.

Active strain and activation models in cardiac electromechanics

We present a model for mechanical activation of the cardiac tissue depending on the evolution of the transmembrane electrical potential and certain gating/ionic variables that are available in most of electrophysiological descriptions of the cardiac membrane. The basic idea consists in adding to the chosen ionic model one ordinary differential equation for the kinetics of the mechanical activation function. A relevant example illustrates the desired properties of the proposed model, such as delayed muscle contraction and correct magnitude of the muscle fibers' shortening.

Aging and motor programming: differential effects according to the selection of action

There are controversial reports about the effect of aging on movement preparation, and it is unclear to which extent cognitive and/or motor related cerebral processes may be affected. This study examines the age effects on electro-cortical oscillatory patterns during various motor programming tasks, in order to assess potential differences according to the mode of action selection. Twenty elderly (EP, 60-84 years) and 20 young (YP, 20-29 years) participants with normal cognition underwent 3 pre-cued response tasks (S1-S2 paradigm). S1 carried either complete information on response side (Full; stimulus-driven motor preparation), no information (None; general motor alertness), or required free response side selection (Free; internally-driven motor preparation). Electroencephalogram (EEG) was recorded using 64 surface electrodes. Alpha (8-12 Hz) desynchronization (ERD)/synchronization (ERS) and motor-related amplitude asymmetries (MRAA) were analyzed during the S1-S2 interval. Reaction times (RTs) to S2 were slower in EP than YP, and in None than in the other 2 tasks. There was an Age x Task interaction due to increased RTs in Free compared to Full in EP only. Central bilateral and midline activation (alpha ERD) was smaller in EP than YP in None. In Full just before S2, readiness to move was reflected by posterior midline inhibition (alpha ERS) in both groups. In Free, such inhibition was present only in YP. Moreover, MRAA showed motor activity lateralization in both groups in Full, but only in YP in Free. The results indicate reduced recruitment of motor regions for motor alertness in the elderly. They further show less efficient cerebral processes subtending free selection of movement in elders, suggesting reduced capacity for internally-driven action with age.

DIFFERENTIAL EXPRESSION OF GLUTARYL-CoA DEHYDROGENASE IN ADULT RAT CNS, PERIPHERAL TISSUES AND DURING EMBRYONIC DEVELOPMENT

Glutaryl-CoA dehydrogenase (GCDH, EC 1.3.99.7) deficiency, known as glutaric acidemia type I, is one of the more common organic acidurias. To investigate the role of this pathway in different organs we studied the tissue-specific expression pattern of rat Gcdh. The open reading frame cDNA of the rat Gcdh gene was cloned from rat brain mRNA by RT-PCR, allowing the synthesis of digoxigenin-labeled in situ hybridization (ISH) riboprobes. Gcdh mRNA expression was analyzed by ISH on cryosections of adult rat brain, kidney, liver, spleen and heart muscle, as well as on E15 and E18 rat embryos. Gcdh was found expressed in the whole rat brain, almost exclusively in neurons. Gcdh was absent from astrocytes but expressed in rare oligodendrocytes. Strong Gcdh expression was found in liver and spleen, where expression appears predominant to lymphatic nodules. In kidney, the highest Gcdh expression is found in the juxtamedullar cortex (but not in glomerula), and at lower levels in medulla. Heart muscle was negative. During embryonic development, Gcdh was found well expressed in liver, intestinal mucosa and skin, as well as at lower levels in CNS. Further studies are ongoing to provide evidence on the presence of the entire pathway in CNS in order to understand the mechanisms leading to neurotoxicity in glutaric aciduria. The high expression of Gcdh in kidney may explain why certain patients with residual enzyme activity are low excretors at the urine metabolite level.

Differential association of drinking motives with alcohol use on weekdays and weekends.

Drinking motives (DM) reflect the reasons why individuals drink alcohol. Weekdays are mainly dedicated to work, whereas weekends are generally associated with spending time with friends during special events or leisure activities; using alcohol on weekdays and weekends may also be related to different DM. This study examined whether DM were differentially associated with drinking volume (DV) on weekdays and weekends. A representative sample of 5,391 young Swiss men completed a questionnaire assessing weekday and weekend DV, as well as their DM, namely, enhancement, social, coping, and conformity motives. Associations of DM with weekday and weekend DV were examined using structural equation models. Each DM was tested individually in a separate model; all associations were positive and generally stronger (except conformity) for weekend rather than for weekday DV. Further specific patterns of association were found when DM were entered into a single model simultaneously. Associations with weekday and with weekend DV were positive for enhancement and coping motives. However, associations were stronger with weekend rather than with weekday DV for enhancement, and stronger with weekday than with weekend DV for coping motives. Associations of social motives were not significant with weekend DV and negative with weekday DV. Conformity motives were negatively associated with weekend DV and positively related to weekday DV. These results suggest that interventions targeting enhancement motives should be particularly effective at decreasing weekend drinking, whereas interventions targeted at coping motives would be particularly effective at reducing alcohol use on weekdays. (PsycINFO Database Record (c) 2014 APA, all rights reserved).

An algebraic operator approach to the analysis of Gerber-Shiu functions

We introduce an algebraic operator framework to study discounted penalty functions in renewal risk models. For inter-arrival and claim size distributions with rational Laplace transform, the usual integral equation is transformed into a boundary value problem, which is solved by symbolic techniques. The factorization of the differential operator can be lifted to the level of boundary value problems, amounting to iteratively solving first-order problems. This leads to an explicit expression for the Gerber-Shiu function in terms of the penalty function.

Comparative genomics of ParaHox clusters of teleost fishes: gene cluster breakup and the retention of gene sets following whole genome duplications.

BACKGROUND: The evolutionary lineage leading to the teleost fish underwent a whole genome duplication termed FSGD or 3R in addition to two prior genome duplications that took place earlier during vertebrate evolution (termed 1R and 2R). Resulting from the FSGD, additional copies of genes are present in fish, compared to tetrapods whose lineage did not experience the 3R genome duplication. Interestingly, we find that ParaHox genes do not differ in number in extant teleost fishes despite their additional genome duplication from the genomic situation in mammals, but they are distributed over twice as many paralogous regions in fish genomes. RESULTS: We determined the DNA sequence of the entire ParaHox C1 paralogon in the East African cichlid fish Astatotilapia burtoni, and compared it to orthologous regions in other vertebrate genomes as well as to the paralogous vertebrate ParaHox D paralogons. Evolutionary relationships among genes from these four chromosomal regions were studied with several phylogenetic algorithms. We provide evidence that the genes of the ParaHox C paralogous cluster are duplicated in teleosts, just as it had been shown previously for the D paralogon genes. Overall, however, synteny and cluster integrity seems to be less conserved in ParaHox gene clusters than in Hox gene clusters. Comparative analyses of non-coding sequences uncovered conserved, possibly co-regulatory elements, which are likely to contain promoter motives of the genes belonging to the ParaHox paralogons. CONCLUSION: There seems to be strong stabilizing selection for gene order as well as gene orientation in the ParaHox C paralogon, since with a few exceptions, only the lengths of the introns and intergenic regions differ between the distantly related species examined. The high degree of evolutionary conservation of this gene cluster's architecture in particular - but possibly clusters of genes more generally - might be linked to the presence of promoter, enhancer or inhibitor motifs that serve to regulate more than just one gene. Therefore, deletions, inversions or relocations of individual genes could destroy the regulation of the clustered genes in this region. The existence of such a regulation network might explain the evolutionary conservation of gene order and orientation over the course of hundreds of millions of years of vertebrate evolution. Another possible explanation for the highly conserved gene order might be the existence of a regulator not located immediately next to its corresponding gene but further away since a relocation or inversion would possibly interrupt this interaction. Different ParaHox clusters were found to have experienced differential gene loss in teleosts. Yet the complete set of these homeobox genes was maintained, albeit distributed over almost twice the number of chromosomes. Selection due to dosage effects and/or stoichiometric disturbance might act more strongly to maintain a modal number of homeobox genes (and possibly transcription factors more generally) per genome, yet permit the accumulation of other (non regulatory) genes associated with these homeobox gene clusters.

Differential damage in the frontal cortex with aging, sporadic and familial Alzheimer's disease.

In order to understand relationships between executive and structural deficits in the frontal cortex of patients within normal aging or Alzheimer's disease, we studied frontal pathological changes in young and old controls compared to cases with sporadic (AD) or familial Alzheimer's disease (FAD). We performed a semi-automatic computer assisted analysis of the distribution of beta-amyloid (Abeta) deposits revealed by Abeta immunostaining as well as of neurofibrillary tangles (NFT) revealed by Gallyas silver staining in Brodman areas 10 (frontal polar), 12 (ventro-infero-median) and 24 (anterior cingular), using tissue samples from 5 FAD, 6 sporadic AD and 10 control brains. We also performed densitometric measurements of glial fibrillary acidic protein, principal compound of intermediate filaments of astrocytes, and of phosphorylated neurofilament H and M epitopes in areas 10 and 24. All regions studied seem almost completely spared in normal old controls, with only the oldest ones exhibiting a weak percentage of beta-amyloid deposit and hardly any NFT. On the contrary, all AD and FAD cases were severely damaged as shown by statistically significant increased percentages of beta-amyloid deposit, as well as by a high number of NFT. FAD cases (all from the same family) had statistically more beta-amyloid and GFAP than sporadic AD cases in both areas 10 and 24 and statistically more NFT only in area 24. The correlation between the percentage of beta-amyloid and the number of NFT was significant only for area 24. Altogether, these data suggest that the frontal cortex can be spared by AD type lesions in normal aging, but is severely damaged in sporadic and still more in familial Alzheimer's disease. The frontal regions appear to be differentially vulnerable, with area 12 having the less amyloid burden, area 24 the less NFT and area 10 having both more amyloid and more NFT. This pattern of damage in frontal regions may represent a strong neuroanatomical support for the deterioration of attention and cognitive capacities as well as for the presence of emotional and behavioral troubles in AD patients.

Further insights on the French WISC-IV factor structure through Bayesian structural equation modeling

The interpretation of the Wechsler Intelligence Scale for Children-Fourth Edition (WISC-IV) is based on a 4-factor model, which is only partially compatible with the mainstream Cattell-Horn-Carroll (CHC) model of intelligence measurement. The structure of cognitive batteries is frequently analyzed via exploratory factor analysis and/or confirmatory factor analysis. With classical confirmatory factor analysis, almost all crossloadings between latent variables and measures are fixed to zero in order to allow the model to be identified. However, inappropriate zero cross-loadings can contribute to poor model fit, distorted factors, and biased factor correlations; most important, they do not necessarily faithfully reflect theory. To deal with these methodological and theoretical limitations, we used a new statistical approach, Bayesian structural equation modeling (BSEM), among a sample of 249 French-speaking Swiss children (8-12 years). With BSEM, zero-fixed cross-loadings between latent variables and measures are replaced by approximate zeros, based on informative, small-variance priors. Results indicated that a direct hierarchical CHC-based model with 5 factors plus a general intelligence factor better represented the structure of the WISC-IV than did the 4-factor structure and the higher order models. Because a direct hierarchical CHC model was more adequate, it was concluded that the general factor should be considered as a breadth rather than a superordinate factor. Because it was possible for us to estimate the influence of each of the latent variables on the 15 subtest scores, BSEM allowed improvement of the understanding of the structure of intelligence tests and the clinical interpretation of the subtest scores.

Assessment of speed of human locomotion using a differential satellite global positioning system.

PURPOSE: The objective was to explore whether a satellite-based navigation system, global positioning system used in differential mode (DGPS), could accurately assess the speed of running in humans. METHODS: A subject was equipped with a portable GPS receptor coupled to a receiver for differential corrections, while running outdoors on a straight asphalt road at 27 different speeds. Actual speed (reference method) was assessed by chronometry. RESULTS: The accuracy of speed prediction had a standard deviation (SD) of 0.08 km x h(-1) for walking, 0.11 km x h(-1) for running, yielding a coefficient of variation (SD/mean) of 1.38% and 0.82%, respectively. There was a highly significant linear relationship between actual and DGPS speed assessment (r2 = 0.999) with little bias in the prediction equation, because the slope of the regression line was close to unity (0.997). CONCLUSION: the DGPS technique appears to be a valid and inconspicuous tool for "on line" monitoring of the speed of displacement of individuals located on any field on earth, for prolonged periods of time and unlimited distance, but only in specific environmental conditions ("open sky"). Furthermore, the accuracy of speed assessment using the differential GPS mode was improved by a factor of 10 as compared to non-differential GPS.

Chronic wound healing by fetal cell therapy may be explained by differential gene profiling observed in fetal versus old skin cells.

Engineering of fetal tissue has a high potential for the treatment of acute and chronic wounds of the skin in humans as these cells have high expansion capacity under simple culture conditions and one organ donation can produce Master Cell Banks which can fabricate over 900 million biological bandages (9 x 12cm). In a Phase 1 clinical safety study, cases are presented for the treatment of therapy resistant leg ulcers. All eight patients, representing 13 ulcers, tolerated multiple treatments with fetal biological bandages showing no negative secondary effects and repair processes similar to that seen in 3rd degree burns. Differential gene profiling using Affymetrix gene chips (analyzing 12,500 genes) were accomplished on these banked fetal dermal skin cells compared to banked dermal skin cells of an aged donor in order to point to potential indicators of wound healing. Families of genes involved in cell adhesion and extracellular matrix, cell cycle, cellular signaling, development and immune response show significant differences in regulation between banked fetal and those from banked old skin cells: with approximately 47.0% of genes over-expressed in fetal fibroblasts. It is perhaps these differences which contribute to efficient tissue repair seen in the clinic with fetal cell therapy.

Determination of extreme floods using a distributed hydrological model

Hydrological models developed for extreme precipitation of PMP type are difficult to calibrate because of the scarcity of available data for these events. This article presents the process and results of calibration for a distributed hydrological model at fine scale developed for the estimation of probable maximal floods in the case of a PMP. This calibration is done on two Swiss catchments for two events of summer storms. The calculation done is concentrated on the estimation of the parameters of the model, divided in two parts. The first is necessary for the computation of flow speeds while the second is required for the determination of the initial and final infiltration capacities for each terrain type. The results, validated with the Nash equation show a good correlation between the simulated and observed flows. We also apply this model on two Romanian catchments, showing the river network and estimated flow.

Differential stress reaction of human colon cells to oleic-acid-stabilized and unstabilized ultrasmall iron oxide nanoparticles.

Therapeutic engineered nanoparticles (NPs), including ultrasmall superparamagnetic iron oxide (USPIO) NPs, may accumulate in the lower digestive tract following ingestion or injection. In order to evaluate the reaction of human colon cells to USPIO NPs, the effects of non-stabilized USPIO NPs (NS-USPIO NPs), oleic-acid-stabilized USPIO NPs (OA-USPIO NPs), and free oleic acid (OA) were compared in human HT29 and CaCo2 colon epithelial cancer cells. First the biophysical characteristics of NS-USPIO NPs and OA-USPIO NPs in water, in cell culture medium supplemented with fetal calf serum, and in cell culture medium preconditioned by HT29 and CaCo₂ cells were determined. Then, stress responses of the cells were evaluated following exposure to NS-USPIO NPs, OA-USPIO NPs, and free OA. No modification of the cytoskeletal actin network was observed. Cell response to stress, including markers of apoptosis and DNA repair, oxidative stress and degradative/autophagic stress, induction of heat shock protein, or lipid metabolism was determined in cells exposed to the two NPs. Induction of an autophagic response was observed in the two cell lines for both NPs but not free OA, while the other stress responses were cell- and NP-specific. The formation of lipid vacuoles/droplets was demonstrated in HT29 and CaCo₂ cells exposed to OA-USPIO NPs but not to NS-USPIO NPs, and to a much lower level in cells exposed to equimolar concentrations of free OA. Therefore, the induction of lipid vacuoles in colon cells exposed to OA utilized as a stabilizer for USPIO NPs is higly amplified compared to free OA, and is not observed in the absence of this lipid in NS-USPIO NPs.