Fast Bayesian parameter estimation for stochastic logistic growth models

The transition density of a stochastic, logistic population growth model with multiplicative intrinsic noise is analytically intractable. Inferring model parameter values by fitting such stochastic differential equation (SDE) models to data therefore requires relatively slow numerical simulation. Where such simulation is prohibitively slow, an alternative is to use model approximations which do have an analytically tractable transition density, enabling fast inference. We introduce two such approximations, with either multiplicative or additive intrinsic noise, each derived from the linear noise approximation of the logistic growth SDE. After Bayesian inference we find that our fast LNA models, using Kalman filter recursion for computation of marginal likelihoods, give similar posterior distributions to slow arbitrarily exact models. We also demonstrate that simulations from our LNA models better describe the characteristics of the stochastic logistic growth models than a related approach. Finally, we demonstrate that our LNA model with additive intrinsic noise and measurement error best describes an example set of longitudinal observations of microbial population size taken from a typical, genome-wide screening experiment.Comment: 24 pages, 5 figures and 2 table

Pif1- and Exo1-dependent nucleases coordinate checkpoint activation following telomere uncapping

In the absence of the telomere capping protein Cdc13, budding yeast telomeres erode, resulting in checkpoint arrest. This study shows that the helicase Pif1, known as a telomerase inhibitor, also has a direct role in the resection of uncapped telomeres, acting in parallel to the nuclease Exo1

Pho4 mediates phosphate acquisition in Candida albicans and is vital for stress resistance and metal homeostasis

ACKNOWLEDGMENTS We thank Karl Kuchler for the C. albicans superoxide dismutase mutants used in this study. This work was funded by a Medical Research Council Doctoral Training Program studentship to M.A.C.I.; Wellcome Trust Grants 089930 to J.Q., 080088 to A.J.P.B., and 097377 to J.Q., A.J.P.B., and L.P.E.; Biotechnology and Biological Sciences Research Council Grants BB/K016393/1 to J.Q. and BB/F00513X/1 and BB/K017365/1 to A.J.P.B.; European Research Council STRIFE Advanced Grant ERC-2009-AdG-249793 to A.J.P.B.; and Wellcome Trust and Royal Society Sir Henry Dale Fellowship 098375/Z/12/Z to K.J.W. and E.T. The funders had no role in study design, data collection, or interpretation or the decision to submit the work for publication.Peer reviewedPublisher PD

Colonyzer: automated quantification of micro-organism growth characteristics on solid agar

<p>Abstract</p> <p>Background</p> <p>High-throughput screens comparing growth rates of arrays of distinct micro-organism cultures on solid agar are useful, rapid methods of quantifying genetic interactions. Growth rate is an informative phenotype which can be estimated by measuring cell densities at one or more times after inoculation. Precise estimates can be made by inoculating cultures onto agar and capturing cell density frequently by plate-scanning or photography, especially throughout the exponential growth phase, and summarising growth with a simple dynamic model (e.g. the logistic growth model). In order to parametrize such a model, a robust image analysis tool capable of capturing a wide range of cell densities from plate photographs is required.</p> <p>Results</p> <p>Colonyzer is a collection of image analysis algorithms for automatic quantification of the size, granularity, colour and location of micro-organism cultures grown on solid agar. Colonyzer is uniquely sensitive to extremely low cell densities photographed after dilute liquid culture inoculation (spotting) due to image segmentation using a mixed Gaussian model for plate-wide thresholding based on pixel intensity. Colonyzer is robust to slight experimental imperfections and corrects for lighting gradients which would otherwise introduce spatial bias to cell density estimates without the need for imaging dummy plates. Colonyzer is general enough to quantify cultures growing in any rectangular array format, either growing after pinning with a dense inoculum or growing with the irregular morphology characteristic of spotted cultures. Colonyzer was developed using the open source packages: Python, RPy and the Python Imaging Library and its source code and documentation are available on SourceForge under GNU General Public License. Colonyzer is adaptable to suit specific requirements: e.g. automatic detection of cultures at irregular locations on streaked plates for robotic picking, or decreasing analysis time by disabling components such as lighting correction or colour measures.</p> <p>Conclusion</p> <p>Colonyzer can automatically quantify culture growth from large batches of captured images of microbial cultures grown during genome-wide scans over the wide range of cell densities observable after highly dilute liquid spot inoculation, as well as after more concentrated pinning inoculation. Colonyzer is open-source, allowing users to assess it, adapt it to particular research requirements and to contribute to its development.</p

Deletion of the gabra2 gene results in hypersensitivity to the acute effects of ethanol but does not alter ethanol self administration

Human genetic studies have suggested that polymorphisms of the GABRA2 gene encoding the GABA(A) α2-subunit are associated with ethanol dependence. Variations in this gene also convey sensitivity to the subjective effects of ethanol, indicating a role in mediating ethanol-related behaviours. We therefore investigated the consequences of deleting the α2-subunit on the ataxic and rewarding properties of ethanol in mice. Ataxic and sedative effects of ethanol were explored in GABA(A) α2-subunit wildtype (WT) and knockout (KO) mice using a Rotarod apparatus, wire hang and the duration of loss of righting reflex. Following training, KO mice showed shorter latencies to fall than WT littermates under ethanol (2 g/kg i.p.) in both Rotarod and wire hang tests. After administration of ethanol (3.5 g/kg i.p.), KO mice took longer to regain the righting reflex than WT mice. To ensure the acute effects are not due to the gabra2 deletion affecting pharmacokinetics, blood ethanol concentrations were measured at 20 minute intervals after acute administration (2 g/kg i.p.), and did not differ between genotypes. To investigate ethanol's rewarding properties, WT and KO mice were trained to lever press to receive increasing concentrations of ethanol on an FR4 schedule of reinforcement. Both WT and KO mice self-administered ethanol at similar rates, with no differences in the numbers of reinforcers earned. These data indicate a protective role for α2-subunits, against the acute sedative and ataxic effects of ethanol. However, no change was observed in ethanol self administration, suggesting the rewarding effects of ethanol remain unchange

Histone methyltransferase Dot1 and Rad9 inhibit single-stranded DNA accumulation at DSBs and uncapped telomeres

Cells respond to DNA double-strand breaks (DSBs) and uncapped telomeres by recruiting checkpoint and repair factors to the site of lesions. Single-stranded DNA (ssDNA) is an important intermediate in the repair of DSBs and is produced also at uncapped telomeres. Here, we provide evidence that binding of the checkpoint protein Rad9, through its Tudor domain, to methylated histone H3-K79 inhibits resection at DSBs and uncapped telomeres. Loss of DOT1 or mutations in RAD9 influence a Rad50-dependent nuclease, leading to more rapid accumulation of ssDNA, and faster activation of the critical checkpoint kinase, Mec1. Moreover, deletion of RAD9 or DOT1 partially bypasses the requirement for CDK1 in DSB resection. Interestingly, Dot1 contributes to checkpoint activation in response to low levels of telomere uncapping but is not essential with high levels of uncapping. We suggest that both Rad9 and histone H3 methylation allow transmission of the damage signal to checkpoint kinases, and keep resection of damaged DNA under control influencing, both positively and negatively, checkpoint cascades and contributing to a tightly controlled response to DNA damage

Quantitative Fitness Analysis Shows That NMD Proteins and Many Other Protein Complexes Suppress or Enhance Distinct Telomere Cap Defects

To better understand telomere biology in budding yeast, we have performed systematic suppressor/enhancer analyses on yeast strains containing a point mutation in the essential telomere capping gene CDC13 (cdc13-1) or containing a null mutation in the DNA damage response and telomere capping gene YKU70 (yku70Δ). We performed Quantitative Fitness Analysis (QFA) on thousands of yeast strains containing mutations affecting telomere-capping proteins in combination with a library of systematic gene deletion mutations. To perform QFA, we typically inoculate 384 separate cultures onto solid agar plates and monitor growth of each culture by photography over time. The data are fitted to a logistic population growth model; and growth parameters, such as maximum growth rate and maximum doubling potential, are deduced. QFA reveals that as many as 5% of systematic gene deletions, affecting numerous functional classes, strongly interact with telomere capping defects. We show that, while Cdc13 and Yku70 perform complementary roles in telomere capping, their genetic interaction profiles differ significantly. At least 19 different classes of functionally or physically related proteins can be identified as interacting with cdc13-1, yku70Δ, or both. Each specific genetic interaction informs the roles of individual gene products in telomere biology. One striking example is with genes of the nonsense-mediated RNA decay (NMD) pathway which, when disabled, suppress the conditional cdc13-1 mutation but enhance the null yku70Δ mutation. We show that the suppressing/enhancing role of the NMD pathway at uncapped telomeres is mediated through the levels of Stn1, an essential telomere capping protein, which interacts with Cdc13 and recruitment of telomerase to telomeres. We show that increased Stn1 levels affect growth of cells with telomere capping defects due to cdc13-1 and yku70Δ. QFA is a sensitive, high-throughput method that will also be useful to understand other aspects of microbial cell biology

Taming the tiger by the tail: modulation of DNA damage responses by telomeres

Telomeres are by definition stable and inert chromosome ends, whereas internal chromosome breaks are potent stimulators of the DNA damage response (DDR). Telomeres do not, as might be expected, exclude DDR proteins from chromosome ends but instead engage with many DDR proteins. However, the most powerful DDRs, those that might induce chromosome fusion or cell-cycle arrest, are inhibited at telomeres. In budding yeast, many DDR proteins that accumulate most rapidly at double strand breaks (DSBs), have important functions in physiological telomere maintenance, whereas DDR proteins that arrive later tend to have less important functions. Considerable diversity in telomere structure has evolved in different organisms and, perhaps reflecting this diversity, different DDR proteins seem to have distinct roles in telomere physiology in different organisms. Drawing principally on studies in simple model organisms such as budding yeast, in which many fundamental aspects of the DDR and telomere biology have been established; current views on how telomeres harness aspects of DDR pathways to maintain telomere stability and permit cell-cycle division are discussed

Zds2p Regulates Swe1p-dependent Polarized Cell Growth in Saccharomyces cerevisiae via a Novel Cdc55p Interaction Domain

A C-terminal region in Zds2p (ZH4) is required for regulation of Swe1p-dependent polarized cell growth and this region is necessary and sufficient for interaction with protein phosphatase 2A regulatory subunit, Cdc55p. Our results indicate that the Zds proteins regulate the Swe1p-dependent G2/M checkpoint in a CDC55-dependent manner