Identification of four novel associations for B-cell acute lymphoblastic leukaemia risk

The evolving role of B cells in malignancies

B cells play diverse roles in different pathological circumstances, such as neoplastic diseases, autoimmune disorders, and neurological maladies. B cells, which are essential elements of the adaptive immune system, demonstrate exceptional functional variety, including the generation of antibodies, the presentation of antigens, and the secretion of cytokines. Within the field of oncology, B cells display a multifaceted nature in the tumor microenvironment, simultaneously manifesting both tumor-promoting and tumor-suppressing characteristics. Studies have found that the existence of tertiary lymphoid structures, which consist of B cells, is linked to better survival rates in different types of cancers. This article examines the involvement of B cells in different types of malignancies, emphasizing their importance in the development of the diseases and their potential as biomarkers. Additionally, the review also examines the crucial role of B cells in autoimmune illnesses and their potential as targets for therapy. The article also analyses the role of B cells in immunization and exploring their potential uses in cancer immunotherapy. This analysis highlights the intricate and occasionally contradictory roles of B cells, underlining the necessity for additional research to clarify their varied actions in various illness scenarios.

Re-envisioning genetic predisposition to childhood and adolescent cancers

Although cancer is rare in children and adolescents, it remains a leading cause of death within this age range, and genetic predisposition is the main known risk factor. Since the discovery of retinoblastoma-predisposing RB1 pathogenic germline variants in 1985, several additional high-penetrance cancer predisposition genes (CPGs) have been identified. Although few clinically recognizable genetic conditions display moderate cancer phenotypes, burden testing has revealed low-to-moderate penetrance CPGs. In addition to germline pathogenic variants in CPGs, postzygotic somatic mosaic CPG pathogenic variants acquired during embryonic development are increasingly recognized as factors that predispose children and adolescents to malignancies. Genome-wide association studies of various childhood and adolescent cancer types have identified some common low-risk cancer susceptibility alleles. Although the clinical utility of polygenic risk scores is currently limited in children and adolescents, polygenic risk scores developed for adults can predict subsequent cancer risks in childhood and adolescent cancer survivors. In this Review, I describe our current knowledge of genetic predisposition to childhood and adolescent cancers. Survival rates in children and adolescents with cancer and CPGs are often poor, necessitating better integration of genomic testing into clinical care to improve cancer prevention, surveillance and therapies.

Overview of the etiology of childhood cancer and future directions

PURPOSE OF REVIEW

We provide an overview of the etiology of childhood cancer, the state of the literature, and highlight some opportunities for future research, including technological advancements that could be applied to etiologic studies of childhood cancer to accelerate our understanding.

RECENT FINDINGS

Risk factors of childhood cancer were summarized based on demographics and perinatal factors, environmental risk factors, and genetic risk factors. Overall, demographics and perinatal factors are the most well studied in relation to childhood cancer. While environmental risk factors have been implicated, more work is needed to pinpoint specific exposures, identify window(s) of susceptibility, and understand mechanisms. With genome-wide association studies (GWAS), genetic risk factors of eight childhood cancers have emerged, and opportunities remain to conduct GWAS for other cancer types and determine whether risk variants are inherited or de novo. Technological advancements that can shed light into the susceptibility of childhood cancer include metabolomics, using primary teeth as an exposure matrix, and long-read sequencing.

SUMMARY

The development of childhood cancer remains largely not well understood. Collaboration to increase sample size to conduct analyses by histology and/or molecular subtype and application of novel technologies will accelerate our understanding of childhood cancer.

The impact of Indigenous American-like ancestry on risk of acute lymphoblastic leukemia in Hispanic/Latino children

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer, with Hispanic/Latino children having a higher incidence of ALL than other racial/ethnic groups. Genetic variants, particularly ones found enriched in Indigenous American (IA)-like ancestry and inherited by Hispanics/Latinos, may contribute to this disparity. In this study, we characterized the impact of IA-like ancestry on overall ALL risk and the frequency and effect size of known risk alleles in a large cohort of self-reported Hispanic/Latino individuals. We also performed genome-wide admixture mapping analysis to identify potentially novel ALL risk loci. We found that global IA ancestry was positively associated with ALL risk, but the association was not significant after adjusting for socio-economic indicators. In a series of local ancestry analyses, we uncovered that at known ALL risk loci, increasing copies of the IA-like haplotype were positively and significantly associated with ALL case-control status. Further, the IA-like haplotype had ~1.33 times the odds of harboring the risk allele compared to non-IA-like haplotypes. We found no evidence of interaction between genotype and ancestry (local or global) in relation to ALL risk. Admixture mapping identified association signals on chromosomes 2 (2q21.2), 7 (7p12.2), 10 (10q21.2), and 15 (15q22.31); however, only the variants at 7p12.2 and 10q21.2 replicated in additional cohorts. Taken together, our results suggest that increased risk of ALL in Hispanic/Latino children may be conferred by higher frequency of risk alleles within IA-like ancestry, which can be leveraged as targets of new precision health strategies and therapeutics.

A phenome-wide association study of polygenic scores for selected childhood cancer: Results from the UK Biobank

The aim of this study was to scan phenotypes in adulthood associated with polygenic risk scores (PRS) for childhood cancers with well-articulated genetic architectures-acute lymphoblastic leukemia (ALL), Ewing sarcoma, and neuroblastoma-to examine genetic pleiotropy. Furthermore, we aimed to determine which SNPs could drive associations. Per-SNP summary statistics were extracted for PRS calculation. Participants with white British ancestry were exclusively included for analyses. SNPs were queried from the UK Biobank genotype imputation data. Records from the cancer registry, death registry, and inpatient diagnoses were abstracted for phenome-wide scans. Firth logistic regression was used to estimate odds ratios (ORs) and 95% confidence intervals (CIs) alongside corresponding p values, adjusting for age at recruitment and sex. A total of 244,332 unrelated white British participants were included. We observed a significant association between ALL-PRS and ALL (OR: 1.20e+24, 95% CI: 9.08e+14-1.60e+33). In addition, we observed a significant association between high-risk neuroblastoma PRS and nonrheumatic aortic valve disorders (OR: 43.9, 95% CI: 7.42-260). There were no significant phenotype associations with Ewing sarcoma and neuroblastoma PRS. Regarding individual SNPs, rs17607816 increased the risk of ALL (OR: 6.40, 95% CI: 3.26-12.57). For high-risk neuroblastoma, rs80059929 elevated the risk of atrioventricular block (OR: 3.04, 95% CI: 1.85-4.99). Our findings suggest that individuals with genetic susceptibility to ALL may face a lifelong risk for developing ALL, along with a genetic pleiotropic association between high-risk neuroblastoma and circulatory diseases.

Multi-step gene set analysis identified HTR3 family genes involving childhood acute lymphoblastic leukemia susceptibility

In our previous conventional genome-wide association study (GWAS), WWOX was a susceptibility gene associated with acute lymphoblastic leukemia (ALL) development. Nowadays, advancements in genetic association analyses promote an in-depth exploration of ALL genomics. We conducted a two-step enrichment analysis at both gene and pathway levels based on ALL GWAS data including 269 cases and 1039 controls of Chinese descent. The following functional prediction and experiments were used to evaluate the genetic biology of candidate variants and genes. The serotonin-activated cation-selective channel complex gene-set was a potential biological pathway involved in ALL occurrence. Of which, individuals carrying the T allele of rs33940208 exhibited a prominent reduced risk of ALL [odds ratio (OR) = 0.71, 95% confidence interval (CI) = 0.53 to 0.96, P = 2.81 × 10-2], whereas those with the A allele of rs6779545 demonstrated an increased risk (OR = 1.23, 95% CI = 1.01 to 1.51, P = 4.11 × 10-2). Notably, the two variants displayed a better prediction capability when combined, that the risk of developing childhood ALL increased by 131% in subjects with 2-4 risk alleles compared to those with 0-1 risk alleles (Ptrend = 2.05 × 10-3). In addition, the T allele of rs33940208 could reduce HTR3A mRNA level, while the A allele of rs6779545 increased HTR3D mRNA expression. In this study, we identified HTR3A rs33940208 and HTR3D rs6779545 as potential susceptibility loci for ALL in Chinese children. Future validation and functional research will elucidate the underlying molecular processes, refining preventive strategies for this disease.

Expression quantitative trait loci influence DNA damage-induced apoptosis in cancer

BACKGROUND

Genomic instability and evading apoptosis are two fundamental hallmarks of cancer and closely linked to DNA damage response (DDR). By analyzing expression quantitative trait loci (eQTL) upon cell stimulation (called exposure eQTL (e2QTL)) it is possible to identify context specific gene regulatory variants and connect them to oncological diseases based on genome-wide association studies (GWAS).

RESULTS

We isolate CD8+ T cells from 461 healthy donors and stimulate them with high doses of 5 different carcinogens to identify regulatory mechanisms of DNA damage-induced apoptosis. Across all stimuli, we find 5,373 genes to be differentially expressed, with 85% to 99% of these genes being suppressed. While upregulated genes are specific to distinct stimuli, downregulated genes are shared across conditions but exhibit enrichment in biological processes depending on the DNA damage type. Analysis of eQTL reveals 654 regulated genes across conditions. Among them, 47 genes are significant e2QTL, representing a fraction of 4% to 5% per stimulus. To unveil disease relevant genetic variants, we compare eQTL and e2QTL with GWAS risk variants. We identify gene regulatory variants for KLF2, PIP4K2A, GPR160, RPS18, ARL17B and XBP1 that represent risk variants for oncological diseases.

CONCLUSION

Our study highlights the relevance of gene regulatory variants influencing DNA damage-induced apoptosis in cancer. The results provide new insights in cellular mechanisms and corresponding genes contributing to inter-individual effects in cancer development.

Concepts in B cell acute lymphoblastic leukemia pathogenesis

B cell acute lymphoblastic leukemia (B-ALL) arises from genetic alterations impacting B cell progenitors, ultimately leading to clinically overt disease. Extensive collaborative efforts in basic and clinical research have significantly improved patient prognoses. Nevertheless, a subset of patients demonstrate resistance to conventional chemotherapeutic approaches and emerging immunotherapeutic interventions. This review highlights the mechanistic underpinnings governing B-ALL transformation. Beginning with exploring normative B cell lymphopoiesis, we delineate the influence of recurrent germline and somatic genetic aberrations on the perturbation of B cell progenitor differentiation and protumorigenic signaling, thereby facilitating the neoplastic transformation underlying B-ALL progression. Additionally, we highlight recent advances in the multifaceted landscape of B-ALL, encompassing metabolic reprogramming, microbiome influences, inflammation, and the discernible impact of socioeconomic and racial disparities on B-ALL transformation and patient survival.

Acute lymphoblastic leukaemia

Acute lymphoblastic leukaemia (ALL) is a haematological malignancy characterized by the uncontrolled proliferation of immature lymphoid cells. Over past decades, significant progress has been made in understanding the biology of ALL, resulting in remarkable improvements in its diagnosis, treatment and monitoring. Since the advent of chemotherapy, ALL has been the platform to test for innovative approaches applicable to cancer in general. For example, the advent of omics medicine has led to a deeper understanding of the molecular and genetic features that underpin ALL. Innovations in genomic profiling techniques have identified specific genetic alterations and mutations that drive ALL, inspiring new therapies. Targeted agents, such as tyrosine kinase inhibitors and immunotherapies, have shown promising results in subgroups of patients while minimizing adverse effects. Furthermore, the development of chimeric antigen receptor T cell therapy represents a breakthrough in ALL treatment, resulting in remarkable responses and potential long-term remissions. Advances are not limited to treatment modalities alone. Measurable residual disease monitoring and ex vivo drug response profiling screening have provided earlier detection of disease relapse and identification of exceptional responders, enabling clinicians to adjust treatment strategies for individual patients. Decades of supportive and prophylactic care have improved the management of treatment-related complications, enhancing the quality of life for patients with ALL.

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.

Shared genetic architecture between autoimmune disorders and B-cell acute lymphoblastic leukemia: insights from large-scale genome-wide cross-trait analysis

BACKGROUND

To study the shared genetic structure between autoimmune diseases and B-cell acute lymphoblastic leukemia (B-ALL) and identify the shared risk loci and genes and genetic mechanisms involved.

METHODS

Based on large-scale genome-wide association study (GWAS) summary-level data sets, we observed genetic overlaps between autoimmune diseases and B-ALL, and cross-trait pleiotropic analysis was performed to detect shared pleiotropic loci and genes. A series of functional annotation and tissue-specific analysis were performed to determine the influence of pleiotropic genes. The heritability enrichment analysis was used to detect crucial immune cells and tissues. Finally, bidirectional Mendelian randomization (MR) methods were utilized to investigate the casual associations.

RESULTS

Our research highlighted shared genetic mechanisms between seven autoimmune disorders and B-ALL. A total of 73 pleiotropic loci were identified at the genome-wide significance level (P < 5 × 10-8), 16 of which had strong evidence of colocalization. We demonstrated that several loci have been previously reported (e.g., 17q21) and discovered some novel loci (e.g., 10p12, 5p13). Further gene-level identified 194 unique pleiotropic genes, for example IKZF1, GATA3, IKZF3, GSDMB, and ORMDL3. Pathway analysis determined the key role of cellular response to cytokine stimulus, B cell activation, and JAK-STAT signaling pathways. SNP-level and gene-level tissue enrichment suggested that crucial role pleiotropic mechanisms involved in the spleen, whole blood, and EBV-transformed lymphocytes. Also, hyprcoloc and stratified LD score regression analyses revealed that B cells at different developmental stages may be involved in mechanisms shared between two different diseases. Finally, two-sample MR analysis determined causal effects of asthma and rheumatoid arthritis on B-ALL.

CONCLUSIONS

Our research proved shared genetic architecture between autoimmune disorders and B-ALL and shed light on the potential mechanism that might involve in.

A noncoding regulatory variant in IKZF1 increases acute lymphoblastic leukemia risk in Hispanic/Latino children

Hispanic/Latino children have the highest risk of acute lymphoblastic leukemia (ALL) in the US compared to other racial/ethnic groups, yet the basis of this remains incompletely understood. Through genetic fine-mapping analyses, we identified a new independent childhood ALL risk signal near IKZF1 in self-reported Hispanic/Latino individuals, but not in non-Hispanic White individuals, with an effect size of ∼1.44 (95% confidence interval = 1.33-1.55) and a risk allele frequency of ∼18% in Hispanic/Latino populations and <0.5% in European populations. This risk allele was positively associated with Indigenous American ancestry, showed evidence of selection in human history, and was associated with reduced IKZF1 expression. We identified a putative causal variant in a downstream enhancer that is most active in pro-B cells and interacts with the IKZF1 promoter. This variant disrupts IKZF1 autoregulation at this enhancer and results in reduced enhancer activity in B cell progenitors. Our study reveals a genetic basis for the increased ALL risk in Hispanic/Latino children.

Epigenomic Characterization of Lymphoid Neoplasms

Lymphoid neoplasms represent a heterogeneous group of disease entities and subtypes with markedly different molecular and clinical features. Beyond genetic alterations, lymphoid tumors also show widespread epigenomic changes. These severely affect the levels and distribution of DNA methylation, histone modifications, chromatin accessibility, and three-dimensional genome interactions. DNA methylation stands out as a tracer of cell identity and memory, as B cell neoplasms show epigenetic imprints of their cellular origin and proliferative history, which can be quantified by an epigenetic mitotic clock. Chromatin-associated marks are informative to uncover altered regulatory regions and transcription factor networks contributing to the development of distinct lymphoid tumors. Tumor-intrinsic epigenetic and genetic aberrations cooperate and interact with microenvironmental cells to shape the transcriptome at different phases of lymphoma evolution, and intraclonal heterogeneity can now be characterized by single-cell profiling. Finally, epigenetics offers multiple clinical applications, including powerful diagnostic and prognostic biomarkers as well as therapeutic targets.

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.

RNA binding protein IGF2BP1 synergizes with ETV6-RUNX1 to drive oncogenic signaling in B-cell Acute Lymphoblastic Leukemia

BACKGROUND

Acute lymphoblastic leukemia (ALL) is the most common pediatric hematological malignancy, with ETV6::RUNX1 being the most prevalent translocation whose exact pathogenesis remains unclear. IGF2BP1 (Insulin-like Growth Factor 2 Binding Protein 1) is an oncofetal RNA binding protein seen to be specifically overexpressed in ETV6::RUNX1 positive B-ALL. In this study, we have studied the mechanistic role of IGF2BP1 in leukemogenesis and its synergism with the ETV6::RUNX1 fusion protein.

METHODS

Gene expression was analyzed from patient bone marrow RNA using Real Time RT-qPCR. Knockout cell lines were created using CRISPR-Cas9 based lentiviral vectors. RNA-Seq and RNA Immunoprecipitation sequencing (RIP-Seq) after IGF2BP1 pulldown were performed using the Illumina platform. Mouse experiments were done by retroviral overexpression of donor HSCs followed by lethal irradiation of recipients using a bone marrow transplant model.

RESULTS

We observed specific overexpression of IGF2BP1 in ETV6::RUNX1 positive patients in an Indian cohort of pediatric ALL (n=167) with a positive correlation with prednisolone resistance. IGF2BP1 expression was essential for tumor cell survival in multiple ETV6::RUNX1 positive B-ALL cell lines. Integrated analysis of transcriptome sequencing after IGF2BP1 knockout and RIP-Seq after IGF2BP1 pulldown in Reh cell line revealed that IGF2BP1 targets encompass multiple pro-oncogenic signalling pathways including TNFα/NFκB and PI3K-Akt pathways. These pathways were also dysregulated in primary ETV6::RUNX1 positive B-ALL patient samples from our center as well as in public B-ALL patient datasets. IGF2BP1 showed binding and stabilization of the ETV6::RUNX1 fusion transcript itself. This positive feedback loop led to constitutive dysregulation of several oncogenic pathways. Enforced co-expression of ETV6::RUNX1 and IGF2BP1 in mouse bone marrow resulted in marrow hypercellularity which was characterized by multi-lineage progenitor expansion and strong Ki67 positivity. This pre-leukemic phenotype confirmed their synergism in-vivo. Clonal expansion of cells overexpressing both ETV6::RUNX1 and IGF2BP1 was clearly observed. These mice also developed splenomegaly indicating extramedullary hematopoiesis.

CONCLUSION

Our data suggest a combined impact of the ETV6::RUNX1 fusion protein and RNA binding protein, IGF2BP1 in activating multiple oncogenic pathways in B-ALL which makes IGF2BP1 and these pathways as attractive therapeutic targets and biomarkers.

Genes of Predisposition to Childhood Beta-Cell Acute Lymphoblastic Leukemia in the Kazakh Population

BACKGROUND

Today, acute lymphoblastic leukemia is one of the most common malignant diseases of the hematopoietic system. The genetic predisposition to ALL is not fully explored in various ethnic populations.

OBJECTIVE

The study aimed to conduct a comparative analysis of the population frequencies of alleles and genotypes of polymorphic gene variants: immune regulation GATA3 (rs3824662); transcription and differentiation of B cells: ARID5B (rs7089424, rs10740055), IKZF1 (rs4132601); differentiation of hematopoietic cells: PIP4K2A (rs7088318); apoptosis: CEBPE (rs2239633), tumor suppressors: CDKN2A (rs3731249), TP53 (rs1042522); carcinogen metabolism: CBR3 (rs1056892), CYP1A1 (rs104894, rs4646903), according to genome-wide association studies analyses associated with the risk of developing pediatric beta-cell acute lymphoblastic leukemia (B-cell ALL), in an ethnically homogeneous population of Kazakhs with studied populations.

METHODS

The genomic database consists of 1800 conditionally healthy persons of Kazakh nationality, genotyped using OmniChip 2.5-8 Illumina chips at the deCODE genetics as part of the InterPregGen 7 project of the European Union (EU) framework program under Grant Agreement No. 282540.

RESULTS

High population frequencies of single nucleotide polymorphism (SNP) minor alleles identified for immune regulation genes - GATA3 rs3824662 - 42.5%; transcription and differentiation of B-cells genes - ARID5B rs7089424 - 33.1% and rs10740055 - 48.5%, which suggests their significant genetic contribution to the risk of development and prognosis of the effectiveness of B-cell ALL therapy in the Kazakh population. The significantly lower population frequency of the minor allele G rs1056892 CBR3 gene - 38.6% in the Kazakhs suggests its significant protective effect in reducing the risk of childhood B-cell ALL and the smaller number of cardiac complications after anthracycline therapy.

CONCLUSION

The obtained results will serve as a basis for developing effective methods for predicting the risk of development, early diagnosis, and effectiveness of treatment of B-cell ALL in children.

Genomic landscape of Down syndrome-associated acute lymphoblastic leukemia

Trisomy 21, the genetic cause of Down syndrome (DS), is the most common congenital chromosomal anomaly. It is associated with a 20-fold increased risk of acute lymphoblastic leukemia (ALL) during childhood and results in distinctive leukemia biology. To comprehensively define the genomic landscape of DS-ALL, we performed whole-genome sequencing and whole-transcriptome sequencing (RNA-Seq) on 295 cases. Our integrated genomic analyses identified 15 molecular subtypes of DS-ALL, with marked enrichment of CRLF2-r, IGH::IGF2BP1, and C/EBP altered (C/EBPalt) subtypes compared with 2257 non-DS-ALL cases. We observed abnormal activation of the CEBPD, CEBPA, and CEBPE genes in 10.5% of DS-ALL cases via a variety of genomic mechanisms, including chromosomal rearrangements and noncoding mutations leading to enhancer hijacking. A total of 42.3% of C/EBP-activated DS-ALL also have concomitant FLT3 point mutations or insertions/deletions, compared with 4.1% in other subtypes. CEBPD overexpression enhanced the differentiation of mouse hematopoietic progenitor cells into pro-B cells in vitro, particularly in a DS genetic background. Notably, recombination-activating gene-mediated somatic genomic abnormalities were common in DS-ALL, accounting for a median of 27.5% of structural alterations, compared with 7.7% in non-DS-ALL. Unsupervised hierarchical clustering analyses of CRLF2-rearranged DS-ALL identified substantial heterogeneity within this group, with the BCR::ABL1-like subset linked to an inferior event-free survival, even after adjusting for known clinical risk factors. These results provide important insights into the biology of DS-ALL and point to opportunities for targeted therapy and treatment individualization.

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.

Significance of CEBPE Gene Promoter Polymorphism (Rs2239630 G > A ) Assessment in Childhood B-cell Acute Lymphoblastic Leukemia

BACKGROUND

A significant association has been reported between CEBPE gene promoter polymorphisms (rs2239630 G > A ) and the incidence of B-cell acute lymphoblastic leukemia (B-ALL). However, no previous study on this issue has been included among the Egyptian cohort of pediatric patients with B-ALL. Therefore, this study was designed to address the associations between CEBPE polymorphisms and susceptibility to B-ALL, as well as its impact on the outcome of B-ALL Egyptian patients with B-ALL.

PATIENTS AND METHODS

In the current study, we evaluated the rs2239630 polymorphism in 225 pediatric patients and 228 controls to assess the association of different rs2239630 genotypes with childhood susceptibility to B-ALL and the impact on the outcome of the patients.

RESULTS

The frequency of the A allele was significantly higher in the cases of B-ALL compared with the control group ( P = 0.004). By analyzing different genotypes for the predictive value of disease development, the GA and AA genotypes have been identified to be the highest among multivariate factors with an odds ratio of 3.330 (95% CI: 1.105-10.035). Likewise, the A allele was significantly associated with the shortest overall survival.

CONCLUSIONS

CEBPE gene promoter polymorphism (rs2239630 G > A ) AA is frequently associated with B-ALL; and has the worst overall survival among the 3 genotypes, followed by the GA and GG genotypes ( P < 0.001).

Hacking hematopoiesis - emerging tools for examining variant effects

Hematopoiesis is a continuous process of blood and immune cell production. It is orchestrated by thousands of gene products that respond to extracellular signals by guiding cell fate decisions to meet the needs of the organism. Although much of our knowledge of this process comes from work in model systems, we have learned a great deal from studies on human genetic variation. Considerable insight has emerged from studies on presumed monogenic blood disorders, which continue to provide key insights into the mechanisms critical for hematopoiesis. Furthermore, the emergence of large-scale biobanks and cohorts has uncovered thousands of genomic loci associated with blood cell traits and diseases. Some of these blood cell trait-associated loci act as modifiers of what were once thought to be monogenic blood diseases. However, most of these loci await functional validation. Here, we discuss the validation bottleneck and emerging methods to more effectively connect variant to function. In particular, we highlight recent innovations in genome editing, which have paved the path forward for high-throughput functional assessment of loci. Finally, we discuss existing barriers to progress, including challenges in manipulating the genomes of primary hematopoietic cells.

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.

Chemotherapy-related hyperbilirubinemia in pediatric acute lymphoblastic leukemia: a genome-wide association study from the AIEOP-BFM ALL study group

BACKGROUND

Characterization of clinical phenotypes in context with tumor and host genomic information can aid in the development of more effective and less toxic risk-adapted and targeted treatment strategies. To analyze the impact of therapy-related hyperbilirubinemia on treatment outcome and to identify contributing genetic risk factors of this well-recognized adverse effect we evaluated serum bilirubin levels in 1547 pediatric patients with acute lymphoblastic leukemia (ALL) and conducted a genome-wide association study (GWAS).

PATIENTS AND METHODS

Patients were treated in multicenter trial AIEOP-BFM ALL 2000 for pediatric ALL. Bilirubin toxicity was graded 0 to 4 according to the Common Toxicity Criteria (CTC) of the National Cancer Institute. In the GWAS discovery cohort, including 650 of the 1547 individuals, genotype frequencies of 745,895 single nucleotide variants were compared between 435 patients with hyperbilirubinemia (CTC grades 1-4) during induction/consolidation treatment and 215 patients without it (grade 0). Replication analyses included 224 patients from the same trial.

RESULTS

Compared to patients with no (grade 0) or moderate hyperbilirubinemia (grades 1-2) during induction/consolidation, patients with grades 3-4 had a poorer 5-year event free survival (76.6 ± 3% versus 87.7 ± 1% for grades 1-2, P = 0.003; 85.2 ± 2% for grade 0, P < 0.001) and a higher cumulative incidence of relapse (15.6 ± 3% versus 9.0 ± 1% for grades 1-2, P = 0.08; 11.1 ± 1% for grade 0, P = 0.007). GWAS identified a strong association of the rs6744284 variant T allele in the UGT1A gene cluster with risk of hyperbilirubinemia (allelic odds ratio (OR) = 2.1, P = 7 × 10^- 8). TT-homozygotes had a 6.5-fold increased risk of hyperbilirubinemia (grades 1-4; 95% confidence interval (CI) = 2.9-14.6, P = 7 × 10^- 6) and a 16.4-fold higher risk of grade 3-4 hyperbilirubinemia (95% CI 6.1-43.8, P = 2 × 10^- 8). Replication analyses confirmed these associations with joint analysis yielding genome-wide significance (allelic OR = 2.1, P = 6 × 10^- 11; 95% CI 1.7-2.7). Moreover, rs6744284 genotypes were strongly linked to the Gilbert's syndrome-associated UGT1A1*28/*37 allele (r² = 0.70), providing functional support for study findings. Of clinical importance, the rs6744284 TT genotype counterbalanced the adverse prognostic impact of high hyperbilirubinemia on therapy outcome.

CONCLUSIONS

Chemotherapy-related hyperbilirubinemia is a prognostic factor for treatment outcome in pediatric ALL and genetic variation in UGT1A aids in predicting the clinical impact of hyperbilirubinemia.

TRIAL REGISTRATION

http://www.

CLINICALTRIALS

gov ; #NCT00430118.

Lymphocytes from B-acute lymphoblastic leukemia patients present differential regulation of the adenosinergic axis depending on risk stratification

PURPOSE

Although risk-stratified chemotherapy regimens improve B-cell acute lymphoblastic leukemia (B-ALL) clinical outcome, relapse occurs in a significant number of cases. The identification of new therapeutic targets as well as prognostic and diagnostic biomarkers can improve B-ALL patients' clinical outcomes. Purinergic signaling is an important pathway in cancer progression, however the expression of ectonucleotidases and their impact on immune cells in B-ALL lacks exploration. We aimed to analyze the expression of ectonucleotidases in B-ALL patients' lymphocyte subpopulations.

METHODS

Peripheral blood samples from 15 patients diagnosed with B-ALL were analyzed. Flow cytometry was used to analyze cellularity, expression level of CD38, CD39, and CD73, and frequency of [Formula: see text], and [Formula: see text] in lymphocyte subpopulations. Plasma was used for cytokines (by CBA kit) and adenine nucleosides/nucleotides detection (by HPLC).

RESULTS

Comparing B-ALL patients to health donors, we observed an increase of CD4 + and CD8 + T-cells. In addition, a decrease in CD38 expression in B and Treg subpopulations and an increase in CD39⁺ CD73⁺  frequency in Breg and CD8+ T-cells. Analyzing cytokines and adenine nucleosides/nucleotides, we found a decrease in TNF, IL-1β, and ADO concentrations, together with an increase in AMP in B-ALL patients' plasma.

CONCLUSION

As immunomodulators, the expression of ectonucleotidases might be associated with activation states, as well as the abundance of different cellular subsets. We observed a pro-tumor immunity expression profile in B-ALL patients at diagnosis, being associated with cell exhaustion and immune evasion in B-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.

Heterozygous BRCA1 and BRCA2 and Mismatch Repair Gene Pathogenic Variants in Children and Adolescents With Cancer

BACKGROUND

Genetic predisposition is has been identified as a cause of cancer, yet little is known about the role of adult cancer predisposition syndromes in childhood cancer. We examined the extent to which heterozygous pathogenic germline variants in BRCA1, BRCA2, PALB2, ATM, CHEK2, MSH2, MSH6, MLH1, and PMS2 contribute to cancer risk in children and adolescents.

METHODS

We conducted a meta-analysis of 11 studies that incorporated comprehensive germline testing for children and adolescents with cancer. ClinVar pathogenic or likely pathogenic variants (PVs) in genes of interest were compared with 2 control groups. Results were validated in a cohort of mainly European patients and controls. We employed the Proxy External Controls Association Test to account for different pipelines.

RESULTS

Among 3975 children and adolescents with cancer, statistically significant associations with cancer risk were observed for PVs in BRCA1 and 2 (26 PVs vs 63 PVs among 27 501 controls, odds ratio = 2.78, 95% confidence interval = 1.69 to 4.45; P < .001) and mismatch repair genes (19 PVs vs 14 PVs among 27 501 controls, odds ratio = 7.33, 95% confidence interval = 3.64 to 14.82; P <.001). Associations were seen in brain and other solid tumors but not in hematologic neoplasms. We confirmed similar findings in 1664 pediatric cancer patients primarily of European descent.

CONCLUSION

These data suggest that heterozygous PVs in BRCA1 and 2 and mismatch repair genes contribute with reduced penetrance to cancer risk in children and adolescents. No changes to predictive genetic testing and surveillance recommendations are required.

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.

Nested epistasis enhancer networks for robust genome regulation

Mammalian genomes possess multiple enhancers spanning an ultralong distance (>megabases) to modulate important genes, yet it is unclear how these enhancers coordinate to achieve this task. Here, we combine multiplexed CRISPRi screening with machine learning to define quantitative enhancer-enhancer interactions. We find that the ultralong distance enhancer network possesses a nested multi-layer architecture that confers functional robustness of gene expression. Experimental characterization reveals that enhancer epistasis is maintained by three-dimensional chromosomal interactions and BRD4 condensation. Machine learning prediction of synergistic enhancers provides an effective strategy to identify non-coding variant pairs associated with pathogenic genes in diseases beyond Genome-Wide Association Studies (GWAS) analysis. Our work unveils nested epistasis enhancer networks, which can better explain enhancer functions within cells and in diseases.

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.

Identification of Genomic Variants Associated with the Risk of Acute Lymphoblastic Leukemia in Native Americans from Brazilian Amazonia

A number of genomic variants related to native American ancestry may be associated with an increased risk of developing Acute Lymphoblastic Leukemia (ALL), which means that Latin American and hispanic populations from the New World may be relatively susceptible to this disease. However, there has not yet been any comprehensive investigation of the variants associated with susceptibility to ALL in traditional Amerindian populations from Brazilian Amazonia. We investigated the exomes of the 18 principal genes associated with susceptibility to ALL in samples of 64 Amerindians from this region, including cancer-free individuals and patients with ALL. We compared the findings with the data on populations representing five continents available in the 1000 Genomes database. The variation in the allele frequencies found between the different groups was evaluated using Fisher’s exact test. The analyses of the exomes of the Brazilian Amerindians identified 125 variants, seven of which were new. The comparison of the allele frequencies between the two Amerindian groups analyzed in the present study (ALL patients vs. cancer-free individuals) identified six variants (rs11515, rs2765997, rs1053454, rs8068981, rs3764342, and rs2304465) that may be associated with susceptibility to ALL. These findings contribute to the identification of genetic variants that represent a potential risk for ALL in Amazonian Amerindian populations and might favor precision oncology measures.

Identification of diagnosis and prognosis gene markers in B-ALL with ETV6-RUNX1 fusion by integrated bioinformatics analysis

B-acute lymphoblastic leukemia (B-ALL) is characterized by clonal expansion of immature B-lymphocytes in the bone marrow, blood, or other tissues. Chromosomal translocations have often been reported in B-ALL, which are important for its prognosis. B-ALL patients with ETV6-RUNX1 fusion have favorable outcomes, but the mechanisms remain to be clarified. In the present study, we crossed the selected WGCNA module genes and differential expression genes to obtain core genes, and random forest algorithm, a type of supervised learning analysis, was conducted to evaluate the importance of those core genes in distinguishing B-ALL samples with ETV6-RUNX2 fusion with extracting 5 genes as gene markers for ETV6-RUNX2 fusion. Moreover, we calculated the immune infiltration profiles and screened out the ETV6-RUNX2 association immune cells using the CIBERSORT algorithm. In conclusion, combined with various solid informatics methods, we depicted the underlying molecular and immune mechanism of ETV6-RUNX2 fusion and providing potential biological targets for diagnosing and treating B-ALL in the future.

Allergies, genetic polymorphisms of Th2 interleukins, and childhood acute lymphoblastic leukemia: The ESTELLE study

CONTEXT

A negative association between a history of allergy and childhood acute lymphoblastic leukemia (ALL) has been reported in previous studies, but remains debated. This work aimed to investigate this association accounting for genetic polymorphisms of the Th2 pathway cytokines (IL4, IL10, IL13, and IL4R).

METHODS

Analyses were based on the French case-control study ESTELLE (2010-2011). The complete sample included 629 ALL cases and 1421 population-based controls frequency-matched on age and gender. The child's medical history was collected through standardized maternal interview. Biological samples were collected, and genotyping data were available for 411 cases and 704 controls of European origin. Odds ratios (OR) were estimated using unconditional regression models adjusted for potential confounders.

RESULTS

In the complete sample, a significant inverse association was observed between ALL and reported history of allergic rhinitis or sinusitis (OR = 0.65 [0.42-0.98]; P = 0.04), but there was no obvious association with allergies overall. There was an interaction between genetic polymorphisms in IL4 and IL4R (P_{interaction =} 0.003), as well as a gene-environment interaction between IL4R-rs1801275 and a reported history of asthma (IOR = 0.23; P_int  = 0.008) and eczema (IOR = 0.47; P_int  = 0.06). We observed no interaction with the candidate polymorphisms in IL4 and IL13.

CONCLUSION

These results suggest that the association between allergic symptoms and childhood ALL could be modified by IL4R-rs1801275, and that this variant could also interact with a functional variant in IL4 gene. Although they warrant confirmation, these results could help understand the pathological mechanisms under the reported inverse association between allergy and childhood ALL.

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.

Antibiotics in early life and childhood pre-B-ALL. Reasons to analyze a possible new piece in the puzzle

Acute lymphoblastic leukemia (ALL) is the most common pediatric cancer with precursor B-cell ALL (pB-ALL) accounting for ~ 85% of the cases. Childhood pB-ALL development is influenced by genetic susceptibility and host immune responses. The role of the intestinal microbiome in leukemogenesis is gaining increasing attention since Vicente-Dueñas' seminal work demonstrated that the gut microbiome is distinct in mice genetically predisposed to ALL and that the alteration of this microbiome by antibiotics is able to trigger pB-ALL in Pax5 heterozygous mice in the absence of infectious stimuli. In this review we provide an overview on novel insights on the role of the microbiome in normal and preleukemic hematopoiesis, inflammation, the effect of dysbiosis on hematopoietic stem cells and the emerging importance of the innate immune responses in the conversion from preleukemic to leukemic state in childhood ALL. Since antibiotics, which represent one of the most widely used medical interventions, alter the gut microbial composition and can cause a state of dysbiosis, this raises exciting epidemiological questions regarding the implications for antibiotic use in early life, especially in infants with a a preleukemic "first hit". Sheading light through a rigorous study on this piece of the puzzle may have broad implications for clinical practice.

Recent Advances in Pediatric Cancer Research

Over the past few years, the field of pediatric cancer has experienced a shift in momentum, and this has led to new and exciting findings that have relevance beyond pediatric malignancies. Here we present the current status of key aspects of pediatric cancer research. We have focused on genetic and epigenetic drivers of disease, cellular origins of different pediatric cancers, disease models, the tumor microenvironment, and cellular immunotherapies.

Impact of fetal expression quantitative trait loci on transcriptome-wide association study of childhood leukemia

Transcriptome-wide association studies increase the yield of loci associated with disease phenotypes by focusing on expression quantitative trait loci (eQTL). The major source of eQTL data for is the Gene and Tissue Expression (GTEx) project, which is comprised entirely of adults, mainly those > 50 years of age at death. Since gene expression levels differ by developmental stage, it is not clear whether eQTLs derived from adult data sources are best suited for use in young-onset diseases such as pediatric cancers. To fill in this knowledge gap, we performed a large-scale eQTL mapping analysis in the GenCord study with newborn samples and compared it with GTEx. Under matched conditions, we found around 80% of the eQTLs in one study can be replicated in the other. However, among all eQTLs identified in GenCord (GTEx), 584 (1045) showed statistically significant differences in effect sizes in GTEx (GenCord). We further investigated how using fetal eQTL data can facilitate the genetic association study of acute lymphoblastic leukemia. GenCord and GTEx identified the same genetic loci with statistical significance; however, the overall association pattern was only weakly correlated. Our paper demonstrates age-differential eQTLs and shows their potential influence on childhood leukemia research.

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.

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.

Infectious triggers and novel therapeutic opportunities in childhood B cell leukaemia

B cell acute lymphoblastic leukaemia (B-ALL) is the most common form of childhood cancer. Although treatment has advanced remarkably in the past 50 years, it still fails in ~20% of patients. Recent studies revealed that more than 5% of healthy newborns carry preleukaemic clones that originate in utero, but only a small percentage of these carriers will progress to overt B-ALL. The drivers of progression are unclear, but B-ALL incidence seems to be increasing in parallel with the adoption of modern lifestyles. Emerging evidence shows that a major driver for the conversion from the preleukaemic state to the B-ALL state is exposure to immune stressors, such as infection. Here, we discuss our current understanding of the environmental triggers and genetic predispositions that may lead to B-ALL, highlighting lessons from epidemiology, the clinic and animal models, and identifying priority areas for future research.

Family history of any cancer for childhood leukemia patients in Sweden

Acute lymphoblastic leukemia (ALL) is the most common childhood leukemia, while the other types, acute myeloid leukemia (AML), chronic lymphocytic leukemia (CLL), and chronic myeloid leukemia (CML) are much rarer. While data on familial risks for childhood ALL have been emerging, such data for the other childhood leukemias are hardly available. We aim to fill in the gap of knowledge by assessing familial clustering of each childhood leukemia with childhood and adult leukemia and with any cancer. We identified 4461 childhood leukemias from the Swedish Cancer Registry and obtained their family members from the Multigeneration Register. Standardized incidence ratios (SIRs) were 3.34 for singleton siblings both diagnosed with ALL before age 20 years and 1.64 for those who had a family member diagnosed with ALL in adult age. Other childhood leukemias showed no familial risk, but childhood ALL risk was increased to 1.40 when adult family members were diagnosed with CLL. Childhood ALL was associated with endometrial cancer, and female ALL patients showed increased risk when family members were diagnosed with testicular cancer, melanoma, and skin squamous cell carcinoma. Childhood CLL was associated with rectal cancer, and childhood AML was associated with pancreatic and bladder cancers. As most of these associations are reported for the first time, there is a need to replicate the findings from independent sources.

Super-enhancer Acquisition Drives FOXC2 Expression in Middle Ear Cholesteatoma

Distinct histone modifications regulate gene expression in certain diseases, but little is known about histone epigenetics in middle ear cholesteatoma. It is known that histone acetylation destabilizes the nucleosome and chromatin structure and induces gene activation. The association of histone acetylation with chronic inflammatory diseases has been indicated in recent studies. In this study, we examined the localization of variously modified histone H3 acetylation at lysine 9, 14, 18, 23, and 27 in paraffin-embedded sections of human middle ear cholesteatoma (cholesteatoma) tissues and the temporal bones of an animal model of cholesteatoma immunohistochemically. As a result, we found that there was a significant increase of the expression levels of H3K27ac both in human cholesteatoma tissues and the animal model. In genetics, super-enhancers are clusters of enhancers that drive the transcription of genes involved in cell identity. Super-enhancers were originally defined using the H3K27ac signal, and then we used H3K27ac chromatin immunoprecipitation followed by sequencing to map the active cis-regulatory landscape in human cholesteatoma. Based on the results, we identified increased H3K27ac signals as super-enhancers of the FOXC2 loci, as well as increased protein of FOXC2 in cholesteatoma. Recent studies have indicated that menin-MLL inhibitor could suppress tumor growth through the control of histone H3 modification. In this study, we demonstrated that the expression of FOXC2 was inhibited by menin-MLL inhibitor in vivo. These findings indicate that FOXC2 expression under histone modifications promoted the pathogenesis of cholesteatoma and suggest that it may be a therapeutic target of cholesteatoma.

Racial and ethnic disparities in pediatric cancer incidence among children and young adults in the United States by single year of age

BACKGROUND

Incidence rates of pediatric cancers in the United States are typically reported in 5-year age groups, obscuring variation by single year of age. Additionally, racial and ethnic variation in incidence is typically presented in broad categories rather than by narrow age ranges.

METHODS

The Surveillance, Epidemiology, and End Results (SEER) 18 data (2000-2017) were examined to calculate frequencies and age-adjusted incidence rates among individuals aged birth to 39 years. Incidence rate ratios (IRRs) and 95% confidence intervals (95% CIs) were estimated as the measure of association for rate comparisons by race and Hispanic origin overall and by single year of age.

RESULTS

Several histologic types showed substantial variation in race/ethnicity-specific and overall rates by single year of age. Overall, Black children and young adults experienced substantially decreased incidence of acute lymphoid leukemia (IRR, 0.52; 95% CI, 0.49-0.55) compared to Whites, and this decreased incidence was strongest at ages 1 through 7 years and 16 through 20 years. Hispanic individuals experienced decreased overall incidence of Hodgkin lymphoma (IRR, 0.50; 95% CI, 0.48-0.52) and astrocytoma (IRR, 0.54; 95% CI, 0.52-0.56) and increased risk of acute lymphoblastic leukemia (IRR, 1.46; 95% CI, 1.42-1.51) compared to non-Hispanic Whites, and the increased risk was strongest at ages 10 through 23 years. Substantial decreased risk across many tumor types was also observed for Asian/Pacific Islanders and American Indian/Alaska Natives.

CONCLUSIONS

Examination of incidence rates for pediatric cancers by narrow age groups may provide insights regarding etiological differences in subgroups. Additionally, variation in age-specific incidence rates by race and ethnicity may enable hypothesis generation on drivers of disparities observed.

Genome-Wide Association Analyses Identify Variants in IRF4 Associated With Acute Myeloid Leukemia and Myelodysplastic Syndrome Susceptibility

The role of common genetic variation in susceptibility to acute myeloid leukemia (AML), and myelodysplastic syndrome (MDS), a group of rare clonal hematologic disorders characterized by dysplastic hematopoiesis and high mortality, remains unclear. We performed AML and MDS genome-wide association studies (GWAS) in the DISCOVeRY-BMT cohorts (2,309 cases and 2,814 controls). Association analysis based on subsets (ASSET) was used to conduct a summary statistics SNP-based analysis of MDS and AML subtypes. For each AML and MDS case and control we used PrediXcan to estimate the component of gene expression determined by their genetic profile and correlate this imputed gene expression level with risk of developing disease in a transcriptome-wide association study (TWAS). ASSET identified an increased risk for de novo AML and MDS (OR = 1.38, 95% CI, 1.26-1.51, Pmeta = 2.8 × 10-12) in patients carrying the T allele at s12203592 in Interferon Regulatory Factor 4 (IRF4), a transcription factor which regulates myeloid and lymphoid hematopoietic differentiation. Our TWAS analyses showed increased IRF4 gene expression is associated with increased risk of de novo AML and MDS (OR = 3.90, 95% CI, 2.36-6.44, Pmeta = 1.0 × 10-7). The identification of IRF4 by both GWAS and TWAS contributes valuable insight on the role of genetic variation in AML and MDS susceptibility.

Family history of early onset acute lymphoblastic leukemia is suggesting genetic associations

Childhood acute lymphoblastic leukemia (ALL) has an origin in the fetal period which may distinguish it from ALL diagnosed later in life. We wanted to test whether familial risks differ in ALL diagnosed in the very early childhood from ALL diagnosed later. The Swedish nation-wide family-cancer data were used until year 2016 to calculate standardized incidence ratios (SIRs) for familial risks in ALL in three diagnostic age-groups: 0–4, 5–34 and 35 + years. Among 1335 ALL patients diagnosed before age 5, familial risks were increased for esophageal (4.78), breast (1.42), prostate (1.40) and connective tissue (2.97) cancers and leukemia (2.51, ALL 7.81). In age-group 5–34 years, rectal (1.73) and endometrial (2.40) cancer, myeloma (2.25) and leukemia (2.00, ALL 4.60) reached statistical significance. In the oldest age-group, the only association was with Hodgkin lymphoma (3.42). Diagnostic ages of family members of ALL patients were significantly lower compared to these cancers in the population for breast, prostate and rectal cancers. The patterns of increased familial cancers suggest that BRCA2 mutations could contribute to associations of ALL with breast and prostate cancers, and mismatch gene PMS2 mutations with rectal and endometrial cancers. Future DNA sequencing data will be a test for these familial predictions.

Age-related differences of genetic susceptibility to patients with acute lymphoblastic leukemia

Inherited predispositions to acute lymphoblastic leukemia have been well investigated in pediatric patients, but studies on adults, particularly Chinese patients, are limited. In this study, we conducted a genome-wide association study in 466 all-age Chinese patients with Acute lymphoblastic leukemia (ALL) and 1,466 non-ALL controls to estimate the impact of age on ALL susceptibility in the Chinese population. Among the 17 reported loci, 8 have been validated in pediatric and 1 in adult patients. The strongest association signal was identified at ARID5B locus and gradually decreased with age, while the signal at GATA3 exhibited the opposite trend and significantly impact on adult patients. With genome-wide approaches, germline variants at 2q14.3 rank as the top inherited predisposition to adult patients (e.g., rs73956024, P = 4.3 × 10^-5) and separate the genetic risk of pediatric vs. adult patients (P = 3.6 × 10^-6), whereas variants at 15q25.3 (e.g., rs11638062) have a similar impact on patients in different age groups (overall P = 2.9 × 10^-7). Our analysis highlights the impact of age on genetic susceptibility to ALL in Chinese patients.

RNA methylation in hematological malignancies and its interactions with other epigenetic modifications

Hematological malignancies are a class of malignant neoplasms attributed to abnormal differentiation of hematopoietic stem cells (HSCs). The systemic involvement, poor prognosis, chemotherapy resistance, and recurrence common in hematological malignancies urge researchers to look for novel treatment targets and mechanisms. In recent years, epigenetic abnormalities have been shown to play a vital role in tumorigenesis and progression in hematological malignancies. In addition to DNA methylation and histone modifications, which are most studied, RNA methylation has become increasingly significant. In this review, we elaborate recent advances in the understanding of RNA modification in the pathogenesis, diagnosis and molecular targeted therapies of hematological malignancies and discuss its intricate interactions with other epigenetic modifications, including DNA methylation, histone modifications and noncoding RNAs.

An immune window of opportunity to prevent childhood B cell leukemia

The prevalence of childhood B cell acute lymphoblastic leukemia (B-ALL) is increasing, particularly in developed countries. There is no clear explanation for this increment, but recent data suggest that, besides genetic predisposition, stress in the immune system (e.g., an infection) might have an important role in B-ALL leukemogenesis. Here, we speculate on how this knowledge might impact B-ALL prevention strategies.

An intact gut microbiome protects genetically predisposed mice against leukemia

The majority of childhood leukemias are precursor B-cell acute lymphoblastic leukemias (pB-ALLs) caused by a combination of prenatal genetic predispositions and oncogenic events occurring after birth. Although genetic predispositions are frequent in children (>1% to 5%), fewer than 1% of genetically predisposed carriers will develop pB-ALL. Although infectious stimuli are believed to play a major role in leukemogenesis, the critical determinants are not well defined. Here, by using murine models of pB-ALL, we show that microbiome disturbances incurred by antibiotic treatment early in life were sufficient to induce leukemia in genetically predisposed mice, even in the absence of infectious stimuli and independent of T cells. By using V4 and full-length 16S ribosomal RNA sequencing of a series of fecal samples, we found that genetic predisposition to pB-ALL (Pax5 heterozygosity or ETV6-RUNX1 fusion) shaped a distinct gut microbiome. Machine learning accurately (96.8%) predicted genetic predisposition using 40 of 3983 amplicon sequence variants as proxies for bacterial species. Transplantation of either wild-type (WT) or Pax5+/- hematopoietic bone marrow cells into WT recipient mice revealed that the microbiome is shaped and determined in a donor genotype-specific manner. Gas chromatography-mass spectrometry (GC-MS) analyses of sera from WT and Pax5+/- mice demonstrated the presence of a genotype-specific distinct metabolomic profile. Taken together, our data indicate that it is a lack of commensal microbiota rather than the presence of specific bacteria that promotes leukemia in genetically predisposed mice. Future large-scale longitudinal studies are required to determine whether targeted microbiome modification in children predisposed to pB-ALL could become a successful prevention strategy.