BMI1 enhancer polymorphism underlies chromosome 10p12.31 association with childhood acute lymphoblastic leukemia

Oncogenic Enhancers in Leukemia

Although the study of leukemogenesis has traditionally focused on protein-coding genes, the role of enhancer dysregulation is becoming increasingly recognized. The advent of high-throughput sequencing, together with a better understanding of enhancer biology, has revealed how various genetic and epigenetic lesions produce oncogenic enhancers that drive transformation. These aberrations include translocations that lead to enhancer hijacking, point mutations that modulate enhancer activity, and copy number alterations that modify enhancer dosage. In this review, we describe these mechanisms in the context of leukemia and discuss potential therapeutic avenues to target these regulatory elements. Significance: Large-scale sequencing projects have uncovered recurrent gene mutations in leukemia, but the picture remains incomplete: some patients harbor no such aberrations, whereas others carry only a few that are insufficient to bring about transformation on their own. One of the missing pieces is enhancer dysfunction, which only recently has emerged as a critical driver of leukemogenesis. Knowledge of the various mechanisms of enhancer dysregulation is thus key for a complete understanding of leukemia and its causes, as well as the development of targeted therapies in the era of precision medicine.

The genetic risk of acute lymphoblastic leukemia and its implications for children of Latin American origin

Acute lymphoblastic leukemia (ALL) is the most common cancer in children, and disproportionately affects children of Hispanic/Latino ethnicity in the United States, who have the highest incidence of disease compared with other racial/ethnic groups. Incidence of childhood ALL is similarly high in several Latin American countries, notably in Mexico, and of concern is the rising incidence of childhood ALL in some Hispanic/Latino populations that may further widen this disparity. Prior studies have implicated common germline genetic variants in the increased risk of ALL among Hispanic/Latino children. In this review, we describe the known disparities in ALL incidence as well as patient outcomes that disproportionately affect Hispanic/Latino children across the Americas, and we focus on the role of genetic variation as well as Indigenous American ancestry in the etiology of these disparities. Finally, we discuss future avenues of research to further our understanding of the causes of the disparities in ALL incidence and outcomes in children of Latin American origin, which will be required for future precision prevention efforts.

Evaluating genomic polygenic risk scores for childhood acute lymphoblastic leukemia in Latinos

The utility of polygenic risk score (PRS) models has not been comprehensively evaluated for childhood acute lymphoblastic leukemia (ALL), the most common type of cancer in children. Previous PRS models for ALL were based on significant loci observed in genome-wide association studies (GWASs), even though genomic PRS models have been shown to improve prediction performance for a number of complex diseases. In the United States, Latino (LAT) children have the highest risk of ALL, but the transferability of PRS models to LAT children has not been studied. In this study, we constructed and evaluated genomic PRS models based on either non-Latino White (NLW) GWAS or a multi-ancestry GWAS. We found that the best PRS models performed similarly between held-out NLW and LAT samples (PseudoR2 = 0.086 ± 0.023 in NLW vs. 0.060 ± 0.020 in LAT), and can be improved for LAT if we performed GWAS in LAT-only (PseudoR2 = 0.116 ± 0.026) or multi-ancestry samples (PseudoR2 = 0.131 ± 0.025). However, the best genomic models currently do not have better prediction accuracy than a conventional model using all known ALL-associated loci in the literature (PseudoR2 = 0.166 ± 0.025), which includes loci from GWAS populations that we could not access to train genomic PRS models. Our results suggest that larger and more inclusive GWASs may be needed for genomic PRS to be useful for ALL. Moreover, the comparable performance between populations may suggest a more oligogenic architecture for ALL, where some large effect loci may be shared between populations. Future PRS models that move away from the infinite causal loci assumption may further improve PRS for ALL.

Treatment outcomes of neuroendoscopic and microscopic trans-sphenoidal pituitary adenomectomies and the effects on hormone levels

BACKGROUND

We aimed to compare the treatment outcomes of neuroendoscopic and microscopic trans-sphenoidal pituitary adenomectomies, as well as the effects on hormone levels and clinical symptoms.

METHODS

A total of 82 patients with pituitary adenomas that were surgically resected from June 2018 to March 2021 were selected and divided into a group receiving neuroendoscopic trans-sphenoidal pituitary adenomectomy (group A, N.=40), and the other group receiving microscopic surgery (group B, N.=42). Surgery-related indices, hormone levels before discharge and alleviation of symptoms 24 weeks after surgery were compared.

RESULTS

Both groups had significantly different degrees of tumor resection (P<0.05). The proportion of cases receiving total adenomectomy in group A significantly exceeded that of group B (P<0.05). The surgical time of group A was significantly longer than that of group B (P<0.05). Group A had significantly shorter mean hospitalization stay than that of group B (P<0.05). The postoperative hormone levels of both groups decreased significantly differently (P<0.05). Before discharge, the hormone recovery rate of group A significantly surpassed that of group B (P<0.05). The hormone levels of cases with prolactinoma, adrenocorticotropic hormone adenoma and growth hormone adenoma in group A dropped more significantly than those of group B did (P<0.05).

CONCLUSIONS

Compared with microscopic surgery, neuroendoscopic trans-sphenoidal pituitary adenomectomy worked more effectively, induced fewer postoperative complications and better promoted the postoperative recovery of hormone levels.

Evaluating Genomic Polygenic Risk Scores for Childhood Acute Lymphoblastic Leukemia in Latinos

The utility of polygenic risk score (PRS) models has not been comprehensively evaluated for childhood acute lymphoblastic leukemia (ALL), the most common type of cancer in children. Previous PRS models for ALL were based on significant loci observed in genome-wide association studies (GWAS), even though genomic PRS models have been shown to improve prediction performance for a number of complex diseases. In the United States, Latino (LAT) children have the highest risk of ALL, but the transferability of PRS models to LAT children has not been studied. In this study we constructed and evaluated genomic PRS models based on either non-Latino white (NLW) GWAS or a multi-ancestry GWAS. We found that the best PRS models performed similarly between held-out NLW and LAT samples (PseudoR 2 = 0.086 ± 0.023 in NLW vs. 0.060 ± 0.020 in LAT), and can be improved for LAT if we performed GWAS in LAT-only (PseudoR 2 = 0.116 ± 0.026) or multi-ancestry samples (PseudoR 2 = 0.131 ± 0.025). However, the best genomic models currently do not have better prediction accuracy than a conventional model using all known ALL-associated loci in the literature (PseudoR 2 = 0.166 ± 0.025), which includes loci from GWAS populations that we could not access to train genomic PRS models. Our results suggest that larger and more inclusive GWAS may be needed for genomic PRS to be useful for ALL. Moreover, the comparable performance between populations may suggest a more oligo-genic architecture for ALL, where some large effect loci may be shared between populations. Future PRS models that move away from the infinite causal loci assumption may further improve PRS for ALL.

Overview on Aneuploidy in Childhood B-Cell Acute Lymphoblastic Leukemia

Recent years have brought significant progress in the treatment of B-cell acute lymphoblastic leukemia (ALL). This was influenced by both the improved schemes of conventionally used therapy, as well as the development of new forms of treatment. As a consequence, 5-year survival rates have increased and now exceed 90% in pediatric patients. For this reason, it would seem that everything has already been explored in the context of ALL. However, delving into its pathogenesis at the molecular level shows that there are many variations that still need to be analyzed in more detail. One of them is aneuploidy, which is among the most common genetic changes in B-cell ALL. It includes both hyperdiploidy and hypodiploidy. Knowledge of the genetic background is important already at the time of diagnosis, because the first of these forms of aneuploidy is characterized by a good prognosis, in contrast to the second, which is in favor of an unfavorable course. In our work, we will focus on summarizing the current state of knowledge on aneuploidy, along with an indication of all the consequences that may be correlated with it in the context of the treatment of patients with B-cell ALL.

Hyperdiploidy: the longest known, most prevalent, and most enigmatic form of acute lymphoblastic leukemia in children

Hyperdiploidy is the largest genetic entity B-cell precursor acute lymphoblastic leukemia in children. The diagnostic hallmark of its two variants that will be discussed in detail herein is a chromosome count between 52 and 67, respectively. The classical HD form consists of heterozygous di-, tri-, and tetrasomies, whereas the nonclassical one (usually viewed as "duplicated hyperhaploid") contains only disomies and tetrasomies. Despite their apparently different clinical behavior, we show that these two sub-forms can in principle be produced by the same chromosomal maldistribution mechanism. Moreover, their respective array, gene expression, and mutation patterns also indicate that they are biologically more similar than hitherto appreciated. Even though in-depth analyses of the genomic intricacies of classical HD leukemias are indispensable for the elucidation of the disease process, the ensuing results play at present surprisingly little role in treatment stratification, a fact that can be attributed to the overall good prognoses and low relapse rates of the concerned patients and, consequently, their excellent treatment outcome. Irrespective of this underutilization, however, the detailed genetic characterization of HD leukemias may, especially in planned treatment reduction trials, eventually become important for further treatment stratification, patient management, and the clinical elucidation of outcome data. It should therefore become an integral part of all upcoming treatment studies.

Investigating DNA methylation as a mediator of genetic risk in childhood acute lymphoblastic leukemia

Genome-wide association studies have identified a growing number of single nucleotide polymorphisms (SNPs) associated with childhood acute lymphoblastic leukemia (ALL), yet the functional roles of most SNPs are unclear. Multiple lines of evidence suggest that epigenetic mechanisms may mediate the impact of heritable genetic variation on phenotypes. Here, we investigated whether DNA methylation mediates the effect of genetic risk loci for childhood ALL. We performed an epigenome-wide association study (EWAS) including 808 childhood ALL cases and 919 controls from California-based studies using neonatal blood DNA. For differentially methylated CpG positions (DMPs), we next conducted association analysis with 23 known ALL risk SNPs followed by causal mediation analyses addressing the significant SNP-DMP pairs. DNA methylation at CpG cg01139861, in the promoter region of IKZF1, mediated the effects of the intronic IKZF1 risk SNP rs78396808, with the average causal mediation effect (ACME) explaining ~30% of the total effect (ACME P = 0.0031). In analyses stratified by self-reported race/ethnicity, the mediation effect was only significant in Latinos, explaining ~41% of the total effect of rs78396808 on ALL risk (ACME P = 0.0037). Conditional analyses confirmed the presence of at least three independent genetic risk loci for childhood ALL at IKZF1, with rs78396808 unique to non-European populations. We also demonstrated that the most significant DMP in the EWAS, CpG cg13344587 at gene ARID5B (P = 8.61 × 10-10), was entirely confounded by the ARID5B ALL risk SNP rs7090445. Our findings provide new insights into the functional pathways of ALL risk SNPs and the DNA methylation differences associated with risk of childhood ALL.

The landscape of GWAS validation; systematic review identifying 309 validated non-coding variants across 130 human diseases

Background

The remarkable growth of genome-wide association studies (GWAS) has created a critical need to experimentally validate the disease-associated variants, 90% of which involve non-coding variants.

Methods

To determine how the field is addressing this urgent need, we performed a comprehensive literature review identifying 36,676 articles. These were reduced to 1454 articles through a set of filters using natural language processing and ontology-based text-mining. This was followed by manual curation and cross-referencing against the GWAS catalog, yielding a final set of 286 articles.

Results

We identified 309 experimentally validated non-coding GWAS variants, regulating 252 genes across 130 human disease traits. These variants covered a variety of regulatory mechanisms. Interestingly, 70% (215/309) acted through cis-regulatory elements, with the remaining through promoters (22%, 70/309) or non-coding RNAs (8%, 24/309). Several validation approaches were utilized in these studies, including gene expression (n = 272), transcription factor binding (n = 175), reporter assays (n = 171), in vivo models (n = 104), genome editing (n = 96) and chromatin interaction (n = 33).

Conclusions

This review of the literature is the first to systematically evaluate the status and the landscape of experimentation being used to validate non-coding GWAS-identified variants. Our results clearly underscore the multifaceted approach needed for experimental validation, have practical implications on variant prioritization and considerations of target gene nomination. While the field has a long way to go to validate the thousands of GWAS associations, we show that progress is being made and provide exemplars of validation studies covering a wide variety of mechanisms, target genes, and disease areas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01216-w.

Genetic determinants of blood-cell traits influence susceptibility to childhood acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer. Despite overlap between genetic risk loci for ALL and hematologic traits, the etiological relevance of dysregulated blood-cell homeostasis remains unclear. We investigated this question in a genome-wide association study (GWAS) of childhood ALL (2,666 affected individuals, 60,272 control individuals) and a multi-trait GWAS of nine blood-cell indices in the UK Biobank. We identified 3,000 blood-cell-trait-associated (p < 5.0 × 10-8) variants, explaining 4.0% to 23.9% of trait variation and including 115 loci associated with blood-cell ratios (LMR, lymphocyte-to-monocyte ratio; NLR, neutrophil-to-lymphocyte ratio; PLR, platelet-to-lymphocyte ratio). ALL susceptibility was genetically correlated with lymphocyte counts (rg = 0.088, p = 4.0 × 10-4) and PLR (rg = -0.072, p = 0.0017). In Mendelian randomization analyses, genetically predicted increase in lymphocyte counts was associated with increased ALL risk (odds ratio [OR] = 1.16, p = 0.031) and strengthened after accounting for other cell types (OR = 1.43, p = 8.8 × 10-4). We observed positive associations with increasing LMR (OR = 1.22, p = 0.0017) and inverse effects for NLR (OR = 0.67, p = 3.1 × 10-4) and PLR (OR = 0.80, p = 0.002). Our study shows that a genetically induced shift toward higher lymphocyte counts, overall and in relation to monocytes, neutrophils, and platelets, confers an increased susceptibility to childhood ALL.

Novel Associations Between METTL3 Gene Polymorphisms and Pediatric Acute Lymphoblastic Leukemia: A Five-Center Case-Control Study

Objective

To reveal the contributing role of METTL3 gene SNPs in pediatric ALL risk.

Patients and Methods

A total of 808 pediatric ALL cases and 1,340 cancer-free controls from five hospitals in South China were recruited. A case-control study by genotyping three SNPs in the METTL3 gene was conducted. Genomic DNA was abstracted from peripheral blood. Three SNPs (rs1263801 C>G, rs1139130 A>G, and rs1061027 A>C) in the METTL3 gene were chosen to be detected by taqman real-time polymerase chain reaction assay.

Results

That rs1263801 C>G, rs1139130 A>G, and rs1061027 A>C polymorphisms were significantly associated with increased pediatric ALL risk was identified. In stratification analyses, it was discovered that rs1263801 CC, rs1061027 AA, and rs1139130 GG carriers were more likely to develop ALL in subgroups of common B-ALL, MLL gene fusion. Rs1263801 CC and rs10610257 AA carriers were more possible to increase the risk of ALL in subgroups of low hyperdiploid, and all of these three SNPs exhibited a trend toward the risk of ALL. All of these three polymorphisms were associated with the primitive/naïve lymphocytes and MRD in marrow after chemotherapy in ALL children. Rs1263801 CC and rs1139130 AA alleles provided a protective effect on MRD ≥0.01% among CCCG-treated children. As for rs1139130, AA alleles provided a protective effect on MRD in marrow ≥0.01% on 33 days and 12 weeks among CCCG-treated children, but provided a risk effect on MRD in the marrow ≥0.01% among SCCLG-treated children. As for rs1263801 CC and rs1139130 AA, these two alleles provided a protective effect on MRD in the marrow ≥0.01% among CCCG-treated children.

Conclusion

In this study, we revealed that METTL3 gene polymorphisms were associated with increased pediatric ALL risk and indicated that METTL3 gene polymorphisms might be a potential biomarker for choosing ALL chemotherapeutics.

MYBL2 Gene Polymorphism Is Associated With Acute Lymphoblastic Leukemia Susceptibility in Children

Purpose

Although MYBL2 had been validated to participate in multiple cancers including leukemia, the role of MYBL2 polymorphisms in acute lymphoblastic leukemia (ALL) was still not clear. In this study, we aimed to evaluate the association between MYBL2 single nucleotide polymorphisms (SNPs) and ALL risk in children.

Methods

A total of 687 pediatric ALL cases and 971 cancer-free controls from two hospitals in South China were recruited. A case-control study by genotyping three SNPs in the MYBL2 gene (rs285162 C>T, rs285207 A>C, and rs2070235 A>G) was conducted. The associations were assessed by odds ratios (ORs) with corresponding 95% confidence intervals (CIs). Subgroup and stratification analyses were conducted to explore the association of rs285207 with ALL risk in terms of age, sex, immunophenotype, risk level, and other clinical characteristics. The false-positive report probability (FPRP) analysis was performed to verify each significant finding. Functional analysis in silico was used to evaluate the probability that rs285207 might influence the regulation of MYBL2 .

Results

Our study demonstrated that rs285207 was related to a decreased ALL risk (adjusted OR = 0.78; 95% CI = 0.63-0.97, P = 0.022) in the dominant model. The associations of rs285207 with ALL risk appeared stronger in patients with pre B ALL (adjusted OR=0.56; 95% CI=0.38-0.84, P=0.004), with normal diploid (adjusted OR=0.73; 95% CI=0.57-0.95, P=0.017), with low risk (adjusted OR=0.68; 95% CI=0.49-0.94, P=0.021), with lower WBC (adjusted OR=0.62; 95% CI=0.43-0.87, P=0.007) or lower platelet level (adjusted OR=0.76; 95% CI=0.59-0.96, P=0.023). With FPRP analysis, the significant association between the rs285207 polymorphism and decreased ALL risk was still noteworthy (FPRP=0.128). Functional analysis showed that IKZF1 bound to DNA motif overlapping rs285207 and had a higher preference for the risk allele A. As for rs285162 C>T and rs2070235 A>G, no significant was found between them and ALL risk.

Conclusion

In this study, we revealed that rs285207 polymorphism decreased the ALL risk in children, and rs285207 might alter the binding to IKZF1, which indicated that the MYBL2 gene polymorphism might be a potential biomarker of childhood ALL.

Enhancing B-Cell Malignancies-On Repurposing Enhancer Activity towards Cancer

B-cell lymphomas and leukemias derive from B cells at various stages of maturation and are the 6th most common cancer-related cause of death. While the role of several oncogenes and tumor suppressors in the pathogenesis of B-cell neoplasms was established, recent research indicated the involvement of non-coding, regulatory sequences. Enhancers are DNA elements controlling gene expression in a cell type- and developmental stage-specific manner. They ensure proper differentiation and maturation of B cells, resulting in production of high affinity antibodies. However, the activity of enhancers can be redirected, setting B cells on the path towards cancer. In this review we discuss different mechanisms through which enhancers are exploited in malignant B cells, from the well-studied translocations juxtaposing oncogenes to immunoglobulin loci, through enhancer dysregulation by sequence variants and mutations, to enhancer hijacking by viruses. We also highlight the potential of therapeutic targeting of enhancers as a direction for future investigation.

Genetic Variants in Transcription Factor Binding Sites in Humans: Triggered by Natural Selection and Triggers of Diseases

Variants of transcription factor binding sites (TFBSs) constitute an important part of the human genome. Current evidence demonstrates close links between nucleotides within TFBSs and gene expression. There are multiple pathways through which genomic sequences located in TFBSs regulate gene expression, and recent genome-wide association studies have shown the biological significance of TFBS variation in human phenotypes. However, numerous challenges remain in the study of TFBS polymorphisms. This article aims to cover the current state of understanding as regards the genomic features of TFBSs and TFBS variants; the mechanisms through which TFBS variants regulate gene expression; the approaches to studying the effects of nucleotide changes that create or disrupt TFBSs; the challenges faced in studies of TFBS sequence variations; the effects of natural selection on collections of TFBSs; in addition to the insights gained from the study of TFBS alleles related to gout, its associated comorbidities (increased body mass index, chronic kidney disease, diabetes, dyslipidemia, coronary artery disease, ischemic heart disease, hypertension, hyperuricemia, osteoporosis, and prostate cancer), and the treatment responses of patients.

Enhancer polymorphisms at the IKZF1 susceptibility locus for acute lymphoblastic leukemia impact B-cell proliferation and differentiation in both Down syndrome and non-Down syndrome genetic backgrounds

Children with Down syndrome have an approximately 10-fold increased risk of developing acute lymphoblastic leukemia and this risk is influenced by inherited genetic variation. Genome-wide association studies have identified IKZF1 as a strong acute lymphoblastic leukemia susceptibility locus in children both with and without Down syndrome, with association signals reported at rs4132601 in non-Down syndrome and rs58923657 in individuals with Down syndrome (r² = 0.98 for these two loci). Expression quantitative trait locus analysis in non-Down syndrome lymphoblastoid cell lines has demonstrated an association between the rs4132601 risk allele and decreased IKZF1 mRNA levels. In this study, we provide further mechanistic evidence linking the region encompassing IKZF1-associated polymorphisms to pro-leukemogenic effects in both human lymphoblastoid cell lines and murine hematopoietic stem cells. CRISPR/Cas9-mediated deletion of the region encompassing the rs17133807 major allele (r² with rs58923657 = 0.97) resulted in significant reduction of IKZF1 mRNA levels in lymphoblastoid cell lines, with a greater effect in Down syndrome versus non-Down syndrome cells. Since rs17133807 is highly conserved in mammals, we also evaluated the orthologous murine locus at rs263378223, in hematopoietic stem cells from the Dp16(1)Yey mouse model of Down syndrome as well as non-Down syndrome control mice. Homozygous deletion of the region encompassing rs263378223 resulted in significantly reduced Ikzf1 mRNA, confirming that this polymorphism maps to a strong murine Ikzf1 enhancer, and resulted in increased B-lymphoid colony growth and decreased B-lineage differentiation. Our results suggest that both the region encompassing rs17133807 and its conserved orthologous mouse locus have functional effects that may mediate increased leukemia susceptibility in both the Down syndrome and non-Down syndrome genetic backgrounds.

Functional genomics of autoimmune diseases

For more than a decade, genome-wide association studies have been applied to autoimmune diseases and have expanded our understanding on the pathogeneses. Genetic risk factors associated with diseases and traits are essentially causative. However, elucidation of the biological mechanism of disease from genetic factors is challenging. In fact, it is difficult to identify the causal variant among multiple variants located on the same haplotype or linkage disequilibrium block and thus the responsible biological genes remain elusive. Recently, multiple studies have revealed that the majority of risk variants locate in the non-coding region of the genome and they are the most likely to regulate gene expression such as quantitative trait loci. Enhancer, promoter and long non-coding RNA appear to be the main target mechanisms of the risk variants. In this review, we discuss functional genetics to challenge these puzzles.

Leveraging Genome and Phenome-Wide Association Studies to Investigate Genetic Risk of Acute Lymphoblastic Leukemia

BACKGROUND

Genome-wide association studies (GWAS) of childhood cancers remain limited, highlighting the need for novel analytic strategies. We describe a hybrid GWAS and phenome-wide association study (PheWAS) approach to uncover genotype-phenotype relationships and candidate risk loci, applying it to acute lymphoblastic leukemia (ALL).

METHODS

PheWAS was performed for 12 ALL SNPs identified by prior GWAS and two control SNP-sets using UK Biobank data. PheWAS-traits significantly associated with ALL SNPs compared with control SNPs were assessed for association with ALL risk (959 cases, 2,624 controls) using polygenic score and Mendelian randomization analyses. Trait-associated SNPs were tested for association with ALL risk in single-SNP analyses, with replication in an independent case-control dataset (1,618 cases, 9,409 controls).

RESULTS

Platelet count was the trait most enriched for association with known ALL risk loci. A polygenic score for platelet count (223 SNPs) was not associated with ALL risk (P = 0.82) and Mendelian randomization did not suggest a causal relationship. However, twelve platelet count-associated SNPs were nominally associated with ALL risk in COG data and three were replicated in UK data (rs10058074, rs210142, rs2836441).

CONCLUSIONS

In our hybrid GWAS-PheWAS approach, we identify pleiotropic genetic variation contributing to ALL risk and platelet count. Three SNPs known to influence platelet count were reproducibly associated with ALL risk, implicating genomic regions containing IRF1, proapoptotic protein BAK1, and ERG in platelet production and leukemogenesis.

IMPACT

Incorporating PheWAS data into association studies can leverage genetic pleiotropy to identify cancer risk loci, highlighting the utility of our novel approach.

Inherited genetic susceptibility to acute lymphoblastic leukemia in Down syndrome

Children with Down syndrome (DS) have a 20-fold increased risk of acute lymphoblastic leukemia (ALL) and distinct somatic features, including CRLF2 rearrangement in ∼50% of cases; however, the role of inherited genetic variation in DS-ALL susceptibility is unknown. We report the first genome-wide association study of DS-ALL, comprising a meta-analysis of 4 independent studies, with 542 DS-ALL cases and 1192 DS controls. We identified 4 susceptibility loci at genome-wide significance: rs58923657 near IKZF1 (odds ratio [OR], 2.02; Pmeta = 5.32 × 10-15), rs3731249 in CDKN2A (OR, 3.63; Pmeta = 3.91 × 10-10), rs7090445 in ARID5B (OR, 1.60; Pmeta = 8.44 × 10-9), and rs3781093 in GATA3 (OR, 1.73; Pmeta = 2.89 × 10-8). We performed DS-ALL vs non-DS ALL case-case analyses, comparing risk allele frequencies at these and other established susceptibility loci (BMI1, PIP4K2A, and CEBPE) and found significant association with DS status for CDKN2A (OR, 1.58; Pmeta = 4.1 × 10-4). This association was maintained in separate regression models, both adjusting for and stratifying on CRLF2 overexpression and other molecular subgroups, indicating an increased penetrance of CDKN2A risk alleles in children with DS. Finally, we investigated functional significance of the IKZF1 risk locus, and demonstrated mapping to a B-cell super-enhancer, and risk allele association with decreased enhancer activity and differential protein binding. IKZF1 knockdown resulted in significantly higher proliferation in DS than non-DS lymphoblastoid cell lines. Our findings demonstrate a higher penetrance of the CDKN2A risk locus in DS and serve as a basis for further biological insights into DS-ALL etiology.

Predisposing germline mutations in high hyperdiploid acute lymphoblastic leukemia in children

High hyperdiploidy (HD) is the most common cytogenetic subtype of childhood acute lymphoblastic leukemia (ALL), and a higher incidence of HD has been reported in ALL patients with congenital cancer syndromes. We assessed the frequency of predisposing germline mutations in 57 HD-ALL patients from the California Childhood Leukemia Study via targeted sequencing of cancer-relevant genes. Three out of 57 patients (5.3%) harbored confirmed germline mutations that were likely causal, in NBN, ETV6, and FLT3, with an additional six patients (10.5%) harboring putative predisposing mutations that were rare in unselected individuals (<0.01% allele frequency in the Exome Aggregation Consortium, ExAC) and predicted functional (scaled CADD score ≥ 20) in known or potential ALL predisposition genes (SH2B3, CREBBP, PMS2, MLL, ABL1, and MYH9). Three additional patients carried rare and predicted damaging germline mutations in GAB2, a known activator of the ERK/MAPK and PI3K/AKT pathways and binding partner of PTPN11-encoded SHP2. The frequency of rare and predicted functional germline GAB2 mutations was significantly higher in our patients (2.6%) than in ExAC (0.28%, P = 4.4 × 10^-3 ), an observation that was replicated in ALL patients from the TARGET project (P = .034). We cloned patient GAB2 mutations and expressed mutant proteins in HEK293 cells and found that frameshift mutation P621fs led to reduced SHP2 binding and ERK1/2 phosphorylation but significantly increased AKT phosphorylation, suggesting possible RAS-independent leukemogenic effects. Our results support a significant contribution of rare, high penetrance germline mutations to HD-ALL etiology, and pinpoint GAB2 as a putative novel ALL predisposition gene.