{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,10,12]],"date-time":"2025-10-12T04:18:05Z","timestamp":1760242685580,"version":"build-2065373602"},"reference-count":48,"publisher":"MDPI AG","issue":"1","license":[{"start":{"date-parts":[[2016,2,4]],"date-time":"2016-02-04T00:00:00Z","timestamp":1454544000000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"name":"Austrian Science Fund (FWF)","award":["W1225-B20"],"award-info":[{"award-number":["W1225-B20"]}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["Computation"],"abstract":"<jats:p>In population genetics, information about evolutionary forces, e.g., mutation, selection and genetic drift, is often inferred from DNA sequence information. Generally, DNA consists of two long strands of nucleotides or sites that pair via the complementary bases cytosine and guanine (C and G), on the one hand, and adenine and thymine (A and T), on the other. With whole genome sequencing, most genomic information stored in the DNA has become available for multiple individuals of one or more populations, at least in humans and model species, such as fruit flies of the genus Drosophila. In a genome-wide sample of L sites for M (haploid) individuals, the state of each site may be made binary, by binning the complementary bases, e.g., C with G to C\/G, and contrasting C\/G to A\/T, to obtain a \u201csite frequency spectrum\u201d (SFS). Two such samples of either a single population from different time-points or two related populations from a single time-point are called joint site frequency spectra (joint SFS). While mathematical models describing the interplay of mutation, drift and selection have been available for more than 80 years, calculation of exact likelihoods from joint SFS is difficult. Sufficient statistics for inference of, e.g., mutation or selection parameters that would make use of all the information in the genomic data are rarely available. Hence, often suites of crude summary statistics are combined in simulation-based computational approaches. In this article, we use a bi-allelic boundary-mutation and drift population genetic model to compute the transition probabilities of joint SFS using orthogonal polynomials. This allows inference of population genetic parameters, such as the mutation rate (scaled by the population size) and the time separating the two samples. We apply this inference method to a population dataset of neutrally-evolving short intronic sites from six DNA sequences of the fruit fly Drosophila melanogaster and the reference sequence of the related species Drosophila sechellia.<\/jats:p>","DOI":"10.3390\/computation4010006","type":"journal-article","created":{"date-parts":[[2016,2,4]],"date-time":"2016-02-04T10:26:22Z","timestamp":1454581582000},"page":"6","update-policy":"https:\/\/doi.org\/10.3390\/mdpi_crossmark_policy","source":"Crossref","is-referenced-by-count":4,"title":["Computation of the Likelihood of Joint Site Frequency Spectra Using Orthogonal Polynomials"],"prefix":"10.3390","volume":"4","author":[{"given":"Claus","family":"Vogl","sequence":"first","affiliation":[{"name":"Institute of Animal Breedings and Genetics, Veterin\u00e4rmedizinische Universit\u00e4t Wien, Veterin\u00e4rplatz 1, A-1210 Vienna, Austria"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-9415-2747","authenticated-orcid":false,"given":"Juraj","family":"Bergman","sequence":"additional","affiliation":[{"name":"Institut f\u00fcr Populationsgenetik, Veterin\u00e4rmedizinische Universit\u00e4t Wien, Veterin\u00e4rplatz 1, A-1210 Vienna, Austria"},{"name":"Vienna Graduate School of Population Genetics, Veterin\u00e4rmedizinische Universit\u00e4t Wien, Veterin\u00e4rplatz 1, A-1210 Vienna, Austria"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"1968","published-online":{"date-parts":[[2016,2,4]]},"reference":[{"key":"ref_1","doi-asserted-by":"crossref","unstructured":"Fisher, R. (1930). The Genetical Theory of Natural Selection, Clarendon Press.","DOI":"10.5962\/bhl.title.27468"},{"key":"ref_2","doi-asserted-by":"crossref","first-page":"97","DOI":"10.1093\/genetics\/16.2.97","article-title":"Evolution in Mendelian populations","volume":"16","author":"Wright","year":"1931","journal-title":"Genetics"},{"key":"ref_3","doi-asserted-by":"crossref","first-page":"19","DOI":"10.1016\/j.tpb.2014.10.002","article-title":"Estimating the scaled mutation rate and mutation bias with site frequency data","volume":"98","author":"Vogl","year":"2014","journal-title":"Theor. Popul. Biol."},{"key":"ref_4","doi-asserted-by":"crossref","first-page":"71","DOI":"10.1016\/j.tpb.2015.10.003","article-title":"Inference of directional selection and mutation parameters assuming equilibrium","volume":"106","author":"Vogl","year":"2015","journal-title":"Theor. Popul. Biol."},{"key":"ref_5","doi-asserted-by":"crossref","first-page":"144","DOI":"10.1073\/pnas.41.3.144","article-title":"Solution of a process of random genetic drift with a continuous model","volume":"41","author":"Kimura","year":"1955","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"ref_6","unstructured":"Griffiths, R., and Span\u00f2, D. (2010). Probability and Mathematical Genetics: Papers in Honour of Sir John Kingman, Cambridge University Press."},{"key":"ref_7","doi-asserted-by":"crossref","first-page":"1117","DOI":"10.1534\/genetics.111.136929","article-title":"A simple method for finding explicit analytic transition densities of diffusion processes with general diploid selection","volume":"190","author":"Song","year":"2012","journal-title":"Genetics"},{"key":"ref_8","doi-asserted-by":"crossref","first-page":"73","DOI":"10.1007\/s12064-012-0170-3","article-title":"An introduction to the mathematical structure of the Wright-Fisher model of population genetics","volume":"132","author":"Tran","year":"2013","journal-title":"Theory Biosci."},{"key":"ref_9","doi-asserted-by":"crossref","first-page":"199","DOI":"10.3390\/computation2040199","article-title":"Computation of the likelihood in biallelic diffusion models using orthogonal polynomials","volume":"2","author":"Vogl","year":"2014","journal-title":"Computation"},{"key":"ref_10","doi-asserted-by":"crossref","first-page":"197","DOI":"10.1016\/j.tpb.2012.01.001","article-title":"The allele-frequency spectrum in a decoupled Moran model with mutation, drift, and directional selection, assuming small mutation rates","volume":"81","author":"Vogl","year":"2012","journal-title":"Theor. Popul. Genet."},{"key":"ref_11","doi-asserted-by":"crossref","first-page":"1226","DOI":"10.1093\/molbev\/msq046","article-title":"On the utility of short intron sequences as a reference for the detection of positive and negative selection in Drosophila","volume":"27","author":"Parsch","year":"2010","journal-title":"Mol. Biol. Evol."},{"key":"ref_12","doi-asserted-by":"crossref","first-page":"1975","DOI":"10.1111\/j.1420-9101.2012.02580.x","article-title":"Unconstrained evolution in short introns?\u2014An analysis of genome-wide polymorphism and divergence data from Drosophila","volume":"25","author":"Clemente","year":"2012","journal-title":"J. Evol. Biol."},{"key":"ref_13","doi-asserted-by":"crossref","first-page":"2582","DOI":"10.1111\/jeb.12003","article-title":"Evidence for complex selection on four-fold degenerate sites in Drosophila melanogaster","volume":"25","author":"Clemente","year":"2012","journal-title":"J. Evol. Biol."},{"key":"ref_14","doi-asserted-by":"crossref","first-page":"1229","DOI":"10.1534\/genetics.115.174664","article-title":"The Drosophila Genome Nexus: A population genomic resource of 623 Drosophila melanogaster genomes, including 197 from a single ancestral range population","volume":"199","author":"Lack","year":"2015","journal-title":"Genetics"},{"key":"ref_15","unstructured":"NCBI Updates of Drosophila Annotations. Available online: http:\/\/www.flybase.org\/."},{"key":"ref_16","doi-asserted-by":"crossref","unstructured":"Carlin, B., and Louis, T. (2000). Bayes and Empirical Bayes Methods for Data Analysis, Chapman and Hall. [2nd ed.].","DOI":"10.1201\/9781420057669"},{"key":"ref_17","doi-asserted-by":"crossref","first-page":"118","DOI":"10.1016\/j.tpb.2010.05.003","article-title":"Sufficiency of the number of segregating sites in the limit under finite-sites mutation","volume":"78","author":"RoyChoudhury","year":"2010","journal-title":"Theor. Popul. Biol."},{"key":"ref_18","doi-asserted-by":"crossref","first-page":"143","DOI":"10.1016\/0040-5809(74)90020-3","article-title":"A note on the sampling theory for infinite alleles and infinite sites models","volume":"6","author":"Ewens","year":"1974","journal-title":"Theor. Popul. Biol."},{"key":"ref_19","doi-asserted-by":"crossref","first-page":"256","DOI":"10.1016\/0040-5809(75)90020-9","article-title":"On the number of segregating sites in genetical models without recombination","volume":"7","author":"Watterson","year":"1975","journal-title":"Theor. Popul. Biol."},{"key":"ref_20","doi-asserted-by":"crossref","first-page":"177","DOI":"10.2307\/3211856","article-title":"Diffusion models in population genetics","volume":"1","author":"Kimura","year":"1964","journal-title":"J. Appl. Probab."},{"key":"ref_21","doi-asserted-by":"crossref","first-page":"893","DOI":"10.1093\/genetics\/61.4.893","article-title":"The number of heterozygous nucleotide sites maintained in a finite population due to steady flux of mutations","volume":"61","author":"Kimura","year":"1969","journal-title":"Genetics"},{"key":"ref_22","doi-asserted-by":"crossref","unstructured":"Chan, A., Jenkins, P., and Song, Y. (2012). Genome-Wide Fine-Scale Recombination Rate Variation in Drosophila melanogaster. PLoS Genet., 8.","DOI":"10.1371\/journal.pgen.1003090"},{"key":"ref_23","doi-asserted-by":"crossref","first-page":"811","DOI":"10.1093\/molbev\/mss222","article-title":"Codon usage bias and effective population sizes on the X chromosome versus the autosomes in Drosophila melanogaster","volume":"30","author":"Campos","year":"2013","journal-title":"Mol. Biol. Evol."},{"key":"ref_24","doi-asserted-by":"crossref","first-page":"1010","DOI":"10.1093\/molbev\/msu056","article-title":"The relation between recombination rate and patterns of molecular evolution and variation in Drosophila melanogaster","volume":"31","author":"Campos","year":"2014","journal-title":"Mol. Biol. Evol."},{"key":"ref_25","doi-asserted-by":"crossref","first-page":"87","DOI":"10.1016\/0040-5809(72)90035-4","article-title":"The sampling theory of selectively neutral alleles","volume":"3","author":"Ewens","year":"1972","journal-title":"Theor. Popul. Biol."},{"key":"ref_26","doi-asserted-by":"crossref","first-page":"1161","DOI":"10.1093\/genetics\/132.4.1161","article-title":"Population genetics of polymorphism and divergence","volume":"132","author":"Sawyer","year":"1992","journal-title":"Genetics"},{"key":"ref_27","doi-asserted-by":"crossref","first-page":"1779","DOI":"10.1093\/genetics\/159.4.1779","article-title":"Directional selection and the site-frequency spectrum","volume":"159","author":"Bustamante","year":"2001","journal-title":"Genetics"},{"key":"ref_28","doi-asserted-by":"crossref","first-page":"91","DOI":"10.1016\/S0040-5809(02)00050-3","article-title":"Maximum likelihood and Bayesian methods for estimating the distribution of selective effects among classes of mutations using DNA polymorphism data","volume":"63","author":"Bustamante","year":"2003","journal-title":"Theor. Popul. Biol."},{"key":"ref_29","doi-asserted-by":"crossref","first-page":"463","DOI":"10.1534\/genetics.103.024745","article-title":"Population genetics of polymorphism and divergence for diploid selection models with arbitrary dominance","volume":"168","author":"Williamson","year":"2004","journal-title":"Genetics"},{"key":"ref_30","unstructured":"Ewens, W. (1979). Mathematical Population Genetics, Springer."},{"key":"ref_31","unstructured":"Abramowitz, M., and Stegun, I. (1970). Handbook of Mathematical Functions, Dover. [9th ed.]."},{"key":"ref_32","doi-asserted-by":"crossref","first-page":"973","DOI":"10.1534\/genetics.113.152017","article-title":"Complete numerical solution of the diffusion equation of random genetic drift","volume":"194","author":"Zhao","year":"2013","journal-title":"Genetics"},{"key":"ref_33","doi-asserted-by":"crossref","unstructured":"Ewens, W. (2004). Mathematical Population Genetics, Springer. [2nd ed.].","DOI":"10.1007\/978-0-387-21822-9"},{"key":"ref_34","doi-asserted-by":"crossref","unstructured":"Gelman, A., Carlin, J., Stern, H., and Rubin, D. (1995). Bayesian Data Analysis, Chapman & Hall.","DOI":"10.1201\/9780429258411"},{"key":"ref_35","doi-asserted-by":"crossref","first-page":"159","DOI":"10.1007\/978-1-4613-0931-4_4","article-title":"Historical biogeography of the Drosophila melanogaster species subgroup","volume":"22","author":"Lachaise","year":"1988","journal-title":"Evol. Biol."},{"key":"ref_36","first-page":"391","article-title":"Molecular phylogeny and divergence times of Drosophilid species","volume":"12","author":"Russo","year":"1995","journal-title":"Mol. Biol. Evol."},{"key":"ref_37","doi-asserted-by":"crossref","first-page":"778","DOI":"10.1093\/molbev\/msn024","article-title":"Divergence times in Caenorhabditis and Drosophila inferred from direct estimates of the neutral mutation rate","volume":"25","author":"Cutter","year":"2008","journal-title":"Mol. Biol. Evol."},{"key":"ref_38","doi-asserted-by":"crossref","first-page":"768","DOI":"10.1093\/bioinformatics\/btk051","article-title":"LAMARC 2.0: Maximum likelihood and Bayesian estimation of population parameters","volume":"15","author":"Kuhner","year":"2006","journal-title":"Bioinformatics"},{"key":"ref_39","doi-asserted-by":"crossref","unstructured":"Gutenkunst, R., Hernandez, R., Williamson, S., and Bustamante, C. (2009). Inferring the Joint Demographic History of Multiple Populations from Multidimensional SNP Frequency Data. PLoS Genet., 5.","DOI":"10.1371\/journal.pgen.1000695"},{"key":"ref_40","doi-asserted-by":"crossref","first-page":"109","DOI":"10.1016\/j.tpb.2006.06.005","article-title":"Non-equilibrium theory of the allele frequency spectrum","volume":"71","author":"Evans","year":"2007","journal-title":"Theor. Popul. Biol."},{"key":"ref_41","doi-asserted-by":"crossref","first-page":"601","DOI":"10.1534\/genetics.115.175265","article-title":"Transition densities and sample frequency spectra of diffusion processes with selection and variable population size","volume":"200","author":"Zivkovic","year":"2015","journal-title":"Genetics"},{"key":"ref_42","doi-asserted-by":"crossref","unstructured":"Hein, J., Schierup, M., and Wiuf, C. (2005). Gene Genealogies, Variation, and Evolution: A Primer in Coalescent Theory, Oxford University Press.","DOI":"10.1093\/oso\/9780198529958.001.0001"},{"key":"ref_43","unstructured":"Wakeley, J. (2009). Coalescent Theory: An Introduction, Roberts and Co."},{"key":"ref_44","doi-asserted-by":"crossref","first-page":"179","DOI":"10.1016\/j.tpb.2011.11.004","article-title":"The joint allele frequency spectrum of multiple populations: A coalescent theory approach","volume":"81","author":"Chen","year":"2012","journal-title":"Theor. Popul. Biol."},{"key":"ref_45","doi-asserted-by":"crossref","first-page":"158","DOI":"10.1111\/ahg.12007","article-title":"Intercoalescence time distribution of incomplete gene genealogies in temporally varying populations and applications in population genetic inference","volume":"77","author":"Chen","year":"2013","journal-title":"Ann. Hum. Genet."},{"key":"ref_46","unstructured":"Kamm, J., Terhorst, J., and Song, Y. Efficient computation of the joint sample frequency spectra for multiple populations. Available online: http:\/\/arxiv.org\/abs\/1503.01133."},{"key":"ref_47","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1016\/j.tpb.2012.10.006","article-title":"An explicit transition density expansion for a multi-allelic Wright-Fisher diffusion with general diploid selection","volume":"83","author":"Wang","year":"2013","journal-title":"Theor. Popul. Biol."},{"key":"ref_48","first-page":"2203","article-title":"A novel method for inferring general diploid selection from time series genetic data","volume":"8","author":"Bhaskar","year":"2014","journal-title":"Ann. Appl. Stat."}],"container-title":["Computation"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.mdpi.com\/2079-3197\/4\/1\/6\/pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,10,11]],"date-time":"2025-10-11T19:18:49Z","timestamp":1760210329000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.mdpi.com\/2079-3197\/4\/1\/6"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2016,2,4]]},"references-count":48,"journal-issue":{"issue":"1","published-online":{"date-parts":[[2016,3]]}},"alternative-id":["computation4010006"],"URL":"https:\/\/doi.org\/10.3390\/computation4010006","relation":{},"ISSN":["2079-3197"],"issn-type":[{"type":"electronic","value":"2079-3197"}],"subject":[],"published":{"date-parts":[[2016,2,4]]}}}