Application of a computational model for complex fluvial ecosystems: The population dynamics of zebra mussel Dreissena polymorpha as a case study

The potential and adaptive flexibility of population dynamic P-systems (PDP) to study population dynamics suggests that they may be suitable for modelling complex fluvial ecosystems, characterized by a composition of dynamic habitats with many variables that interact simultaneously. Using as a model a reservoir occupied by the zebra mussel Dreissena polymorpha, we designed a computational model based on P systems to study the population dynamics of larvae, in order to evaluate management actions to control or eradicate this invasive species. The population dynamics of this species was simulated under different scenarios ranging from the absence of water flow change to a weekly variation with different flow rates, to the actual hydrodynamic situation of an intermediate flow rate. Our results show that PDP models can be very useful tools to model complex, partially desynchronized, processes that work in parallel. This allows the study of complex hydroecological processes such as the one presented, where reproductive cycles, temperature and water dynamics are involved in the desynchronization of the population dynamics both, within areas and among them. The results obtained may be useful in the management of other reservoirs with similar hydrodynamic situations in which the presence of this invasive species has been documented.

Habitat-specific demography and source-sink dynamics in a population of Siberian jays

P>1. There are a number of models describing population structure, many of which have the capacity to incorporate spatial habitat effects. One such model is the source-sink model, that describes a system where some habitats have a natality that is higher than mortality (source) and others have a mortality that exceeds natality (sink). A source can be maintained in the absence of migration, whereas a sink will go extinct. 2. However, the interaction between population dynamics and habitat quality is complex, and concerns have been raised about the validity of published empirical studies addressing source-sink dynamics. In particular, some of these studies fail to provide data on survival, a significant component in disentangling a sink from a low quality source. Moreover, failing to account for a density-dependent increase in mortality, or decrease in fecundity, can result in a territory being falsely assigned as a sink, when in fact, this density-dependent suppression only decreases the population size to a lower level, hence indicating a 'pseudo-sink'. 3. In this study, we investigate a long-term data set for key components of territory-specific demography (mortality and reproduction) and their relationship to habitat characteristics in the territorial, group-living Siberian jay (Perisoreus infaustus). We also assess territory-specific population growth rates (r), to test whether spatial population dynamics are consistent with the ideas of source-sink dynamics. 4. Although average mortality did not differ between sexes, habitat-specific mortality did. Female mortality was higher in older forests, a pattern not observed in males. Male mortality only increased with an increasing amount of open areas. Moreover, reproductive success was higher further away from human settlement, indicating a strong effect of human-associated nest predators. 5. Averaged over all years, 76% of the territories were sources. These territories generally consisted of less open areas, and were located further away from human settlement. 6. The source-sink model provides a tool for modelling demography in distinct habitat patches of different quality, which can aid in identifying key habitats within the landscape, and thus, reduce the risk of implementing unsound management decisions.

Self-organised transients in a neural mass model of epileptogenic tissue dynamics

Stimulation of human epileptic tissue can induce rhythmic, self-terminating responses on the EEG or ECoG. These responses play a potentially important role in localising tissue involved in the generation of seizure activity, yet the underlying mechanisms are unknown. However, in vitro evidence suggests that self-terminating oscillations in nervous tissue are underpinned by non-trivial spatio-temporal dynamics in an excitable medium. In this study, we investigate this hypothesis in spatial extensions to a neural mass model for epileptiform dynamics. We demonstrate that spatial extensions to this model in one and two dimensions display propagating travelling waves but also more complex transient dynamics in response to local perturbations. The neural mass formulation with local excitatory and inhibitory circuits, allows the direct incorporation of spatially distributed, functional heterogeneities into the model. We show that such heterogeneities can lead to prolonged reverberating responses to a single pulse perturbation, depending upon the location at which the stimulus is delivered. This leads to the hypothesis that prolonged rhythmic responses to local stimulation in epileptogenic tissue result from repeated self-excitation of regions of tissue with diminished inhibitory capabilities. Combined with previous models of the dynamics of focal seizures this macroscopic framework is a first step towards an explicit spatial formulation of the concept of the epileptogenic zone. Ultimately, an improved understanding of the pathophysiologic mechanisms of the epileptogenic zone will help to improve diagnostic and therapeutic measures for treating epilepsy.

A novel computer simulation method for simulating the multiscale transduction dynamics of signal proteins

Signal proteins are able to adapt their response to a change in the environment, governing in this way a broad variety of important cellular processes in living systems. While conventional molecular-dynamics (MD) techniques can be used to explore the early signaling pathway of these protein systems at atomistic resolution, the high computational costs limit their usefulness for the elucidation of the multiscale transduction dynamics of most signaling processes, occurring on experimental timescales. To cope with the problem, we present in this paper a novel multiscale-modeling method, based on a combination of the kinetic Monte-Carlo- and MD-technique, and demonstrate its suitability for investigating the signaling behavior of the photoswitch light-oxygen-voltage-2-Jα domain from Avena Sativa (AsLOV2-Jα) and an AsLOV2-Jα-regulated photoactivable Rac1-GTPase (PA-Rac1), recently employed to control the motility of cancer cells through light stimulus. More specifically, we show that their signaling pathways begin with a residual re-arrangement and subsequent H-bond formation of amino acids near to the flavin-mononucleotide chromophore, causing a coupling between β-strands and subsequent detachment of a peripheral α-helix from the AsLOV2-domain. In the case of the PA-Rac1 system we find that this latter process induces the release of the AsLOV2-inhibitor from the switchII-activation site of the GTPase, enabling signal activation through effector-protein binding. These applications demonstrate that our approach reliably reproduces the signaling pathways of complex signal proteins, ranging from nanoseconds up to seconds at affordable computational costs.

Learning real-world stimuli in a neural network with spike-driven synaptic dynamics

We present a model of spike-driven synaptic plasticity inspired by experimental observations and motivated by the desire to build an electronic hardware device that can learn to classify complex stimuli in a semisupervised fashion. During training, patterns of activity are sequentially imposed on the input neurons, and an additional instructor signal drives the output neurons toward the desired activity. The network is made of integrate-and-fire neurons with constant leak and a floor. The synapses are bistable, and they are modified by the arrival of presynaptic spikes. The sign of the change is determined by both the depolarization and the state of a variable that integrates the postsynaptic action potentials. Following the training phase, the instructor signal is removed, and the output neurons are driven purely by the activity of the input neurons weighted by the plastic synapses. In the absence of stimulation, the synapses preserve their internal state indefinitely. Memories are also very robust to the disruptive action of spontaneous activity. A network of 2000 input neurons is shown to be able to classify correctly a large number (thousands) of highly overlapping patterns (300 classes of preprocessed Latex characters, 30 patterns per class, and a subset of the NIST characters data set) and to generalize with performances that are better than or comparable to those of artificial neural networks. Finally we show that the synaptic dynamics is compatible with many of the experimental observations on the induction of long-term modifications (spike-timing-dependent plasticity and its dependence on both the postsynaptic depolarization and the frequency of pre- and postsynaptic neurons).

UPLC-ESI-TOFMS-based metabolomics and gene expression dynamics inspector self-organizing metabolomic maps as tools for understanding the cellular response to ionizing radiation

Global transcriptomic and proteomic profiling platforms have yielded important insights into the complex response to ionizing radiation (IR). Nonetheless, little is known about the ways in which small cellular metabolite concentrations change in response to IR. Here, a metabolomics approach using ultraperformance liquid chromatography coupled with electrospray time-of-flight mass spectrometry was used to profile, over time, the hydrophilic metabolome of TK6 cells exposed to IR doses ranging from 0.5 to 8.0 Gy. Multivariate data analysis of the positive ions revealed dose- and time-dependent clustering of the irradiated cells and identified certain constituents of the water-soluble metabolome as being significantly depleted as early as 1 h after IR. Tandem mass spectrometry was used to confirm metabolite identity. Many of the depleted metabolites are associated with oxidative stress and DNA repair pathways. Included are reduced glutathione, adenosine monophosphate, nicotinamide adenine dinucleotide, and spermine. Similar measurements were performed with a transformed fibroblast cell line, BJ, and it was found that a subset of the identified TK6 metabolites were effective in IR dose discrimination. The GEDI (Gene Expression Dynamics Inspector) algorithm, which is based on self-organizing maps, was used to visualize dynamic global changes in the TK6 metabolome that resulted from IR. It revealed dose-dependent clustering of ions sharing the same trends in concentration change across radiation doses. "Radiation metabolomics," the application of metabolomic analysis to the field of radiobiology, promises to increase our understanding of cellular responses to stressors such as radiation.

Diversification in the South American Pampas: the genetic and morphological variation of the widespread Petunia axillaris complex (Solanaceae)

Understanding the spatiotemporal distribution of genetic variation and the ways in which this distribution is connected to the ecological context of natural populations is fundamental for understanding the nature and mode of intraspecific and, ultimately, interspecific differentiation. The Petunia axillaris complex is endemic to the grasslands of southern South America and includes three subspecies: P.a.axillaris, P.a.parodii and P.a.subandina. These subspecies are traditionally delimited based on both geography and floral morphology, although the latter is highly variable. Here, we determined the patterns of genetic (nuclear and cpDNA), morphological and ecological (bioclimatic) variation of a large number of P.axillaris populations and found that they are mostly coincident with subspecies delimitation. The nuclear data suggest that the subspecies are likely independent evolutionary units, and their morphological differences may be associated with local adaptations to diverse climatic and/or edaphic conditions and population isolation. The demographic dynamics over time estimated by skyline plot analyses showed different patterns for each subspecies in the last 100000years, which is compatible with a divergence time between 35000 and 107000years ago between P.a.axillaris and P.a.parodii, as estimated with the IMa program. Coalescent simulation tests using Approximate Bayesian Computation do not support previous suggestions of extensive gene flow between P.a.axillaris and P.a.parodii in their contact zone.

Recruitment of EB1, a master regulator of microtubule dynamics, to the surface of the Theileria annulata schizont

The apicomplexan parasite Theileria annulata transforms infected host cells, inducing uncontrolled proliferation and clonal expansion of the parasitized cell population. Shortly after sporozoite entry into the target cell, the surrounding host cell membrane is dissolved and an array of host cell microtubules (MTs) surrounds the parasite, which develops into the transforming schizont. The latter does not egress to invade and transform other cells. Instead, it remains tethered to host cell MTs and, during mitosis and cytokinesis, engages the cell's astral and central spindle MTs to secure its distribution between the two daughter cells. The molecular mechanism by which the schizont recruits and stabilizes host cell MTs is not known. MT minus ends are mostly anchored in the MT organizing center, while the plus ends explore the cellular space, switching constantly between phases of growth and shrinkage (called dynamic instability). Assuming the plus ends of growing MTs provide the first point of contact with the parasite, we focused on the complex protein machinery associated with these structures. We now report how the schizont recruits end-binding protein 1 (EB1), a central component of the MT plus end protein interaction network and key regulator of host cell MT dynamics. Using a range of in vitro experiments, we demonstrate that T. annulata p104, a polymorphic antigen expressed on the schizont surface, functions as a genuine EB1-binding protein and can recruit EB1 in the absence of any other parasite proteins. Binding strictly depends on a consensus SxIP motif located in a highly disordered C-terminal region of p104. We further show that parasite interaction with host cell EB1 is cell cycle regulated. This is the first description of a pathogen-encoded protein to interact with EB1 via a bona-fide SxIP motif. Our findings provide important new insight into the mode of interaction between Theileria and the host cell cytoskeleton.

The role of zinc dynamics in growth hormone secretion

Human growth hormone (GH) causes a variety of physiological and metabolic effects in humans and plays a pivotal role in postnatal growth. In somatotroph cells of the anterior pituitary, GH is stored in concentrated forms in secretory granules to be rapidly released upon GH-releasing hormone stimulation. During the process of secretory granule biogenesis, self-association of GH occurs in the compartments of the early secretory pathway (endoplasmic reticulum and Golgi complex). Since this process is greatly facilitated by the presence of zinc ions, it is of importance to understand the potential role of zinc transporters that participate in the fine-tuning of zinc homeostasis and dynamics, particularly in the early secretory pathway. Thus, the role of zinc transporters in supplying the secretory pathway with the sufficient amount of zinc required for the biogenesis of GH-containing secretory granules is essential for normal secretion. This report, illustrated by a clinical case report on transient neonatal zinc deficiency, focuses on the role of zinc in GH storage in the secretory granules and highlights the role of specific zinc transporters in the early secretory pathway.

Sediment dynamics in the subaquatic channel of the Rhone delta (Lake Geneva, France/Switzerland)

bstract With its smaller size, well-known boundary conditions, and the availability of detailed bathymetric data, Lake Geneva’s subaquatic canyon in the Rhone Delta is an excellent analogue to understand sedimentary pro- cesses in deep-water submarine channels. A multidisciplinary research effort was undertaken to unravel the sediment dynamics in the active canyon. This approach included innovative coring using the Russian MIR sub- mersibles, in situ geotechnical tests, and geophysical, sedimentological, geochemical and radiometric analysis techniques. The canyon floor/levee complex is character- ized by a classic turbiditic system with frequent spillover events. Sedimentary evolution in the active canyon is controlled by a complex interplay between erosion and sedimentation processes. In situ profiling of sediment strength in the upper layer was tested using a dynamic penetrometer and suggests that erosion is the governing mechanism in the proximal canyon floor while sedimen- tation dominates in the levee structure. Sedimentation rates progressively decrease down-channel along the levee structure, with accumulation exceeding 2.6 cm/year in the proximal levee. A decrease in the frequency of turbidites upwards along the canyon wall suggests a progressive confinement of the flow through time. The multi-proxy methodology has also enabled a qualitative slope-stability assessment in the levee structure. The rapid sediment loading, slope undercutting and over-steepening, and increased pore pressure due to high methane concentrations hint at a potential instability of the proximal levees. Fur- thermore, discrete sandy intervals show very high methane concentrations and low shear strength and thus could cor- respond to potentially weak layers prone to scarp failures.

Inverse modelling of the 14C bomb pulse in stalagmites to constrain the dynamics of soil carbon cycling at selected European cave sites

The decomposition of soil organic matter (SOM) is temperature dependent, but its response to a future warmer climate remains equivocal. Enhanced rates of decomposition of SOM under increased global temperatures might cause higher CO2 emissions to the atmosphere, and could therefore constitute a strong positive feedback. The magnitude of this feedback however remains poorly understood, primarily because of the difficulty in quantifying the temperature sensitivity of stored, recalcitrant carbon that comprises the bulk (>90%) of SOM in most soils. In this study we investigated the effects of climatic conditions on soil carbon dynamics using the attenuation of the 14C ‘bomb’ pulse as recorded in selected modern European speleothems. These new data were combined with published results to further examine soil carbon dynamics, and to explore the sensitivity of labile and recalcitrant organic matter decomposition to different climatic conditions. Temporal changes in 14C activity inferred from each speleothem was modelled using a three pool soil carbon inverse model (applying a Monte Carlo method) to constrain soil carbon turnover rates at each site. Speleothems from sites that are characterised by semi-arid conditions, sparse vegetation, thin soil cover and high mean annual air temperatures (MAATs), exhibit weak attenuation of atmospheric 14C ‘bomb’ peak (a low damping effect, D in the range: 55–77%) and low modelled mean respired carbon ages (MRCA), indicating that decomposition is dominated by young, recently fixed soil carbon. By contrast, humid and high MAAT sites that are characterised by a thick soil cover and dense, well developed vegetation, display the highest damping effect (D = c. 90%), and the highest MRCA values (in the range from 350 ± 126 years to 571 ± 128 years). This suggests that carbon incorporated into these stalagmites originates predominantly from decomposition of old, recalcitrant organic matter. SOM turnover rates cannot be ascribed to a single climate variable, e.g. (MAAT) but instead reflect a complex interplay of climate (e.g. MAAT and moisture budget) and vegetation development.

Dynamics of the last four glacial terminations recorded in Lake Van, Turkey

A well-dated suite of Lake Van climate-proxy data covering the last 360 ka documents environmental changes over 4 glacial/interglacial cycles in Eastern Anatolia, Turkey. The picture of cold and dry glacials and warm and wet interglacials emerging from pollen, organic carbon, authigenic carbonate content, elemental profiling by XRF and lithological analyses is inconsistent with classical interpretation of ox- ygen isotopic composition of carbonates pointing to a more complex pattern in Lake Van region. Detailed analysis of glacial terminations allows for the constraining of a depositional model explaining different patterns observed in all the proxies. We hypothesize that variations in relative contribution of rainfall, snowmelt and glacier meltwater recharging the basin have a very important role for all sedimentary processes in Lake Van. Lake level of glacial Lake Van, predominantly fed by snowmelt, was low, the water column was oxic, and carbonates precipitating in the epilimnion recorded the light isotopic signature of inflow. During terminations, increasing rainfall and significant supply of mountain glaciers' meltwater contributed to lake level rise. Increased rainfall enhanced density gradients in the water column, and hindered mixing leading to development of bottom-water anoxia. Carbonates precipitating during terminations show large fluctuations in their isotopic composition. Full interglacial conditions in Lake Van are characterized by high or slowly falling lake level. Rainfall and snowmelt feed the lake but due to re-established mixing, the isotopic composition of authigenic carbonates is heavier and closer to that of evaporation-influenced lake water than that of runoff representing snowmelt and atmospheric precipitation.