ABCtoolbox: a versatile toolkit for approximate Bayesian computations

Background The estimation of demographic parameters from genetic data often requires the computation of likelihoods. However, the likelihood function is computationally intractable for many realistic evolutionary models, and the use of Bayesian inference has therefore been limited to very simple models. The situation changed recently with the advent of Approximate Bayesian Computation (ABC) algorithms allowing one to obtain parameter posterior distributions based on simulations not requiring likelihood computations. Results Here we present ABCtoolbox, a series of open source programs to perform Approximate Bayesian Computations (ABC). It implements various ABC algorithms including rejection sampling, MCMC without likelihood, a Particle-based sampler and ABC-GLM. ABCtoolbox is bundled with, but not limited to, a program that allows parameter inference in a population genetics context and the simultaneous use of different types of markers with different ploidy levels. In addition, ABCtoolbox can also interact with most simulation and summary statistics computation programs. The usability of the ABCtoolbox is demonstrated by inferring the evolutionary history of two evolutionary lineages of Microtus arvalis. Using nuclear microsatellites and mitochondrial sequence data in the same estimation procedure enabled us to infer sex-specific population sizes and migration rates and to find that males show smaller population sizes but much higher levels of migration than females. Conclusion ABCtoolbox allows a user to perform all the necessary steps of a full ABC analysis, from parameter sampling from prior distributions, data simulations, computation of summary statistics, estimation of posterior distributions, model choice, validation of the estimation procedure, and visualization of the results.

A new model construction by making a detour via intuitionistic theories I: Operational set theory without choice is Π1-equivalent to KP

We introduce a version of operational set theory, OST−, without a choice operation, which has a machinery for Δ0Δ0 separation based on truth functions and the separation operator, and a new kind of applicative set theory, so-called weak explicit set theory WEST, based on Gödel operations. We show that both the theories and Kripke–Platek set theory KPKP with infinity are pairwise Π1Π1 equivalent. We also show analogous assertions for subtheories with ∈-induction restricted in various ways and for supertheories extended by powerset, beta, limit and Mahlo operations. Whereas the upper bound is given by a refinement of inductive definition in KPKP, the lower bound is by a combination, in a specific way, of realisability, (intuitionistic) forcing and negative interpretations. Thus, despite interpretability between classical theories, we make “a detour via intuitionistic theories”. The combined interpretation, seen as a model construction in the sense of Visser's miniature model theory, is a new way of construction for classical theories and could be said the third kind of model construction ever used which is non-trivial on the logical connective level, after generic extension à la Cohen and Krivine's classical realisability model.

The relation of secondary student's career choice readiness to a six-phase model of career decision-making

Based on common aspects of recent models of career decision-making (CDM) a sixphase model of CDM for secondary students is presented and empirically evaluated. The study tested the hypothesis that students who are in later phases possess more career choice readiness and consider different numbers of career alternatives. 266 Swiss secondary students completed measures tapping phase of CDM, career choice readiness, and number of considered career options. Career choice readiness showed an increase with phase of CDM. Later phases were generally associated with a larger increase in career choice readiness. Number of considered career options showed a curve-linear development with fewer options considered at the beginning and at the end of the process. Male students showed a larger variability in their distribution among the process with more male than female students in the first and last phase of the process. Implications for theory and practice are presented.

The choice of anesthesia influences oxidative energy metabolism and tissue survival in critically ischemic murine skin

In surgical animal studies anesthesia is used regularly. Several reports in the literature demonstrate respiratory and cardiovascular side effects of anesthesiologic agents. The aim of this study was to compare two frequently used anesthesia cocktails (ketamine/xylazine [KX] versus medetomidine/climazolam/fentanyl [MCF]) in skin flap mouse models. Systemic blood values, local metabolic parameters, and surgical outcome should be analyzed in critical ischemic skin flap models. Systemic hypoxia was found in the animals undergoing KX anesthesia compared with normoxia in the MCF group (sO(2): 89.2% +/- 2.4% versus 98.5% +/- 1.2%, P < 0.01). Analysis of tissue metabolism revealed impaired anaerobic oxygen metabolism and increased cellular damage in critical ischemic flap tissue under KX anesthesia (lactate/pyruvate ratio: KX 349.86 +/- 282.38 versus MCF 64.53 +/- 18.63; P < 0.01 and glycerol: KX 333.50 +/- 83.91 micromol/L versus MCF 195.83 +/- 29.49 micromol/L; P < 0.01). After 6 d, different rates of flap tissue necrosis could be detected (MCF 57% +/- 6% versus KX 68% +/- 6%, P < 0.01). In summary we want to point out that the type of anesthesia, the animal model and the goal of the study have to be well correlated. Comparing the effects of KX and MCF anesthesia in mice on surgical outcome was a novel aspect of our study.

Femur Specific Polyaffine Model to Regularize the Log-Domain Demons Registration

Osteoarticular allograft transplantation is a popular treatment method in wide surgical resections with large defects. For this reason hospitals are building bone data banks. Performing the optimal allograft selection on bone banks is crucial to the surgical outcome and patient recovery. However, current approaches are very time consuming hindering an efficient selection. We present an automatic method based on registration of femur bones to overcome this limitation. We introduce a new regularization term for the log-domain demons algorithm. This term replaces the standard Gaussian smoothing with a femur specific polyaffine model. The polyaffine femur model is constructed with two affine (femoral head and condyles) and one rigid (shaft) transformation. Our main contribution in this paper is to show that the demons algorithm can be improved in specific cases with an appropriate model. We are not trying to find the most optimal polyaffine model of the femur, but the simplest model with a minimal number of parameters. There is no need to optimize for different number of regions, boundaries and choice of weights, since this fine tuning will be done automatically by a final demons relaxation step with Gaussian smoothing. The newly developed synthesis approach provides a clear anatomically motivated modeling contribution through the specific three component transformation model, and clearly shows a performance improvement (in terms of anatomical meaningful correspondences) on 146 CT images of femurs compared to a standard multiresolution demons. In addition, this simple model improves the robustness of the demons while preserving its accuracy. The ground truth are manual measurements performed by medical experts.

A genetically explicit model of speciation by sensory drive within a continuous population in aquatic environments

BACKGROUND: The sensory drive hypothesis predicts that divergent sensory adaptation in different habitats may lead to premating isolation upon secondary contact of populations. Speciation by sensory drive has traditionally been treated as a special case of speciation as a byproduct of adaptation to divergent environments in geographically isolated populations. However, if habitats are heterogeneous, local adaptation in the sensory systems may cause the emergence of reproductively isolated species from a single unstructured population. In polychromatic fishes, visual sensitivity might become adapted to local ambient light regimes and the sensitivity might influence female preferences for male nuptial color. In this paper, we investigate the possibility of speciation by sensory drive as a byproduct of divergent visual adaptation within a single initially unstructured population. We use models based on explicit genetic mechanisms for color vision and nuptial coloration. RESULTS: We show that in simulations in which the adaptive evolution of visual pigments and color perception are explicitly modeled, sensory drive can promote speciation along a short selection gradient within a continuous habitat and population. We assumed that color perception evolves to adapt to the modal light environment that individuals experience and that females prefer to mate with males whose nuptial color they are most sensitive to. In our simulations color perception depends on the absorption spectra of an individual's visual pigments. Speciation occurred most frequently when the steepness of the environmental light gradient was intermediate and dispersal distance of offspring was relatively small. In addition, our results predict that mutations that cause large shifts in the wavelength of peak absorption promote speciation, whereas we did not observe speciation when peak absorption evolved by stepwise mutations with small effect. CONCLUSION: The results suggest that speciation can occur where environmental gradients create divergent selection on sensory modalities that are used in mate choice. Evidence for such gradients exists from several animal groups, and from freshwater and marine fishes in particular. The probability of speciation in a continuous population under such conditions may then critically depend on the genetic architecture of perceptual adaptation and female mate choice.

Variation in habitat choice and delayed reproduction: Adaptive queuing strategies or individual quality differences?

In most species, some individuals delay reproduction or occupy inferior breeding positions. The queue hypothesis tries to explain both patterns by proposing that individuals strategically delay breeding (queue) to acquire better breeding or social positions. In 1995, Ens, Weissing, and Drent addressed evolutionarily stable queuing strategies in situations with habitat heterogeneity. However, their model did not consider the non - mutually exclusive individual quality hypothesis, which suggests that some individuals delay breeding or occupy inferior breeding positions because they are poor competitors. Here we extend their model with individual differences in competitive abilities, which are probably plentiful in nature. We show that including even the smallest competitive asymmetries will result in individuals using queuing strategies completely different from those in models that assume equal competitors. Subsequently, we investigate how well our models can explain settlement patterns in the wild, using a long-term study on oystercatchers. This long-lived shorebird exhibits strong variation in age of first reproduction and territory quality. We show that only models that include competitive asymmetries can explain why oystercatchers' settlement patterns depend on natal origin. We conclude that predictions from queuing models are very sensitive to assumptions about competitive asymmetries, while detecting such differences in the wild is often problematic.

The antioxidant glutathione in the fish cell lines EPC and BCF-2: response to model pro-oxidants as measured by three different fluorescent dyes

Reduced glutathione (GSH) protects cells against injury by oxidative stress and maintains a range of vital functions. In vitro cell cultures have been used as experimental models to study the role of GSH in chemical toxicity in mammals; however, this approach has been rarely used with fish cells to date. The present study aimed to evaluate sensitivity and specificity of three fluorescent dyes for measuring pro-oxidant-induced changes of GSH contents in fish cell lines: monochlorobimane (mBCl), 5-chloromethylfluorescein diacetate (CMFDA) and 7-amino-4-chloromethylcoumarin (CMAC-blue). Two cell lines were studied, the EPC line established from a skin tumour of carp Cyprinus carpio, and BF-2 cells established from fins of bluegill sunfish Lepomis macrochirus. The cells were exposed for 6 and 24 h to low cytotoxic concentrations of pro-oxidants including hydrogen peroxide, paraquat (PQ), copper and the GSH synthesis inhibitor, L-buthionine-SR-sulfoximine (BSO). The results indicate moderate differences in the GSH response between EPC and BF-2 cells, but distinct differences in the magnitude of the GSH response for the four pro-oxidants. Further, the choice of GSH dye can critically affect the results, with CMFDA appearing to be less specific for GSH than mBCl and CMAC-blue.

Carbon dioxide and climate impulse response functions for the computation of greenhouse gas metrics: a multi-model analysis

The responses of carbon dioxide (CO2) and other climate variables to an emission pulse of CO2 into the atmosphere are often used to compute the Global Warming Potential (GWP) and Global Temperature change Potential (GTP), to characterize the response timescales of Earth System models, and to build reduced-form models. In this carbon cycle-climate model intercomparison project, which spans the full model hierarchy, we quantify responses to emission pulses of different magnitudes injected under different conditions. The CO2 response shows the known rapid decline in the first few decades followed by a millennium-scale tail. For a 100 Gt-C emission pulse added to a constant CO2 concentration of 389 ppm, 25 ± 9% is still found in the atmosphere after 1000 yr; the ocean has absorbed 59 ± 12% and the land the remainder (16 ± 14%). The response in global mean surface air temperature is an increase by 0.20 ± 0.12 °C within the first twenty years; thereafter and until year 1000, temperature decreases only slightly, whereas ocean heat content and sea level continue to rise. Our best estimate for the Absolute Global Warming Potential, given by the time-integrated response in CO2 at year 100 multiplied by its radiative efficiency, is 92.5 × 10−15 yr W m−2 per kg-CO2. This value very likely (5 to 95% confidence) lies within the range of (68 to 117) × 10−15 yr W m−2 per kg-CO2. Estimates for time-integrated response in CO2 published in the IPCC First, Second, and Fourth Assessment and our multi-model best estimate all agree within 15% during the first 100 yr. The integrated CO2 response, normalized by the pulse size, is lower for pre-industrial conditions, compared to present day, and lower for smaller pulses than larger pulses. In contrast, the response in temperature, sea level and ocean heat content is less sensitive to these choices. Although, choices in pulse size, background concentration, and model lead to uncertainties, the most important and subjective choice to determine AGWP of CO2 and GWP is the time horizon.

Regional climate model simulations for Europe at 6 and 0.2 k BP: sensitivity to changes in anthropogenic deforestation

This study aims to evaluate the direct effects of anthropogenic deforestation on simulated climate at two contrasting periods in the Holocene, ~6 and ~0.2 k BP in Europe. We apply We apply the Rossby Centre regional climate model RCA3, a regional climate model with 50 km spatial resolution, for both time periods, considering three alternative descriptions of the past vegetation: (i) potential natural vegetation (V) simulated by the dynamic vegetation model LPJ-GUESS, (ii) potential vegetation with anthropogenic land use (deforestation) from the HYDE3.1 (History Database of the Global Environment) scenario (V + H3.1), and (iii) potential vegetation with anthropogenic land use from the KK10 scenario (V + KK10). The climate model results show that the simulated effects of deforestation depend on both local/regional climate and vegetation characteristics. At ~6 k BP the extent of simulated deforestation in Europe is generally small, but there are areas where deforestation is large enough to produce significant differences in summer temperatures of 0.5–1 °C. At ~0.2 k BP, extensive deforestation, particularly according to the KK10 model, leads to significant temperature differences in large parts of Europe in both winter and summer. In winter, deforestation leads to lower temperatures because of the differences in albedo between forested and unforested areas, particularly in the snow-covered regions. In summer, deforestation leads to higher temperatures in central and eastern Europe because evapotranspiration from unforested areas is lower than from forests. Summer evaporation is already limited in the southernmost parts of Europe under potential vegetation conditions and, therefore, cannot become much lower. Accordingly, the albedo effect dominates in southern Europe also in summer, which implies that deforestation causes a decrease in temperatures. Differences in summer temperature due to deforestation range from −1 °C in south-western Europe to +1 °C in eastern Europe. The choice of anthropogenic land-cover scenario has a significant influence on the simulated climate, but uncertainties in palaeoclimate proxy data for the two time periods do not allow for a definitive discrimination among climate model results.

Reconstructing the subsurface ocean decadal variability using surface nudging in a perfect model framework

Initialising the ocean internal variability for decadal predictability studies is a new area of research and a variety of ad hoc methods are currently proposed. In this study, we explore how nudging with sea surface temperature (SST) and salinity (SSS) can reconstruct the threedimensional variability of the ocean in a perfect model framework. This approach builds on the hypothesis that oceanic processes themselves will transport the surface information into the ocean interior as seen in ocean-only simulations. Five nudged simulations are designed to reconstruct a 150 years ‘‘target’’ simulation, defined as a portion of a long control simulation. The nudged simulations differ by the variables restored to, SST or SST + SSS, and by the area where the nudging is applied. The strength of the heat flux feedback is diagnosed from observations and the restoring coefficients for SSS use the same time-scale. We observed that this choice prevents spurious convection at high latitudes and near sea-ice border when nudging both SST and SSS. In the tropics, nudging the SST is enough to reconstruct the tropical atmosphere circulation and the associated dynamical and thermodynamical impacts on the underlying ocean. In the tropical Pacific Ocean, the profiles for temperature show a significant correlation from the surface down to 2,000 m, due to dynamical adjustment of the isopycnals. At mid-tohigh latitudes, SSS nudging is required to reconstruct both the temperature and the salinity below the seasonal thermocline. This is particularly true in the North Atlantic where adding SSS nudging enables to reconstruct the deep convection regions of the target. By initiating a previously documented 20-year cycle of the model, the SST + SSS nudging is also able to reproduce most of the AMOC variations, a key source of decadal predictability. Reconstruction at depth does not significantly improve with amount of time spent nudging and the efficiency of the surface nudging rather depends on the period/events considered. The joint SST + SSS nudging applied verywhere is the most efficient approach. It ensures that the right water masses are formed at the right surface density, the subsequent circulation, subduction and deep convection further transporting them at depth. The results of this study underline the potential key role of SSS for decadal predictability and further make the case for sustained largescale observations of this field.

A conceptual model of compensation/decompensation in lumbar segmental instability

Lumbar spinal instability (LSI) is a common spinal disorder and can be associated with substantial disability. The concept of defining clinically relevant classifications of disease or 'target condition' is used in diagnostic research. Applying this concept to LSI we hypothesize that a set of clinical and radiological criteria can be developed to identify patients with this target condition who are at high risk of 'irreversible' decompensated LSI for whom surgery becomes the treatment of choice. In LSI, structural deterioration of the lumbar disc initiates a degenerative cascade of segmental instability. Over time, radiographic signs become visible: traction spurs, facet joint degeneration, misalignment, stenosis, olisthesis and de novo scoliosis. Ligaments, joint capsules, local and distant musculature are the functional elements of the lumbar motion segment. Influenced by non-functional factors, these functional elements allow a compensation of degeneration of the motion segment. Compensation may happen on each step of the degenerative cascade but cannot reverse it. However, compensation of LSI may lead to an alleviation or resolution of clinical symptoms. In return, the target condition of decompensation of LSI may cause the new occurrence of symptoms and pain. Functional compensation and decompensation are subject to numerous factors that can change which makes estimation of an individual's long-term prognosis difficult. Compensation and decompensation may influence radiographic signs of degeneration, e.g. the degree of misalignment and segmental angulation caused by LSI is influenced by the tonus of the local musculature. This conceptual model of compensation/decompensation may help solve the debate on functional and psychosocial factors that influence low back pain and to establish a new definition of non-specific low back pain. Individual differences of identical structural disorders could be explained by compensated or decompensated LSI leading to changes in clinical symptoms and pain. Future spine surgery will have to carefully define and measure functional aspects of LSI, e.g. to identify a point of no return where multidisciplinary interventions do not allow a re-compensation and surgery becomes the treatment of choice.