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.

Association of immunoglobulin E levels with glioma risk and survival.

Background: Previous epidemiological studies have reported an association of serum immunoglobulin E (IgE) levels with reduced glioma risk, but the association between IgE and glioma prognosis is not well characterized. This study aimed to examine how sex, tumor subtype, and IgE class modulate the association of serum IgE levels with glioma risk and survival. Methods: We conducted a case-control study using participants from the University of California, San Francisco Adult Glioma Study (1997-2010). Serum IgE levels for total, respiratory and food allergy were measured in adults diagnosed with glioma (n=1,696) and cancer-free controls (n=1,135) matched based on age, sex, and race/ethnicity. Logistic regression was adjusted for patient demographics to assess the association between IgE levels and glioma risk. Multivariable Cox regression adjusted for patient-specific and tumor-specific factors compared survival between the elevated and normal IgE groups. Results: Elevated total IgE was associated with reduced risk of IDH wildtype (OR=0.65, 95% CI: 0.54-0.78) and IDH mutant glioma (OR=0.65, 95% CI: 0.50-0.85). In multivariable Cox regression, elevated respiratory IgE was associated with improved survival for IDH wildtype glioma (HR=0.78, 95% CI: 0.67-0.91). The reduction in mortality risk was more pronounced in females (HR=0.71, 95% CI: 0.53-0.96) than in males (HR=0.80, 95% CI: 0.66-0.97), with improvements in median survival of 6.2 months (P<.001) and 1.6 months (P=0.003), respectively. Conclusion: Elevated serum IgE was associated with improved prognosis for IDH wildtype glioma, with a more pronounced protective effect in females. These results suggest a possible sexual dimorphism and antitumor activity of IgE-mediated immune responses.

Genetic predisposition to altered blood cell homeostasis is associated with glioma risk and survival.

Glioma is a highly fatal brain tumor comprised of molecular subtypes with distinct clinical trajectories. Observational studies have suggested that variability in immune response may play a role in glioma etiology. However, their findings have been inconsistent and susceptible to reverse causation due to treatment effects and the immunosuppressive nature of glioma. We applied genetic variants associated (p<5×10-8) with blood cell traits to a meta-analysis of 3418 glioma cases and 8156 controls. Genetically predicted increase in the platelet to lymphocyte ratio (PLR) was associated with an increased risk of glioma (odds ratio (OR)=1.25, p=0.005), especially in IDH-mutant (IDHmut OR=1.38, p=0.007) and IDHmut 1p/19q non-codeleted (IDHmut-noncodel OR=1.53, p=0.004) tumors. However, reduced glioma risk was observed for higher counts of lymphocytes (IDHmut-noncodel OR=0.70, p=0.004) and neutrophils (IDHmut OR=0.69, p=0.019; IDHmut-noncodel OR=0.60, p=0.009), which may reflect genetic predisposition to enhanced immune-surveillance. In contrast to susceptibility, there was no association with survival in IDHmut-noncodel; however, in IDHmut 1p/19q co-deleted tumors, we observed higher mortality with increasing genetically predicted counts of lymphocytes (hazard ratio (HR)=1.65, 95% CI: 1.24-2.20), neutrophils (HR=1.49, 1.13-1.97), and eosinophils (HR=1.59, 1.18-2.14). Polygenic scores for blood cell traits were also associated with tumor immune microenvironment features, with heterogeneity by IDH status observed for 17 signatures related to interferon signaling, PD-1 expression, and T-cell/Cytotoxic responses. In summary, we identified novel, immune-mediated susceptibility mechanisms for glioma with potential disease management implications.

Antibodies to varicella-zoster virus and three other herpesviruses and survival in adults with glioma.

BACKGROUND: Lifetime exposure to the varicella-zoster virus (VZV) has been consistently inversely associated with glioma risk, however, the relationship of VZV with survival in adults with glioma has not been investigated. In this study, we analyzed the survival of adults with glioma in relation to their antibody measurements to 4 common herpes viral infections, including VZV, measured post-diagnosis. METHODS: We analyzed IgG antibody measurements to VZV, cytomegalovirus (CMV), herpes simplex virus 1/2 (HSV), and Epstein-Barr virus (EBV) collected from 1378 adults with glioma diagnosed between 1991 and 2010. Blood was obtained a median of 3 months after surgery. Associations of patient IgG levels with overall survival were estimated using Cox models adjusted for age, sex, self-reported race, surgery type, dexamethasone usage at blood draw, and tumor grade. Models were stratified by recruitment series and meta-analyzed to account for time-dependent treatment effects. RESULTS: VZV antibody seropositivity was associated with improved survival outcomes in adults with glioma (Hazard ratio, HR = 0.70, 95% Confidence Interval 0.54-0.90, P = .006). Amongst cases who were seropositive for VZV antibodies, survival was significantly improved for those above the 25th percentile of continuous reactivity measurements versus those below (HR = 0.76, 0.66-0.88, P = .0003). Antibody seropositivity to EBV was separately associated with improved survival (HR = 0.71, 0.53-0.96, P = .028). Antibody positivity to 2 other common viruses (CMV, HSV) was not associated with altered survival. CONCLUSIONS: Low levels of VZV or EBV antibodies are associated with poorer survival outcomes for adults with glioma. Differential immune response rather than viral exposure may explain these findings.

Enrichment of a neutrophil-like monocyte transcriptional state in glioblastoma myeloid suppressor cells.

Glioblastomas (GBM) are lethal central nervous system cancers associated with tumor and systemic immunosuppression. Heterogeneous monocyte myeloid-derived suppressor cells (M-MDSC) are implicated in the altered immune response in GBM, but M-MDSC ontogeny and definitive phenotypic markers are unknown. Using single-cell transcriptomics, we revealed heterogeneity in blood M-MDSC from GBM subjects and an enrichment in a transcriptional state reminiscent of neutrophil-like monocytes (NeuMo), a newly described pathway of monopoiesis in mice. Human NeuMo gene expression and Neu-like deconvolution fraction algorithms were created to quantitate the enrichment of this transcriptional state in GBM subjects. NeuMo populations were also observed in M-MDSCs from lung and head and neck cancer subjects. Dexamethasone (DEX) and prednisone exposures increased the usage of Neu-like states, which were inversely associated with tumor purity and survival in isocitrate dehydrogenase wildtype (IDH WT) gliomas. Anti-inflammatory ZC3HA12/Regnase-1 transcripts were highly correlated with NeuMo expression in tumors and in blood M-MDSC from GBM, lung, and head and neck cancer subjects. Additional novel transcripts of immune-modulating proteins were identified. Collectively, these findings provide a framework for understanding the heterogeneity of M-MDSCs in GBM as cells with different clonal histories and may reshape approaches to study and therapeutically target these cells.

Evaluation of cross-platform compatibility of a DNA methylation-based glucocorticoid response biomarker.

BACKGROUND: Identifying blood-based DNA methylation patterns is a minimally invasive way to detect biomarkers in predicting age, characteristics of certain diseases and conditions, as well as responses to immunotherapies. As microarray platforms continue to evolve and increase the scope of CpGs measured, new discoveries based on the most recent platform version and how they compare to available data from the previous versions of the platform are unknown. The neutrophil dexamethasone methylation index (NDMI 850) is a blood-based DNA methylation biomarker built on the Illumina MethylationEPIC (850K) array that measures epigenetic responses to dexamethasone (DEX), a synthetic glucocorticoid often administered for inflammation. Here, we compare the NDMI 850 to one we built using data from the Illumina Methylation 450K (NDMI 450). RESULTS: The NDMI 450 consisted of 22 loci, 15 of which were present on the NDMI 850. In adult whole blood samples, the linear composite scores from NDMI 450 and NDMI 850 were highly correlated and had equivalent predictive accuracy for detecting DEX exposure among adult glioma patients and non-glioma adult controls. However, the NDMI 450 scores of newborn cord blood were significantly lower than NDMI 850 in samples measured with both assays. CONCLUSIONS: We developed an algorithm that reproduces the DNA methylation glucocorticoid response score using 450K data, increasing the accessibility for researchers to assess this biomarker in archived or publicly available datasets that use the 450K version of the Illumina BeadChip array. However, the NDMI850 and NDMI450 do not give similar results in cord blood, and due to data availability limitations, results from sample types of newborn cord blood should be interpreted with care.

DNA methylation as a pharmacodynamic marker of glucocorticoid response and glioma survival.

Assessing individual responses to glucocorticoid drug therapies that compromise immune status and affect survival outcomes in neuro-oncology is a great challenge. Here we introduce a blood-based neutrophil dexamethasone methylation index (NDMI) that provides a measure of the epigenetic response of subjects to dexamethasone. This marker outperforms conventional approaches based on leukocyte composition as a marker of glucocorticoid response. The NDMI is associated with low CD4 T cells and the accumulation of monocytic myeloid-derived suppressor cells and also serves as prognostic factor in glioma survival. In a non-glioma population, the NDMI increases with a history of prednisone use. Therefore, it may also be informative in other conditions where glucocorticoids are employed. We conclude that DNA methylation remodeling within the peripheral immune compartment is a rich source of clinically relevant markers of glucocorticoid response.

Accelerated epigenetic aging in newborns with Down syndrome.

Accelerated aging is a hallmark of Down syndrome (DS), with adults experiencing early-onset Alzheimer's disease and premature aging of the skin, hair, and immune and endocrine systems. Accelerated epigenetic aging has been found in the blood and brain tissue of adults with DS but when premature aging in DS begins remains unknown. We investigated whether accelerated aging in DS is already detectable in blood at birth. We assessed the association between age acceleration and DS using five epigenetic clocks in 346 newborns with DS and 567 newborns without DS using Illumina MethylationEPIC DNA methylation array data. We compared two epigenetic aging clocks (DNAmSkinBloodClock and pan-tissue DNAmAge) and three epigenetic gestational age clocks (Haftorn, Knight, and Bohlin) between DS and non-DS newborns using linear regression adjusting for observed age, sex, batch, deconvoluted blood cell proportions, and genetic ancestry. Targeted sequencing of GATA1 was performed in a subset of 184 newborns with DS to identify somatic mutations associated with transient abnormal myelopoiesis. DS was significantly associated with increased DNAmSkinBloodClock (effect estimate = 0.2442, p < 0.0001), with an epigenetic age acceleration of 244 days in newborns with DS after adjusting for potential confounding factors (95% confidence interval: 196-292 days). We also found evidence of epigenetic age acceleration associated with somatic GATA1 mutations among newborns with DS (p = 0.015). DS was not associated with epigenetic gestational age acceleration. We demonstrate that accelerated epigenetic aging in the blood of DS patients begins prenatally, with implications for the pathophysiology of immunosenescence and other aging-related traits in DS.

Pleiotropic MLLT10 variation confers risk of meningioma and estrogen-mediated cancers.

Background: Risk of tumors of the breast, ovary, and meninges has been associated with hormonal factors and with one another. Genome-wide association studies (GWAS) identified a meningioma risk locus on 10p12 near previous GWAS hits for breast and ovarian cancers, raising the possibility of genetic pleiotropy. Methods: We performed imputation-based fine-mapping in three case-control datasets of meningioma (927 cases, 790 controls), female breast cancer (28 108 cases, 22 209 controls), and ovarian cancer (25 509 cases, 40 941 controls). Analyses were stratified by sex (meningioma), estrogen receptor (ER) status (breast), and histotype (ovarian), then combined using subset-based meta-analysis in ASSET. Lead variants were assessed for association with additional traits in UK Biobank to identify potential effect-mediators. Results: Two-sided subset-based meta-analysis identified rs7084454, an expression quantitative trait locus (eQTL) near the MLLT10 promoter, as lead variant (5.7 × 10-14). The minor allele was associated with increased risk of meningioma in females (odds ratio (OR) = 1.42, 95% Confidence Interval (95%CI):1.20-1.69), but not males (OR = 1.19, 95%CI: 0.91-1.57). It was positively associated with ovarian (OR = 1.09, 95%CI:1.06-1.12) and ER+ breast (OR = 1.05, 95%CI: 1.02-1.08) cancers, and negatively associated with ER- breast cancer (OR = 0.91, 95%CI: 0.86-0.96). It was also associated with several adiposity traits (P < 5.0 × 10-8), but adjusting for body mass index did not attenuate its association with meningioma. MLLT10 and ESR1 expression were positively correlated in normal meninges (P = .058) and meningioma tumors (P = .0065). Conclusions: We identify a MLLT10 eQTL positively associated with risk of female meningioma, ER+ breast cancer, ovarian cancer, and obesity, and implicate a potential estrogenic mechanism underlying this pleiotropy.

Pre-surgery immune profiles of adult glioma patients.

INTRODUCTION: Although immunosuppression is a known characteristic of glioma, no previous large studies have reported peripheral blood immune cell profiles prior to patient surgery and chemoradiation. This report describes blood immune cell characteristics and associated variables prior to surgery among typical glioma patients seen at a large University practice. METHODS: We analyzed pre-surgery blood samples from 139 glioma patients diagnosed with a new or recurrent grade II/III glioma (LrGG, n = 64) or new glioblastoma (GBM, n = 75) and 454 control participants without glioma. Relative cell fractions of CD4, CD8, B-cells, Natural Killer cells, monocytes, and neutrophils, were estimated via a validated deconvolution algorithm from blood DNA methylation measures from Illumina EPIC arrays. RESULTS: Dexamethasone use at time of blood draw varied by glioma type being highest among patients with IDH wild-type (wt) GBM (75%) and lowest for those with oligodendroglioma (14%). Compared to controls, glioma patients showed statistically significant lower cell fractions for all immune cell subsets except for neutrophils which were higher (all p-values < 0.001), in part because of the higher prevalence of dexamethasone use at time of blood draw for IDHwt GBM. Patients who were taking dexamethasone were more likely to have a low CD4 count (< 200, < 500), increased neutrophils, low absolute lymphocyte counts, higher total cell count and higher NLR. CONCLUSION: We show that pre-surgery blood immune profiles vary by glioma subtype, age, and more critically, by use of dexamethasone. Our results highlight the importance of considering dexamethasone exposures in all studies of immune profiles and of obtaining immune measures prior to use of dexamethasone, if possible.

The immunogenetics of viral antigen response is associated with subtype-specific glioma risk and survival.

Glioma is a highly fatal cancer with prognostically significant molecular subtypes and few known risk factors. Multiple studies have implicated infections in glioma susceptibility, but evidence remains inconsistent. Genetic variants in the human leukocyte antigen (HLA) region modulate host response to infection and have been linked to glioma risk. In this study, we leveraged genetic predictors of antibody response to 12 viral antigens to investigate the relationship with glioma risk and survival. Genetic reactivity scores (GRSs) for each antigen were derived from genome-wide-significant (p < 5 × 10-8) variants associated with immunoglobulin G antibody response in the UK Biobank cohort. We conducted parallel analyses of glioma risk and survival for each GRS and HLA alleles imputed at two-field resolution by using data from 3,418 glioma-affected individuals subtyped by somatic mutations and 8,156 controls. Genetic reactivity scores to Epstein-Barr virus (EBV) ZEBRA and EBNA antigens and Merkel cell polyomavirus (MCV) VP1 antigen were associated with glioma risk and survival (Bonferroni-corrected p < 0.01). GRSZEBRA and GRSMCV were associated in opposite directions with risk of IDH wild-type gliomas (ORZEBRA = 0.91, p = 0.0099/ORMCV = 1.11, p = 0.0054). GRSEBNA was associated with both increased risk for IDH mutated gliomas (OR = 1.09, p = 0.040) and improved survival (HR = 0.86, p = 0.010). HLA-DQA1∗03:01 was significantly associated with decreased risk of glioma overall (OR = 0.85, p = 3.96 × 10-4) after multiple testing adjustment. This systematic investigation of the role of genetic determinants of viral antigen reactivity in glioma risk and survival provides insight into complex immunogenomic mechanisms of glioma pathogenesis. These results may inform applications of antiviral-based therapies in glioma treatment.

Interaction between maternal killer immunoglobulin-like receptors and offspring HLAs and susceptibility of childhood ALL.

Acute lymphoblastic leukemia (ALL) in children is associated with a distinct neonatal cytokine profile. The basis of this neonatal immune phenotype is unknown but potentially related to maternal-fetal immune receptor interactions. We conducted a case-control study of 226 case child-mother pairs and 404 control child-mother pairs to evaluate the role of interaction between HLA genotypes in the offspring and maternal killer immunoglobulin-like receptor (KIR) genotypes in the etiology of childhood ALL, while considering potential mediation by neonatal cytokines and the immune-modulating enzyme arginase-II (ARG-II). We observed different associations between offspring HLA-maternal KIR activating profiles and the risk of ALL in different predicted genetic ancestry groups. For instance, in Latino subjects who experience the highest risk of childhood leukemia, activating profiles were significantly associated with a lower risk of childhood ALL (odds ratio [OR] = 0.59; 95% confidence interval [CI], 0.49-0.71) and a higher level of ARG-II at birth (coefficient = 0.13; 95% CI, 0.04-0.22). HLA-KIR activating profiles were also associated with a lower risk of ALL in non-Latino Asians (OR = 0.63; 95% CI, 0.38-1.01), although they had a lower tumor necrosis factor-α level (coefficient = -0.27; 95% CI, -0.49 to -0.06). Among non-Latino White subjects, no significant association was observed between offspring HLA-maternal KIR interaction and ALL risk or cytokine levels. The current study reports the association between offspring HLA-maternal KIR interaction and the development of childhood ALL with variation by predicted genetic ancestry. We also observed some associations between activating profiles and immune factors related to cytokine control; however, cytokines did not demonstrate causal mediation of the activating profiles on ALL risk.

Enhanced cell deconvolution of peripheral blood using DNA methylation for high-resolution immune profiling.

DNA methylation microarrays can be employed to interrogate cell-type composition in complex tissues. Here, we expand reference-based deconvolution of blood DNA methylation to include 12 leukocyte subtypes (neutrophils, eosinophils, basophils, monocytes, naïve and memory B cells, naïve and memory CD4 + and CD8 + T cells, natural killer, and T regulatory cells). Including derived variables, our method provides 56 immune profile variables. The IDOL (IDentifying Optimal Libraries) algorithm was used to identify libraries for deconvolution of DNA methylation data for current and previous platforms. The accuracy of deconvolution estimates obtained using our enhanced libraries was validated using artificial mixtures and whole-blood DNA methylation with known cellular composition from flow cytometry. We applied our libraries to deconvolve cancer, aging, and autoimmune disease datasets. In conclusion, these libraries enable a detailed representation of immune-cell profiles in blood using only DNA and facilitate a standardized, thorough investigation of immune profiles in human health and disease.

Interactions of Age and Blood Immune Factors and Noninvasive Prediction of Glioma Survival.

BACKGROUND: Tumor-based classification of human glioma portends patient prognosis, but considerable unexplained survival variability remains. Host factors (eg, age) also strongly influence survival times, partly reflecting a compromised immune system. How blood epigenetic measures of immune characteristics and age augment molecular classifications in glioma survival has not been investigated. We assess the prognostic impact of immune cell fractions and epigenetic age in archived blood across glioma molecular subtypes for the first time. METHODS: We evaluated immune cell fractions and epigenetic age in archived blood from the University of California San Francisco Adult Glioma Study, which included a training set of 197 patients with IDH-wild type, 1p19q intact, TERT wild type (IDH/1p19q/TERT-WT) glioma, an evaluation set of 350 patients with other subtypes of glioma, and 454 patients without glioma. RESULTS: IDH/1p19q/TERT-WT patients had lower lymphocyte fractions (CD4+ T, CD8+ T, natural killer, and B cells) and higher neutrophil fractions than people without glioma. Recursive partitioning analysis delineated 4 statistically significantly different survival groups for patients with IDH/1p19q/TERT-WT based on an interaction between chronological age and 2 blood immune factors, CD4+ T cells, and neutrophils. Median overall survival ranged from 0.76 years (95% confidence interval = 0.55-0.99) for the worst survival group (n = 28) to 9.72 years (95% confidence interval = 6.18 to not available) for the best (n = 33). The recursive partitioning analysis also statistically significantly delineated 4 risk groups in patients with other glioma subtypes. CONCLUSIONS: The delineation of different survival groups in the training and evaluation sets based on an interaction between chronological age and blood immune characteristics suggests that common host immune factors among different glioma types may affect survival. The ability of DNA methylation-based markers of immune status to capture diverse, clinically relevant information may facilitate noninvasive, personalized patient evaluation in the neuro-oncology clinic.

The genome-wide impact of trisomy 21 on DNA methylation and its implications for hematopoiesis.

Down syndrome is associated with genome-wide perturbation of gene expression, which may be mediated by epigenetic changes. We perform an epigenome-wide association study on neonatal bloodspots comparing 196 newborns with Down syndrome and 439 newborns without Down syndrome, adjusting for cell-type heterogeneity, which identifies 652 epigenome-wide significant CpGs (P < 7.67 × 10-8) and 1,052 differentially methylated regions. Differential methylation at promoter/enhancer regions correlates with gene expression changes in Down syndrome versus non-Down syndrome fetal liver hematopoietic stem/progenitor cells (P < 0.0001). The top two differentially methylated regions overlap RUNX1 and FLI1, both important regulators of megakaryopoiesis and hematopoietic development, with significant hypermethylation at promoter regions of these two genes. Excluding Down syndrome newborns harboring preleukemic GATA1 mutations (N = 30), identified by targeted sequencing, has minimal impact on the epigenome-wide association study results. Down syndrome has profound, genome-wide effects on DNA methylation in hematopoietic cells in early life, which may contribute to the high frequency of hematological problems, including leukemia, in children with Down syndrome.

Genetic predisposition to longer telomere length and risk of childhood, adolescent and adult-onset ependymoma.

Ependymoma is the third most common brain tumor in children, with well-described molecular characterization but poorly understood underlying germline risk factors. To investigate whether genetic predisposition to longer telomere length influences ependymoma risk, we utilized case-control data from three studies: a population-based pediatric and adolescent ependymoma case-control sample from California (153 cases, 696 controls), a hospital-based pediatric posterior fossa type A (EPN-PF-A) ependymoma case-control study from Toronto's Hospital for Sick Children and the Children's Hospital of Philadelphia (83 cases, 332 controls), and a multicenter adult-onset ependymoma case-control dataset nested within the Glioma International Case-Control Consortium (GICC) (103 cases, 3287 controls). In the California case-control sample, a polygenic score for longer telomere length was significantly associated with increased risk of ependymoma diagnosed at ages 12-19 (P = 4.0 × 10-3), but not with ependymoma in children under 12 years of age (P = 0.94). Mendelian randomization supported this observation, identifying a significant association between genetic predisposition to longer telomere length and increased risk of adolescent-onset ependymoma (ORPRS = 1.67; 95% CI 1.18-2.37; P = 3.97 × 10-3) and adult-onset ependymoma (PMR-Egger = 0.042), but not with risk of ependymoma diagnosed before age 12 (OR = 1.12; 95% CI 0.94-1.34; P = 0.21), nor with EPN-PF-A (PMR-Egger = 0.59). These findings complement emerging literature suggesting that augmented telomere maintenance is important in ependymoma pathogenesis and progression, and that longer telomere length is a risk factor for diverse nervous system malignancies.

European genetic ancestry associated with risk of childhood ependymoma.

BACKGROUND: Ependymoma is a histologically defined central nervous system tumor most commonly occurring in childhood. Population-level incidence differences by race/ethnicity are observed, with individuals of European ancestry at highest risk. We aimed to determine whether extent of European genetic ancestry is associated with ependymoma risk in US populations. METHODS: In a multi-ethnic study of Californian children (327 cases, 1970 controls), we estimated the proportions of European, African, and Native American ancestry among recently admixed Hispanic and African American subjects and estimated European admixture among non-Hispanic white subjects using genome-wide data. We tested whether genome-wide ancestry differences were associated with ependymoma risk and performed admixture mapping to identify associations with local ancestry. We also evaluated race/ethnicity-stratified ependymoma incidence data from the Central Brain Tumor Registry of the United States (CBTRUS). RESULTS: CBTRUS data revealed that African American and Native American children have 33% and 36%, respectively, reduced incidence of ependymoma compared with non-Hispanic whites. In genetic analyses, a 20% increase in European ancestry was associated with a 1.31-fold higher odds of ependymoma among self-reported Hispanics and African Americans (95% CI: 1.08-1.59, Pmeta = 6.7 × 10-3). Additionally, eastern European ancestral substructure was associated with increased ependymoma risk in non-Hispanic whites (P = 0.030) and in Hispanics (P = 0.043). Admixture mapping revealed a peak at 20p13 associated with increased local European ancestry, and targeted fine-mapping identified a lead variant at rs6039499 near RSPO4 (odds ratio = 1.99; 95% CI: 1.45-2.73; P = 2.2 × 10-5) but which was not validated in an independent set of posterior fossa type A patients. CONCLUSIONS: Interethnic differences in ependymoma risk are recapitulated in the genomic ancestry of ependymoma patients, implicating regions to target in future association studies.