Cellular noise

Last updated

Cellular noise is random variability in quantities arising in cellular biology. For example, cells which are genetically identical, even within the same tissue, are often observed to have different expression levels of proteins, different sizes and structures. [1] [2] These apparently random differences can have important biological and medical consequences. [3]

Contents

Cellular noise was originally, and is still often, examined in the context of gene expression levels – either the concentration or copy number of the products of genes within and between cells. As gene expression levels are responsible for many fundamental properties in cellular biology, including cells' physical appearance, behaviour in response to stimuli, and ability to process information and control internal processes, the presence of noise in gene expression has profound implications for many processes in cellular biology.

Definitions

The most frequent quantitative definition of noise is the coefficient of variation:[ citation needed ]

where is the noise in a quantity , is the mean value of and is the standard deviation of . This measure is dimensionless, allowing a relative comparison of the importance of noise, without necessitating knowledge of the absolute mean.

Other quantities often used for mathematical convenience are the Fano factor:

and the normalized variance:

Experimental measurement

The first experimental account and analysis of gene expression noise in prokaryotes is from Becskei & Serrano [4] and from Alexander van Oudenaarden's lab. [5] The first experimental account and analysis of gene expression noise in eukaryotes is from James J. Collins's lab. [6]

Intrinsic and extrinsic noise

A schematic illustration of a dual reporter study. Each data point corresponds to a measurement of the expression level of two identically regulated genes in a single cell: the scatter reflects measurements of a population of cells. Extrinsic noise is characterised by expression levels of both genes covarying between cells, intrinsic by internal differences. Extrinsic and intrinsic noise in cellular biology.svg
A schematic illustration of a dual reporter study. Each data point corresponds to a measurement of the expression level of two identically regulated genes in a single cell: the scatter reflects measurements of a population of cells. Extrinsic noise is characterised by expression levels of both genes covarying between cells, intrinsic by internal differences.

Cellular noise is often investigated in the framework of intrinsic and extrinsic noise. Intrinsic noise refers to variation in identically regulated quantities within a single cell: for example, the intra-cell variation in expression levels of two identically controlled genes. Extrinsic noise refers to variation in identically regulated quantities between different cells: for example, the cell-to-cell variation in expression of a given gene.

Intrinsic and extrinsic noise levels are often compared in dual reporter studies, in which the expression levels of two identically regulated genes (often fluorescent reporters like GFP and YFP) are plotted for each cell in a population. [7]

An issue with the general depiction of extrinsic noise as a spread along the main diagonal in dual-reporter studies is the assumption that extrinsic factors cause positive expression correlations between the two reporters. In fact, when the two reporters compete for binding of a low-copy regulator, the two reporters become anomalously anticorrelated, and the spread is perpendicular to the main diagonal. In fact, any deviation of the dual-reporter scatter plot from circular symmetry indicates extrinsic noise. Information theory offers a way to avoid this anomaly. [8]

Sources

Note: These lists are illustrative, not exhaustive, and identification of noise sources is an active and expanding area of research.

Intrinsic noise
Extrinsic noise

Note that extrinsic noise can affect levels and types of intrinsic noise: [19] for example, extrinsic differences in the mitochondrial content of cells lead, through differences in ATP levels, to some cells transcribing faster than others, affecting the rates of gene expression and the magnitude of intrinsic noise across the population. [17]

Effects

Note: These lists are illustrative, not exhaustive, and identification of noise effects is an active and expanding area of research.

Analysis

A canonical model for stochastic gene expression, known as the two-state or telegraph model. DNA flips between "inactive" and "active" states (involving, for example, chromatin remodelling and transcription factor binding). Active DNA is transcribed to produce mRNA which is translated to produce protein, both of which are degraded. All processes are Poissonian with given rates. Modelling stochastic gene expression.svg
A canonical model for stochastic gene expression, known as the two-state or telegraph model. DNA flips between "inactive" and "active" states (involving, for example, chromatin remodelling and transcription factor binding). Active DNA is transcribed to produce mRNA which is translated to produce protein, both of which are degraded. All processes are Poissonian with given rates.

As many quantities of cell biological interest are present in discrete copy number within the cell (single DNAs, dozens of mRNAs, hundreds of proteins), tools from discrete stochastic mathematics are often used to analyse and model cellular noise. [31] [32] In particular, master equation treatments – where the probabilities of observing a system in a state at time are linked through ODEs – have proved particularly fruitful. A canonical model for noise gene expression, where the processes of DNA activation, transcription and translation are all represented as Poisson processes with given rates, gives a master equation which may be solved exactly (with generating functions) under various assumptions or approximated with stochastic tools like Van Kampen's system size expansion.

Numerically, the Gillespie algorithm or stochastic simulation algorithm is often used to create realisations of stochastic cellular processes, from which statistics can be calculated.

The problem of inferring the values of parameters in stochastic models (parametric inference) for biological processes, which are typically characterised by sparse and noisy experimental data, is an active field of research, with methods including Bayesian MCMC and approximate Bayesian computation proving adaptable and robust. [33] Regarding the two-state model, a moment-based method was described for parameters inference from mRNAs distributions. [30]

Related Research Articles

<span class="mw-page-title-main">Apoptosis</span> Programmed cell death in multicellular organisms

Apoptosis is a form of programmed cell death that occurs in multicellular organisms and in some eukaryotic, single-celled microorganisms such as yeast. Biochemical events lead to characteristic cell changes (morphology) and death. These changes include blebbing, cell shrinkage, nuclear fragmentation, chromatin condensation, DNA fragmentation, and mRNA decay. The average adult human loses between 50 and 70 billion cells each day due to apoptosis. For an average human child between eight and fourteen years old, approximately twenty to thirty billion cells die per day.

<span class="mw-page-title-main">Bistability</span> Quality of a system having two stable equilibrium states

In a dynamical system, bistability means the system has two stable equilibrium states. A bistable structure can be resting in either of two states. An example of a mechanical device which is bistable is a light switch. The switch lever is designed to rest in the "on" or "off" position, but not between the two. Bistable behavior can occur in mechanical linkages, electronic circuits, nonlinear optical systems, chemical reactions, and physiological and biological systems.

<span class="mw-page-title-main">Cellular differentiation</span> Developmental biology

Cellular differentiation is the process in which a stem cell changes from one type to a differentiated one. Usually, the cell changes to a more specialized type. Differentiation happens multiple times during the development of a multicellular organism as it changes from a simple zygote to a complex system of tissues and cell types. Differentiation continues in adulthood as adult stem cells divide and create fully differentiated daughter cells during tissue repair and during normal cell turnover. Some differentiation occurs in response to antigen exposure. Differentiation dramatically changes a cell's size, shape, membrane potential, metabolic activity, and responsiveness to signals. These changes are largely due to highly controlled modifications in gene expression and are the study of epigenetics. With a few exceptions, cellular differentiation almost never involves a change in the DNA sequence itself. However, metabolic composition does get altered quite dramatically where stem cells are characterized by abundant metabolites with highly unsaturated structures whose levels decrease upon differentiation. Thus, different cells can have very different physical characteristics despite having the same genome.

<span class="mw-page-title-main">Gene expression</span> Conversion of a genes sequence into a mature gene product or products

Gene expression is the process by which information from a gene is used in the synthesis of a functional gene product that enables it to produce end products, proteins or non-coding RNA, and ultimately affect a phenotype. These products are often proteins, but in non-protein-coding genes such as transfer RNA (tRNA) and small nuclear RNA (snRNA), the product is a functional non-coding RNA. Gene expression is summarized in the central dogma of molecular biology first formulated by Francis Crick in 1958, further developed in his 1970 article, and expanded by the subsequent discoveries of reverse transcription and RNA replication.

<span class="mw-page-title-main">Gene regulatory network</span> Collection of molecular regulators

A generegulatory network (GRN) is a collection of molecular regulators that interact with each other and with other substances in the cell to govern the gene expression levels of mRNA and proteins which, in turn, determine the function of the cell. GRN also play a central role in morphogenesis, the creation of body structures, which in turn is central to evolutionary developmental biology (evo-devo).

A circadian clock, or circadian oscillator, is a biochemical oscillator that cycles with a stable phase and is synchronized with solar time.

Intrinsic immunity refers to a set of cellular-based anti-viral defense mechanisms, notably genetically encoded proteins which specifically target eukaryotic retroviruses. Unlike adaptive and innate immunity effectors, intrinsic immune proteins are usually expressed at a constant level, allowing a viral infection to be halted quickly. Intrinsic antiviral immunity refers to a form of innate immunity that directly restricts viral replication and assembly, thereby rendering a cell non-permissive to a specific class or species of viruses. Intrinsic immunity is conferred by restriction factors preexisting in certain cell types, although these factors can be further induced by virus infection. Intrinsic viral restriction factors recognize specific viral components, but unlike other pattern recognition receptors that inhibit viral infection indirectly by inducing interferons and other antiviral molecules, intrinsic antiviral factors block viral replication immediately and directly.

Developmental noise or stochastic noise is a concept within developmental biology in which the observable characteristics or traits (phenotype) varies between individuals even though both individuals share the same genetic code (genotypes) and the other environmental factors are completely the same. Factors that influence the effect include stochastic, or randomized, gene expression and other cellular noise.

<span class="mw-page-title-main">Alexander van Oudenaarden</span> Dutch biophysicist and systems biologist

Alexander van Oudenaarden is a Dutch biophysicist and systems biologist. He is a leading researcher in stem cell biology, specialising in single cell techniques. In 2012 he started as director of the Hubrecht Institute and was awarded three times an ERC Advanced Grant, in 2012, 2017, and 2022. He was awarded the Spinoza Prize in 2017.

Epigenomics is the study of the complete set of epigenetic modifications on the genetic material of a cell, known as the epigenome. The field is analogous to genomics and proteomics, which are the study of the genome and proteome of a cell. Epigenetic modifications are reversible modifications on a cell's DNA or histones that affect gene expression without altering the DNA sequence. Epigenomic maintenance is a continuous process and plays an important role in stability of eukaryotic genomes by taking part in crucial biological mechanisms like DNA repair. Plant flavones are said to be inhibiting epigenomic marks that cause cancers. Two of the most characterized epigenetic modifications are DNA methylation and histone modification. Epigenetic modifications play an important role in gene expression and regulation, and are involved in numerous cellular processes such as in differentiation/development and tumorigenesis. The study of epigenetics on a global level has been made possible only recently through the adaptation of genomic high-throughput assays.

Transcriptional bursting, also known as transcriptional pulsing, is a fundamental property of genes in which transcription from DNA to RNA can occur in "bursts" or "pulses", which has been observed in diverse organisms, from bacteria to mammals.

Transcriptional noise is a primary cause of the variability (noise) in gene expression occurring between cells in isogenic populations. A proposed source of transcriptional noise is transcriptional bursting although other sources of heterogeneity, such as unequal separation of cell contents at mitosis are also likely to contribute considerably. Bursting transcription, as opposed to simple probabilistic models of transcription, reflects multiple states of gene activity, with fluctuations between states separated by irregular intervals, generating uneven protein expression between cells. Noise in gene expression can have tremendous consequences on cell behaviour, and must be mitigated or integrated. In certain contexts, such as establishment of viral latency, the survival of microbes in rapidly changing stressful environments, or several types of scattered differentiation, the variability may be essential. Variability also impacts upon the effectiveness of clinical treatment, with resistance of bacteria and yeast to antibiotics demonstrably caused by non-genetic differences. Variability in gene expression may also contribute to resistance of sub-populations of cancer cells to chemotherapy and appears to be a barrier to curing HIV.

<span class="mw-page-title-main">Robustness (evolution)</span> Persistence of a biological trait under uncertain conditions

In evolutionary biology, robustness of a biological system is the persistence of a certain characteristic or trait in a system under perturbations or conditions of uncertainty. Robustness in development is known as canalization. According to the kind of perturbation involved, robustness can be classified as mutational, environmental, recombinational, or behavioral robustness etc. Robustness is achieved through the combination of many genetic and molecular mechanisms and can evolve by either direct or indirect selection. Several model systems have been developed to experimentally study robustness and its evolutionary consequences.

Johan Paulsson is a Swedish mathematician and systems biologist at Harvard Medical School. He is a leading researcher in systems biology and stochastic processes, specializing in stochasticity in gene networks and plasmid reproduction.

Within bioinformatics, intrinsic Noise Analyzer (iNA) is an open source software for studying reaction kinetics in living cells. The software analyzes mathematical models of intracellular reaction kinetics such as gene expression, regulatory networks or signaling pathways to quantify concentration fluctuations due to the random nature of chemical reactions.

Single-cell sequencing examines the nucleic acid sequence information from individual cells with optimized next-generation sequencing technologies, providing a higher resolution of cellular differences and a better understanding of the function of an individual cell in the context of its microenvironment. For example, in cancer, sequencing the DNA of individual cells can give information about mutations carried by small populations of cells. In development, sequencing the RNAs expressed by individual cells can give insight into the existence and behavior of different cell types. In microbial systems, a population of the same species can appear genetically clonal. Still, single-cell sequencing of RNA or epigenetic modifications can reveal cell-to-cell variability that may help populations rapidly adapt to survive in changing environments.

<span class="mw-page-title-main">Single-cell variability</span>

In cell biology, single-cell variability occurs when individual cells in an otherwise similar population differ in shape, size, position in the cell cycle, or molecular-level characteristics. Such differences can be detected using modern single-cell analysis techniques. Investigation of variability within a population of cells contributes to understanding of developmental and pathological processes,

State switching is a fundamental physiological process in which a cell/organism undergoes spontaneous, and potentially reversible, transitions between different phenotypes. Thus, the ability to switch states/phenotypes is a key feature of development and normal function of cells within most multicellular organisms that enables the cell to respond to various intrinsic and extrinsic cues and stimuli in a concerted fashion enabling them to ‘make’ appropriate cellular decisions. Although state switching is essential for normal functioning, the repertoire of phenotypes in a normal cell is albeit limited.

Single-cell transcriptomics examines the gene expression level of individual cells in a given population by simultaneously measuring the RNA concentration of hundreds to thousands of genes. Single-cell transcriptomics makes it possible to unravel heterogeneous cell populations, reconstruct cellular developmental pathways, and model transcriptional dynamics — all previously masked in bulk RNA sequencing.

Monoallelic gene expression (MAE) is the phenomenon of the gene expression, when only one of the two gene copies (alleles) is actively expressed (transcribed), while the other is silent. Diploid organisms bear two homologous copies of each chromosome (one from each parent), a gene can be expressed from both chromosomes (biallelic expression) or from only one (monoallelic expression). MAE can be Random monoallelic expression (RME) or Constitutive monoallelic expression (constitutive). Constitutive monoallelic expression occurs from the same specific allele throughout the whole organism or tissue, as a result of genomic imprinting. RME is a broader class of monoallelic expression, which is defined by random allelic choice in somatic cells, so that different cells of the multi-cellular organism express different alleles.

References

  1. Kaern, M.; Elston, T.R.; Blake, W.J. & Collins, J.J. (2005). "Stochasticity in gene expression: from theories to phenotypes". Nat. Rev. Genet. 6 (6): 451–464. doi:10.1038/nrg1615. PMID   15883588. S2CID   1028111.
  2. Maheshri N, O'Shea EK (2007). "Living with noisy genes: how cells function reliably with inherent variability in gene expression". Annu. Rev. Biophys. Biomol. Struct. 36: 413–434. doi:10.1146/annurev.biophys.36.040306.132705. PMID   17477840.
  3. Johnston, I. G. (2012). "The chaos within: exploring noise in cellular biology". Significance. 19 (4): 17–21. arXiv: 1208.2250 . Bibcode:2012arXiv1208.2250J. doi:10.1111/j.1740-9713.2012.00586.x. S2CID   16368991.
  4. Becskei, Attila; Serrano, Luis (2000). "Engineering Stability in Gene Networks by Autoregulation" (PDF). Nature. 405 (6786): 590–593. Bibcode:2000Natur.405..590B. doi:10.1038/35014651. PMID   10850721. S2CID   4407358.
  5. Ozbudak, Ertugrul M; Thattai, Mukund; Kurtser, Iren; Grossman, Alan D; van Oudenaarden, Alexander (2002). "Regulation of Noise in the Expression of a Single Gene". Nature Genetics. 31 (1): 69–73. doi: 10.1038/ng869 . PMID   11967532.
  6. Blake, William J; Kærn, Mads; Cantor, Charles R; Collins, James J (2003). "Noise in Eukaryotic Gene Expression". Nature. 422 (6932): 633–637. Bibcode:2003Natur.422..633B. doi:10.1038/nature01546. PMID   12687005. S2CID   4347106.
  7. Elowitz, M.B.; Levine, A.J.; Siggia, E.D.; Swain, P.S. (2002). "Stochastic gene expression in a single cell". Science. 297 (5584): 1183–6. Bibcode:2002Sci...297.1183E. doi:10.1126/science.1070919. PMID   12183631. S2CID   10845628.
  8. 1 2 Stamatakis, Michail; Adams, Rhys M; Balázsi, Gábor (2011). "A Common Repressor Pool Results in Indeterminacy of Extrinsic Noise". Chaos. 21 (4): 047523–047523–12. Bibcode:2011Chaos..21d7523S. doi:10.1063/1.3658618. PMC   3258287 . PMID   22225397.
  9. Morelli, M.J.; Allen, R.J. & ten Wolde, P.R. (2011). "Effects of macromolecular crowding on genetic networks". Biophys. J. 101 (12): 2882–2891. Bibcode:2011BpJ...101.2882M. doi:10.1016/j.bpj.2011.10.053. PMC   3244068 . PMID   22208186.
  10. Komorowski M, Miekisz J, Stumpf M PH (2013). "Decomposing Noise in Biochemical Signaling Systems Highlights the Role of Protein Degradation". Biophys. J. 104 (8): 1783–1793. Bibcode:2013BpJ...104.1783K. doi:10.1016/j.bpj.2013.02.027. PMC   3627874 . PMID   23601325.
  11. Jetka T, Charzynska A, Gambin A, Stumpf M PH, Komorowski M (2013). "StochDecomp - Matlab package for noise decomposition in stochastic biochemical systems". Bioinformatics. 30 (1): 137–138. arXiv: 1308.3103 . Bibcode:2013arXiv1308.3103J. doi:10.1093/bioinformatics/btt631. PMID   24191070.
  12. Thomas, Philipp (2019-01-24). "Intrinsic and extrinsic noise of gene expression in lineage trees". Scientific Reports. 9 (1): 474. Bibcode:2019NatSR...9..474T. doi:10.1038/s41598-018-35927-x. ISSN   2045-2322. PMC   6345792 . PMID   30679440.
  13. Newman JR, Ghaemmaghami S, Ihmels J, Breslow DK, Noble M, DeRisi JL, Weissman JS (2006). "Single-cell proteomic analysis of S. cerevisiae reveals the architecture of biological noise". Nature. 441 (7095): 840–846. Bibcode:2006Natur.441..840N. doi:10.1038/nature04785. PMID   16699522. S2CID   4300851.
  14. Weiße, Andrea Y.; Vincent Danos; Terradot, Guillaume; Thomas, Philipp (2018-10-30). "Sources, propagation and consequences of stochasticity in cellular growth". Nature Communications. 9 (1): 4528. Bibcode:2018NatCo...9.4528T. doi:10.1038/s41467-018-06912-9. ISSN   2041-1723. PMC   6207721 . PMID   30375377.
  15. Athale, C.A.; Chaudhari, H. (2011). "Population length variability and nucleoid numbers in Escherichia coli". Bioinformatics. 27 (21): 2944–2998. doi:10.1093/bioinformatics/btr501. PMID   21930671.
  16. das Neves RP, Jones NS, Andreu L, Gupta R, Enver T, Iborra FJ (2010). "Connecting Variability in Global Transcription Rate to Mitochondrial Variability". PLOS Biol. 8 (12): e1000560. doi: 10.1371/journal.pbio.1000560 . PMC   3001896 . PMID   21179497.
  17. 1 2 3 4 Johnston IG, Gaal B, das Neves RP, Enver T, Iborra FJ, Jones NS (2012). "Mitochondrial Variability as a Source of Extrinsic Cellular Noise". PLOS Comput. Biol. 8 (3): e1002416. arXiv: 1107.4499 . Bibcode:2012PLSCB...8E2416J. doi: 10.1371/journal.pcbi.1002416 . PMC   3297557 . PMID   22412363.
  18. Huh, D.; Paulsson, J. (2011). "Random partitioning of molecules at cell division". Proc. Natl. Acad. Sci. USA. 108 (36): 15004–15009. Bibcode:2011PNAS..10815004H. doi: 10.1073/pnas.1013171108 . PMC   3169110 . PMID   21873252.
  19. Shahrezaei, V. & Swain, P.S. (2008). "Analytical distributions for stochastic gene expression". Proc. Natl. Acad. Sci. USA. 105 (45): 17256–17261. arXiv: 0812.3344 . Bibcode:2008PNAS..10517256S. doi: 10.1073/pnas.0803850105 . PMC   2582303 . PMID   18988743.
  20. 1 2 Lestas, I.; Vinnicombe, G.; Paulsson, J. (2010). "Fundamental limits on the suppression of molecular fluctuations". Nature. 467 (7312): 174–8. Bibcode:2010Natur.467..174L. doi:10.1038/nature09333. PMC   2996232 . PMID   20829788.
  21. Fraser D, Kaern M (2009). "A chance at survival: gene expression noise and phenotypic diversification strategies". Mol. Microbiol. 71 (6): 1333–1340. doi: 10.1111/j.1365-2958.2009.06605.x . PMID   19220745.
  22. T. Enver; C.M. Heyworth & T.M. Dexter (1998). "Do Stem Cells Play Dice?". Blood. 92 (2): 348–51, discussion 352. doi:10.1182/blood.V92.2.348. PMID   9657728.
  23. Balázsi, Gábor; van Oudenaarden, Alexander; Collins, James J (2011). "Cellular Decision Making and Biological Noise: From Microbes to Mammals". Cell. 144 (6): 910–925. doi:10.1016/j.cell.2011.01.030. PMC   3068611 . PMID   21414483.
  24. Blake, William J; Balázsi, Gábor; Kohanski, Michael A; Isaacs, Farren J; Murphy, Kevin F; Kuang, Yina; Cantor, Charles R; Walt, David R; Collins, James J (2006). "Phenotypic Consequences of Promoter-Mediated Transcriptional Noise". Molecular Cell. 24 (6): 853–865. doi: 10.1016/j.molcel.2006.11.003 . PMID   17189188.
  25. Farquhar, Kevin F; Charlebois, Daniel A; Szenk, Mariola; Cohen, Joseph; Nevozhay, Dmitry; Balázsi, Gábor (2019). "Role of Network-Mediated Stochasticity in Mammalian Drug Resistance". Nature Communications. 10 (1): 2766. doi: 10.1038/s41467-019-10330-w . PMC   6591227 . PMID   31235692.
  26. Spencer, S.L.; Gaudet, S.; Albeck, J.G.; Burke, J.M.; Sorger, P.K. (2009). "Non-genetic origins of cell-to-cell variability in TRAIL-induced apoptosis". Nature. 459 (7245): 428–432. Bibcode:2009Natur.459..428S. doi:10.1038/nature08012. PMC   2858974 . PMID   19363473.
  27. Liao D, Estévez-Salmerón L, Tlsty TD (2012). "Conceptualizing a tool to optimize therapy based on dynamic heterogeneity". Phys. Biol. 9 (6): 065005. Bibcode:2012PhBio...9f5005L. doi:10.1088/1478-3975/9/6/065005. PMC   3618714 . PMID   23197078.
  28. Batada NN, Hurst LD (2007). "Evolution of chromosome organization driven by selection for reduced gene expression noise". Nature Genetics. 39 (8): 945–9. doi:10.1038/ng2071. PMID   17660811. S2CID   19546863.
  29. Perkins, T.J. & Swain, P.S. (2009). "Strategies for cellular decision-making". Mol. Syst. Biol. 5 (236): 326. doi:10.1038/msb.2009.83. PMC   2795477 . PMID   19920811.
  30. 1 2 Peccoud, J. & Ycart, B. (1995). "Markovian Modelling of Gene Product Synthesis". Theoretical Population Biology. 48 (2): 222–234. doi:10.1006/tpbi.1995.1027.
  31. Paulsson, J. (2005). "Models of stochastic gene expression". Phys. Life Rev. 2 (2): 157–175. Bibcode:2005PhLRv...2..157P. doi:10.1016/j.plrev.2005.03.003.
  32. Wilkinson, D.J. (2009). "Stochastic modelling for quantitative description of heterogeneous biological systems". Nat. Rev. Genet. 10 (2): 122–133. doi:10.1038/nrg2509. PMID   19139763. S2CID   14731499.
  33. Sunnåker, Mikael; Busetto, Alberto Giovanni; Numminen, Elina; Corander, Jukka; Foll, Matthieu; Dessimoz, Christophe (2013). "Approximate Bayesian Computation". PLOS Computational Biology. 9 (1): e1002803. Bibcode:2013PLSCB...9E2803S. doi: 10.1371/journal.pcbi.1002803 . PMC   3547661 . PMID   23341757.