A stochastic Farris transform for genetic data under the multispecies coalescent with applications to data requirements
Name
285_2022_1731_ReferencePDF.pdf
Size
770.94 KB
Format
Unknown
Checksum (MD5)
8b0a11a179fb6f2bc6366b103c64f958
Author(s)
Dasarathy, Gautam
Mossel, Elchanan
Nowak, Robert
Roch, Sebastien
Date Issued
April 2022
Journal
Journal of Mathematical Biology
Publisher
Springer Berlin Heidelberg
Citation
Journal of Mathematical Biology. 2022 Apr 08;84(5):36
Version
Author's final manuscript
Abstract
Abstract
Species tree estimation faces many significant hurdles. Chief among them is that the trees describing the ancestral lineages of each individual gene—the gene trees—often differ from the species tree. The multispecies coalescent is commonly used to model this gene tree discordance, at least when it is believed to arise from incomplete lineage sorting, a population-genetic effect. Another significant challenge in this area is that molecular sequences associated to each gene typically provide limited information about the gene trees themselves. While the modeling of sequence evolution by single-site substitutions is well-studied, few species tree reconstruction methods with theoretical guarantees actually address this latter issue. Instead, a standard—but unsatisfactory—assumption is that gene trees are perfectly reconstructed before being fed into a so-called summary method. Hence much remains to be done in the development of inference methodologies that rigorously account for gene tree estimation error—or completely avoid gene tree estimation in the first place. In previous work, a data requirement trade-off was derived between the number of loci m needed for an accurate reconstruction and the length of the locus sequences k. It was shown that to reconstruct an internal branch of length f, one needs m to be of the order of
$$1/[f^{2} \sqrt{k}]$$
1
/
[
f
2
k
]
. That previous result was obtained under the restrictive assumption that mutation rates as well as population sizes are constant across the species phylogeny. Here we further generalize this result beyond this assumption. Our main contribution is a novel reduction to the molecular clock case under the multispecies coalescent, which we refer to as a stochastic Farris transform. As a corollary, we also obtain a new identifiability result of independent interest: for any species tree with
$$n \ge 3$$
n
≥
3
species, the rooted topology of the species tree can be identified from the distribution of its unrooted weighted gene trees even in the absence of a molecular clock.
MIT Department
Massachusetts Institute of Technology. Department of Mathematics
Massachusetts Institute of Technology. Institute for Data, Systems, and Society
Terms of Use
Creative Commons Attribution-Noncommercial-Share Alike
Persistent DSpace Link
DOI of Published Version
https://doi.org/10.1007/s00285-022-01731-5
