How to identify additional SNPs on EPICv2

Question

[Basti]

Hi everyone,

I realized that many CpGs had a 3-distribution DNA methylation profile on EPIC v2 chips, suggesting that they are SNPs.

I've already used dropLociWithSnps with the most recent annotation, so there must be additional unidentified SNPs. To your knowledge, are there any tools that would enable me to eliminate these probes?

I thought of developing an algorithm that would identify this distribution in 3 groups at 0%, 50% and 100%, but there might be a more precise tool based on the identification of these SNPs according to their genomic positions.

Thank you for your suggestions.

Answer 1

[Papyrus]

There is actually an algorithm, you can try the gaphunter function from the minfi package. I typically play with the threshold and outCutoff parameters because they depend on your dataset size, etc. but have found the function to really detect these types of SNP-like probes.

[Basti]

Thank you very much, I'm definitely going to test this !

Answer 2

[Joe]

Our recent work provides annotation files for unambiguous probes on the Illumina 450K, EPIC/850K, and EPICv2.0/935K methylation arrays. Additionally, we offer T2T-CHM13-based annotation information, including genomic features, SNPs with a minor allele frequency > 0.05, repetitive elements, and epigenetic marks. These files are available for download at https://github.com/functionalepigenomics/Illumina_Infinium_HumanMethylation_BeadChips_Annotation.

[Basti]

Thank you very much for this valuable new annotation. However, I still observe unambiguous CpGs that are not annotated to any SNPs in your work but that still show similar pattern in my dataset. I will try to check with other EPICv2 datasets, but I'm convinced there is a SNP there.