Genetic and regulatory mechanism of susceptibility to high-hyperdiploid acute lymphoblastic leukaemia at 10p21.2

Discovery of Cis-Regulatory Mechanisms via Non-Coding Mutations in Acute Lymphoblastic Leukemia

The non-coding genome, constituting 98% of human DNA, remains largely unexplored, yet holds potential for identifying new biomarkers and therapeutic targets in acute lymphoblastic leukemia (ALL). In this study, we conducted a systematic analysis of recurrent somatic non-coding single nucleotide variants (SNVs) in pediatric B-cell precursor (BCP) ALL. We leveraged whole genome sequencing (WGS) data from 345 pediatric BCP ALL cases, representing all major genetic subtypes and identified 346 mutational hotspots that harbored somatic SNVs in at least three cases. Through the integration of paired RNA sequencing along with published ChIP-seq and ATAC-seq data, we found 128 non-coding hotspots associated with differentially expressed genes nearby, which were enriched for cis-regulatory elements, demonstrating the effectiveness of multi-omics integration in distinguishing pathogenic mutations from passengers. We identified one mutational hotspot that was associated with increased expression of the leukemia-associated gene NRAS in three primary ALLs. Micro-C analysis in the leukemia cell line demonstrated interactions between the hotspot region and NRAS regulatory elements. Dual luciferase assays indicated that the mutations disrupted regulatory interactions and CRISPR-mediated deletion of the region significantly upregulated NRAS, confirming the hypothesized regulatory link. Altogether, we provide new insights into the functional roles of non-coding mutations in 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.

ARID5B is a negative modulator of IL-6 production in rheumatoid arthritis synovial fibroblasts

Recent single-cell RNA-sequencing analysis of rheumatoid arthritis (RA) synovial tissues revealed the heterogeneity of RA synovial fibroblasts (SFs) with distinct functions such as high IL-6 production. The molecular mechanisms responsible for high IL-6 production will become a promising drug target of RASFs to treat RA. In this study, we performed siRNA screening of 65 transcription factors (TFs) differentially expressed among RASF subsets to identify TFs involved in IL-6 production. The siRNA screening identified 7 TFs including ARID5B, a RA risk gene, that affected IL-6 production. Both long and short isoforms of ARID5B were expressed and negatively regulated by TNF-α in RASFs. The siRNA knockdown and lentiviral overexpression of long and short isoforms of ARID5B revealed that the long isoform suppressed IL-6 production stimulated with TNF-α. eQTL analysis using 58 SFs demonstrated that RA risk allele, rs10821944, in intron 4 of the ARID5B gene had a trend of eQTL effects to the expression of long isoform of ARID5B in SFs treated with TNF-α. ARID5B was found to be a negative modulator of IL-6 production in RASFs. The RA risk allele of ARID5B intron may cause high IL-6 production, suggesting that ARID5B will become a promising drug target to treat RA.

Constitutional and acquired genetic variants in ARID5B in pediatric B-cell precursor acute lymphoblastic leukemia

Constitutional polymorphisms in ARID5B are associated with an increased risk of developing high hyperdiploid (HeH; 51-67 chromosomes) pediatric B-cell precursor acute lymphoblastic leukemia (BCP ALL). Here, we investigated constitutional and somatic ARID5B variants in 1335 BCP ALL cases from five different cohorts, with a particular focus on HeH cases. In 353 HeH ALL that were heterozygous for risk alleles and trisomic for chromosome 10, where ARID5B is located, a significantly higher proportion of risk allele duplication was seen for the SNPs rs7090445 (p = 0.009), rs7089424 (p = 0.005), rs7073837 (p = 0.03), and rs10740055 (p = 0.04). Somatic ARID5B deletions were seen in 16/1335 cases (1.2%), being more common in HeH than in other genetic subtypes (2.2% vs. 0.4%; p = 0.002). The expression of ARID5B in HeH cases with genomic deletions was reduced, consistent with a functional role in leukemogenesis. Whole-genome sequencing and RNA-sequencing in HeH revealed additional somatic events involving ARID5B, resulting in a total frequency of 3.6% of HeH cases displaying a somatic ARID5B aberration. Overall, our results show that both constitutional and somatic events in ARID5B are involved in the leukemogenesis of pediatric BCP ALL, particularly in the HeH subtype.

Investigation of inherited noncoding genetic variation impacting the pharmacogenomics of childhood acute lymphoblastic leukemia treatment

Defining genetic factors impacting chemotherapy failure can help to better predict response and identify drug resistance mechanisms. However, there is limited understanding of the contribution of inherited noncoding genetic variation on inter-individual differences in chemotherapy response in childhood acute lymphoblastic leukemia (ALL). Here we map inherited noncoding variants associated with treatment outcome and/or chemotherapeutic drug resistance to ALL cis-regulatory elements and investigate their gene regulatory potential and target gene connectivity using massively parallel reporter assays and three-dimensional chromatin looping assays, respectively. We identify 54 variants with transcriptional effects and high-confidence gene connectivity. Additionally, functional interrogation of the top variant, rs1247117, reveals changes in chromatin accessibility, PU.1 binding affinity and gene expression, and deletion of the genomic interval containing rs1247117 sensitizes cells to vincristine. Together, these data demonstrate that noncoding regulatory variants associated with diverse pharmacological traits harbor significant effects on allele-specific transcriptional activity and impact sensitivity to antileukemic agents.

Identification of an alternative short ARID5B isoform associated with B-ALL survival

Utilizing RNA sequence (RNA-Seq) splice junction data from a cohort of 1841 B-cell acute lymphoblastic leukemia (B-ALL) patients we define transcriptionally distinct isoforms of ARID5B, a risk-associated gene identified in genome wide association studies (GWAS), which associate with disease survival. Short (S) and long (L) ARID5B transcripts, which differ in an encoded BAH-like chromatin interaction domain, show remarkable correlation to the isoform splicing pattern. Testing of the ARID5B proximal promoter of the S & L isoforms indicated that both are functionally independent in luciferase reporter assays. Increased short isoform expression is associated with decreased event-free and overall survival. The abundance of short and long transcripts strongly correlates to B-ALL prognostic stratification, where B-ALL subtypes with poor outcomes express a higher proportion of the S-isoform. These data demonstrate that the analysis of independent promoters and alternative splicing events are essential for improved risk stratification and a more complete understanding of disease pathology.

Genome-wide assessment of genetic risk loci for childhood acute lymphoblastic leukemia in Japanese patients

Not available.

Chromatin accessibility landscape of relapsed pediatric B-lineage acute lymphoblastic leukemia

For around half of the pediatric B-lineage acute lymphoblastic leukemia (B-ALL) patients, the molecular mechanism of relapse remains unclear. To fill this gap in knowledge, here we characterize the chromatin accessibility landscape in pediatric relapsed B-ALL. We observe rewired accessible chromatin regions (ACRs) associated with transcription dysregulation in leukemia cells as compared with normal B-cell progenitors. We show that over a quarter of the ACRs in B-ALL are in quiescent regions with high heterogeneity among B-ALLs. We identify subtype-specific and allele-imbalanced chromatin accessibility by integrating multi-omics data. By characterizing the differential ACRs between diagnosis and relapse in B-ALL, we identify alterations in chromatin accessibility during drug treatment. Further analysis of ACRs associated with relapse free survival leads to the identification of a subgroup of B-ALL which show early relapse. These data provide an advanced and integrative portrait of the importance of chromatin accessibility alterations in tumorigenesis and drug responses.

ARID5B regulates fatty acid metabolism and proliferation at the Pre-B cell stage during B cell development

During B cell development in bone marrow, large precursor B cells (large Pre-B cells) proliferate rapidly, exit the cell cycle, and differentiate into non-proliferative (quiescent) small Pre-B cells. Dysregulation of this process may result in the failure to produce functional B cells and pose a risk of leukemic transformation. Here, we report that AT rich interacting domain 5B (ARID5B), a B cell acute lymphoblastic leukemia (B-ALL) risk gene, regulates B cell development at the Pre-B stage. In both mice and humans, we observed a significant upregulation of ARID5B expression that initiates at the Pre-B stage and is maintained throughout later stages of B cell development. In mice, deletion of Arid5b in vivo and ex vivo exhibited a significant reduction in the proportion of immature B cells but an increase in large and small Pre-B cells. Arid5b inhibition ex vivo also led to an increase in proliferation of both Pre-B cell populations. Metabolic studies in mouse and human bone marrow revealed that fatty acid uptake peaked in proliferative B cells then decreased during non-proliferative stages. We showed that Arid5b ablation enhanced fatty acid uptake and oxidation in Pre-B cells. Furthermore, decreased ARID5B expression was observed in tumor cells from B-ALL patients when compared to B cells from non-leukemic individuals. In B-ALL patients, ARID5B expression below the median was associated with decreased survival particularly in subtypes originating from Pre-B cells. Collectively, our data indicated that Arid5b regulates fatty acid metabolism and proliferation of Pre-B cells in mice, and reduced expression of ARID5B in humans is a risk factor for B cell leukemia.

Functional investigation of inherited noncoding genetic variation impacting the pharmacogenomics of childhood acute lymphoblastic leukemia treatment

Although acute lymphoblastic leukemia (ALL) is the most common childhood cancer, there is limited understanding of the contribution of inherited genetic variation on inter-individual differences in chemotherapy response. Defining genetic factors impacting therapy failure can help better predict response and identify drug resistance mechanisms. We therefore mapped inherited noncoding variants associated with chemotherapeutic drug resistance and/or treatment outcome to ALL cis-regulatory elements and investigated their gene regulatory potential and genomic connectivity using massively parallel reporter assays and promoter capture Hi-C, respectively. We identified 53 variants with reproducible allele-specific effects on transcription and high-confidence gene targets. Subsequent functional interrogation of the top variant (rs1247117) determined that it disrupted a PU.1 consensus motif and PU.1 binding affinity. Importantly, deletion of the genomic interval containing rs1247117 sensitized ALL cells to vincristine. Together, these data demonstrate that noncoding regulatory variation associated with diverse pharmacological traits harbor significant effects on allele-specific transcriptional activity and impact sensitivity to chemotherapeutic agents in ALL.

Association of two ARID5B gene variant single nucleotide polymorphisms with acute lymphoblastic leukemia in the Egyptian population

BACKGROUND

ARID5B SNPs have been linked to ALL in many research studies in which it was identified as a risk factor. From this context, we had great interest to investigate the relationship between ARID5B rs4948488 and ARID5B rs2893881 genotypes and ALL susceptibility and relapse in this study.

MATERIALS AND METHODS

Peripheral blood mononuclear cells were analyzed for ARID5B rs4948488 and rs2893881 gene polymorphisms by real-time quantitative polymerase chain reaction in 80 ALL patients and 80 controls.

RESULTS

Our results showed that the C/C genotype of ARID5B rs4948488 and A/G genotype and G-allele of rs2893881 were linked to higher ALL incidence. Regarding the relapse of ALL, rs4948488 C/C genotype and C-alleles were significantly associated with relapse of ALL. Meanwhile, rs4948488 C/C genotype and rs2893881 A/A genotype and A-allele are associated with T-ALL, while rs2893881 A/G genotype and G-allele are associated with B-ALL.

CONCLUSION

The results of our study suggested that ARID5B rs4948488 and rs2893881 SNPs might be used risk factors for genetic susceptibility for B-ALL and T-ALL, and that ARID5B s4948488 is related to relapse in ALL patients.<br />.

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.

Molecular Mechanisms of ARID5B-Mediated Genetic Susceptibility to Acute Lymphoblastic Leukemia

BACKGROUND

There is growing evidence for the inherited basis of susceptibility to childhood acute lymphoblastic leukemia (ALL). Genome-wide association studies have identified non-coding ALL risk variants at the ARID5B gene locus, but their exact functional effects and the molecular mechanism linking ARID5B to B-cell ALL leukemogenesis remain largely unknown.

METHODS

We performed targeted sequencing of ARID5B in germline DNA of 5008 children with ALL. Variants were evaluated for association with ALL susceptibility using 3644 patients from the UK10K cohort as non-ALL controls, under an additive model. Cis-regulatory elements in ARID5B were systematically identified using dCas9-KRAB-mediated enhancer interference system enhancer screen in ALL cells. Disruption of transcription factor binding by ARID5B variant was predicted informatically and then confirmed using chromatin immunoprecipitation and coimmunoprecipitation. ARID5B variant association with hematological traits was examined using UK Biobank dataset. All statistical tests were 2-sided.

RESULTS

We identified 54 common variants in ARID5B statistically significantly associated with leukemia risk, all of which were noncoding. Six cis-regulatory elements at the ARID5B locus were discovered using CRISPR-based high-throughput enhancer screening. Strikingly, the top ALL risk variant (rs7090445, P = 5.57 × 10-45) is located precisely within the strongest enhancer element, which is also distally tethered to the ARID5B promoter. The variant allele disrupts the MEF2C binding motif sequence, resulting in reduced MEF2C affinity and decreased local chromosome accessibility. MEF2C influences ARID5B expression in ALL, likely via a transcription factor complex with RUNX1. Using the UK Biobank dataset (n = 349 861), we showed that rs7090445 was also associated with lymphocyte percentage and count in the general population (P = 8.6 × 10-22 and 2.1 × 10-18, respectively).

CONCLUSIONS

Our results indicate that ALL risk variants in ARID5B function by modulating cis-regulatory elements at this locus.

ARID5B influences B-cell development and function in mouse

There is growing evidence for an inherited basis of susceptibility to childhood acute lymphoblastic leukemia. Genomewide association studies by us and others have identified non-coding acute lymphoblastic leukemia risk variants at the ARID5B gene locus, but the molecular mechanisms linking ARID5B to normal and malignant hematopoiesis remain largely unknown. Using a Vav1-driven transgenic mouse model, we characterized the role of Arid5b in hematopoiesis in vivo. Arid5b overexpression resulted in a dramatic reduction in the proportion of circulating B cells, immature, and mature Bcell fractions in the peripheral blood and the bone marrow, and also a decrease of follicular B cells in the spleen. There were significant defects in B-cell activation upon Arid5b overexpression in vitro with hyperactivation of B-cell receptor signaling at baseline. In addition, increased mitochondrial oxygen consumption rate of naïve or stimulated B cells of Arid5b-overexpressing mice was observed, compared to the rate of wild-type counterparts. Taken together, our results indicate that ARID5B may play an important role in B-cell development and function.

Impact of Variants in the ATIC and ARID5B Genes on Therapeutic Failure with Imatinib in Patients with Chronic Myeloid Leukemia

Chronic myeloid leukemia (CML) is a myeloproliferative neoplasm derived from the balanced reciprocal translocation of chromosomes 9 and 22 t (9q34 and 22q11), which leads to the formation of the Philadelphia chromosome and fusion of the BCR-ABL genes. The first-line treatment for CML is imatinib, a tyrosine kinase inhibitor that acts on the BCR-ABL protein. However, even though it is a target-specific drug, about 25% of patients do not respond to this treatment. The resistance mechanisms involved in this process have been investigated and studies have shown that germinal alterations can influence this mechanism. The aim of this work was to investigate 32 polymorphisms in 24 genes of carcinogenic pathway to verify the influence of these genetic variants on the response to treatment with imatinib. Our results demonstrated that individuals with the recessive GG genotype for the rs2372536 variant in the ATIC gene are approximately three times more likely to experience treatment failure with imatinib (p = 0.045, HR = 2.726, 95% CI = 0.9986–7.441), as well as individuals with the TT genotype for the rs10821936 variant in the ARID5B gene, who also have a higher risk for treatment failure with imatinib over time (p = 0.02, HR = 0.4053, IC 95% = 0.1802–0.911). In conclusion, we show that variants in the ATIC and ARIDB5 gene, never screened in previous studies, could potentially influence the therapeutic response to imatinib in patients treated for CML.

Toward prevention of childhood ALL by early-life immune training

B-cell precursor acute lymphoblastic leukemia (BCP-ALL) is the most common form of childhood cancer. Chemotherapy is associated with life-long health sequelae and fails in ∼20% of cases. Thus, prevention of leukemia would be preferable to treatment. Childhood leukemia frequently starts before birth, during fetal hematopoiesis. A first genetic hit (eg, the ETV6-RUNX1 gene fusion) leads to the expansion of preleukemic B-cell clones, which are detectable in healthy newborn cord blood (up to 5%). These preleukemic clones give rise to clinically overt leukemia in only ∼0.2% of carriers. Experimental evidence suggests that a major driver of conversion from the preleukemic to the leukemic state is exposure to immune challenges. Novel insights have shed light on immune host responses and how they shape the complex interplay between (1) inherited or acquired genetic predispositions, (2) exposure to infection, and (3) abnormal cytokine release from immunologically untrained cells. Here, we integrate the recently emerging concept of "trained immunity" into existing models of childhood BCP-ALL and suggest future avenues toward leukemia prevention.

Chronic lymphocytic leukemia (CLL) risk is mediated by multiple enhancer variants within CLL risk loci

Chronic lymphocytic leukemia (CLL) is the most common adult leukemia in Western countries. It has a strong genetic basis, showing a ~ 8-fold increased risk of CLL in first-degree relatives. Genome-wide association studies (GWAS) have identified 41 risk variants across 41 loci. However, for a majority of the loci, the functional variants and the mechanisms underlying their causal roles remain undefined. Here, we examined the genetic and epigenetic features associated with 12 index variants, along with any correlated (r2 ≥ 0.5) variants, at the CLL risk loci located outside of gene promoters. Based on publicly available ChIP-seq and chromatin accessibility data as well as our own ChIP-seq data from CLL patients, we identified six candidate functional variants at six loci and at least two candidate functional variants at each of the remaining six loci. The functional variants are predominantly located within enhancers or super-enhancers, including bi-directionally transcribed enhancers, which are often restricted to immune cell types. Furthermore, we found that, at 78% of the functional variants, the alternative alleles altered the transcription factor binding motifs or histone modifications, indicating the involvement of these variants in the change of local chromatin state. Finally, the enhancers carrying functional variants physically interacted with genes enriched in the type I interferon signaling pathway, apoptosis, or TP53 network that are known to play key roles in CLL. These results support the regulatory roles for inherited noncoding variants in the pathogenesis of CLL.

Somatic Sex: On the Origin of Neoplasms With Chromosome Counts in Uneven Ploidy Ranges

Stable aneuploid genomes with nonrandom numerical changes in uneven ploidy ranges define distinct subsets of hematologic malignancies and solid tumors. The idea put forward herein suggests that they emerge from interactions between diploid mitotic and G0/G1 cells, which can in a single step produce all combinations of mono-, di-, tri-, tetra- and pentasomic paternal/maternal homologue configurations that define such genomes. A nanotube-mediated influx of interphase cell cytoplasm into mitotic cells would thus be responsible for the critical nondisjunction and segregation errors by physically impeding the proper formation of the cell division machinery, whereas only a complete cell fusion can simultaneously generate pentasomies, uniparental trisomies as well as biclonal hypo- and hyperdiploid cell populations. The term "somatic sex" was devised to accentuate the similarities between germ cell and somatic cell fusions. A somatic cell fusion, in particular, recapitulates many processes that are also instrumental in the formation of an abnormal zygote that involves a diploid oocyte and a haploid sperm, which then may further develop into a digynic triploid embryo. Despite their somehow deceptive differences and consequences, the resemblance of these two routes may go far beyond of what has hitherto been appreciated. Based on the arguments put forward herein, I propose that embryonic malignancies of mesenchymal origin with these particular types of aneuploidies can thus be viewed as the kind of flawed somatic equivalent of a digynic triploid embryo.

Identification of consistent post-translational regulatory triplets related to oncogenic and tumour suppressive modulators in childhood acute lymphoblastic leukemia

Background

Acute lymphoblastic leukemia (ALL) is the most common type of childhood cancer. It can be caused by mutations that turn on oncogenes or turn off tumour suppressor genes. For instance, changes in certain genes including Rb and p53 are common in ALL cells. Oncogenes and TSGs may serve as a modulator gene to regulate the gene expression level via their respective target genes. To investigate the regulatory relationship between oncogenes, tumour suppressor genes and transcription factors at the post translational level in childhood ALL, we performed an integrative network analysis on the gene regulation in the post-translational level for childhood ALL based on many publicly available cancer gene expression data including TARGET and GEO database.

Methods

We collected 259 childhood ALL-related genes from the latest online leukemia database, Leukemia Gene Literature Database. These 259 genes were selected from a comprehensive systematic literature with experimental evidences. The identified and curated genes were also associated with patient survival cases and we incorporated this pediatric ALL-related gene list into our analysis. We extracted the known human TFs from the TRRUST database. Among 259 childhood ALL-related genes, 101 unique regulators were mapped to the list of oncogene and tumour suppressor genes (TSGs) from the ONGene and the TSGene databases, and these included 74 TSGs, 62 oncogenes and 46 TF genes.

Results

The resulted regulation was presented as a hierarchical regulatory network with transcription factors (TFs) as intermediate regulators connecting the top modulators (oncogene and TSGs) to the common target genes. Cross-validation was applied to the results from the TARGET dataset by identifying the consistent regulatory motifs based on three independent ALL expression datasets. A three-layer regulatory network of consistent positive modulators in childhood ALL was constructed in which 74 modulators (40 oncogenes, 34 TSGs) are considered as the most important regulators. The middle layer and the bottom layer contain 34 TFs and 176 target genes, respectively. Oncogenes mostly participated in positive regulation of gene expression and the transcription process of RNA II polymerase, while TSGs were mainly involved in the negative regulation of gene expression. In addition, the oncogene-specific targets were enriched with regulators of the MAPK cascade while tumour suppressor-specific targets were associated with cell death.

Conclusion

The results revealed that oncogenes and TSGs possess a different functional regulatory pattern with regard to not only their biological functions but also their specific target genes in childhood ALL cancer progression. Taken together, our findings could contribute to a better understanding of the important regulatory mechanisms and this method could be used to analyse the targeted genes at the post-translational level in childhood ALL through integrative network analysis.

Germline predisposition to haematopoietic malignancies

Once thought to be exceedingly rare, the advent of next-generation sequencing has revealed a plethora of germline predisposition disorders that confer risk for haematopoietic malignancies (HMs). These syndromes are now recognized to be much more common than previously thought. The recognition of a germline susceptibility risk allele in an individual impacts the clinical management and health surveillance strategies in the index patient and relatives who share the causative DNA variant. Challenges to accurate clinical testing include a lack of familiarity in many health care providers, the requirement for DNA samples that reasonably approximate the germline state, and a lack of standardization among diagnostic platforms as to which genes are sequenced and their capabilities in detecting the full range of variant types that confer risk. Current knowledge gaps include a comprehensive understanding of all predisposition genes; whether scenarios exist in which an allogeneic stem cell transplant using donor haematopoietic stem cells with deleterious variants is permissive; and effective means of delivering genetic counseling and results disclosure for these conditions. We are hopeful that comprehensive germline genetic testing, universal germline testing for all patients with an HM, universal germline testing for allogeneic haematopoietic stem cell donors, and the development of preventive strategies to delay or even prevent malignancies will be available in the near future. These factors will likely contribute to improved health outcomes for at-risk individuals and their family members.

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.

Association between apurinic/apyrimidinic endonuclease 1 rs1760944 T>G polymorphism and susceptibility of cancer: a meta-analysis involving 21764 subjects

BACKGROUND

Previous case-control studies have suggested that apurinic/apyrimidinic endonuclease 1 (APE1) rs1760944 T>G polymorphism may be associated with cancer risk. Here, we carried out an updated meta-analysis to focus on the correlation between APE1 rs1760944 T>G locus and the risk of cancer.

METHODS

We used the crude odds ratios (ORs) with their 95% confidence intervals (CIs) to evaluate the possible relationship between the APE1 rs1760944 T>G polymorphism and cancer risk. Heterogeneity, publication bias and sensitivity analysis were also harnessed to check the potential bias of the present study.

RESULTS

Twenty-three independent studies involving 10166 cancer cases and 11598 controls were eligible for this pooled analysis. We found that APE1 rs1760944 T>G polymorphism decreased the risk of cancer in four genetic models (G vs. T: OR, 0.87; 95% CI, 0.83-0.92; P<0.001; GG vs. TT: OR, 0.77; 95% CI, 0.69-0.86; P<0.001; GG/TG vs. TT: OR, 0.83; 95% CI, 0.77-0.89, P<0.001 and GG vs. TT/TG: OR, 0.85; 95% CI, 0.80-0.92, P<0.001). Results of subgroup analyses also demonstrated that this single-nucleotide polymorphism (SNP) modified the risk among lung cancer, breast cancer, osteosarcoma, and Asians. Evidence of publication bias was found in the present study. When we treated the publication bias with 'trim-and-fill' method, the adjusted ORs and CIs were not significantly changed.

CONCLUSION

In conclusion, current evidence highlights that the APE1 rs1760944 T>G polymorphism is a protective factor for cancer susceptibility. In the future, case-control studies with detailed risk factors are needed to confirm or refute our findings.

Association of relapse-linked ARID5B single nucleotide polymorphisms with drug resistance in B-cell precursor acute lymphoblastic leukemia cell lines

Background

The genetic variants of the ARID5B gene have recently been reported to be associated with disease susceptibility and treatment outcome in childhood acute lymphoblastic leukemia (ALL). However, few studies have explored the association of ARID5B with sensitivities to chemotherapeutic agents.

Methods

We genotyped susceptibility-linked rs7923074 and rs10821936 as well as relapse-linked rs4948488, rs2893881, and rs6479778 of ARDI5B by direct sequencing of polymerase chain reaction (PCR) products in 72 B-cell precursor-ALL (BCP-ALL) cell lines established from Japanese patients. We also quantified their ARID5B expression levels by real-time reverse transcription PCR, and determined their 50% inhibitory concentration (IC50) values by alamarBlue assays in nine representative chemotherapeutic agents used for ALL treatment.

Results

No significant associations were observed in genotypes of the susceptibility-linked single nucleotide polymorphisms (SNPs) and the relapsed-linked SNPs with ARID5B gene expression levels. Of note, IC50 values of vincristine (VCR) (median IC50: 39.6 ng/ml) in 12 cell lines with homozygous genotype of risk allele (C) in the relapse-linked rs4948488 were significantly higher (p = 0.031 in Mann–Whitney U test) than those (1.04 ng/ml) in 60 cell lines with heterozygous or homozygous genotypes of the non-risk allele (T). Furthermore, the IC50 values of mafosfamide [Maf; active metabolite of cyclophosphamide (CY)] and cytarabine (AraC) tended to be associated with the genotype of rs4948488. Similar associations were observed in genotypes of the relapse-linked rs2893881 and rs6479778, but not in those of the susceptibility-linked rs7923074 and rs10821936. In addition, the IC50 values of methotrexate (MTX) were significantly higher (p = 0.023) in 36 cell lines with lower ARID5B gene expression (median IC50: 37.1 ng/ml) than those in the other 36 cell lines with higher expression (16.9 ng/ml).

Conclusion

These observations in 72 BCP-ALL cell lines suggested that the risk allele of the relapse-linked SNPs of ARID5B may be involved in a higher relapse rate because of resistance to chemotherapeutic agents such as VCR, CY, and AraC. In addition, lower ARID5B gene expression may be associated with MTX resistance.

Validations of Top and Novel Susceptibility Variants in All-Age Chinese Patients With Acute Lymphoblastic Leukemia

Through genome-wide association studies (GWAS), multiple inherited predispositions to acute lymphoblastic leukemia (ALL) have been identified in children. Most recently, a novel susceptibility locus at ERG was localized, exhibiting Hispanic-specific manner. In this study, we conducted a replication study to in all-age Chinese patients (N = 451), not only validating the novel ERG locus, but also systematically determining the impact of age on association status of the top GWAS signals. We found that single nucleotide polymorphisms at ARID5B, IKZF1, CEBPE, PIP4K2A were only significantly associated with ALL susceptibility in childhood patients with no BCR-ABL fusion, while GATA3 signal exhibited its significance in adults no matter carrying BCR-ABL fusion or not. Moreover, the novel ERG SNP can be validated in pediatric patients without both BCR-ABL and ETV6-RUNX1 fusion. Our finding suggests the modifying effects of age on genetic predisposition to ALL, and highlights the impact of ERG SNP in Chinese patients.

The Role of ARID5B in Acute Lymphoblastic Leukemia and Beyond

Acute lymphoblastic leukemia (ALL) is the most common malignancy in children with distinct characteristics among different subtypes. Although the etiology of ALL has not been fully unveiled, initiation of ALL has been demonstrated to partly depend on genetic factors. As indicated by several genome wide association studies (GWASs) and candidate gene analyses, ARID5B, a member of AT-rich interactive domain (ARID) protein family, is associated with the occurrence and prognosis of ALL. However, the mechanisms by which ARID5B genotype impact on the susceptibility and treatment outcome remain vague. In this review, we outline developments in the understanding of ARID5B in the susceptibility of ALL and its therapeutic perspectives, and summarize the underlying mechanisms based on the limited functional studies, hoping to illustrate the possible mechanisms of ARID5B impact and highlight the potential treatment regimens.

Novel phenotypes observed in patients with ETV6-linked leukaemia/familial thrombocytopenia syndrome and a biallelic ARID5B risk allele as leukaemogenic cofactor

Background. The phenotypes of patients with the recently discovered, dominant, ETV6-linked leukaemia predisposition and familial thrombocytopenia syndrome are variable, and the exact mechanism of leukaemogenesis remains unclear. Patients and Methods. Here, we present novel clinical and laboratory phenotypes of seven individuals from three families with ETV6 germline mutations and a refined genetic analysis of one child with additional high-hyperdiploid acute lymphoblastic leukaemia (HD-ALL), aiming to elucidate second oncogenic hits. Results. Four individuals from two pedigrees harboured one novel or one previously described variant in the central domain of ETV6 (c.592C>T, p.Gln198* or c.641C>T, p.Pro241Leu, respectively). Neutropenia was an accompanying feature in one of these families that also harboured a variant in RUNX1 (c.1098_1103dup, p.Ile366_Gly367dup), while in the other, an autism-spectrum disorder was observed. In the third family, the index patient suffered from HD-ALL and life-threatening pulmonary mucor mycosis, and had a positive family history of 'immune' thrombocytopenia. Genetic analyses revealed a novel heterozygous mutation in the ETS domain of ETV6 (c.1136T>C, p.Leu379Pro) along with absence of heterozygosity of chromosome (10)(q21.2q21.3), yielding a biallelic leukaemia risk allele in ARID5B (rs7090445-C). The neutrophil function was normal in all individuals tested, and the platelet immune histochemistry of all three pedigrees showed delta-storage-pool defect-like features and cytoskeletal defects. Conclusions. Our clinical observations and results of high-resolution genetic analyses extend the spectrum of possible phenotypes cosegregating with ETV6 germline mutations. Further, we propose ARID5B as potential leukaemogenic cofactor in patients with ETV6-linked leukaemia predisposition and familial thrombocytopenia syndrome.

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.

Variants in ARID5B gene are associated with the development of acute lymphoblastic leukemia in Mexican children

A high impact of ARID5B SNPs on acute lymphoblastic leukemia (ALL) susceptibility has been described in Hispanic children; therefore, it is relevant to know if they influence the high incidence of childhood-ALL in Mexicans. Seven SNPs (rs10821936, rs10994982, rs7089424, rs2393732, rs2393782, rs2893881, rs4948488) of ARID5B were analyzed in 384 controls and 298 ALL children using genomic DNA and TaqMan probes. The SNPs were analyzed for deviation of Hardy-Weinberg equilibrium; Fisher's exact test was used to compare the genotypic and allelic frequencies between controls and patients. The association between SNPs and ALL susceptibility was calculated, and haplotype and ancestry analyses were conducted. All SNPs were associated with ALL, pre-B ALL, and hyperdiploid-ALL susceptibility (p < 0.05). No association with T-ALL and gene fusions was found (p > 0.05). The seven SNPs were associated with risk of pre-B ALL in younger children; however, rs2393732, rs2393782, rs2893881, and rs4948488 were not associated with susceptibility in older children and adolescents. The CAG haplotype (rs10821936, rs10994982, rs7089424) was strongly associated with ALL risk in our population (p < 0.00001). The frequency of all risk alleles in our ALL, pre-B, and hyperdiploid-ALL patients was higher than that in Hispanic children reported. This is the first report showing the association between rs2393732, rs2393782, and rs4948488 with pre-B hyperdiploid-ALL children. The G allele at rs2893881 confers major risk for pre-B hyperdiploid-ALL in Mexican (OR, 2.29) than in Hispanic children (OR, 1.71). The genetic background of our population could influence the susceptibility to ALL and explain its high incidence in Mexico.

Analysis of 153 115 patients with hematological malignancies refines the spectrum of familial risk

Estimating familial cancer risks is clinically important in being able to discriminate between individuals in the population at differing risk for malignancy. To gain insight into the familial risk for the different hematological malignancies and their possible inter-relationship, we analyzed data on more than 16 million individuals from the Swedish Family-Cancer Database. After identifying 153 115 patients diagnosed with a primary hematological malignancy, we quantified familial relative risks (FRRs) by calculating standardized incident ratios (SIRs) in 391 131 of their first-degree relatives. The majority of hematological malignancies showed increased FRRs for the same tumor type, with the highest FRRs being observed for mixed cellularity Hodgkin lymphoma (SIR, 16.7), lymphoplasmacytic lymphoma (SIR, 15.8), and mantle cell lymphoma (SIR, 13.3). There was evidence for pleiotropic relationships; notably, chronic lymphocytic leukemia was associated with an elevated familial risk for other B-cell tumors and myeloproliferative neoplasms. Collectively, these data provide evidence for shared etiological factors for many hematological malignancies and provide information for identifying individuals at increased risk, as well as informing future gene discovery initiatives.

Microarray testing as an efficient tool to redefine hyperdiploid paediatric B-cell precursor acute lymphoblastic leukaemia patients

The aim of our study was to characterize genetic alterations in a cohort of paediatric patients with B-cell progenitors (BCP-ALL) and a hyperdiploid karyotype. In our study, we analysed 55 childhood hyperdiploid BCP-ALL patients using single nucleotide polymorphism (SNP) microarray testing. The group consisted mostly of patients with the modal number of chromosomes between 54 and 58 (34 cases). Within this group, Trisomy 4 and Trisomy 10 (30 cases) were the most frequent cases. Additionally, a total of 93 structural abnormalities mainly affecting chromosomes 1, 6, 9, 12, and 17 as well as 68 copy number alterations (CNAs) were identified. The microarray testing revealed a loss of ETV6, IKZF1, CDKN2A/CDKN2B, PAX5, and RB1. Moreover, chromosomal abnormalities resulting in the loss of heterozygosity (LOH) were also observed. Currently, patients with hyperdiploidy constitute a genetically heterogeneous group, and therefore, it is insufficient to rely only on banding cytogenetic analysis for the identification of hyperdiploid karyotype. Microarray testing has been proven an effective and satisfactory tool for the analysis of molecular karyotypes and to redefine the prognostic criteria in hyperdiploid patients.

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.

Ionising radiation and childhood leukaemia revisited

Increased incidences of childhood acute leukaemia were noted among survivors of the atomic bombings of Hiroshima and Nagasaki. In Western societies, Childhood Acute Lymphoblastic Leukaemia has a distinct epidemiology peaking at 3 years old. Exposure to ionising radiation is an established hazard but it is difficult to gauge the precise risk of less than 100 mSv. Since 1983 significant leukaemia incidences have been reported among families residing near nuclear installations. The target cells (naïve neonatal lymphocytes) get exposed to multiple xenobiotic challenges and undergo extraordinary proliferation and physiological somatic genetic change. Population movements and ionising radiation are considered taking account of updated understanding of radiation biology, cancer cytogenetics and immunological diversity. Double Strand Breaks in DNA arise through metabolic generation of Reactive Oxygen Species, and nearly always are repaired; but mis-repairs can be oncogenic. Recombinant Activating Gene enzymes in rapidly dividing perinatal pre-B lymphocytes being primed for antibody diversity are targeted to Signal Sequences in the Immunoglobulin genes. off target pseudo-sequences may allow RAG enzymes to create autosomal DSBs which, when mis-repaired, become translocated oncogenes. Immunogens acting by chance at crucial stages may facilitate this. In such circumstances, oncogenic DSBs from ionising radiation are less likely to be significant.

Comparative features and outcomes between paediatric T-cell and B-cell acute lymphoblastic leukaemia

Contemporary paediatric clinical trials have improved 5-year event-free survival above 85% and 5-year overall survival above 90% in B-cell acute lymphoblastic leukaemia (ALL) in many study groups, whilst outcomes for T-cell ALL are still lagging behind by 5-10% in most studies. Several factors have contributed to this discrepant outcome. First, patients with T-cell ALL are generally older than those with B-cell ALL and, therefore, have poorer tolerance to chemotherapy, especially dexamethasone and asparaginase, and have increased risk of extramedullary relapse. Second, a higher proportion of patients with B-cell ALL have favourable genetic subtypes (eg, ETV6-RUNX1 and high hyperdiploidy), which confer a superior outcome compared with favourable subtypes of T-cell ALL. Third, T-cell ALL blasts are generally more resistant to conventional chemotherapeutic drugs than are B-cell ALL blasts. Finally, patients with B-cell ALL are more amendable to available targeted therapies, such as Philadelphia chromosome-positive and some Philadelphia chromosome-like ALL cases to ABL-class tyrosine kinase inhibitors, and CD19-positive and CD22-postive B-cell ALL cases to a variety of immunotherapies. Several novel treatments under investigation might narrow the gap in survival between T-cell ALL and B-cell ALL, although novel treatment options for T-cell ALL are limited.

Primary Immunodeficiency and Cancer Predisposition Revisited: Embedding Two Closely Related Concepts Into an Integrative Conceptual Framework

Common understanding suggests that the normal function of a "healthy" immune system safe-guards and protects against the development of malignancies, whereas a genetically impaired one might increase the likelihood of their manifestation. This view is primarily based on and apparently supported by an increased incidence of such diseases in patients with specific forms of immunodeficiencies that are caused by high penetrant gene defects. As I will review and discuss herein, such constellations merely represent the tip of an iceberg. The overall situation is by far more varied and complex, especially if one takes into account the growing difficulties to define what actually constitutes an immunodeficiency and what defines a cancer predisposition. The enormous advances in genome sequencing, in bioinformatic analyses and in the functional in vitro and in vivo assessment of novel findings together with the availability of large databases provide us with a wealth of information that steadily increases the number of sequence variants that concur with clinically more or less recognizable immunological problems and their consequences. Since many of the newly identified hard-core defects are exceedingly rare, their tumor predisposing effect is difficult to ascertain. The analyses of large data sets, on the other hand, continuously supply us with low penetrant variants that, at least in statistical terms, are clearly tumor predisposing, although their specific relevance for the respective carriers still needs to be carefully assessed on an individual basis. Finally, defects and variants that affect the same gene families and pathways in both a constitutional and somatic setting underscore the fact that immunodeficiencies and cancer predisposition can be viewed as two closely related errors of development. Depending on the particular genetic and/or environmental context as well as the respective stage of development, the same changes can have either a neutral, predisposing and, in some instances, even a protective effect. To understand the interaction between the immune system, be it "normal" or "deficient" and tumor predisposition and development on a systemic level, one therefore needs to focus on the structure and dynamic functional organization of the entire immune system rather than on its isolated individual components alone.

Insights of Acute Lymphoblastic Leukemia with Development of Genomic Investigation

Treatment outcomes for acute lymphoblastic leukemia (ALL), especially pediatric ALL, have greatly improved due to the risk-adapted therapy. Combination of drug development, clinical practice, as well as basic genetic researches has brought the survival rate of ALL from less than 10% to more than 90% today, not only increasing the treatment efficacy but also limiting adverse drug reactions (ADRs). In this review, we summarized the landscape identification of ALL genetic alterations, which provided the opportunity to increase the survival rate and especially minimize the relapse risk of ALL, and highlighted the importance of the development of new technologies of genomic investigation for translational medicine.

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

Genome-wide association studies of childhood acute lymphoblastic leukemia (ALL) have identified regions of association at PIP4K2A and upstream of BMI1 at chromosome 10p12.31-12.2. The contribution of both loci to ALL risk and underlying functional variants remain to be elucidated. We carried out single nucleotide polymorphism (SNP) imputation across chromosome 10p12.31-12.2 in Latino and non-Latino white ALL cases and controls from two independent California childhood leukemia studies, and additional Genetic Epidemiology Research on Aging study controls. Ethnicity-stratified association analyses were performed using logistic regression, with meta-analysis including 3,133 cases (1,949 Latino, 1,184 non-Latino white) and 12,135 controls (8,584 Latino, 3,551 non-Latino white). SNP associations were identified at both BMI1 and PIP4K2A. After adjusting for the lead PIP4K2A SNP, genome-wide significant associations remained at BMI1, and vice-versa (pmeta < 10-10 ), supporting independent effects. Lead SNPs differed by ethnicity at both peaks. We sought functional variants in tight linkage disequilibrium with both the lead Latino SNP among Admixed Americans and lead non-Latino white SNP among Europeans. This pinpointed rs11591377 (pmeta = 2.1 x 10-10 ) upstream of BMI1, residing within a hematopoietic stem cell enhancer of BMI1, and which showed significant preferential binding of the risk allele to MYBL2 (p = 1.73 x 10-5 ) and p300 (p = 1.55 x 10-3 ) transcription factors using binomial tests on ChIP-Seq data from a SNP heterozygote. At PIP4K2A, we identified rs4748812 (pmeta = 1.3 x 10-15 ), which alters a RUNX1 binding motif and demonstrated chromosomal looping to the PIP4K2A promoter. Fine-mapping chromosome 10p12 in a multi-ethnic ALL GWAS confirmed independent associations and identified putative functional variants upstream of BMI1 and at PIP4K2A.

A causal mechanism for childhood acute lymphoblastic leukaemia

In this Review, I present evidence supporting a multifactorial causation of childhood acute lymphoblastic leukaemia (ALL), a major subtype of paediatric cancer. ALL evolves in two discrete steps. First, in utero initiation by fusion gene formation or hyperdiploidy generates a covert, pre-leukaemic clone. Second, in a small fraction of these cases, the postnatal acquisition of secondary genetic changes (primarily V(D)J recombination-activating protein (RAG) and activation-induced cytidine deaminase (AID)-driven copy number alterations in the case of ETS translocation variant 6 (ETV6)-runt-related transcription factor 1 (RUNX1)+ ALL) drives conversion to overt leukaemia. Epidemiological and modelling studies endorse a dual role for common infections. Microbial exposures earlier in life are protective but, in their absence, later infections trigger the critical secondary mutations. Risk is further modified by inherited genetics, chance and, probably, diet. Childhood ALL can be viewed as a paradoxical consequence of progress in modern societies, where behavioural changes have restrained early microbial exposure. This engenders an evolutionary mismatch between historical adaptations of the immune system and contemporary lifestyles. Childhood ALL may be a preventable cancer.

Genetics and epigenetics of pediatric leukemia in the era of precision medicine

Pediatric leukemia represents a heterogeneous group of diseases characterized by germline and somatic mutations that manifest within the context of disturbances in the epigenetic machinery and genetic regulation. Advances in genomic medicine have allowed finer resolution of genetic and epigenetic strategies that can be effectively used to risk-stratify patients and identify novel targets for therapy. This review discusses the genetic and epigenetic mechanisms of leukemogenesis, particularly as it relates to acute lymphocytic leukemias, the mechanisms of epigenetic control of leukemogenesis, namely DNA methylation, histone modifications, microRNAs, and LINE-1 retroelements, and highlights opportunities for precision medicine therapeutics in further guiding disease management. Future efforts to broaden the integration of advances in genomic and epigenomic science into the practice of pediatric oncology will not only identify novel therapeutic strategies to improve clinical outcomes but also improve the quality of life for this unique patient population. Recent findings in precision therapeutics of acute lymphocytic leukemias over the past three years, along with some provocative areas of epigenetics research, are reviewed here.

Genetic predisposition to B-cell acute lymphoblastic leukemia at 14q11.2 is mediated by a CEBPE promoter polymorphism

Acute lymphoblastic leukaemia (ALL) is the most common paediatric malignancy. Genome-wide association studies have shown variation at 14q11.2 influences ALL risk. We sought to decipher causal variant(s) at 14q11.2 and the mechanism of tumorigenesis. We show rs2239630 G>A resides in the promoter of the CCAT enhancer-binding protein epsilon (CEBPE) gene. The rs2239630-A risk allele is associated with increased promotor activity and CEBPE expression. Depletion of CEBPE in ALL cells reduces cell growth, correspondingly CEBPE binds to the promoters of electron transport and energy generation genes. RNA-seq in CEBPE depleted cells demonstrates CEBPE regulates the expression of genes involved in B-cell development (IL7R), apoptosis (BCL2), and methotrexate resistance (RASS4L). CEBPE regulated genes significantly overlapped in CEBPE depleted cells, ALL blasts and IGH-CEBPE translocated ALL. This suggests CEBPE regulates a similar set of genes in each, consistent with a common biological mechanism of leukemogenesis for rs2239630 associated and CEBPE translocated ALL. Finally, we map IGH-CEBPE translocation breakpoints in two cases, implicating RAG recombinase activity in their formation.

Genome-wide association study identifies susceptibility loci for B-cell childhood acute lymphoblastic leukemia

Genome-wide association studies (GWAS) have advanced our understanding of susceptibility to B-cell precursor acute lymphoblastic leukemia (BCP-ALL); however, much of the heritable risk remains unidentified. Here, we perform a GWAS and conduct a meta-analysis with two existing GWAS, totaling 2442 cases and 14,609 controls. We identify risk loci for BCP-ALL at 8q24.21 (rs28665337, P = 3.86 × 10-9, odds ratio (OR) = 1.34) and for ETV6-RUNX1 fusion-positive BCP-ALL at 2q22.3 (rs17481869, P = 3.20 × 10-8, OR = 2.14). Our findings provide further insights into genetic susceptibility to ALL and its biology.

GWAS in childhood acute lymphoblastic leukemia reveals novel genetic associations at chromosomes 17q12 and 8q24.21

Childhood acute lymphoblastic leukemia (ALL) (age 0-14 years) is 20% more common in Latino Americans than non-Latino whites. We conduct a genome-wide association study in a large sample of 3263 Californian children with ALL (including 1949 of Latino heritage) and 3506 controls matched on month and year of birth, sex, and ethnicity, and an additional 12,471 controls from the Kaiser Resource for Genetic Epidemiology Research on Aging Cohort. Replication of the strongest genetic associations is performed in two independent datasets from the Children's Oncology Group and the California Childhood Leukemia Study. Here we identify new risk loci on 17q12 near IKZF3/ZPBP2/GSDMB/ORMDL3, a locus encompassing a transcription factor important for lymphocyte development (IKZF3), and at an 8q24 region known for structural contacts with the MYC oncogene. These new risk loci may impact gene expression via local (four 17q12 genes) or long-range (8q24) interactions, affecting function of well-characterized hematopoietic and growth-regulation pathways.

Regional evaluation of childhood acute lymphoblastic leukemia genetic susceptibility loci among Japanese

Genome-wide association studies (GWAS) performed mostly in populations of European and Hispanic ancestry have confirmed an inherited genetic basis for childhood acute lymphoblastic leukemia (ALL), but these associations are less clear in other races/ethnicities. DNA samples from ALL patients (aged 0–19 years) previously enrolled onto a Tokyo Children’s Cancer Study Group trial were collected during 2013–2015, and underwent single nucleotide polymorphism (SNP) microarray genotyping resulting in 527 B-cell ALL for analysis. Cases and control data for 3,882 samples from the Nagahama Study Group and Aichi Cancer Center Study were combined, and association analyses across 10 previous GWAS-identified regions were performed after targeted SNP imputation. Linkage disequilibrium (LD) patterns in Japanese and other populations were evaluated using the varLD score based on 1000 Genomes data. Risk associations for ARID5B (rs10821936, OR = 1.84, P = 6 × 10⁻¹⁷) and PIP4K2A (rs7088318, OR = 0.76, P = 2 × 10⁻⁴) directly transferred to Japanese, and the IKZF1 association was detected by an alternate SNP (rs1451367, OR = 1.52, P = 2 × 10⁻⁶). Marked regional LD differences between Japanese and Europeans was observed for most of the remaining loci for which associations did not transfer, including CEBPE, CDKN2A, CDKN2B, and ELK3. This study represents a first step towards characterizing the role of genetic susceptibility in childhood ALL risk in Japanese.

ARID5B as a critical downstream target of the TAL1 complex that activates the oncogenic transcriptional program and promotes T-cell leukemogenesis

Leong et al. identified ARID5B as a collaborating oncogenic factor involved in the transcriptional program in T-ALL. ARID5B positively regulates the expression of TAL1 and its regulatory partners and also activates the expression of the oncogene MYC.

The oncogenic transcription factor TAL1/SCL induces an aberrant transcriptional program in T-cell acute lymphoblastic leukemia (T-ALL) cells. However, the critical factors that are directly activated by TAL1 and contribute to T-ALL pathogenesis are largely unknown. Here, we identified AT-rich interactive domain 5B (ARID5B) as a collaborating oncogenic factor involved in the transcriptional program in T-ALL. ARID5B expression is down-regulated at the double-negative 2–4 stages in normal thymocytes, while it is induced by the TAL1 complex in human T-ALL cells. The enhancer located 135 kb upstream of the ARID5B gene locus is activated under a superenhancer in T-ALL cells but not in normal T cells. Notably, ARID5B-bound regions are associated predominantly with active transcription. ARID5B and TAL1 frequently co-occupy target genes and coordinately control their expression. ARID5B positively regulates the expression of TAL1 and its regulatory partners. ARID5B also activates the expression of the oncogene MYC. Importantly, ARID5B is required for the survival and growth of T-ALL cells, and forced expression of ARID5B in immature thymocytes results in thymus retention, differentiation arrest, radioresistance, and tumor formation in zebrafish. Our results indicate that ARID5B reinforces the oncogenic transcriptional program by positively regulating the TAL1-induced regulatory circuit and MYC in T-ALL, thereby contributing to T-cell leukemogenesis.

Genome-wide association studies of cancer: current insights and future perspectives

Genome-wide association studies (GWAS) provide an agnostic approach for investigating the genetic basis of complex diseases. In oncology, GWAS of nearly all common malignancies have been performed, and over 450 genetic variants associated with increased risks have been identified. As well as revealing novel pathways important in carcinogenesis, these studies have shown that common genetic variation contributes substantially to the heritable risk of many common cancers. The clinical application of GWAS is starting to provide opportunities for drug discovery and repositioning as well as for cancer prevention. However, deciphering the functional and biological basis of associations is challenging and is in part a barrier to fully unlocking the potential of GWAS.

Epigenetics in pediatric acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most common malignancy in children. ALL arises from the malignant transformation of progenitor B- and T-cells in the bone marrow into leukemic cells, but the mechanisms underlying this transformation are not well understood. Recent technical advances and decreasing costs of methods for high-throughput DNA sequencing and SNP genotyping have stimulated systematic studies of the epigenetic changes in leukemic cells from pediatric ALL patients. The results emerging from these studies are increasing our understanding of the epigenetic component of leukemogenesis and have demonstrated the potential of DNA methylation as a biomarker for lineage and subtype classification, prognostication, and disease progression in ALL. In this review, we provide a concise examination of the epigenetic studies in ALL, with a focus on DNA methylation and mutations perturbing genes involved in chromatin modification, and discuss the future role of epigenetic analyses in research and clinical management of ALL.