TY - JOUR
T1 - Weight of the evidence of genetic investigations of ancestry informative markers
AU - Tvedebrink, Torben
AU - Eriksen, Poul Svante
AU - Mogensen, Helle Smidt
AU - Morling, Niels
N1 - Copyright © 2017 Elsevier Inc. All rights reserved.
PY - 2018/3
Y1 - 2018/3
N2 - Ancestry-informative markers (AIMs) are markers that give information about the ancestry of individuals. They are used in forensic genetics for predicting the geographic origin of the investigated individual in crime and identification cases. In the exploration of the genogeographic origin of an AIMs profile, the likelihoods of the AIMs profile in various populations may be calculated. However, there may not be an appropriate reference population in the database. The fact that the likelihood ratio (LR) of one population compared to that of another population is large does not imply that any of the populations is relevant. To handle this phenomena, we derived a likelihood ratio test (LRT) that is a measure of absolute concordance between an AIMs profile and a population rather than a relative measure of the AIMs profile's likelihood in two populations. The LRT is similar to a Fisher's exact test. By aggregating over markers, the central limit theorem suggests that the resulting quantity is approximately normally distributed. If only a few markers are genotyped or if the majority of the markers are fixed in a given population, the approximation may fail. We overcome this using importance sampling and show how exponential tilting results in an efficient proposal distribution. By simulations and published AIMs profiles, we demonstrate the applicability of the derived methodology. For the genotyped AIMs, the LRT approach achieves the nominal levels of rejection when tested on data from five major continental regions.
AB - Ancestry-informative markers (AIMs) are markers that give information about the ancestry of individuals. They are used in forensic genetics for predicting the geographic origin of the investigated individual in crime and identification cases. In the exploration of the genogeographic origin of an AIMs profile, the likelihoods of the AIMs profile in various populations may be calculated. However, there may not be an appropriate reference population in the database. The fact that the likelihood ratio (LR) of one population compared to that of another population is large does not imply that any of the populations is relevant. To handle this phenomena, we derived a likelihood ratio test (LRT) that is a measure of absolute concordance between an AIMs profile and a population rather than a relative measure of the AIMs profile's likelihood in two populations. The LRT is similar to a Fisher's exact test. By aggregating over markers, the central limit theorem suggests that the resulting quantity is approximately normally distributed. If only a few markers are genotyped or if the majority of the markers are fixed in a given population, the approximation may fail. We overcome this using importance sampling and show how exponential tilting results in an efficient proposal distribution. By simulations and published AIMs profiles, we demonstrate the applicability of the derived methodology. For the genotyped AIMs, the LRT approach achieves the nominal levels of rejection when tested on data from five major continental regions.
KW - Journal Article
U2 - 10.1016/j.tpb.2017.12.004
DO - 10.1016/j.tpb.2017.12.004
M3 - Journal article
C2 - 29278682
SN - 0040-5809
VL - 120
SP - 1
EP - 10
JO - Theoretical Population Biology
JF - Theoretical Population Biology
ER -