Missense SNP rs3731249 Explains the CDKN2A Association with Childhood ALL and Shows Risk Allele Selection in Tumors with Somatic CDKN2A Alterations

Little is known about the interaction between constitutive and somatic genetic variation in carcinogenesis. Genome-wide association studies (GWAS) of childhood acute lymphoblastic leukemia (ALL) have identified 6 genetic associations, including the single nucleotide polymorphism (SNP) rs3731217, located between the CDKN2A and CDKN2B genes. Somatic deletions of this chromosome 9p21.3 region are common in cancer, and occur in ~30% of childhood ALL.

We replicated the rs3731217 association with childhood ALL in Hispanics (P = 0.021) via a GWAS of 323 childhood ALL cases and 454 controls, using the Illumina OmniExpress SNP array. SNP imputation analysis across chromosome 9p21.3 revealed a stronger association at SNP rs3731249 (P = 6.4x10^-4; OR = 2.77), a low-frequency missense variant in exon 2 (p.Ala148Thr) of CDKN2A (MAF_controls = 2.1%). These associations were replicated in an additional 378 ALL cases and 536 controls using Taqman assays, where the missense variant again had a greater magnitude of effect (P_rs3731249 = 6.5x10^-3; OR _rs3731249 = 2.00 versus P_rs3731217 = 0.029; OR _rs3731217= 1.38). We assessed the effect of rs3731249 on risk of other cancers with known chr9p21.3 GWAS associations, but observed no association between this variant and either adult glioblastoma or melanoma risk, suggesting an ALL-specific role for rs3731249.

We hypothesized that cases heterozygous for rs3731249 may show allelic imbalance, with preferential loss of the wildtype (WT) allele when hemizygous chr9p21.3 deletions occur in leukemic cells. To test this, we developed a novel method termed “SMART-ddPCR” (Somatic Mutation Allelic Ratio Test using Droplet Digital PCR). Concentration of risk and WT alleles was measured by ddPCR using the rs3731249 Taqman assay in diagnostic bone marrow (i.e. tumor) DNA from 35 known heterozygotes. Allelic imbalance due to somatic CDKN2A alterations was identified in 17 cases, of which 14 had a higher risk:WT allele ratio (P = 0.006, 1-sided binomial significance test), thus demonstrating preferential loss of WT rs3731249 and retention of the missense allele. Subclonal heterogeneity in CDKN2A copy number was also evident, suggesting that CDKN2A alterations are secondary events in leukemogenesis.

Preferential allelic imbalance towards the rs3731249 risk allele in ALL tumor DNA provides strong evidence of an important role for this missense variant in leukemogenesis. We show, for the first time, a direct relationship between constitutive and somatic genetic variants underlying ALL development, suggesting that such constitutive variants can provide an additional substrate for selection during tumor evolution.