Simulation of Drosophila circadian oscillations, mutations, and light responses by a model with VRI, PDP-1, and CLK.

A model of Drosophila circadian rhythm generation was developed to represent feedback loops based on transcriptional regulation of per, Clk (dclock), Pdp-1, and vri (vrille). The model postulates that histone acetylation kinetics make transcriptional activation a nonlinear function of [CLK]. Such a nonlinearity is essential to simulate robust circadian oscillations of transcription in our model and in previous models. Simulations suggest that two positive feedback loops involving Clk are not essential for oscillations, because oscillations of [PER] were preserved when Clk, vri, or Pdp-1 expression was fixed. However, eliminating positive feedback by fixing vri expression altered the oscillation period. Eliminating the negative feedback loop in which PER represses per expression abolished oscillations. Simulations of per or Clk null mutations, of per overexpression, and of vri, Clk, or Pdp-1 heterozygous null mutations altered model behavior in ways similar to experimental data. The model simulated a photic phase-response curve resembling experimental curves, and oscillations entrained to simulated light-dark cycles. Temperature compensation of oscillation period could be simulated if temperature elevation slowed PER nuclear entry or PER phosphorylation. The model makes experimental predictions, some of which could be tested in transgenic Drosophila.

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.

The nonsteroidal anti-inflammatory drug indomethacin induces heterogeneity in lipid membranes: potential implication for its diverse biological action.

BACKGROUND: The nonsteroidal anti-inflammatory drug (NSAID), indomethacin (Indo), has a large number of divergent biological effects, the molecular mechanism(s) for which have yet to be fully elucidated. Interestingly, Indo is highly amphiphilic and associates strongly with lipid membranes, which influence localization, structure and function of membrane-associating proteins and actively regulate cell signaling events. Thus, it is possible that Indo regulates diverse cell functions by altering micro-environments within the membrane. Here we explored the effect of Indo on the nature of the segregated domains in a mixed model membrane composed of dipalmitoyl phosphatidyl-choline (di16:0 PC, or DPPC) and dioleoyl phosphatidyl-choline (di18:1 PC or DOPC) and cholesterol that mimics biomembranes. METHODOLOGY/PRINCIPAL FINDINGS: Using a series of fluorescent probes in a fluorescence resonance energy transfer (FRET) study, we found that Indo induced separation between gel domains and fluid domains in the mixed model membrane, possibly by enhancing the formation of gel-phase domains. This effect originated from the ability of Indo to specifically target the ordered domains in the mixed membrane. These findings were further confirmed by measuring the ability of Indo to affect the fluidity-dependent fluorescence quenching and the level of detergent resistance of membranes. CONCLUSION/SIGNIFICANCE: Because the tested lipids are the main lipid constituents in cell membranes, the observed formation of gel phase domains induced by Indo potentially occurs in biomembranes. This marked Indo-induced change in phase behavior potentially alters membrane protein functions, which contribute to the wide variety of biological activities of Indo and other NSAIDs.

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.

Model-Based Estimates of Calanus Finmarchicus Abundance in the Gulf of Maine

Ocean observing systems and satellites routinely collect a wealth of information on physical conditions in the ocean. With few exceptions, such as chlorophyll concentrations, information on biological properties is harder to measure autonomously. Here, we present a system to produce estimates of the distribution and abundance of the copepod Calanus finmarchicus in the Gulf of Maine. Our system uses satellite-based measurements of sea surface temperature and chlorophyll concentration to determine the developmental and reproductive rates of C. finmarchicus. The rate information then drives a population dynamics model of C. finmarchicus that is embedded in a 2-dimensional circulation field. The first generation of this system produces realistic information on interannual variability in C. finmarchicus distribution and abundance during the winter and spring. The model can also be used to identify key drivers of interannual variability in C. finmarchicus. Experiments with the model suggest that changes in initial conditions are overwhelmed by variability in growth rates after approximately 50 d. Temperature has the largest effect on growth rate. Elevated chlorophyll during the late winter can lead to increased C. finmarchicus abundance during the spring, but the effect of variations in chlorophyll concentrations is secondary to the other inputs. Our system could be used to provide real-time estimates or even forecasts of C. finmarchicus distribution. These estimates could then be used to support management of copepod predators such as herring and right whales.

Decadal Variability in the Northeast Pacific in a Physical-Ecosystem Model: Role of Mixed Layer Depth and Trophic Interactions

A basin-wide interdecadal change in both the physical state and the ecology of the North Pacific occurred near the end of 1976. Here we use a physical-ecosystem model to examine whether changes in the physical environment associated with the 1976-1977 transition influenced the lower trophic levels of the food web and if so by what means. The physical component is an ocean general circulation model, while the biological component contains 10 compartments: two phytoplankton, two zooplankton, two detritus pools, nitrate, ammonium, silicate, and carbon dioxide. The model is forced with observed atmospheric fields during 1960-1999. During spring, there is a similar to 40% reduction in plankton biomass in all four plankton groups during 1977-1988 relative to 1970-1976 in the central Gulf of Alaska (GOA). The epoch difference in plankton appears to be controlled by the mixed layer depth. Enhanced Ekman pumping after 1976 caused the halocline to shoal, and thus the mixed layer depth, which extends to the top of the halocline in late winter, did not penetrate as deep in the central GOA. As a result, more phytoplankton remained in the euphotic zone, and phytoplankton biomass began to increase earlier in the year after the 1976 transition. Zooplankton biomass also increased, but then grazing pressure led to a strong decrease in phytoplankton by April followed by a drop in zooplankton by May: Essentially, the mean seasonal cycle of plankton biomass was shifted earlier in the year. As the seasonal cycle progressed, the difference in plankton concentrations between epochs reversed sign again, leading to slightly greater zooplankton biomass during summer in the later epoch.

Bilateral uterine vessel ligation as a model of intrauterine growth restriction in mice

BACKGROUND Intrauterine growth restriction (IUGR) occurs in up to 10% of pregnancies and is considered as a major risk to develop various diseases in adulthood, such as cardiovascular diseases, insulin resistance, hypertension or end stage kidney disease. Several IUGR models have been developed in order to understand the biological processes linked to fetal growth retardation, most of them being rat or mouse models and nutritional models. In order to reproduce altered placental flow, surgical models have also been developed, and among them bilateral uterine ligation has been frequently used. Nevertheless, this model has never been developed in the mouse, although murine tools display multiple advantages for biological research. The aim of this work was therefore to develop a mouse model of bilateral uterine ligation as a surgical model of IUGR. RESULTS In this report, we describe the set up and experimental data obtained from three different protocols (P1, P2, P3) of bilateral uterine vessel ligation in the mouse. Ligation was either performed at the cervical end of each uterine horn (P1) or at the central part of each uterine horn (P2 and P3). Time of surgery was E16 (P1), E17 (P2) or E16.5 (P3). Mortality, maternal weight and abortion parameters were recorded, as well as placentas weights, fetal resorption, viability, fetal weight and size. Results showed that P1 in test animals led to IUGR but was also accompanied with high mortality rate of mothers (50%), low viability of fetuses (8%) and high resorption rate (25%). P2 and P3 improved most of these parameters (decreased mortality and improved pregnancy outcomes; improved fetal viability to 90% and 27%, respectively) nevertheless P2 was not associated to IUGR contrary to P3. Thus P3 experimental conditions enable IUGR with better pregnancy and fetuses outcomes parameters that allow its use in experimental studies. CONCLUSIONS Our results show that bilateral uterine artery ligation according to the protocol we have developed and validated can be used as a surgical mouse model of IUGR.

Molecular Consequences of the SERPINH1/HSP47 Mutation in the Dachshund Natural Model of Osteogenesis Imperfecta.

Osteogenesis imperfecta (OI) is a heritable connective tissue disease characterized by bone fragility and increased risk of fractures. Up to now, mutations in at least 18 genes have been associated with dominant and recessive forms of OI that affect the production or post-translational processing of procollagen or alter bone homeostasis. Among those, SERPINH1 encoding heat shock protein 47 (HSP47), a chaperone exclusive for collagen folding in the ER, was identified to cause a severe form of OI in dachshunds (L326P) as well as in humans (one single case with a L78P mutation). To elucidate the disease mechanism underlying OI in the dog model, we applied a range of biochemical assays to mutant and control skin fibroblasts as well as on bone samples. These experiments revealed that type I collagen synthesized by mutant cells had decreased electrophoretic mobility. Procollagen was retained intracellularly with concomitant dilation of ER cisternae and activation of the ER stress response markers GRP78 and phospho-eIF2α, thus suggesting a defect in procollagen processing. In line with the migration shift detected on SDS-PAGE of cell culture collagen, extracts of bone collagen from the OI dog showed a similar mobility shift, and on tandem mass spectrometry, the chains were post-translationally overmodified. The bone collagen had a higher content of pyridinoline than control dog bone. We conclude that the SERPINH1 mutation in this naturally occurring model of OI impairs how HSP47 acts as a chaperone in the ER. This results in abnormal post-translational modification and cross-linking of the bone collagen.

Intensification of tropical Pacific biological productivity due to volcanic eruptions

Major volcanic eruptions generate widespread ocean cooling, which reduces upper ocean stratification. This effect has the potential to increase nutrient delivery into the euphotic zone and boost biological productivity. Using externally forced last millennium simulations of three climate/Earth System models (Model for Interdisciplinary Research On Climate (MIROC), Community Earth System Model (CESM), and LOch-Vecode-Ecbilt-CLio-agIsm Model (LOVECLIM)), we test the hypothesis that large volcanic eruptions intensify nutrient-driven export production. It is found that strong volcanic radiative forcing enhances the likelihood of eastern Pacific El Niño-like warming in CESM and LOVECLIM. This leads to an initial reduction of nutrients and export production in the eastern equatorial Pacific. However, this initial response reverses after about 3 years in association with La Niña cooling. The resulting delayed enhancement of biological production resembles the multiyear response in MIROC. The model simulations show that volcanic impacts on tropical Pacific dynamics and biogeochemistry persist for several years, thus providing a new source for potential multiyear ecosystem predictability.

Effects of Surface Forcing on Interannual Variability of the Fall Phytoplankton Bloom in the Gulf of Maine Revealed Using a Process-Oriented Model

SeaWiFS (Sea-viewing Wide Field-of-view Sensor) chlorophyll data revealed strong interannual variability in fall phytoplankton dynamics in the Gulf of Maine, with 3 general features in any one year: (1) rapid chlorophyll increases in response to storm events in fall; (2) gradual chlorophyll increases in response to seasonal wind-and cooling-induced mixing that gradually deepens the mixed layer; and (3) the absence of any observable fall bloom. We applied a mixed-layer box model and a 1-dimensional physical-biological numerical model to examine the influence of physical forcing (surface wind, heat flux, and freshening) on the mixed-layer dynamics and its impact on the entrainment of deep-water nutrients and thus on the appearance of fall bloom. The model results suggest that during early fall, the surface mixed-layer depth is controlled by both wind-and cooling-induced mixing. Strong interannual variability in mixed-layer depth has a direct impact on short-and long-term vertical nutrient fluxes and thus the fall bloom. Phytoplankton concentrations over time are sensitive to initial pre-bloom profiles of nutrients. The strength of the initial stratification can affect the modeled phytoplankton concentration, while the timing of intermittent freshening events is related to the significant interannual variability of fall blooms.

Physical and Biological Controls on the Latitudinal Asymmetry of Surface Nutrients and pCO(2) in the Central and Eastern Equatorial Pacific

Surface nutrients and dissolved inorganic carbon (DIC) in the central (CEP) and eastern equatorial Pacific (EEP) show much higher concentrations to the south than to the north of the equator. In this study, the physical and biological controls on this asymmetry are investigated using a coupled physical-biogeochemical model. Two numerical experiments are conducted to examine the effects of asymmetrical photosynthetic efficiency (a) due to asymmetrical iron supply about the equator. The experiment with asymmetrical photosynthesis produces improved results as compared with historical observations. A nitrate budget analysis suggests that in the EEP the divergence of upwelling waters controls the surface nitrate asymmetry with additional contribution from the South Equatorial Current (SEC) carrying nutrient-rich Peru upwelling water. The changes of a affect the surface nitrate distribution but not the overall asymmetry. The SEC further carries excess nitrate to the west and thus extends the asymmetry in the east to the CEP. In the CEP, however, stronger northward than southward transport tends to reduce the nitrate asymmetry, while the asymmetrical photosynthesis would help to maintain it. Similar processes also control the distributions of surface silicate and DIC in the equatorial Pacific, which is also affected by the air-sea CO(2) exchange. The asymmetrical photosynthesis influences the distribution of surface DIC, pCO(2), and the air-sea CO(2) flux, by redistributing about 20% CO(2) flux from the north to the south of the equator. Owing to the adjustment of air-sea CO(2) flux, however, the net surface DIC change is smaller than the direct change associated with primary production.

Emotions, Not Just Decision-Making Processes, Are Critical to an Evolutionary Model of Human Behavior

An evolutionary model of human behavior should privilege emotions: essential, phylogenetically ancient behaviors that learning and decision making only subserve. Infants and non-mammals lack advanced cognitive powers but still survive. Decision making is only a means to emotional ends, which organize and prioritize behavior. The emotion of pride/shame, or dominance striving, bridges the social and biological sciences via internalization of cultural norms.

Validation of a Spatially Continuous EDEN Water-Surface Model for the Everglades, Florida

The Everglades Depth Estimation Network (EDEN) is an integrated network of realtime water-level monitoring, ground-elevation modeling, and water-surface modeling that provides scientists and managers with current (2000-present), online water-stage and water-depth information for the entire freshwater portion of the Greater Everglades. Continuous daily spatial interpolations of the EDEN network stage data are presented on grid with 400-square-meter spacing. EDEN offers a consistent and documented dataset that can be used by scientists and managers to: (1) guide large-scale field operations, (2) integrate hydrologic and ecological responses, and (3) support biological and ecological assessments that measure ecosystem responses to the implementation of the Comprehensive Everglades Restoration Plan (CERP) (U.S. Army Corps of Engineers, 1999). The target users are biologists and ecologists examining trophic level responses to hydrodynamic changes in the Everglades. The first objective of this report is to validate the spatially continuous EDEN water-surface model for the Everglades, Florida developed by Pearlstine et al. (2007) by using an independent field-measured data-set. The second objective is to demonstrate two applications of the EDEN water-surface model: to estimate site-specific ground elevation by using the validated EDEN water-surface model and observed water depth data; and to create water-depth hydrographs for tree islands. We found that there are no statistically significant differences between model-predicted and field-observed water-stage data in both southern Water Conservation Area (WCA) 3A and WCA 3B. Tree island elevations were derived by subtracting field water-depth measurements from the predicted EDEN water-surface. Water-depth hydrographs were then computed by subtracting tree island elevations from the EDEN water stage. Overall, the model is reliable by a root mean square error (RMSE) of 3.31 cm. By region, the RMSE is 2.49 cm and 7.77 cm in WCA 3A and 3B, respectively. This new landscape-scale hydrological model has wide applications for ongoing research and management efforts that are vital to restoration of the Florida Everglades. The accurate, high-resolution hydrological data, generated over broad spatial and temporal scales by the EDEN model, provides a previously missing key to understanding the habitat requirements and linkages among native and invasive populations, including fish, wildlife, wading birds, and plants. The EDEN model is a powerful tool that could be adapted for other ecosystem-scale restoration and management programs worldwide.

High-resolution microarray analysis of chromosome 20q in human colon cancer metastasis model systems

Amplification of human chromosome 20q DNA is the most frequently occurring chromosomal abnormality detected in sporadic colorectal carcinomas and shows significant correlation with liver metastases. Through comprehensive high-resolution microarray comparative genomic hybridization and microarray gene expression profiling, we have characterized chromosome 20q amplicon genes associated with human colorectal cancer metastasis in two in vitro metastasis model systems. The results revealed increasing complexity of the 20q genomic profile from the primary tumor-derived cell lines to the lymph node and liver metastasis derived cell lines. Expression analysis of chromosome 20q revealed a subset of over expressed genes residing within the regions of genomic copy number gain in all the tumor cell lines, suggesting these are Chromosome 20q copy number responsive genes. Bases on their preferential expression levels in the model system cell lines and known biological function, four of the over expressed genes mapping to the common intervals of genomic copy gain were considered the most promising candidate colorectal metastasis-associated genes. Validation of genomic copy number and expression array data was carried out on these genes, with one gene, DNMT3B, standing out as expressed at a relatively higher levels in the metastasis-derived cell lines compared with their primary-derived counterparts in both the models systems analyzed. The data provide evidence for the role of chromosome 20q genes with low copy gain and elevated expression in the clonal evolution of metastatic cells and suggests that such genes may serve as early biomarkers of metastatic potential. The data also support the utility of the combined microarray comparative genomic hybridization and expression array analysis for identifying copy number responsive genes in areas of low DNA copy gain in cancer cells. ^

Concomitant Mycobacterium tuberculosis infection and type 2 diabetes: A descriptive model with global health implications

The type 2 diabetes (diabetes) pandemic is recognized as a threat to tuberculosis (TB) control worldwide. This secondary data analysis project estimated the contribution of diabetes to TB in a binational community on the Texas-Mexico border where both diseases occur. Newly-diagnosed TB patients > 20 years of age were prospectively enrolled at Texas-Mexico border clinics between January 2006 and November 2008. Upon enrollment, information regarding social, demographic, and medical risks for TB was collected at interview, including self-reported diabetes. In addition, self-reported diabetes was supported by blood-confirmation according to guidelines published by the American Diabetes Association (ADA). For this project, data was compared to existing statistics for TB incidence and diabetes prevalence from the corresponding general populations of each study site to estimate the relative and attributable risks of diabetes to TB. In concordance with historical sociodemographic data provided for TB patients with self-reported diabetes, our TB patients with diabetes also lacked the risk factors traditionally associated with TB (alcohol abuse, drug abuse, history of incarceration, and HIV infection); instead, the majority of our TB patients with diabetes were characterized by overweight/obesity, chronic hyperglycemia, and older median age. In addition, diabetes prevalence among our TB patients was significantly higher than in the corresponding general populations. Findings of this study will help accurately characterize TB patients with diabetes, thus aiding in the timely recognition and diagnosis of TB in a population not traditionally viewed as at-risk. We provide epidemiological and biological evidence that diabetes continues to be an increasingly important risk factor for TB.^