Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2017 May 5;7(5):1569-1575.
doi: 10.1534/g3.117.039693.

asymptoticMK: A Web-Based Tool for the Asymptotic McDonald-Kreitman Test

Affiliations

asymptoticMK: A Web-Based Tool for the Asymptotic McDonald-Kreitman Test

Benjamin C Haller et al. G3 (Bethesda). .

Abstract

The McDonald-Kreitman (MK) test is a widely used method for quantifying the role of positive selection in molecular evolution. One key shortcoming of this test lies in its sensitivity to the presence of slightly deleterious mutations, which can severely bias its estimates. An asymptotic version of the MK test was recently introduced that addresses this problem by evaluating polymorphism levels for different mutation frequencies separately, and then extrapolating a function fitted to that data. Here, we present asymptoticMK, a web-based implementation of this asymptotic MK test. Our web service provides a simple R-based interface into which the user can upload the required data (polymorphism and divergence data for the genomic test region and a neutrally evolving reference region). The web service then analyzes the data and provides plots of the test results. This service is free to use, open-source, and available at http://benhaller.com/messerlab/asymptoticMK.html We provide results from simulations to illustrate the performance and robustness of the asymptoticMK test under a wide range of model parameters.

Keywords: molecular evolution; positive selection; web service.

PubMed Disclaimer

Figures

Figure 1
Figure 1
A screenshot of the web page for asymptoticMK. After entering values for d and d0, choosing an input file with binned values for x, p, and p0, and choosing the x interval to fit, the user can click the Submit button and asymptoticMK will provide its results in a new browser window or tab.
Figure 2
Figure 2
Results from asymptoticMK for three test data sets. (A) Normalized site frequency spectrum (SFS) for the Drosophila data set used in Messer and Petrov (2013). Points show normalized binned polymorphism frequencies for the neutral region (black) and the test region (red). (B) Result of asymptoticMK’s analysis of that data set. The two vertical blue lines show the limits of the frequency cutoff interval used for fitting. Points indicate binned values of α(x), estimated according to Equation 2; points are gray if they are outside the cutoff interval (and thus not used in fitting). The solid red curve shows the fitted αfit(x) (here, exponential). The dashed red line shows the estimate of αasymptotic, obtained from the fitted function according to Equation 3. The gray band indicates the 95% C.I. around this αasymptotic estimate. The dotted gray line shows the estimate of αoriginal, obtained from the original (nonasymptotic) McDonald–Kreitman (MK) test, for comparison (also calculated using only the data within the cutoff interval). (C) and (D) show corresponding results from one SLiM simulation run, and (E) and (F) show results from another SLiM simulation run; in each case, the first panel shows the result of an automated fit using asymptoticMK, whereas the second shows the improvement after hand tailoring of the fit (see Results and Discussion). Note that in all four cases, the linear fit was deemed more appropriate by asymptoticMK. The solid green horizontal lines, finally, show the true value of α in the simulation runs for comparison.

Similar articles

Cited by

References

    1. Andolfatto P., 2005. Adaptive evolution of non-coding DNA in Drosophila. Nature 437: 1149–1152. - PubMed
    1. Boyko A. R., Williamson S. H., Indap A. R., Degenhardt J. D., Hernandez R. D., et al. , 2008. Assessing the evolutionary impact of amino acid mutations in the human genome. PLoS Genet. 4: e1000083. - PMC - PubMed
    1. Bustamante C. D., Fledel-Alon A., Williamson S., Nielsen R., Hubisz M. T., et al. , 2005. Natural selection on protein-coding genes in the human genome. Nature 437: 1153–1157. - PubMed
    1. Charlesworth J., Eyre-Walker A., 2008. The McDonald-Kreitman test and slightly deleterious mutations. Mol. Biol. Evol. 25: 1007–1015. - PubMed
    1. Clark A. G., Eisen M. B., Smith D. R., Bergman C. M., Oliver B., et al. , 2007. Evolution of genes and genomes on the Drosophila phylogeny. Nature 450: 203–218. - PubMed

Publication types

LinkOut - more resources