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
Comment
. 2023 Dec 21;6(1):1295.
doi: 10.1038/s42003-023-05658-5.

Genetic differentiation at probe SNPs leads to spurious results in meQTL discovery

Affiliations
Comment

Genetic differentiation at probe SNPs leads to spurious results in meQTL discovery

Gillian L Meeks et al. Commun Biol. .
No abstract available

PubMed Disclaimer

Conflict of interest statement

The authors declare no competing interests.

Figures

Fig. 1
Fig. 1. SNPs in Li et al.’s meQTL CpG probes are more differentiated than SNPs in non-meQTL CpG probes between 1000 Genomes EUR and AFR super-populations.
a Fst distributions between meQTL (red) and non-meQTL (gray) probe SNPs are significantly different (two-tailed t-test p-value = 2.2e−16, 95% confidence interval for the difference in mean Fst [0.20–0.23]). The boxplots show the median (horizontal line), interquartile range (box), and 1.5x interquartile range (whiskers). b Probe SNP frequencies in AFR vs. EUR of the same probe SNPs (red points, Li et al.’s meQTL; gray points, non-meQTL) further highlight that they tend to be significantly differentiated, i.e., they tend to deviate from the dashed black line that indicates identical EUR and AFR allele frequencies.
Fig. 2
Fig. 2. Probe SNPs bias effect size estimates.
a MeQTL involving CpGs with probe SNPs (blue) have larger effect sizes than those without (yellow) (95% confidence interval for difference in mean effect size [0.28, 0.45], two-tailed t-test p-value = 1.44e−15). b A cubic regression of the distance between the probe SNP and target CpG explains 42% of the variation in meQTL effect size for the 351 CpG probes with SNPs (F-test p-value over a quadratic model = 0.0006). c LA-specific meQTL involving CpGs with probe SNPs (blue) have larger differences in effect sizes (Δ) between European and African local ancestry than those without (yellow) (95% confidence interval for difference in mean Δ effect size [0.24, 0.70], two-tailed t-test p-value = 0.0002). d A quadratic regression of the distance between the probe SNP and target CpG explains 25% of the variation in Δ effect sizes for the 58 CpG probes with SNPs (F-test p-value over a linear model = 0.014). The boxplots in a and c show the median (horizontal line), interquartile range (box), and 1.5× interquartile range (whiskers). The ribbons in b and d show 2 standard errors above and below the model predictions.

Comment in

Comment on

References

    1. Carja O, et al. Worldwide patterns of human epigenetic variation. Nat. Ecol. Evol. 2017;1:1577–1583. doi: 10.1038/s41559-017-0299-z. - DOI - PMC - PubMed
    1. Fraser HB, Lam LL, Neumann SM, Kobor MS. Population-specificity of human DNA methylation. Genome Biol. 2012;13:R8. doi: 10.1186/gb-2012-13-2-r8. - DOI - PMC - PubMed
    1. Heyn H, et al. DNA methylation contributes to natural human variation. Genome Res. 2013;23:1363–1372. doi: 10.1101/gr.154187.112. - DOI - PMC - PubMed
    1. Rahmani E, et al. Genome-wide methylation data mirror ancestry information. Epigenet. Chromatin. 2017;10:1. doi: 10.1186/s13072-016-0108-y. - DOI - PMC - PubMed
    1. Rawlik K, et al. Evidence of epigenetic admixture in the Colombian population. Hum. Mol. Genet. 2017;26:501–508. - PMC - PubMed