Non-linear regression models for Approximate Bayesian Computation

Approximate Bayesian inference on the basis of summary statistics is well-suited to complex problems for which the likelihood is either mathematically or computationally intractable. However the methods that use rejection suffer from the curse of dimensionality when the number of summary statistics is increased. Here we propose a machine-learning approach to the estimation of the posterior density by introducing two innovations. The new method fits a nonlinear conditional heteroscedastic regression of the parameter on the summary statistics, and then adaptively improves estimation using importance sampling. The new algorithm is compared to the state-of-the-art approximate Bayesian methods, and achieves considerable reduction of the computational burden in two examples of inference in statistical genetics and in a queueing model.Comment: 4 figures; version 3 minor changes; to appear in Statistics and Computin

An Evolutionary Reduction Principle for Mutation Rates at Multiple Loci

A model of mutation rate evolution for multiple loci under arbitrary selection is analyzed. Results are obtained using techniques from Karlin (1982) that overcome the weak selection constraints needed for tractability in prior studies of multilocus event models. A multivariate form of the reduction principle is found: reduction results at individual loci combine topologically to produce a surface of mutation rate alterations that are neutral for a new modifier allele. New mutation rates survive if and only if they fall below this surface - a generalization of the hyperplane found by Zhivotovsky et al. (1994) for a multilocus recombination modifier. Increases in mutation rates at some loci may evolve if compensated for by decreases at other loci. The strength of selection on the modifier scales in proportion to the number of germline cell divisions, and increases with the number of loci affected. Loci that do not make a difference to marginal fitnesses at equilibrium are not subject to the reduction principle, and under fine tuning of mutation rates would be expected to have higher mutation rates than loci in mutation-selection balance. Other results include the nonexistence of 'viability analogous, Hardy-Weinberg' modifier polymorphisms under multiplicative mutation, and the sufficiency of average transmission rates to encapsulate the effect of modifier polymorphisms on the transmission of loci under selection. A conjecture is offered regarding situations, like recombination in the presence of mutation, that exhibit departures from the reduction principle. Constraints for tractability are: tight linkage of all loci, initial fixation at the modifier locus, and mutation distributions comprising transition probabilities of reversible Markov chains.Comment: v3: Final corrections. v2: Revised title, reworked and expanded introductory and discussion sections, added corollaries, new results on modifier polymorphisms, minor corrections. 49 pages, 64 reference

Caspase-8 binding to cardiolipin in giant unilamellar vesicles provides a functional docking platform for bid

Caspase-8 is involved in death receptor-mediated apoptosis in type II cells, the proapoptotic programme of which is triggered by truncated Bid. Indeed, caspase-8 and Bid are the known intermediates of this signalling pathway. Cardiolipin has been shown to provide an anchor and an essential activating platform for caspase-8 at the mitochondrial membrane surface. Destabilisation of this platform alters receptor-mediated apoptosis in diseases such as Barth Syndrome, which is characterised by the presence of immature cardiolipin which does not allow caspase-8 binding. We used a simplified in vitro system that mimics contact sites and/or cardiolipin-enriched microdomains at the outer mitochondrial surface in which the platform consisting of caspase-8, Bid and cardiolipin was reconstituted in giant unilamellar vesicles. We analysed these vesicles by flow cytometry and confirm previous results that demonstrate the requirement for intact mature cardiolipin for caspase-8 activation and Bid binding and cleavage. We also used confocal microscopy to visualise the rupture of the vesicles and their revesiculation at smaller sizes due to alteration of the curvature following caspase-8 and Bid binding. Biophysical approaches, including Laurdan fluorescence and rupture/tension measurements, were used to determine the ability of these three components (cardiolipin, caspase-8 and Bid) to fulfil the minimal requirements for the formation and function of the platform at the mitochondrial membrane. Our results shed light on the active functional role of cardiolipin, bridging the gap between death receptors and mitochondria

5-Fluorouracil signaling through a calcium-calmodulin-dependent pathway is required for p53 activation and apoptosis in colon carcinoma cells

5-Fluorouracil (5-FU) is an anti-metabolite that is in clinical use for treatment of several cancers. In cells, it is converted into three distinct fluoro-based nucleotide analogs, which interfere with DNA synthesis and repair, leading to genome impairment and, eventually, apoptotic cell death. Current knowledge states that in certain cell types, 5-FU-induced stress is signaling through a p53-dependent induction of tumor necrosis factor-receptor oligomerization required for death-inducing signaling complex formation and caspase-8 activation. Here we establish a role of calcium (Ca 2+) as a messenger for p53 activation in response to 5-FU. Using a combination of pharmacological and genetic approaches, we show that treatment of colon carcinoma cells stimulates entry of extracellular Ca 2+ through long lasting-type plasma membrane channels, which further directs posttranslational phosphorylation of at least three p53 serine residues (S15, S33 and S37) by means of calmodulin (CaM) activity. Obstructing this pathway by the Ca 2+ -chelator BAPTA (1,2-bis(o-aminophenoxy)ethane- N,N,N',N'-tetraacetic acid) or by inhibitors of CaM efficiently reduces 5-FU-induced caspase activities and subsequent cell death. Moreover, ectopic expression of p53 S15A in HCT116 p53 -/- cells confirmed the importance of a Ca 2+ -CaM-p53 axis in 5-FU-induced extrinsic apoptosis. The fact that a widely used therapeutic drug, such as 5-FU, is operating via this pathway could provide new therapeutic intervention points, or specify new combinatorial treatment regimes. © 2013 Macmillan Publishers Limited.Swedish Science Foundation; Swedish and Stockholm Cancer Societies; Swedish Childhood Cancer Foundation; EC-FP-6 (Chemores); EC-FP-7 (APO-SYS

Caspase-2 is upregulated after sciatic nerve transection and its inhibition protects dorsal root ganglion neurons from Apoptosis after serum withdrawal

Sciatic nerve (SN) transection-induced apoptosis of dorsal root ganglion neurons (DRGN) is one factor determining the efficacy of peripheral axonal regeneration and the return of sensation. Here, we tested the hypothesis that caspase-2(CASP2) orchestrates apoptosis of axotomised DRGN both in vivo and in vitro by disrupting the local neurotrophic supply to DRGN. We observed significantly elevated levels of cleaved CASP2 (C-CASP2), compared to cleaved caspase-3 (C-CASP3), within TUNEL+DRGN and DRG glia (satellite and Schwann cells) after SN transection. A serum withdrawal cell culture model, which induced 40% apoptotic death in DRGN and 60% in glia, was used to model DRGN loss after neurotrophic factor withdrawal. Elevated C-CASP2 and TUNEL were observed in both DRGN and DRG glia, with C-CASP2 localisation shifting from the cytosol to the nucleus, a required step for induction of direct CASP2-mediated apoptosis. Furthermore, siRNAmediated downregulation of CASP2 protected 50% of DRGN from apoptosis after serum withdrawal, while downregulation of CASP3 had no effect on DRGN or DRG glia survival. We conclude that CASP2 orchestrates the death of SN-axotomised DRGN directly and also indirectly through loss of DRG glia and their local neurotrophic factor support. Accordingly, inhibiting CASP2 expression is a potential therapy for improving both the SN regeneration response and peripheral sensory recovery

Differential Dynamics of Transposable Elements during Long-Term Diploidization of Nicotiana Section Repandae (Solanaceae) Allopolyploid Genomes

PubMed ID: 23185607This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Ancient Migratory Events in the Middle East: New Clues from the Y-Chromosome Variation of Modern Iranians

Knowledge of high resolution Y-chromosome haplogroup diversification within Iran provides important geographic context regarding the spread and compartmentalization of male lineages in the Middle East and southwestern Asia. At present, the Iranian population is characterized by an extraordinary mix of different ethnic groups speaking a variety of Indo-Iranian, Semitic and Turkic languages. Despite these features, only few studies have investigated the multiethnic components of the Iranian gene pool. In this survey 938 Iranian male DNAs belonging to 15 ethnic groups from 14 Iranian provinces were analyzed for 84 Y-chromosome biallelic markers and 10 STRs. The results show an autochthonous but non-homogeneous ancient background mainly composed by J2a sub-clades with different external contributions. The phylogeography of the main haplogroups allowed identifying post-glacial and Neolithic expansions toward western Eurasia but also recent movements towards the Iranian region from western Eurasia (R1b-L23), Central Asia (Q-M25), Asia Minor (J2a-M92) and southern Mesopotamia (J1-Page08). In spite of the presence of important geographic barriers (Zagros and Alborz mountain ranges, and the Dasht-e Kavir and Dash-e Lut deserts) which may have limited gene flow, AMOVA analysis revealed that language, in addition to geography, has played an important role in shaping the nowadays Iranian gene pool. Overall, this study provides a portrait of the Y-chromosomal variation in Iran, useful for depicting a more comprehensive history of the peoples of this area as well as for reconstructing ancient migration routes. In addition, our results evidence the important role of the Iranian plateau as source and recipient of gene flow between culturally and genetically distinct population

Recent acquisition of Helicobacter pylori by Baka Pygmies

Both anatomically modern humans and the gastric pathogen Helicobacter pylori originated in Africa, and both species have been associated for at least 100,000 years. Seven geographically distinct H. pylori populations exist, three of which are indigenous to Africa: hpAfrica1, hpAfrica2, and hpNEAfrica. The oldest and most divergent population, hpAfrica2, evolved within San hunter-gatherers, who represent one of the deepest branches of the human population tree. Anticipating the presence of ancient H. pylori lineages within all hunter-gatherer populations, we investigated the prevalence and population structure of H. pylori within Baka Pygmies in Cameroon. Gastric biopsies were obtained by esophagogastroduodenoscopy from 77 Baka from two geographically separated populations, and from 101 non-Baka individuals from neighboring agriculturalist populations, and subsequently cultured for H. pylori. Unexpectedly, Baka Pygmies showed a significantly lower H. pylori infection rate (20.8%) than non-Baka (80.2%). We generated multilocus haplotypes for each H. pylori isolate by DNA sequencing, but were not able to identify Baka-specific lineages, and most isolates in our sample were assigned to hpNEAfrica or hpAfrica1. The population hpNEAfrica, a marker for the expansion of the Nilo-Saharan language family, was divided into East African and Central West African subpopulations. Similarly, a new hpAfrica1 subpopulation, identified mainly among Cameroonians, supports eastern and western expansions of Bantu languages. An age-structured transmission model shows that the low H. pylori prevalence among Baka Pygmies is achievable within the timeframe of a few hundred years and suggests that demographic factors such as small population size and unusually low life expectancy can lead to the eradication of H. pylori from individual human populations. The Baka were thus either H. pylori-free or lost their ancient lineages during past demographic fluctuations. Using coalescent simulations and phylogenetic inference, we show that Baka almost certainly acquired their extant H. pylori through secondary contact with their agriculturalist neighbors

Heterozygosity increases microsatellite mutation rate, linking it to demographic history

Abstract Background Biochemical experiments in yeast suggest a possible mechanism that would cause heterozygous sites to mutate faster than equivalent homozygous sites. If such a process operates, it could undermine a key assumption at the core of population genetic theory, namely that mutation rate and population size are indpendent, because population expansion would increase heterozygosity that in turn would increase mutation rate. Here we test this hypothesis using both direct counting of microsatellite mutations in human pedigrees and an analysis of the relationship between microsatellite length and patterns of demographically-induced variation in heterozygosity. Results We find that microsatellite alleles of any given length are more likely to mutate when their homologue is unusually different in length. Furthermore, microsatellite lengths in human populations do not vary randomly, but instead exhibit highly predictable trends with both distance from Africa, a surrogate measure of genome-wide heterozygosity, and modern population size. This predictability remains even after statistically controlling for non-independence due to shared ancestry among populations. Conclusion Our results reveal patterns that are unexpected under classical population genetic theory, where no mechanism exists capable of linking allele length to extrinsic variables such as geography or population size. However, the predictability of microsatellite length is consistent with heterozygote instability and suggest that this has an important impact on microsatellite evolution. Whether similar processes impact on single nucleotide polymorphisms remains unclear.</p

Essential versus accessory aspects of cell death: recommendations of the NCCD 2015

Cells exposed to extreme physicochemical or mechanical stimuli die in an uncontrollable manner, as a result of their immediate structural breakdown. Such an unavoidable variant of cellular demise is generally referred to as ‘accidental cell death’ (ACD). In most settings, however, cell death is initiated by a genetically encoded apparatus, correlating with the fact that its course can be altered by pharmacologic or genetic interventions. ‘Regulated cell death’ (RCD) can occur as part of physiologic programs or can be activated once adaptive responses to perturbations of the extracellular or intracellular microenvironment fail. The biochemical phenomena that accompany RCD may be harnessed to classify it into a few subtypes, which often (but not always) exhibit stereotyped morphologic features. Nonetheless, efficiently inhibiting the processes that are commonly thought to cause RCD, such as the activation of executioner caspases in the course of apoptosis, does not exert true cytoprotective effects in the mammalian system, but simply alters the kinetics of cellular demise as it shifts its morphologic and biochemical correlates. Conversely, bona fide cytoprotection can be achieved by inhibiting the transduction of lethal signals in the early phases of the process, when adaptive responses are still operational. Thus, the mechanisms that truly execute RCD may be less understood, less inhibitable and perhaps more homogeneous than previously thought. Here, the Nomenclature Committee on Cell Death formulates a set of recommendations to help scientists and researchers to discriminate between essential and accessory aspects of cell death