A reduced model clarifies the role of feedback loops and time delays in the Drosophila circadian oscillator.

Although several detailed models of molecular processes essential for circadian oscillations have been developed, their complexity makes intuitive understanding of the oscillation mechanism difficult. The goal of the present study was to reduce a previously developed, detailed model to a minimal representation of the transcriptional regulation essential for circadian rhythmicity in Drosophila. The reduced model contains only two differential equations, each with time delays. A negative feedback loop is included, in which PER protein represses per transcription by binding the dCLOCK transcription factor. A positive feedback loop is also included, in which dCLOCK indirectly enhances its own formation. The model simulated circadian oscillations, light entrainment, and a phase-response curve with qualitative similarities to experiment. Time delays were found to be essential for simulation of circadian oscillations with this model. To examine the robustness of the simplified model to fluctuations in molecule numbers, a stochastic variant was constructed. Robust circadian oscillations and entrainment to light pulses were simulated with fewer than 80 molecules of each gene product present on average. Circadian oscillations persisted when the positive feedback loop was removed. Moreover, elimination of positive feedback did not decrease the robustness of oscillations to stochastic fluctuations or to variations in parameter values. Such reduced models can aid understanding of the oscillation mechanisms in Drosophila and in other organisms in which feedback regulation of transcription may play an important role.

Last interglacial temperature evolution — a model inter-comparison

There is a growing number of proxy-based reconstructions detailing the climatic changes that occurred during the last interglacial period (LIG). This period is of special interest, because large parts of the globe were characterized by a warmer-than-present-day climate, making this period an interesting test bed for climate models in light of projected global warming. However, mainly because synchronizing the different palaeoclimatic records is difficult, there is no consensus on a global picture of LIG temperature changes. Here we present the first model inter-comparison of transient simulations covering the LIG period. By comparing the different simulations, we aim at investigating the common signal in the LIG temperature evolution, investigating the main driving forces behind it and at listing the climate feedbacks which cause the most apparent inter-model differences. The model inter-comparison shows a robust Northern Hemisphere July temperature evolution characterized by a maximum between 130–125 ka BP with temperatures 0.3 to 5.3 K above present day. A Southern Hemisphere July temperature maximum, −1.3 to 2.5 K at around 128 ka BP, is only found when changes in the greenhouse gas concentrations are included. The robustness of simulated January temperatures is large in the Southern Hemisphere and the mid-latitudes of the Northern Hemisphere. For these regions maximum January temperature anomalies of respectively −1 to 1.2 K and −0.8 to 2.1 K are simulated for the period after 121 ka BP. In both hemispheres these temperature maxima are in line with the maximum in local summer insolation. In a number of specific regions, a common temperature evolution is not found amongst the models. We show that this is related to feedbacks within the climate system which largely determine the simulated LIG temperature evolution in these regions. Firstly, in the Arctic region, changes in the summer sea-ice cover control the evolution of LIG winter temperatures. Secondly, for the Atlantic region, the Southern Ocean and the North Pacific, possible changes in the characteristics of the Atlantic meridional overturning circulation are crucial. Thirdly, the presence of remnant continental ice from the preceding glacial has shown to be important when determining the timing of maximum LIG warmth in the Northern Hemisphere. Finally, the results reveal that changes in the monsoon regime exert a strong control on the evolution of LIG temperatures over parts of Africa and India. By listing these inter-model differences, we provide a starting point for future proxy-data studies and the sensitivity experiments needed to constrain the climate simulations and to further enhance our understanding of the temperature evolution of the LIG period.

Dynamics of a minimal model of interlocked positive and negative feedback loops of transcriptional regulation by cAMP-response element binding proteins.

cAMP-response element binding (CREB) proteins are involved in transcriptional regulation in a number of cellular processes (e.g., neural plasticity and circadian rhythms). The CREB family contains activators and repressors that may interact through positive and negative feedback loops. These loops can be generated by auto- and cross-regulation of expression of CREB proteins, via CRE elements in or near their genes. Experiments suggest that such feedback loops may operate in several systems (e.g., Aplysia and rat). To understand the functional implications of such feedback loops, which are interlocked via cross-regulation of transcription, a minimal model with a positive and negative loop was developed and investigated using bifurcation analysis. Bifurcation analysis revealed diverse nonlinear dynamics (e.g., bistability and oscillations). The stability of steady states or oscillations could be changed by time delays in the synthesis of the activator (CREB1) or the repressor (CREB2). Investigation of stochastic fluctuations due to small numbers of molecules of CREB1 and CREB2 revealed a bimodal distribution of CREB molecules in the bistability region. The robustness of the stable HIGH and LOW states of CREB expression to stochastic noise differs, and a critical number of molecules was required to sustain the HIGH state for days or longer. Increasing positive feedback or decreasing negative feedback also increased the lifetime of the HIGH state, and persistence of this state may correlate with long-term memory formation. A critical number of molecules was also required to sustain robust oscillations of CREB expression. If a steady state was near a deterministic Hopf bifurcation point, stochastic resonance could induce oscillations. This comparative analysis of deterministic and stochastic dynamics not only provides insights into the possible dynamics of CREB regulatory motifs, but also demonstrates a framework for understanding other regulatory processes with similar network architecture.

A model of late long-term potentiation simulates aspects of memory maintenance.

Late long-term potentiation (L-LTP) denotes long-lasting strengthening of synapses between neurons. L-LTP appears essential for the formation of long-term memory, with memories at least partly encoded by patterns of strengthened synapses. How memories are preserved for months or years, despite molecular turnover, is not well understood. Ongoing recurrent neuronal activity, during memory recall or during sleep, has been hypothesized to preferentially potentiate strong synapses, preserving memories. This hypothesis has not been evaluated in the context of a mathematical model representing ongoing activity and biochemical pathways important for L-LTP. In this study, ongoing activity was incorporated into two such models - a reduced model that represents some of the essential biochemical processes, and a more detailed published model. The reduced model represents synaptic tagging and gene induction simply and intuitively, and the detailed model adds activation of essential kinases by Ca(2+). Ongoing activity was modeled as continual brief elevations of Ca(2+). In each model, two stable states of synaptic strength/weight resulted. Positive feedback between synaptic weight and the amplitude of ongoing Ca(2+) transients underlies this bistability. A tetanic or theta-burst stimulus switches a model synapse from a low basal weight to a high weight that is stabilized by ongoing activity. Bistability was robust to parameter variations in both models. Simulations illustrated that prolonged periods of decreased activity reset synaptic strengths to low values, suggesting a plausible forgetting mechanism. However, episodic activity with shorter inactive intervals maintained strong synapses. Both models support experimental predictions. Tests of these predictions are expected to further understanding of how neuronal activity is coupled to maintenance of synaptic strength. Further investigations that examine the dynamics of activity and synaptic maintenance can be expected to help in understanding how memories are preserved for up to a lifetime in animals including humans.

Paleoclimate data-model comparison and the role of climate forcings over the past 1500 years

The past 1500 years provide a valuable opportunity to study the response of the climate system to external forcings. However, the integration of paleoclimate proxies with climate modeling is critical to improving the understanding of climate dynamics. In this paper, a climate system model and proxy records are therefore used to study the role of natural and anthropogenic forcings in driving the global climate. The inverse and forward approaches to paleoclimate data–model comparison are applied, and sources of uncertainty are identified and discussed. In the first of two case studies, the climate model simulations are compared with multiproxy temperature reconstructions. Robust solar and volcanic signals are detected in Southern Hemisphere temperatures, with a possible volcanic signal detected in the Northern Hemisphere. The anthropogenic signal dominates during the industrial period. It is also found that seasonal and geographical biases may cause multiproxy reconstructions to overestimate the magnitude of the long-term preindustrial cooling trend. In the second case study, the model simulations are compared with a coral δ18O record from the central Pacific Ocean. It is found that greenhouse gases, solar irradiance, and volcanic eruptions all influence the mean state of the central Pacific, but there is no evidence that natural or anthropogenic forcings have any systematic impact on El Niño–Southern Oscillation. The proxy climate relationship is found to change over time, challenging the assumption of stationarity that underlies the interpretation of paleoclimate proxies. These case studies demonstrate the value of paleoclimate data–model comparison but also highlight the limitations of current techniques and demonstrate the need to develop alternative approaches.

Understanding the performance of the FLake model over two African Great Lakes

The ability of the one-dimensional lake model FLake to represent the mixolimnion temperatures for tropical conditions was tested for three locations in East Africa: Lake Kivu and Lake Tanganyika's northern and southern basins. Meteorological observations from surrounding automatic weather stations were corrected and used to drive FLake, whereas a comprehensive set of water temperature profiles served to evaluate the model at each site. Careful forcing data correction and model configuration made it possible to reproduce the observed mixed layer seasonality at Lake Kivu and Lake Tanganyika (northern and southern basins), with correct representation of both the mixed layer depth and water temperatures. At Lake Kivu, mixolimnion temperatures predicted by FLake were found to be sensitive both to minimal variations in the external parameters and to small changes in the meteorological driving data, in particular wind velocity. In each case, small modifications may lead to a regime switch, from the correctly represented seasonal mixed layer deepening to either completely mixed or permanently stratified conditions from similar to 10 m downwards. In contrast, model temperatures were found to be robust close to the surface, with acceptable predictions of near-surface water temperatures even when the seasonal mixing regime is not reproduced. FLake can thus be a suitable tool to parameterise tropical lake water surface temperatures within atmospheric prediction models. Finally, FLake was used to attribute the seasonal mixing cycle at Lake Kivu to variations in the near-surface meteorological conditions. It was found that the annual mixing down to 60m during the main dry season is primarily due to enhanced lake evaporation and secondarily to the decreased incoming long wave radiation, both causing a significant heat loss from the lake surface and associated mixolimnion cooling.

Combined systemic and local morpholino treatment rescues the phenotype of the SMA Delta 7 mouse model

Spinal muscular atrophy (SMA) is a childhood fatal motor neuron disease caused by mutations in the Survival Motor Neuron 1 (SMN1) gene, currently without effective treatment. One possible therapeutic approach is the use of antisense oligonucleotides (ASOs) to redirect the splicing of a paralogous gene, SMN2, to increase the production of functional SMN protein. A range of ASOs with different chemical properties is suitable for these applications, including a morpholino (MO) variant, which has a particularly excellent safety, and efficacy profile. We used a 25- nt MO oligomer sequence against the ISS-N1 region of SMN2 (HSMN2Ex7D(-10-34)) with superior efficacy to previously described sequences also in transgenic SMA Δ7 mice. The combined local and systemic administration of MO (bare or conjugated to octa-guanidine) is necessary to increase full-length SMN expression, leading to robust neuropathological features improvement and survival rescue. Additionally, several snRNA levels that are dysregulated in SMA mice could be restored by MO treatment. These results demonstrate that MO therapy is efficacious and can result in phenotypic rescue. These data provide important insights for the development of therapeutic strategies in SMA patients.

Effect of Combined Systemic and Local Morpholino Treatment on the Spinal Muscular Atrophy Δ7 Mouse Model Phenotype

BACKGROUND: Spinal muscular atrophy (SMA) is a fatal motor neuron disease of childhood that is caused by mutations in the SMN1 gene. Currently, no effective treatment is available. One possible therapeutic approach is the use of antisense oligos (ASOs) to redirect the splicing of the paralogous gene SMN2, thus increasing functional SMN protein production. Various ASOs with different chemical properties are suitable for these applications, including a morpholino oligomer (MO) variant with a particularly excellent safety and efficacy profile. OBJECTIVE: We investigated a 25-nt MO sequence targeting the negative intronic splicing silencer (ISS-N1) 10 to 34 region. METHODS: We administered a 25-nt MO sequence against the ISS-N1 region of SMN2 (HSMN2Ex7D[-10-34]) in the SMAΔ7 mouse model and evaluated the effect and neuropathologic phenotype. We tested different concentrations (from 2 to 24 nM) and delivery protocols (intracerebroventricular injection, systemic injection, or both). We evaluated the treatment efficacy regarding SMN levels, survival, neuromuscular phenotype, and neuropathologic features. RESULTS: We found that a 25-nt MO sequence against the ISS-N1 region of SMN2 (HSMN2Ex7D[-10-34]) exhibited superior efficacy in transgenic SMAΔ7 mice compared with previously described sequences. In our experiments, the combination of local and systemic administration of MO (bare or conjugated to octaguanidine) was the most effective approach for increasing full-length SMN expression, leading to robust improvement in neuropathologic features and survival. Moreover, we found that several small nuclear RNAs were deregulated in SMA mice and that their levels were restored by MO treatment. CONCLUSION: These results indicate that MO-mediated SMA therapy is efficacious and can result in phenotypic rescue, providing important insights for further development of ASO-based therapeutic strategies in SMA patients.

A thermodynamic model for di-trioctahedral chlorite from experimental and natural data in the system MgO-FeO-Al2O3-SiO2-H2O. Applications to P-T sections and geothermometry

We present a new thermodynamic activity-composition model for di-trioctahedral chlorite in the system FeO–MgO–Al2O3–SiO2–H2O that is based on the Holland–Powell internally consistent thermodynamic data set. The model is formulated in terms of four linearly independent end-members, which are amesite, clinochlore, daphnite and sudoite. These account for the most important crystal-chemical substitutions in chlorite, the Fe–Mg, Tschermak and di-trioctahedral substitution. The ideal part of end-member activities is modeled with a mixing-on-site formalism, and non-ideality is described by a macroscopic symmetric (regular) formalism. The symmetric interaction parameters were calibrated using a set of 271 published chlorite analyses for which robust independent temperature estimates are available. In addition, adjustment of the standard state thermodynamic properties of sudoite was required to accurately reproduce experimental brackets involving sudoite. This new model was tested by calculating representative P–T sections for metasediments at low temperatures (<400 °C), in particular sudoite and chlorite bearing metapelites from Crete. Comparison between the calculated mineral assemblages and field data shows that the new model is able to predict the coexistence of chlorite and sudoite at low metamorphic temperatures. The predicted lower limit of the chloritoid stability field is also in better agreement with petrological observations. For practical applications to metamorphic and hydrothermal environments, two new semi-empirical chlorite geothermometers named Chl(1) and Chl(2) were calibrated based on the chlorite + quartz + water equilibrium (2 clinochlore + 3 sudoite = 4 amesite + 4 H2O + 7 quartz). The Chl(1) thermometer requires knowledge of the (Fe3+/ΣFe) ratio in chlorite and predicts correct temperatures for a range of redox conditions. The Chl(2) geothermometer which assumes that all iron in chlorite is ferrous has been applied to partially recrystallized detrital chlorite from the Zone houillère in the French Western Alps.

Waddlia chondrophila induces systemic infection, organ pathology, and elicits Th1-associated humoral immunity in a murine model of genital infection

Waddlia chondrophila is a known bovine abortigenic Chlamydia-related bacterium that has been associated with adverse pregnancy outcomes in human. However, there is a lack of knowledge regarding how W. chondrophila infection spreads, its ability to elicit an immune response and induce pathology. A murine model of genital infection was developed to investigate the pathogenicity and immune response associated with a W. chondrophila infection. Genital inoculation of the bacterial agent resulted in a dose-dependent infection that spread to lumbar lymph nodes and successively to spleen and liver. Bacterial-induced pathology peaked on day 14, characterized by leukocyte infiltration (uterine horn, liver, and spleen), necrosis (liver) and extramedullary hematopoiesis (spleen). Immunohistochemistry demonstrated the presence of a large number of W. chondrophila in the spleen on day 14. Robust IgG titers were detected by day 14 and remained high until day 52. IgG isotypes consisted of high IgG2a, moderate IgG3 and no detectable IgG1, indicating a Th1-associated immune response. This study provides the first evidence that W. chondrophila genital infection is capable of inducing a systemic infection that spreads to major organs, induces uterus, spleen, and liver pathology and elicits a Th1-skewed humoral response. This new animal model will help our understanding of the mechanisms related to intracellular bacteria-induced miscarriages, the most frequent complication of pregnancy that affects one in four women.

A randomized trial of the effects of the noble gases helium and argon on neuroprotection in a rodent cardiac arrest model.

BACKGROUND The noble gas xenon is considered as a neuroprotective agent, but availability of the gas is limited. Studies on neuroprotection with the abundant noble gases helium and argon demonstrated mixed results, and data regarding neuroprotection after cardiac arrest are scant. We tested the hypothesis that administration of 50% helium or 50% argon for 24 h after resuscitation from cardiac arrest improves clinical and histological outcome in our 8 min rat cardiac arrest model. METHODS Forty animals had cardiac arrest induced with intravenous potassium/esmolol and were randomized to post-resuscitation ventilation with either helium/oxygen, argon/oxygen or air/oxygen for 24 h. Eight additional animals without cardiac arrest served as reference, these animals were not randomized and not included into the statistical analysis. Primary outcome was assessment of neuronal damage in histology of the region I of hippocampus proper (CA1) from those animals surviving until day 5. Secondary outcome was evaluation of neurobehavior by daily testing of a Neurodeficit Score (NDS), the Tape Removal Test (TRT), a simple vertical pole test (VPT) and the Open Field Test (OFT). Because of the non-parametric distribution of the data, the histological assessments were compared with the Kruskal-Wallis test. Treatment effect in repeated measured assessments was estimated with a linear regression with clustered robust standard errors (SE), where normality is less important. RESULTS Twenty-nine out of 40 rats survived until day 5 with significant initial deficits in neurobehavioral, but rapid improvement within all groups randomized to cardiac arrest. There were no statistical significant differences between groups neither in the histological nor in neurobehavioral assessment. CONCLUSIONS The replacement of air with either helium or argon in a 50:50 air/oxygen mixture for 24 h did not improve histological or clinical outcome in rats subjected to 8 min of cardiac arrest.

Bayesian estimation of multilevel item response model with missing data

This study proposed a novel statistical method that modeled the multiple outcomes and missing data process jointly using item response theory. This method follows the "intent-to-treat" principle in clinical trials and accounts for the correlation between outcomes and missing data process. This method may provide a good solution to chronic mental disorder study. ^ The simulation study demonstrated that if the true model is the proposed model with moderate or strong correlation, ignoring the within correlation may lead to overestimate of the treatment effect and result in more type I error than specified level. Even if the within correlation is small, the performance of proposed model is as good as naïve response model. Thus, the proposed model is robust for different correlation settings if the data is generated by the proposed model.^

STUDYING AGGREGATE FORMATION BY AMYOTROPHIC LATERAL SCLEROSIS-ASSOCIATED MUTANT SOD1 PROTEIN IN DROSOPHILA MODEL

A common pathological hallmark of most neurodegenerative disorders is the presence of protein aggregates in the brain. Understanding the regulation of aggregate formation is thus important for elucidating disease pathogenic mechanisms and finding effective preventive avenues and cures. Amyotrophic Lateral Sclerosis (ALS), also known as Lou Gehrig’s disease, is a selective neurodegenerative disorder predominantly affecting motor neurons. The majority of ALS cases are sporadic, however, mutations in superoxide dismutase 1 (SOD1) are responsible for about 20% of familial ALS (fALS). Mutated SOD1 proteins are prone to misfold and form protein aggregates, thus representing a good candidate for studying aggregate formation. The long-term goal of this project is to identify regulators of aggregate formation by mutant SOD1 and other ALS-associated disease proteins. The specific aim of this thesis project is to assess the possibility of using the well-established Drosophila model system to study aggregation by human SOD1 (hSOD1) mutants. To this end, using wild type and the three mutant hSOD1 (A4V, G85R and G93A) most commonly found among fALS, I have generated 16 different SOD1 constructs containing either eGFP or mCherry in-frame fluorescent reporters, established and tested both cell- and animal-based Drosophila hSOD1 models. The experimental strategy allows for clear visualization of ectopic hSOD1 expression as well as versatile co-expression schemes to fully investigate protein aggregation specifically by mutant hSOD1. I have performed pilot cell-transfection experiments and verified induced expression of hSOD1 proteins. Using several tissue- or cell type-specific Gal4 lines, I have confirmed the proper expression of hSOD1 from established transgenic fly lines. Interestingly, in both Drosophila S2 cells and different fly tissues including the eye and motor neurons, robust aggregate formation by either wild type or mutant hSOD1 proteins was not observed. These preliminary observations suggest that Drosophila might not be a good experimental organism to study aggregation and toxicity of mutant hSOD1 protein. Nevertheless this preliminary conclusion implies the potential existence of a potent protective mechanism against mutant hSOD1 aggregation and toxicity in Drosophila. Thus, results from my SOD1-ALS project in Drosophila will help future studies on how to best employ this classic model organism to study ALS and other human brain degenerative diseases.

An individual based model spatial distribution, production and biomass for Calanus finmarchicus in the Norwegian Sea

The copepod Calanus finmarchicus is the dominant species of the meso-zooplankton in the Norwegian Sea, and constitutes an important link between the phytoplankton and the higher trophic levels in the Norwegian Sea food chain. An individualbased model for C. finmarchicus, based on super-individuals and evolving traits for behaviour, stages, etc., is two-way coupled to the NORWegian ECOlogical Model system (NORWECOM). One year of modelled C. finmarchicus spatial distribution, production and biomass are found to represent observations reasonably well. High C. finmarchicus abundance is found along the Norwegian shelf-break in the early summer, while the overwintering population is found along the slope and in the deeper Norwegian Sea basins. The timing of the spring bloom is generally later than in the observations. Annual Norwegian Sea production is found to be 29 million tonnes of carbon and a production to biomass (P/B) ratio of 4.3 emerges. Sensitivity tests show that the modelling system is robust to initial values of behavioural traits and with regards to the number of super-individuals simulated given that this is above about 50,000 individuals. Experiments with the model system indicate that it provides a valuable tool for studies of ecosystem responses to causative forces such as prey density or overwintering population size. For example, introducing C. finmarchicus food limitations reduces the stock dramatically, but on the other hand, a reduced stock may rebuild in one year under normal conditions. The NetCDF file contains model grid coordinates and bottom topography.

Lagrangian back-tracing results for Kohnen station; output of time dependent numerical experiment with a 3D nested Antarctic ice sheet model

A nested ice flow model was developed for eastern Dronning Maud Land to assist with the dating and interpretation of the EDML deep ice core. The model consists of a high-resolution higher-order ice dynamic flow model that was nested into a comprehensive 3-D thermomechanical model of the whole Antarctic ice sheet. As the drill site is on a flank position the calculations specifically take into account the effects of horizontal advection as deeper ice in the core originated from higher inland. First the regional velocity field and ice sheet geometry is obtained from a forward experiment over the last 8 glacial cycles. The result is subsequently employed in a Lagrangian backtracing algorithm to provide particle paths back to their time and place of deposition. The procedure directly yields the depth-age distribution, surface conditions at particle origin, and a suite of relevant parameters such as initial annual layer thickness. This paper discusses the method and the main results of the experiment, including the ice core chronology, the non-climatic corrections needed to extract the climatic part of the signal, and the thinning function. The focus is on the upper 89% of the ice core (appr. 170 kyears) as the dating below that is increasingly less robust owing to the unknown value of the geothermal heat flux. It is found that the temperature biases resulting from variations of surface elevation are up to half of the magnitude of the climatic changes themselves.