The Sec1/Munc18 protein Vps45 regulates cellular levels of its SNARE binding partners Tlg2 and Snc2 in Saccharomyces cerevisiae

Intracellular membrane trafficking pathways must be tightly regulated to ensure proper functioning of all eukaryotic cells. Central to membrane trafficking is the formation of specific SNARE (soluble N-ethylmeleimide-sensitive factor attachment protein receptor) complexes between proteins on opposing lipid bilayers. The Sec1/Munc18 (SM) family of proteins play an essential role in SNARE-mediated membrane fusion, and like the SNAREs are conserved through evolution from yeast to humans. The SM protein Vps45 is required for the formation of yeast endosomal SNARE complexes and is thus essential for traffic through the endosomal system. Here we report that, in addition to its role in regulating SNARE complex assembly, Vps45 regulates cellular levels of its SNARE binding partners: the syntaxin Tlg2 and the v-SNARE Snc2: Cells lacking Vps45 have reduced cellular levels of Tlg2 and Snc2; and elevation of Vps45 levels results in concomitant increases in the levels of both Tlg2 and Snc2. As well as regulating traffic through the endosomal system, the Snc v-SNAREs are also required for exocytosis. Unlike most vps mutants, cells lacking Vps45 display multiple growth phenotypes. Here we report that these can be reversed by selectively restoring Snc2 levels in vps45 mutant cells. Our data indicate that as well as functioning as part of the machinery that controls SNARE complex assembly, Vps45 also plays a key role in determining the levels of its cognate SNARE proteins; another key factor in regulation of membrane traffic

Construction and Random Generation of Hypergraphs with Prescribed Degree and Dimension Sequences

We propose algorithms for construction and random generation of hypergraphs without loops and with prescribed degree and dimension sequences. The objective is to provide a starting point for as well as an alternative to Markov chain Monte Carlo approaches. Our algorithms leverage the transposition of properties and algorithms devised for matrices constituted of zeros and ones with prescribed row- and column-sums to hypergraphs. The construction algorithm extends the applicability of Markov chain Monte Carlo approaches when the initial hypergraph is not provided. The random generation algorithm allows the development of a self-normalised importance sampling estimator for hypergraph properties such as the average clustering coefficient.We prove the correctness of the proposed algorithms. We also prove that the random generation algorithm generates any hypergraph following the prescribed degree and dimension sequences with a non-zero probability. We empirically and comparatively evaluate the effectiveness and efficiency of the random generation algorithm. Experiments show that the random generation algorithm provides stable and accurate estimates of average clustering coefficient, and also demonstrates a better effective sample size in comparison with the Markov chain Monte Carlo approaches.Comment: 21 pages, 3 figure

Inferring signalling networks from longitudinal data using sampling based approaches in the R-package 'ddepn'

<p>Abstract</p> <p>Background</p> <p>Network inference from high-throughput data has become an important means of current analysis of biological systems. For instance, in cancer research, the functional relationships of cancer related proteins, summarised into signalling networks are of central interest for the identification of pathways that influence tumour development. Cancer cell lines can be used as model systems to study the cellular response to drug treatments in a time-resolved way. Based on these kind of data, modelling approaches for the signalling relationships are needed, that allow to generate hypotheses on potential interference points in the networks.</p> <p>Results</p> <p>We present the R-package 'ddepn' that implements our recent approach on network reconstruction from longitudinal data generated after external perturbation of network components. We extend our approach by two novel methods: a Markov Chain Monte Carlo method for sampling network structures with two edge types (activation and inhibition) and an extension of a prior model that penalises deviances from a given reference network while incorporating these two types of edges. Further, as alternative prior we include a model that learns signalling networks with the scale-free property.</p> <p>Conclusions</p> <p>The package 'ddepn' is freely available on R-Forge and CRAN <url>http://ddepn.r-forge.r-project.org</url>, <url>http://cran.r-project.org</url>. It allows to conveniently perform network inference from longitudinal high-throughput data using two different sampling based network structure search algorithms.</p

A Bayesian method for calculating real-time quantitative PCR calibration curves using absolute plasmid DNA standards

<p>Abstract</p> <p>Background</p> <p>In real-time quantitative PCR studies using absolute plasmid DNA standards, a calibration curve is developed to estimate an unknown DNA concentration. However, potential differences in the amplification performance of plasmid DNA compared to genomic DNA standards are often ignored in calibration calculations and in some cases impossible to characterize. A flexible statistical method that can account for uncertainty between plasmid and genomic DNA targets, replicate testing, and experiment-to-experiment variability is needed to estimate calibration curve parameters such as intercept and slope. Here we report the use of a Bayesian approach to generate calibration curves for the enumeration of target DNA from genomic DNA samples using absolute plasmid DNA standards.</p> <p>Results</p> <p>Instead of the two traditional methods (classical and inverse), a Monte Carlo Markov Chain (MCMC) estimation was used to generate single, master, and modified calibration curves. The mean and the percentiles of the posterior distribution were used as point and interval estimates of unknown parameters such as intercepts, slopes and DNA concentrations. The software WinBUGS was used to perform all simulations and to generate the posterior distributions of all the unknown parameters of interest.</p> <p>Conclusion</p> <p>The Bayesian approach defined in this study allowed for the estimation of DNA concentrations from environmental samples using absolute standard curves generated by real-time qPCR. The approach accounted for uncertainty from multiple sources such as experiment-to-experiment variation, variability between replicate measurements, as well as uncertainty introduced when employing calibration curves generated from absolute plasmid DNA standards.</p

Conservation of pattern as a tool for inference on spatial snapshots in ecological data

As climate change and other anthropogenic factors increase the uncertainty of vegetation ecosystem persistence, the ability to rapidly assess their dynamics is paramount. Vegetation and sessile communities form a variety of striking regular spatial patterns such as stripes, spots and labyrinths, that have been used as indicators of ecosystem current state, through qualitative analysis of simple models. Here we describe a new method for rigorous quantitative estimation of biological parameters from a single spatial snapshot. We formulate a synthetic likelihood through consideration of the expected change in the correlation structure of the spatial pattern. This then allows Bayesian inference to be performed on the model parameters, which includes providing parameter uncertainty. The method was validated against simulated data and then applied to real data in the form of aerial photographs of seagrass banding. The inferred parameters were found to be able to reproduce similar patterns to those observed and able to detect strength of spatial competition, competition-induced mortality and the local range of reproduction. This technique points to a way of performing rapid inference of spatial competition and ecological stability from a single spatial snapshots of sessile communities

Time scales and exponential trends to equilibrium: Gaussian model problems

We review results on the exponential convergence of multi- dimensional Ornstein-Uhlenbeck processes and discuss related notions of characteristic timescales with concrete model systems. We focus, on the one hand, on exit time distributions and provide ecplicit expressions for the exponential rate of the distribution in the small noise limit. On the other hand, we consider relaxation timescales of the process to its equi- librium measured in terms of relative entropy and discuss the connection with exit probabilities. Along these lines, we study examples which il- lustrate specific properties of the relaxation and discuss the possibility of deriving a simulation-based, empirical definition of slow and fast de- grees of freedom which builds upon a partitioning of the relative entropy functional in conjuction with the observed relaxation behaviour

Bayesian profiling of molecular signatures to predict event times

BACKGROUND: It is of particular interest to identify cancer-specific molecular signatures for early diagnosis, monitoring effects of treatment and predicting patient survival time. Molecular information about patients is usually generated from high throughput technologies such as microarray and mass spectrometry. Statistically, we are challenged by the large number of candidates but only a small number of patients in the study, and the right-censored clinical data further complicate the analysis. RESULTS: We present a two-stage procedure to profile molecular signatures for survival outcomes. Firstly, we group closely-related molecular features into linkage clusters, each portraying either similar or opposite functions and playing similar roles in prognosis; secondly, a Bayesian approach is developed to rank the centroids of these linkage clusters and provide a list of the main molecular features closely related to the outcome of interest. A simulation study showed the superior performance of our approach. When it was applied to data on diffuse large B-cell lymphoma (DLBCL), we were able to identify some new candidate signatures for disease prognosis. CONCLUSION: This multivariate approach provides researchers with a more reliable list of molecular features profiled in terms of their prognostic relationship to the event times, and generates dependable information for subsequent identification of prognostic molecular signatures through either biological procedures or further data analysis

An empirical Bayesian approach for model-based inference of cellular signaling networks

Background A common challenge in systems biology is to infer mechanistic descriptions of biological process given limited observations of a biological system. Mathematical models are frequently used to represent a belief about the causal relationships among proteins within a signaling network. Bayesian methods provide an attractive framework for inferring the validity of those beliefs in the context of the available data. However, efficient sampling of high-dimensional parameter space and appropriate convergence criteria provide barriers for implementing an empirical Bayesian approach. The objective of this study was to apply an Adaptive Markov chain Monte Carlo technique to a typical study of cellular signaling pathways. Results As an illustrative example, a kinetic model for the early signaling events associated with the epidermal growth factor (EGF) signaling network was calibrated against dynamic measurements observed in primary rat hepatocytes. A convergence criterion, based upon the Gelman-Rubin potential scale reduction factor, was applied to the model predictions. The posterior distributions of the parameters exhibited complicated structure, including significant covariance between specific parameters and a broad range of variance among the parameters. The model predictions, in contrast, were narrowly distributed and were used to identify areas of agreement among a collection of experimental studies. Conclusion In summary, an empirical Bayesian approach was developed for inferring the confidence that one can place in a particular model that describes signal transduction mechanisms and for inferring inconsistencies in experimental measurements

Dissection of a QTL Hotspot on Mouse Distal Chromosome 1 that Modulates Neurobehavioral Phenotypes and Gene Expression

A remarkably diverse set of traits maps to a region on mouse distal chromosome 1 (Chr 1) that corresponds to human Chr 1q21–q23. This region is highly enriched in quantitative trait loci (QTLs) that control neural and behavioral phenotypes, including motor behavior, escape latency, emotionality, seizure susceptibility (Szs1), and responses to ethanol, caffeine, pentobarbital, and haloperidol. This region also controls the expression of a remarkably large number of genes, including genes that are associated with some of the classical traits that map to distal Chr 1 (e.g., seizure susceptibility). Here, we ask whether this QTL-rich region on Chr 1 (Qrr1) consists of a single master locus or a mixture of linked, but functionally unrelated, QTLs. To answer this question and to evaluate candidate genes, we generated and analyzed several gene expression, haplotype, and sequence datasets. We exploited six complementary mouse crosses, and combed through 18 expression datasets to determine class membership of genes modulated by Qrr1. Qrr1 can be broadly divided into a proximal part (Qrr1p) and a distal part (Qrr1d), each associated with the expression of distinct subsets of genes. Qrr1d controls RNA metabolism and protein synthesis, including the expression of ∼20 aminoacyl-tRNA synthetases. Qrr1d contains a tRNA cluster, and this is a functionally pertinent candidate for the tRNA synthetases. Rgs7 and Fmn2 are other strong candidates in Qrr1d. FMN2 protein has pronounced expression in neurons, including in the dendrites, and deletion of Fmn2 had a strong effect on the expression of few genes modulated by Qrr1d. Our analysis revealed a highly complex gene expression regulatory interval in Qrr1, composed of multiple loci modulating the expression of functionally cognate sets of genes

A Bayesian hierarchical approach for spatial analysis of climate model bias in multi-model ensembles

Coupled atmosphere–ocean general circulation models are key tools to investigate climate dynamics and the climatic response to external forcings, to predict climate evolution and to generate future climate projections. Current general circulation models are, however, undisputedly affected by substantial systematic errors in their outputs compared to observations. The assessment of these so-called biases, both individually and collectively, is crucial for the models’ evaluation prior to their predictive use. We present a Bayesian hierarchical model for a unified assessment of spatially referenced climate model biases in a multi-model framework. A key feature of our approach is that the model quantifies an overall common bias that is obtained by synthesizing bias across the different climate models in the ensemble, further determining the contribution of each model to the overall bias. Moreover, we determine model-specific individual bias components by characterizing them as non-stationary spatial fields. The approach is illustrated based on the case of near-surface air temperature bias in the tropical Atlantic and bordering regions from a multi-model ensemble of historical simulations from the fifth phase of the Coupled Model Intercomparison Project. The results demonstrate the improved quantification of the bias and interpretative advantages allowed by the posterior distributions derived from the proposed Bayesian hierarchical framework, whose generality favors its broader application within climate model assessment