Early-life tobacco exposure is causally implicated in aberrant RAG-mediated recombination in childhood acute lymphoblastic leukemia.

Acute lymphoblastic leukemia (ALL) is the most common cancer in children, yet few environmental risk factors have been identified. We previously found an association between early-life tobacco smoke exposure and frequency of somatic deletions of 8 leukemia driver genes among childhood ALL patients in the California Childhood Leukemia Study. To expand analysis genome-wide and examine potential mechanisms, we conducted tumor whole-genome sequencing in 35 ALL patients, including 18 with high prenatal tobacco exposure and 17 with low exposure as determined by established epigenetic biomarkers. High tobacco exposure patients had significantly more structural variants (P < .001) and deletions (P = .001) genome-wide than low exposure patients. Investigation of off-target RAG recombination revealed that 41% of deletions in the high tobacco exposure patients were putatively RAG-mediated (full RAG motif identified at one or both breakpoints) compared with only 21% in the low exposure group (P = .001). In a multilevel model, deletions in high tobacco exposure patients were 2.44-fold (95% CI:1.13-5.38) more likely to be putatively RAG-mediated than deletions in low exposure patients. No point mutational signatures were associated with prenatal tobacco exposure. Our findings suggest that early-life tobacco smoke exposure may promote leukemogenesis by driving development of somatic deletions in pre-leukemic lymphocytes via off-target RAG recombination.

Genetic variation near GRB10 associated with bone growth and osteosarcoma risk in canine and human populations.

BACKGROUND: Canine and human osteosarcoma are similar in clinical presentation and tumor genomics. Giant breed dogs experience elevated osteosarcoma incidence, and taller stature remains a consistent risk factor for human osteosarcoma. Whether evolutionarily conserved genes contribute to both human and canine osteosarcoma predisposition merits evaluation. METHODS: A multi-center sample of childhood osteosarcoma patients and controls underwent genome-wide genotyping and imputation. Ancestry-adjusted SNP associations were calculated within each dataset using logistic regression, then meta-analyzed across the three datasets, totaling 1091 patients and 3026 controls. Ten regions previously associated with canine osteosarcoma risk were mapped to the human genome, spanning ∼6 Mb. We prioritized association testing of 5985 human SNPs mapping to candidate osteosarcoma risk regions detected in Irish wolfhounds, the largest dog breed studied. Secondary analyses explored 6289 additional human SNPs mapping to candidate osteosarcoma risk regions identified in Rottweilers and greyhounds. RESULTS: Fourteen SNPs were associated with human osteosarcoma risk after adjustment for multiple comparisons, all within a 42 kb region of human Chromosome 7p12.1. The lead variant was rs17454681 (OR=1.25, 95 %CI: 1.12-1.39; P=4.1×10-5), and independent risk variants were not observed in conditional analyses. While the associated region spanned 2.1 Mb and contained eight genes in Irish wolfhounds, associations were localized to a 50-fold smaller region of the human genome and strongly implicate GRB10 (growth factor receptor-bound protein 10) in canine and human osteosarcoma predisposition. PheWAS analysis in UK Biobank data identified noteworthy associations of the rs17454681 risk allele with varied measures of height and pubertal timing. CONCLUSIONS: Our comparative oncology analysis identified a novel human osteosarcoma risk allele near GRB10, a growth inhibitor that suppresses activated receptor tyrosine kinases including IGF1R, PDGFRB, and EGFR. Epidemiologists may benefit from leveraging cross-species comparisons to identify haplotypes in highly susceptible but genetically homogenous populations of domesticated animals, then fine-mapping these associations in diverse human populations.

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.

In Utero Origins of Acute Leukemia in Children.

Acute leukemias, mainly consisting of acute lymphoblastic leukemia (ALL) and acute myeloid leukemia (AML), comprise a major diagnostic group among hematologic cancers. Due to the early age at onset of ALL, particularly, it has long been suspected that acute leukemias of childhood may have an in utero origin. This supposition has motivated many investigations seeking direct proof of prenatal leukemogenesis, in particular, twin and "backtracking studies". The suspected in utero origin has also focused on gestation as a critical window of risk, resulting in a rich literature on prenatal risk factors for pediatric acute leukemias. In this narrative review, we recount the circumstantial and direct evidence for an in utero origin of childhood acute leukemias.

Evaluation of the genetic basis of familial-associated early-onset hematologic cancers in an ancestral/ethnically diverse population.

Genetic predisposition to hematologic malignancies has historically been addressed utilizing patients recruited from clinical trials and pedigrees constructed at major treatment centers. Such efforts leave unexplored the genetic basis of variations in risk by race/ethnic group shown in population-based surveillance data where cancer registration, compulsory by law, delivers universal enrollment. To address this, we performed exome sequencing on DNA isolated from newborn bloodspots derived from sibling pairs with early-onset cancers across California in which at least one of the siblings developed a hematologic cancer, using unbiased recruitment from the full state population. We identified pathogenic/likely pathogenic (P/ LP) variants among 1,172 selected cancer genes that were private or present at low allele frequencies in reference populations. Within 64 subjects from 32 families, we found 9 LP variants shared between siblings, and an additional 7 such variants in singleton children (not shared with their sibling). In 8 of the shared cases, the ancestral origin of the local haplotype that carries P/LP variants matched the dominant global ancestry of study participant families. This was the case for Latino sibling pairs on FLG and CBLB, non-Latino White sibling pairs in TP53 and NOD2, and a shared GATA2 variant for a non-Latino Black sibling pair. A new inherited mutation in HABP2 was identified in a sibling pair, one with diffuse large B-cell lymphoma and the other with neuroblastoma. Overall, the profile of P/LP germline variants across ancestral/ethnic groups suggests that rare alleles contributing to hematologic diseases originate within their race/ethnic origin parental populations, demonstrating the value of this discovery process in diverse, population-based registries.

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.

Gene-Environment Analyses Reveal Novel Genetic Candidates with Prenatal Tobacco Exposure in Relation to Risk for Childhood Acute Lymphoblastic Leukemia.

BACKGROUND: Associations between maternal tobacco exposure during pregnancy and childhood acute lymphoblastic leukemia (ALL) have yielded mixed results. This may be due to biases in self-reported smoking or other differences in individual-level risk factors. We utilized a biological marker of maternal tobacco exposure to evaluate the association between maternal tobacco exposure during pregnancy, genetics, and subsequent childhood ALL risk in two large population-based studies of childhood ALL in California. METHODS: Maternal exposure to tobacco smoke was assessed with a validated methylation marker (cg05575921) of the aryl hydrocarbon receptor repressor (AHRR) gene in newborn dried blood spots. We adjusted for sex, birthweight, gestational age, mode of delivery, year of birth, AHRR quantitative trait locus (mQTL) rs77111113, and a polygenetic risk score for childhood ALL. We additionally adjusted for principal components in a gene-environment interaction testing method that incorporates gene-only and environment-only effects along with interactions. RESULTS: AHRR hypomethylation overall was not associated with childhood ALL. In gene-environment interaction testing, several genetic variants displayed significant interaction with AHRR hypomethylation and childhood ALL. CONCLUSIONS: Our results suggest that novel candidates in PTPRK and DPP6 may play a role in tobacco-related leukemogenesis. Further research is necessary to better understand the effects of tobacco and these variants on childhood ALL risk. IMPACT: Despite the lack of an overall "main effect," tobacco exposure during pregnancy affects childhood ALL risk depending on specific genetic variants.

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.

Multi-ancestry genome-wide association study of 4069 children with glioma identifies 9p21.3 risk locus.

BACKGROUND: Although recent sequencing studies have revealed that 10% of childhood gliomas are caused by rare germline mutations, the role of common variants is undetermined and no genome-wide significant risk loci for pediatric central nervous system tumors have been identified to date. METHODS: Meta-analysis of 3 population-based genome-wide association studies comprising 4069 children with glioma and 8778 controls of multiple genetic ancestries. Replication was performed in a separate case-control cohort. Quantitative trait loci analyses and a transcriptome-wide association study were conducted to assess possible links with brain tissue expression across 18 628 genes. RESULTS: Common variants in CDKN2B-AS1 at 9p21.3 were significantly associated with astrocytoma, the most common subtype of glioma in children (rs573687, P-value of 6.974e-10, OR 1.273, 95% CI 1.179-1.374). The association was driven by low-grade astrocytoma (P-value of 3.815e-9) and exhibited unidirectional effects across all 6 genetic ancestries. For glioma overall, the association approached genome-wide significance (rs3731239, P-value of 5.411e-8), while no significant association was observed for high-grade tumors. Predicted decreased brain tissue expression of CDKN2B was significantly associated with astrocytoma (P-value of 8.090e-8). CONCLUSIONS: In this population-based genome-wide association study meta-analysis, we identify and replicate 9p21.3 (CDKN2B-AS1) as a risk locus for childhood astrocytoma, thereby establishing the first genome-wide significant evidence of common variant predisposition in pediatric neuro-oncology. We furthermore provide a functional basis for the association by showing a possible link to decreased brain tissue CDKN2B expression and substantiate that genetic susceptibility differs between low- and high-grade astrocytoma.

Associations between early-life and in utero infections and cytomegalovirus-positive acute lymphoblastic leukemia in children.

Childhood infections and cytomegalovirus (CMV) are associated with pediatric acute lymphoblastic leukemia (ALL). CMV dysregulates the host immune system and alters the immune response to subsequent antigenic exposures. We suspect that this immune dysregulation contributes to increased numbers of symptomatic infections in childhood allowing for expansion of pre-leukemic clones. We explored the association between childhood infections, maternal infections during pregnancy and CMV-positive ALL. Using a droplet digital PCR assay, we screened diagnostic ALL bone marrow samples from the California Childhood Leukemia Study (1995-2015) for the presence of CMV DNA identifying CMV-positive and CMV-negative cases. We performed a case-only analysis (n = 524) comparing the number and types of childhood infections and maternal infections during pregnancy between CMV-positive and CMV-negative ALL cases using logistic regression. With increasing numbers of infections in the first 12 months of life, children were more likely to classify to the highest tertile of CMV DNA in the bone marrow at diagnosis (OR: 1.04, 95% CI: 1.01-1.08). Specifically, those reporting cough or flu in the first 12 months were more likely to be CMV-positive at ALL diagnosis (OR: 2.15, 95% CI: 1.06-4.37 and OR: 2.06, 95% CI: 1.17-3.63 respectively). Furthermore, those with a history of maternal infection during pregnancy were more likely to be CMV-positive (OR: 2.12, 95% CI: 1.24-3.62). We hypothesize that children with underlying immune dysregulation develop more symptomatic infections in childhood and ultimately CMV-positive ALL; this underlying immune dysregulation may be due to early immune system alterations via CMV exposure (in utero or early infancy) proposing a potential link between CMV and ALL etiology.

Clonal hematopoiesis and risk of prostate cancer in large samples of European ancestry men.

Little is known regarding the potential relationship between clonal hematopoiesis (CH) of indeterminate potential (CHIP), which is the expansion of hematopoietic stem cells with somatic mutations, and risk of prostate cancer, the fifth leading cause of cancer death of men worldwide. We evaluated the association of age-related CHIP with overall and aggressive prostate cancer risk in two large whole-exome sequencing studies of 75 047 European ancestry men, including 7663 prostate cancer cases, 2770 of which had aggressive disease, and 3266 men carrying CHIP variants. We found that CHIP, defined by over 50 CHIP genes individually and in aggregate, was not significantly associated with overall (aggregate HR = 0.93, 95% CI = 0.76-1.13, P = 0.46) or aggressive (aggregate OR = 1.14, 95% CI = 0.92-1.41, P = 0.22) prostate cancer risk. CHIP was weakly associated with genetic risk of overall prostate cancer, measured using a polygenic risk score (OR = 1.05 per unit increase, 95% CI = 1.01-1.10, P = 0.01). CHIP was not significantly associated with carrying pathogenic/likely pathogenic/deleterious variants in DNA repair genes, which have previously been found to be associated with aggressive prostate cancer. While findings from this study suggest that CHIP is likely not a risk factor for prostate cancer, it will be important to investigate other types of CH in association with prostate cancer risk.

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

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

DNA methylation at birth in monozygotic twins discordant for pediatric acute lymphoblastic leukemia.

Aberrant DNA methylation constitutes a key feature of pediatric acute lymphoblastic leukemia at diagnosis, however its role as a predisposing or early contributor to leukemia development remains unknown. Here, we evaluate DNA methylation at birth in 41 leukemia-discordant monozygotic twin pairs using the Illumina EPIC array on archived neonatal blood spots to identify epigenetic variation associated with development of pediatric acute lymphoblastic leukemia, independent of genetic influence. Through conditional logistic regression we identify 240 significant probes and 10 regions associated with the discordant onset of leukemia. We identify a significant negative coefficient bias, indicating DNA hypomethylation in cases, across the array and enhanced in open sea, shelf/shore, and gene body regions compared to promoter and CpG island regions. Here, we show an association between global DNA hypomethylation and future development of pediatric acute lymphoblastic leukemia across disease-discordant genetically identical twins, implying DNA hypomethylation may contribute more generally to leukemia risk.

Periconceptional folate intake influences DNA methylation at birth based on dietary source in an analysis of pediatric acute lymphoblastic leukemia cases and controls.

BACKGROUND: Periconceptional folate intake is associated with the establishment of DNA methylation in offspring; however, variations in this relation by food sources compared with folic acid supplements are not described. Also, maternal folate intake is associated with decreased risk of pediatric acute lymphoblastic leukemia (ALL), but the mechanism is not known. OBJECTIVES: We evaluated the relation between periconceptional folate intake by source and DNA methylation at birth in a cohort of pediatric ALL cases and controls in an epigenome-wide association study. METHODS: Genome-wide DNA methylation status obtained from archived neonatal blood spots from pediatric ALL cases (n = 189) and controls (n = 205) in the California Childhood Leukemia Study (CCLS) from 1995-2008 was compared with periconceptional folate from total, food, and supplemental sources using multivariable linear regression. Further stratification was performed by income, education, ethnicity, and total folate intake. We evaluated variable DNA methylation response to periconceptional folate by ALL case status through an interaction term. RESULTS: Two significant differentially methylated probes (DMPs) were associated with food and supplemental periconceptional folate intake in all subjects (n = 394). The top differentially methylated region at the promoter region of DUSP22(dual specificity phosphatase 22) demonstrated DNA hypermethylation in ALL cases but not in controls in response to total and food folate intake. We further identified 8 interaction term DMPs with variable DNA methylation response to folate intake by ALL case status. Further stratification of the cohort by education and ethnicity revealed a substantially higher number of DMPs associated with supplemental folic acid intake in Hispanic subjects with lower income and educational level. CONCLUSIONS: We identified modest associations between periconceptional folate intake and DNA methylation differing by source, including variation by ALL case status. Hispanic subjects of lower income and education appear uniquely responsive to periconceptional folate supplementation.

Localized variation in ancestral admixture identifies pilocytic astrocytoma risk loci among Latino children.

BACKGROUND: Pilocytic astrocytoma (PA) is the most common pediatric brain tumor. PA has at least a 50% higher incidence in populations of European ancestry compared to other ancestral groups, which may be due in part to genetic differences. METHODS: We first compared the global proportions of European, African, and Amerindian ancestries in 301 PA cases and 1185 controls of self-identified Latino ethnicity from the California Biobank. We then conducted admixture mapping analysis to assess PA risk with local ancestry. RESULTS: We found PA cases had a significantly higher proportion of global European ancestry than controls (case median = 0.55, control median = 0.51, P value = 3.5x10-3). Admixture mapping identified 13 SNPs in the 6q14.3 region (SNX14) contributing to risk, as well as three other peaks approaching significance on chromosomes 7, 10 and 13. Downstream fine mapping in these regions revealed several SNPs potentially contributing to childhood PA risk. CONCLUSIONS: There is a significant difference in genomic ancestry associated with Latino PA risk and several genomic loci potentially mediating this risk.

Ras-mutant cancers are sensitive to small molecule inhibition of V-type ATPases in mice.

Mutations in Ras family proteins are implicated in 33% of human cancers, but direct pharmacological inhibition of Ras mutants remains challenging. As an alternative to direct inhibition, we screened for sensitivities in Ras-mutant cells and discovered 249C as a Ras-mutant selective cytotoxic agent with nanomolar potency against a spectrum of Ras-mutant cancers. 249C binds to vacuolar (V)-ATPase with nanomolar affinity and inhibits its activity, preventing lysosomal acidification and inhibiting autophagy and macropinocytosis pathways that several Ras-driven cancers rely on for survival. Unexpectedly, potency of 249C varies with the identity of the Ras driver mutation, with the highest potency for KRASG13D and G12V both in vitro and in vivo, highlighting a mutant-specific dependence on macropinocytosis and lysosomal pH. Indeed, 249C potently inhibits tumor growth without adverse side effects in mouse xenografts of KRAS-driven lung and colon cancers. A comparison of isogenic SW48 xenografts with different KRAS mutations confirmed that KRASG13D/+ (followed by G12V/+) mutations are especially sensitive to 249C treatment. These data establish proof-of-concept for targeting V-ATPase in cancers driven by specific KRAS mutations such as KRASG13D and G12V.

Genetic Risk of Second Primary Cancer in Breast Cancer Survivors: The Multiethnic Cohort Study.

Women who have had breast cancer in the past are at increased risk of developing a second primary cancer (SPC), including second primary breast cancer (SPBC) or a second primary non-breast cancer (SPNBC). In the Multiethnic Cohort (MEC) Study, we conducted a prospective cohort analysis in 3,223 female breast cancer survivors from five racial/ethnic populations (White, African American, Japanese American, Latino, and Native Hawaiian) to assess the association of rare pathogenic variants (PV) in 37 known cancer predisposition genes with risk of SPC. A total of 719 (22.3%) women developed SPC, of which, 323 (10.0%) were SPBC. Germline PVs in BRCA1 (HR, 2.28; 95% CI, 1.11-4.65) and ERCC2 (HR, 3.51; 95% CI, 1.29-9.54) were significantly enriched in women with SPC. In the subtype analysis for SPBC, a significant association of ERCC2 PVs (HR, 5.09; 95% CI, 1.58-16.4) and a suggestive association of BRCA2 PVs (HR, 2.24; 95% CI, 0.91-5.55) were observed. There was also a higher risk of SPNBC in carriers of BRCA1 PVs (HR, 2.98; 95% CI, 1.21-7.36). These results provide evidence that germline PVs in BRCA1, BRCA2, and ERCC2 contribute to the development of SPC in breast cancer survivors. These findings also suggest that compromised DNA repair mechanisms could be a predisposition factor for SPC in patients with breast cancer, supporting the need for closer monitoring of SPC in women carrying PVs in these genes. SIGNIFICANCE: This multiethnic study links germline pathogenic variants in BRCA1, BRCA2, and ERCC2 to the development of second primary cancer in breast cancer survivors, providing biological insights and biomarkers to guide patient monitoring.

Variant to function mapping at single-cell resolution through network propagation.

Genome-wide association studies in combination with single-cell genomic atlases can provide insights into the mechanisms of disease-causal genetic variation. However, identification of disease-relevant or trait-relevant cell types, states and trajectories is often hampered by sparsity and noise, particularly in the analysis of single-cell epigenomic data. To overcome these challenges, we present SCAVENGE, a computational algorithm that uses network propagation to map causal variants to their relevant cellular context at single-cell resolution. We demonstrate how SCAVENGE can help identify key biological mechanisms underlying human genetic variation, applying the method to blood traits at distinct stages of human hematopoiesis, to monocyte subsets that increase the risk for severe Coronavirus Disease 2019 (COVID-19) and to intermediate lymphocyte developmental states that predispose to acute leukemia. Our approach not only provides a framework for enabling variant-to-function insights at single-cell resolution but also suggests a more general strategy for maximizing the inferences that can be made using single-cell genomic data.