The Community for Technology Leaders
Green Image
Issue No. 05 - Sept.-Oct. (2013 vol. 10)
ISSN: 1545-5963
pp: 1
Yujin Chung , Dept. of Stat., Univ. of Wisconsin, Madison, WI, USA
Nicole T. Perna , Dept. of Genetics, Univ. of Wisconsin, Madison, WI, USA
Cecile Ane , Dept. of Stat., Univ. of Wisconsin, Madison, WI, USA
Ancestral recombination events can cause the underlying genealogy of a site to vary along the genome. We consider Bayesian models to simultaneously detect recombination breakpoints in very long sequence alignments and estimate the phylogenetic tree of each block between breakpoints. The models we consider use a dissimilarity measure between trees in their prior distribution to favor similar trees at neighboring loci. We show empirical evidence in Enterobacteria that neighboring genomic regions have similar trees. The main hurdle to using such models is the need to properly calculate the normalizing function for the prior probabilities on trees. In this work, we quantify the impact of approximating this normalizing function as done in biomc2, a hierarchical Bayesian method to detect recombination based on distance between tree topologies. We then derive an algorithm to calculate the normalizing function exactly, for a Gibbs distribution based on the Robinson-Foulds (RF) distance between gene trees at neighboring loci. At the core is the calculation of the joint distribution of the shape of a random tree and its RF distance to a fixed tree. We also propose fast approximations to the normalizing function, which are shown to be very accurate with little impact on the Bayesian inference.
Bayes methods, Bioinformatics, Genomics, Phylogeny, Radio frequency, Biological system modeling
Yujin Chung, Nicole T. Perna, Cecile Ane, "Computing the Joint Distribution of Tree Shape and Tree Distance for Gene Tree Inference and Recombination Detection", IEEE/ACM Transactions on Computational Biology and Bioinformatics, vol. 10, no. , pp. 1, Sept.-Oct. 2013, doi:10.1109/TCBB.2013.109
165 ms
(Ver 3.3 (11022016))