期刊
MOLECULAR ECOLOGY RESOURCES
卷 21, 期 4, 页码 1359-1368出版社
WILEY
DOI: 10.1111/1755-0998.13326
关键词
bioinfomatics; phyloinfomatics; genomics; proteomics; molecular evolution; population genetics – empirical; software
资金
- Division of Environmental Biology [DEB-1754439]
- National Institute of General Medical Sciences [1R35GM133481-01]
Population genetic analyses often use summary statistics such as pi and d(XY) to describe genetic variation patterns. However, biases can be introduced due to missing data, which can lead to underestimation of genetic diversity within and between populations. In this study, a user-friendly UNIX command line utility, pixy, is introduced to address this issue and provide unbiased estimates of pi and d(XY) regardless of the amount or form of missing data.
Population genetic analyses often use summary statistics to describe patterns of genetic variation and provide insight into evolutionary processes. Among the most fundamental of these summary statistics are pi and d(XY), which are used to describe genetic diversity within and between populations, respectively. Here, we address a widespread issue in pi and d(XY) calculation: systematic bias generated by missing data of various types. Many popular methods for calculating pi and d(XY) operate on data encoded in the variant call format (VCF), which condenses genetic data by omitting invariant sites. When calculating pi and d(XY) using a VCF, it is often implicitly assumed that missing genotypes (including those at sites not represented in the VCF) are homozygous for the reference allele. Here, we show how this assumption can result in substantial downward bias in estimates of pi and d(XY) that is directly proportional to the amount of missing data. We discuss the pervasive nature and importance of this problem in population genetics, and introduce a user-friendly UNIX command line utility, pixy, that solves this problem via an algorithm that generates unbiased estimates of pi and d(XY) in the face of missing data. We compare pixy to existing methods using both simulated and empirical data, and show that pixy alone produces unbiased estimates of pi and d(XY) regardless of the form or amount of missing data. In summary, our software solves a long-standing problem in applied population genetics and highlights the importance of properly accounting for missing data in population genetic analyses.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据