Noncoding genetic variation in GATA3 increases acute lymphoblastic leukemia risk through local and global changes in chromatin conformation

Inherited noncoding genetic variants confer significant disease susceptibility to childhood acute lymphoblastic leukemia (ALL) but the molecular processes linking germline polymorphisms with somatic lesions in this cancer are poorly understood. Through targeted sequencing in 5,008 patients, we identified a key regulatory germline variant in GATA3 associated with Philadelphia chromosome-like ALL (Ph-like ALL). Using CRISPR-Cas9 editing and samples from patients with Ph-like ALL, we showed that this variant activated a strong enhancer that upregulated GATA3 transcription. This, in turn, reshaped global chromatin accessibility and three-dimensional genome organization, including regions proximal to the ALL oncogene CRLF2. Finally, we showed that GATA3 directly regulated CRLF2 and potentiated the JAK-STAT oncogenic effects during leukemogenesis. Taken together, we provide evidence for a distinct mechanism by which a germline noncoding variant contributes to oncogene activation, epigenetic regulation and three-dimensional genome reprogramming.

Abstract 2118: Non-coding germline GATA3 variants alter chromatin topology and contribute to pathogenesis of acute lymphoblastic leukemia

Inherited non-coding genetic variants confer significant disease susceptibility in many cancers. However, the molecular processes of by which germline variants contribute to somatic lesions are poorly understood. We performed targeted sequencing in 5,008 patients and identified a key regulatory germline variant in GATA3 strongly associated with Philadelphia chromosome-like acute lymphoblastic leukemia (Ph-like ALL). By creating an isogenic cellular model with CRISPR-Cas9 system, we showed that this variant activated a strong enhancer that significantly upregulated GATA3 transcription, which in turn reshaped the global chromatin accessibility and 3D genome organization. Remarkably, this genotype switch induced a chromatin loop between the CRLF2 oncogene and a distal enhancer, similar to the somatically acquired super-enhancer hijacking event in patients. GATA3 genotype-related alterations in transcriptional control and 3D chromatin organization were further validated in Ph-like ALL patients. Finally, we showed that GATA3 directly regulates CRLF2 and potentiates the oncogenic effects of JAK-STAT signaling in leukemogenesis. Altogether, our results provide evidence for a novel mechanism by which a germline non-coding variant contributes to oncogene activation epigenetic regulation and 3D genome reprogramming.

Citation Format: Hongbo Yang, Hui Zhang, Yu Luan, Tingting Liu, Kathryn Roberts, Mao-xiang Qian, Bo Zhang, Wenjian Yang, Virginia Perez-Andreu, Jie Xu, Sriranga lyyanki, Da Kuang, Shalini Reshmi, Julie Gastier-Foster, Colton Smith, Ching-Hon Pui, William Evans, Stephen Hunger, Stephen Hunger, Leonidas Platanias, Mary Relling, Charles Mullighan, Mignon Loh, Feng Yue, Jun Yang. Non-coding germline GATA3 variants alter chromatin topology and contribute to pathogenesis of acute lymphoblastic leukemia [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2118.

Novel susceptibility variants at the ERG locus for childhood acute lymphoblastic leukemia in Hispanics

Acute lymphoblastic leukemia (ALL) is the most common malignancy in children. Characterized by high levels of Native American ancestry, Hispanics are disproportionally affected by this cancer with high incidence and inferior survival. However, the genetic basis for this disparity remains poorly understood because of a paucity of genome-wide investigation of ALL in Hispanics. Performing a genome-wide association study (GWAS) in 940 Hispanic children with ALL and 681 ancestry-matched non-ALL controls, we identified a novel susceptibility locus in the ERG gene (rs2836365; P = 3.76 × 10-8; odds ratio [OR] = 1.56), with independent validation (P = .01; OR = 1.43). Imputation analyses pointed to a single causal variant driving the association signal at this locus overlapping with putative regulatory DNA elements. The effect size of the ERG risk variant rose with increasing Native American genetic ancestry. The ERG risk genotype was underrepresented in ALL with the ETV6-RUNX1 fusion (P < .0005) but enriched in the TCF3-PBX1 subtype (P < .05). Interestingly, ALL cases with germline ERG risk alleles were significantly less likely to have somatic ERG deletion (P < .05). Our results provide novel insights into genetic predisposition to ALL and its contribution to racial disparity in this cancer.

Creating a genotype-based dosing algorithm for acenocoumarol steady dose

Acenocoumarol is a commonly prescribed anticoagulant drug for the prophylaxis and treatment of venous and arterial thromboembolic disorders in several countries. In counterpart of warfarin, there is scarce information about pharmacogenetic algorithms for steady acenocoumarol dose estimation. The aim of this study was to develop an algorithm of prediction for acenocoumarol.The algorithm was created using the data from 973 retrospectively selected anticoagulated patients and was validated in a second independent cohort adding up to 2,683 patients. The best regression model to predict stable dosage in the Primary Cohort included clinical factors (age and body mass index, BSA) and genetic variants (VKORC1, CYP2C9* and CYP4F2 polymorphisms) and explained up to 50% of stable dose. In the validation study the clinical algorithm yielded an adjusted R2=0.15 (estimation´s standard error=4.5) and the genetic approach improved the dose forecast up to 30% (estimation´s standard error=4.6). Again, the best model combined clinical and genetic factors (R2= 0.48; estimation´s standard error=4) which provided the best results of doses estimates within 20% of the real dose in patients taking lower (≤7mg/week) or higher (≥25mg/week) acenocoumarol doses. In conclusion, we developed a prediction algorithm using clinical data and three polymorphisms in VKORC1, CYP2C9* and CYP4F2 genes that provided a steady acenocoumarol dose for about 50% of patients in the Validation Cohort. Such algorithm was especially useful to patients who need higher or lower acenocoumarol doses, those patients with higher time required until their stabilisation and are more prone to suffer a treatment derived complication.

Family-based exome-wide assessment of maternal genetic effects on susceptibility to childhood B-cell acute lymphoblastic leukemia in hispanics

BACKGROUND

Children of Hispanic ancestry have a higher incidence of acute lymphoblastic leukemia (ALL) compared with other ethnic groups, but to the authors' knowledge, the genetic basis for these racial disparities remain incompletely understood. Genome-wide association studies of childhood ALL to date have focused on inherited genetic effects; however, maternal genetic effects (the role of the maternal genotype on phenotype development in the offspring) also may play a role in ALL susceptibility.

METHODS

The authors conducted a family-based exome-wide association study of maternal genetic effects among Hispanics with childhood B-cell ALL using the Illumina Infinium HumanExome BeadChip. A discovery cohort of 312 Guatemalan and Hispanic American families and an independent replication cohort of 152 Hispanic American families were used.

RESULTS

Three maternal single-nucleotide polymorphisms (SNPs) approached the study threshold for significance after correction for multiple testing (P<1.0 × 10^-6 ): MTL5 rs12365708 (testis expressed metallothionein-like protein [tesmin]) (relative risk [RR], 2.62; 95% confidence interval [95% CI], 1.61-4.27 [P = 1.8 × 10^-5 ]); ALKBH1 rs6494 (AlkB homolog 1, histone H2A dioxygenase) (RR, 3.77; 95% CI, 1.84-7.74 [P = 3.7 × 10^-5 ]); and NEUROG3 rs4536103 (neurogenin 3) (RR, 1.75; 95% CI, 1.30-2.37 [P = 1.2 × 10^-4 ]). Although effect sizes were similar, these SNPs were not nominally significant in the replication cohort in the current study. In a meta-analysis comprised of the discovery cohort and the replication cohort, these SNPs were still not found to be statistically significant after correction for multiple comparisons (rs12365708: pooled RR, 2.27 [95% CI, 1.48-3.50], P = 1.99 × 10^-4 ; rs6494: pooled RR, 2.31 [95% CI, 1.38-3.85], P = .001; and rs4536103: pooled RR, 1.67 [95% CI, 1.29-2.16] P = 9.23 × 10^-5 ).

CONCLUSIONS

In what to the authors' knowledge is the first family-based based exome-wide association study to investigate maternal genotype effects associated with childhood ALL, the results did not implicate a strong role of maternal genotype on disease risk among Hispanics; however, 3 maternal SNPs were identified that may play a modest role in susceptibility. Cancer 2016;122:3697-704. © 2016 American Cancer Society.

NUDT15 polymorphisms alter thiopurine metabolism and hematopoietic toxicity

Widely used as anticancer and immunosuppressive agents, thiopurines have narrow therapeutic indices owing to frequent toxicities, partly explained by TPMT genetic polymorphisms. Recent studies identified germline NUDT15 variation as another critical determinant of thiopurine intolerance, but the underlying molecular mechanisms and the clinical implications of this pharmacogenetic association remain unknown. In 270 children enrolled in clinical trials for acute lymphoblastic leukemia in Guatemala, Singapore and Japan, we identified four NUDT15 coding variants (p.Arg139Cys, p.Arg139His, p.Val18Ile and p.Val18_Val19insGlyVal) that resulted in 74.4-100% loss of nucleotide diphosphatase activity. Loss-of-function NUDT15 diplotypes were consistently associated with thiopurine intolerance across the three cohorts (P = 0.021, 2.1 × 10(-5) and 0.0054, respectively; meta-analysis P = 4.45 × 10(-8), allelic effect size = -11.5). Mechanistically, NUDT15 inactivated thiopurine metabolites and decreased thiopurine cytotoxicity in vitro, and patients with defective NUDT15 alleles showed excessive levels of thiopurine active metabolites and toxicity. Taken together, these results indicate that a comprehensive pharmacogenetic model integrating NUDT15 variants may inform personalized thiopurine therapy.

Inherited coding variants at the CDKN2A locus influence susceptibility to acute lymphoblastic leukaemia in children

There is increasing evidence from genome-wide association studies for a strong inherited genetic basis of susceptibility to acute lymphoblastic leukaemia (ALL) in children, yet the effects of protein-coding variants on ALL risk have not been systematically evaluated. Here we show a missense variant in CDKN2A associated with the development of ALL at genome-wide significance (rs3731249, P=9.4 × 10⁻²³, odds ratio=2.23). Functional studies indicate that this hypomorphic variant results in reduced tumour suppressor function of p16^INK4A, increases the susceptibility to leukaemic transformation of haematopoietic progenitor cells, and is preferentially retained in ALL tumour cells. Resequencing the CDKN2A–CDKN2B locus in 2,407 childhood ALL cases reveals 19 additional putative functional germline variants. These results provide direct functional evidence for the influence of inherited genetic variation on ALL risk, highlighting the important and complex roles of CDKN2A–CDKN2B tumour suppressors in leukaemogenesis.

Genome-wide association studies indicate a strong genetic susceptibility to acute lymphoblastic leukaemia in children, though the effect on protein-coding genes is not fully understood. Here Xu and Zhang et al. identify a missense variant in CDKN2A which reduces tumour suppression.

Novel associations of VKORC1 variants with higher acenocoumarol requirements

BACKGROUND

Algorithms combining both clinical and genetic data have been developed to improve oral anticoagulant therapy. Three polymorphisms in two genes, VKORC1 and CYP2C9, are the main coumarin dose determinants and no additional polymorphisms of any relevant pharmacogenetic importance have been identified.

OBJECTIVES

To identify new genetic variations in VKORC1 with relevance for oral anticoagulant therapy.

METHODS AND RESULTS

3949 consecutive patients taking acenocoumarol were genotyped for the VKORC1 rs9923231 and CY2C9* polymorphisms. Of these, 145 patients with a dose outside the expected range for the genetic profile determined by these polymorphisms were selected. Clinical factors explained the phenotype in 88 patients. In the remaining 57 patients, all with higher doses than expected, we sequenced the VKORC1 gene and genetic changes were identified in 14 patients. Four patients carried VKORC1 variants previously related to high coumarin doses (L128R, N = 1 and D36Y, N = 3).Three polymorphisms were also detected: rs17878544 (N = 5), rs55894764 (N = 4) and rs7200749 (N = 2) which was in linkage disequilibrium with rs17878544. Finally, 2 patients had lost the rs9923231/rs9934438 linkage. The prevalence of these variations was higher in these patients than in the whole sample. Multivariate linear regression analysis revealed that only D36Y and rs55894764 variants significantly affect the dose, although the improvement in the prediction model is small (from 39% to 40%).

CONCLUSION

Our strategy identified novel associations of VKORC1 variants with higher acenocoumarol doses albeit with a low effect size. Further studies are necessary to test their influence on the VKORC1 function and the cost/benefit of their inclusion in pharmacogenetic algorithms.

Novel susceptibility variants at 10p12.31-12.2 for childhood acute lymphoblastic leukemia in ethnically diverse populations

BACKGROUND

Acute lymphoblastic leukemia (ALL) is the most common cancer in children and the incidence of ALL varies by ethnicity. Although accumulating evidence indicates inherited predisposition to ALL, the genetic basis of ALL susceptibility in diverse ancestry has not been comprehensively examined.

METHODS

We performed a multiethnic genome-wide association study in 1605 children with ALL and 6661 control subjects after adjusting for population structure, with validation in three replication series of 845 case subjects and 4316 control subjects. Association was tested by two-sided logistic regression.

RESULTS

A novel ALL susceptibility locus at 10p12.31-12.2 (BMI1-PIP4K2A, rs7088318, P = 1.1 × 10(-11)) was identified in the genome-wide association study, with independent replication in European Americans, African Americans, and Hispanic Americans (P = .001, .009, and .04, respectively). Association was also validated at four known ALL susceptibility loci: ARID5B, IKZF1, CEBPE, and CDKN2A/2B. Associations at ARID5B, IKZF1, and BMI1-PIP4K2A variants were consistent across ethnicity, with multiple independent signals at IKZF1 and BMI1-PIP4K2A loci. The frequency of ARID5B and BMI1-PIP4K2A variants differed by ethnicity, in parallel with ethnic differences in ALL incidence. Suggestive evidence for modifying effects of age on genetic predisposition to ALL was also observed. ARID5B, IKZF1, CEBPE, and BMI1-PIP4K2A variants cumulatively conferred strong predisposition to ALL, with children carrying six to eight copies of risk alleles at a ninefold (95% confidence interval = 6.9 to 11.8) higher ALL risk relative to those carrying zero to one risk allele at these four single nucleotide polymorphisms.

CONCLUSIONS

These findings indicate strong associations between inherited genetic variation and ALL susceptibility in children and shed new light on ALL molecular etiology in diverse ancestry.