Bayesian estimation of the number of inversions in the history of two chromosomes

Thomas L. York, Richard Durrett, Rasmus Nielsen

Research output: Contribution to journalJournal articleResearchpeer-review

20 Citations (Scopus)

Abstract

We present a Bayesian approach to the problem of inferring the history of inversions separating homologous chromosomes from two different species. The method is based on Markov Chain Monte Carlo (MCMC) and takes full advantage of all the information from marker order. We apply the method both to simulated data and to two real data sets. For the simulated data, we show that the MCMC method provides accurate estimates of the true posterior distributions and in the analysis of the real data we show that the most likely number of inversions in some cases is considerably larger than estimates obtained based on the parsimony inferred number of inversions. Indeed, in the case of the Drosophila repleta-D. melanogaster comparison, the lower boundary of a 95% highest posterior density credible interval for the number of inversions is considerably larger than the most parsimonious number of inversions.

Original languageEnglish
JournalJournal of Computational Biology
Volume9
Issue number6
Pages (from-to)805-818
Number of pages14
ISSN1066-5277
DOIs
Publication statusPublished - 1 Dec 2002
Externally publishedYes

Keywords

  • Bayesian estimation
  • Breakpoint graph
  • Inversions
  • MCMC

Cite this