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.

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.

Inherited genetic susceptibility to acute lymphoblastic leukemia in Down syndrome.

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

Burden of illness in not adequately controlled chronic hypoparathyroidism: Findings from a 13-country patient and caregiver survey.

OBJECTIVE: To address knowledge gaps regarding burdens associated with not adequately controlled chronic hypoparathyroidism. DESIGN: Global patient and caregiver survey. STUDY POPULATIONS: Patients with chronic hypoparathyroidism not adequately controlled on conventional therapy and their caregivers. MEASUREMENTS: Health-related quality of life (HRQoL) and health status were evaluated using the 36-item Short Form version 2 (SF-36 v2.0) and Five-Level EuroQoL 5 Dimensions (EQ-5D-5L) instruments, respectively. Hypoparathyroidism-associated symptoms were assessed by a disease-specific Hypoparathyroidism Symptom Diary and caregiver burden via the Modified Caregiver Strain Index (MCSI). RESULTS: Data were obtained from 398 patients and 207 caregivers. Patients' self-rated hypoparathyroidism-related symptom severity was none (3%), mild (32%), moderate (53%) or severe (12%). Per the Hypoparathyroidism Symptom Diary, patients reported moderate, severe or very severe symptoms of physical fatigue (73%), muscle cramps (55%), heaviness in limbs (55%) and tingling (51%) over a 7-day recall period. Impacts (rated 'somewhat' or 'very much') were reported by 84% of patients for ability to exercise, 78% for sleep, 75% for ability to work and 63% for family relationships. Inverse relationships were observed between patient self-rated overall symptom severity and HRQoL and health status assessment scores-the greater the symptom severity, the lower the SF-36 and EQ-5D-5L scores. Caregiver burden increased with patient self-rated symptom severity: none, 1.7 MCSI; mild, 5.4 MCSI; moderate, 9.5 MCSI; and severe, 12.5 MCSI. CONCLUSION: Patients with not adequately controlled hypoparathyroidism reported substantial symptoms and impacts. Greater patient symptom severity was associated with decreased patient HRQoL and health status assessments and increased caregiver burden.

Predisposing germline mutations in high hyperdiploid acute lymphoblastic leukemia in children.

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

Longer genotypically-estimated leukocyte telomere length is associated with increased meningioma risk.

PURPOSE: Telomere length-associated SNPs have been associated with incidence and survival rates for malignant brain tumors such as glioma. Here, we study the influence of genetically determined lymphocyte telomere length (LTL) by comparing telomerase associated SNPs between the most common non-malignant brain tumor, i.e. meningioma, and healthy controls. METHODS/PATIENTS: One thousand fifty-three (1053) surgically treated meningioma patients and 4437 controls of Western European ancestry were included. Germline DNA was genotyped for 8 SNPs previously significantly associated with LTL. Genotypically-estimated LTL was then calculated by summing each SNP's genotypically-specified telomere length increase in base pairs (bp) for each person. Odds ratios for genotypically-estimated LTL in meningioma cases and controls were evaluated using logistic regression with the first two ancestral principal components and sex as covariates. RESULTS: Three out of the eight evaluated LTL SNPs were significantly associated with increased meningioma risk (rs10936599: OR 1.14, 95% CI 1.01-1.28, rs2736100: OR 1.13, 95% CI 1.03-1.25, rs9420907: OR 1.22, 95% CI 1.07-1.39). Only rs9420907 remained significant after correction for multiple testing. Average genotypically-estimated LTL was significantly longer for those with meningioma compared to controls [mean cases: 560.2 bp (standard error (SE): 4.05 bp), mean controls: 541.5 bp (SE: 2.02 bp), logistic regression p value = 2.13 × 10-5]. CONCLUSION: Increased genotypically-estimated LTL was significantly associated with increased meningioma risk. A role for telomere length in the pathophysiology of meningioma is novel, and could lead to new insights on the etiology of meningioma.

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

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

Age-specific genome-wide association study in glioblastoma identifies increased proportion of 'lower grade glioma'-like features associated with younger age.

Glioblastoma (GBM) is the most common malignant brain tumor in the United States. Incidence of GBM increases with age, and younger age-at-diagnosis is significantly associated with improved prognosis. While the relationship between candidate GBM risk SNPs and age-at-diagnosis has been explored, genome-wide association studies (GWAS) have not previously been stratified by age. Potential age-specific genetic effects were assessed in autosomal SNPs for GBM patients using data from four previous GWAS. Using age distribution tertiles (18-53, 54-64, 65+) datasets were analyzed using age-stratified logistic regression to generate p values, odds ratios (OR), and 95% confidence intervals (95%CI), and then combined using meta-analysis. There were 4,512 total GBM cases, and 10,582 controls used for analysis. Significant associations were detected at two previously identified SNPs in 7p11.2 (rs723527 [p54-63 = 1.50x10-9 , OR54-63 = 1.28, 95%CI54-63 = 1.18-1.39; p64+ = 2.14x10-11 , OR64+ = 1.32, 95%CI64+ = 1.21-1.43] and rs11979158 [p54-63 = 6.13x10-8 , OR54-63 = 1.35, 95%CI54-63 = 1.21-1.50; p64+ = 2.18x10-10 , OR64+ = 1.42, 95%CI64+ = 1.27-1.58]) but only in persons >54. There was also a significant association at the previously identified lower grade glioma (LGG) risk locus at 8q24.21 (rs55705857) in persons ages 18-53 (p18-53 = 9.30 × 10-11 , OR18-53 = 1.76, 95%CI18-53 = 1.49-2.10). Within The Cancer Genome Atlas (TCGA) there was higher prevalence of 'LGG'-like tumor characteristics in GBM samples in those 18-53, with IDH1/2 mutation frequency of 15%, as compared to 2.1% [54-63] and 0.8% [64+] (p = 0.0005). Age-specific differences in cancer susceptibility can provide important clues to etiology. The association of a SNP known to confer risk for IDH1/2 mutant glioma and higher prevalence of IDH1/2 mutation within younger individuals 18-53 suggests that more younger individuals may present initially with 'secondary glioblastoma.'

Genetic determinants of childhood and adult height associated with osteosarcoma risk.

BACKGROUND: Although increased height has been associated with osteosarcoma risk in previous epidemiologic studies, to the authors' knowledge the relative contribution of stature during different developmental timepoints remains unclear. Furthermore, the question of how genetic determinants of height impact osteosarcoma etiology remains unexplored. Genetic variants associated with stature in previous genome-wide association studies may be biomarkers of osteosarcoma risk. METHODS: The authors tested the associations between osteosarcoma risk and polygenic scores for adult height (416 variants), childhood height (6 variants), and birth length (5 variants) in 864 osteosarcoma cases and 1879 controls of European ancestry. RESULTS: Each standard deviation increase in the polygenic score for adult height, corresponding to a 1.7-cm increase in stature, was found to be associated with a 1.10-fold increase in the risk of osteosarcoma (95% confidence interval [95% CI], 1.01-1.19; P =.027). Each standard deviation increase in the polygenic score for childhood height, corresponding to a 0.5-cm increase in stature, was associated with a 1.10-fold increase in the risk of osteosarcoma (95% CI, 1.01-1.20; P =.023). The polygenic score for birth length was not found to be associated with osteosarcoma risk (P =.11). When adult and childhood height scores were modeled together, they were found to be independently associated with osteosarcoma risk (P =.037 and P = .043, respectively). An expression quantitative trait locus for cartilage intermediate layer protein 2 (CILP2), rs8103992, was significantly associated with osteosarcoma risk after adjustment for multiple comparisons (odds ratio, 1.35; 95% CI, 1.16-1.56 [P = 7.93×10-5 and Padjusted =.034]). CONCLUSIONS: A genetic propensity for taller adult and childhood height attainments contributed independently to osteosarcoma risk in the current study data. These results suggest that the biological pathways affecting normal bone growth may be involved in osteosarcoma etiology.

Glioma Groups Based on 1p/19q, IDH, and TERT Promoter Mutations in Tumors.

BACKGROUND: The prediction of clinical behavior, response to therapy, and outcome of infiltrative glioma is challenging. On the basis of previous studies of tumor biology, we defined five glioma molecular groups with the use of three alterations: mutations in the TERT promoter, mutations in IDH, and codeletion of chromosome arms 1p and 19q (1p/19q codeletion). We tested the hypothesis that within groups based on these features, tumors would have similar clinical variables, acquired somatic alterations, and germline variants. METHODS: We scored tumors as negative or positive for each of these markers in 1087 gliomas and compared acquired alterations and patient characteristics among the five primary molecular groups. Using 11,590 controls, we assessed associations between these groups and known glioma germline variants. RESULTS: Among 615 grade II or III gliomas, 29% had all three alterations (i.e., were triple-positive), 5% had TERT and IDH mutations, 45% had only IDH mutations, 7% were triple-negative, and 10% had only TERT mutations; 5% had other combinations. Among 472 grade IV gliomas, less than 1% were triple-positive, 2% had TERT and IDH mutations, 7% had only IDH mutations, 17% were triple-negative, and 74% had only TERT mutations. The mean age at diagnosis was lowest (37 years) among patients who had gliomas with only IDH mutations and was highest (59 years) among patients who had gliomas with only TERT mutations. The molecular groups were independently associated with overall survival among patients with grade II or III gliomas but not among patients with grade IV gliomas. The molecular groups were associated with specific germline variants. CONCLUSIONS: Gliomas were classified into five principal groups on the basis of three tumor markers. The groups had different ages at onset, overall survival, and associations with germline variants, which implies that they are characterized by distinct mechanisms of pathogenesis. (Funded by the National Institutes of Health and others.).

Germline polymorphisms in myeloid-associated genes are not associated with survival in glioma patients.

Immune cells of myeloid origin, including microglia, macrophages, and myeloid-derived suppressor cells adopt immunosuppressive phenotypes that support gliomagenesis. Here, we tested an a priori hypothesis that single nucleotide polymorphisms (SNPs) in genes related to glioma-associated myeloid cell regulation and function are also associated with patient survival after glioma diagnosis. Subjects for this study were 992 glioma patients treated at The University of Texas MD Anderson Cancer Center in Houston, Texas between 1992 and 2008. Haplotype-tagging SNPs in 91 myeloid-associated genes were analyzed for association with survival by Cox regression. Individual SNP- and gene-based tests were performed separately in glioblastoma (WHO grade IV, n = 511) and lower-grade glioma (WHO grade II-III, n = 481) groups. After adjustment for multiple testing, no myeloid-associated gene variants were significantly associated with survival in glioblastoma. Two SNPs, rs147960238 in CD163 (p = 2.2 × 10-5) and rs17138945 in MET (p = 5.6 × 10-5) were significantly associated with survival of patients with lower-grade glioma. However, these associations were not confirmed in an independent analysis of 563 lower-grade glioma cases from the University of California at San Francisco Adult Glioma Study (p = 0.65 and p = 0.41, respectively). The results of this study do not support a role for inherited polymorphisms in myeloid-associated genes in affecting survival of patients diagnosed with glioblastoma or lower-grade glioma.

Genomic characterization of chronic lymphocytic leukemia (CLL) in radiation-exposed Chornobyl cleanup workers.

BACKGROUND: Chronic lymphocytic leukemia (CLL) was the predominant leukemia in a recent study of Chornobyl cleanup workers from Ukraine exposed to radiation (UR-CLL). Radiation risks of CLL significantly increased with increasing bone marrow radiation doses. Current analysis aimed to clarify whether the increased risks were due to radiation or to genetic mutations in the Ukrainian population. METHODS: A detailed characterization of the genomic landscape was performed in a unique sample of 16 UR-CLL patients and age- and sex-matched unexposed general population Ukrainian-CLL (UN-CLL) and Western-CLL (W-CLL) patients (n = 28 and 100, respectively). RESULTS: Mutations in telomere-maintenance pathway genes POT1 and ATM were more frequent in UR-CLL compared to UN-CLL and W-CLL (both p < 0.05). No significant enrichment in copy-number abnormalities at del13q14, del11q, del17p or trisomy12 was identified in UR-CLL compared to other groups. Type of work performed in the Chornobyl zone, age at exposure and at diagnosis, calendar time, and Rai stage were significant predictors of total genetic lesions (all p < 0.05). Tumor telomere length was significantly longer in UR-CLL than in UN-CLL (p = 0.009) and was associated with the POT1 mutation and survival. CONCLUSIONS: No significant enrichment in copy-number abnormalities at CLL-associated genes was identified in UR-CLL compared to other groups. The novel associations between radiation exposure, telomere maintenance and CLL prognosis identified in this unique case series provide suggestive, though limited data and merit further investigation.

Adult infiltrating gliomas with WHO 2016 integrated diagnosis: additional prognostic roles of ATRX and TERT.

The "integrated diagnosis" for infiltrating gliomas in the 2016 revised World Health Organization (WHO) classification of tumors of the central nervous system requires assessment of the tumor for IDH mutations and 1p/19q codeletion. Since TERT promoter mutations and ATRX alterations have been shown to be associated with prognosis, we analyzed whether these tumor markers provide additional prognostic information within each of the five WHO 2016 categories. We used data for 1206 patients from the UCSF Adult Glioma Study, the Mayo Clinic and The Cancer Genome Atlas (TCGA) with infiltrative glioma, grades II-IV for whom tumor status for IDH, 1p/19q codeletion, ATRX, and TERT had been determined. All cases were assigned to one of 5 groups following the WHO 2016 diagnostic criteria based on their morphologic features, and IDH and 1p/19q codeletion status. These groups are: (1) Oligodendroglioma, IDH-mutant and 1p/19q-codeleted; (2) Astrocytoma, IDH-mutant; (3) Glioblastoma, IDH-mutant; (4) Glioblastoma, IDH-wildtype; and (5) Astrocytoma, IDH-wildtype. Within each group, we used univariate and multivariate Cox proportional hazards models to assess associations of overall survival with patient age at diagnosis, grade, and ATRX alteration status and/or TERT promoter mutation status. Among Group 1 IDH-mutant 1p/19q-codeleted oligodendrogliomas, the TERT-WT group had significantly worse overall survival than the TERT-MUT group (HR: 2.72, 95% CI 1.05-7.04, p = 0.04). In both Group 2, IDH-mutant astrocytomas and Group 3, IDH-mutant glioblastomas, neither TERT mutations nor ATRX alterations were significantly associated with survival. Among Group 4, IDH-wildtype glioblastomas, ATRX alterations were associated with favorable outcomes (HR: 0.36, 95% CI 0.17-0.81, p = 0.01). Among Group 5, IDH-wildtype astrocytomas, the TERT-WT group had significantly better overall survival than the TERT-MUT group (HR: 0.48, 95% CI 0.27-0.87), p = 0.02). Thus, we present evidence that in certain WHO 2016 diagnostic groups, testing for TERT promoter mutations or ATRX alterations may provide additional useful prognostic information.

Association of common genetic variation in the protein C pathway genes with clinical outcomes in acute respiratory distress syndrome.

BACKGROUND: Altered plasma levels of protein C, thrombomodulin, and the endothelial protein C receptor are associated with poor clinical outcomes in patients with acute respiratory distress syndrome (ARDS). We hypothesized that common variants in these genes would be associated with mortality as well as ventilator-free and organ failure-free days in patients with ARDS. METHODS: We genotyped linkage disequilibrium-based tag single-nucleotide polymorphisms in the ProteinC, Thrombomodulin and Endothelial Protein C Reptor Genes among 320 self-identified white patients of European ancestry from the ARDS Network Fluid and Catheter Treatment Trial. We then tested their association with mortality as well as ventilator-free and organ-failure free days. RESULTS: The GG genotype of rs1042580 (p = 0.02) and CC genotype of rs3716123 (p = 0.002), both in the thrombomodulin gene, and GC/CC genotypes of rs9574 (p = 0.04) in the endothelial protein C receptor gene were independently associated with increased mortality. An additive effect on mortality (p < 0.001), ventilator-free days (p = 0.01), and organ failure-free days was observed with combinations of these high-risk genotypes. This association was independent of age, severity of illness, presence or absence of sepsis, and treatment allocation. CONCLUSIONS: Genetic variants in thrombomodulin and endothelial protein C receptor genes are additively associated with mortality in ARDS. These findings suggest that genetic differences may be at least partially responsible for the observed associations between dysregulated coagulation and poor outcomes in ARDS.

Analysis of 60 reported glioma risk SNPs replicates published GWAS findings but fails to replicate associations from published candidate-gene studies.

Genomewide association studies (GWAS) and candidate-gene studies have implicated single-nucleotide polymorphisms (SNPs) in at least 45 different genes as putative glioma risk factors. Attempts to validate these associations have yielded variable results and few genetic risk factors have been consistently replicated. We conducted a case-control study of Caucasian glioma cases and controls from the University of California San Francisco (810 cases, 512 controls) and the Mayo Clinic (852 cases, 789 controls) in an attempt to replicate previously reported genetic risk factors for glioma. Sixty SNPs selected from the literature (eight from GWAS and 52 from candidate-gene studies) were successfully genotyped on an Illumina custom genotyping panel. Eight SNPs in/near seven different genes (TERT, EGFR, CCDC26, CDKN2A, PHLDB1, RTEL1, TP53) were significantly associated with glioma risk in the combined dataset (P < 0.05), with all associations in the same direction as in previous reports. Several SNP associations showed considerable differences across histologic subtype. All eight successfully replicated associations were first identified by GWAS, although none of the putative risk SNPs from candidate-gene studies was associated in the full case-control sample (all P values > 0.05). Although several confirmed associations are located near genes long known to be involved in gliomagenesis (e.g., EGFR, CDKN2A, TP53), these associations were first discovered by the GWAS approach and are in noncoding regions. These results highlight that the deficiencies of the candidate-gene approach lay in selecting both appropriate genes and relevant SNPs within these genes.

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 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.

Smoking and genetic risk variation across populations of European, Asian, and African American ancestry--a meta-analysis of chromosome 15q25.

Recent meta-analyses of European ancestry subjects show strong evidence for association between smoking quantity and multiple genetic variants on chromosome 15q25. This meta-analysis extends the examination of association between distinct genes in the CHRNA5-CHRNA3-CHRNB4 region and smoking quantity to Asian and African American populations to confirm and refine specific reported associations. Association results for a dichotomized cigarettes smoked per day phenotype in 27 datasets (European ancestry (N = 14,786), Asian (N = 6,889), and African American (N = 10,912) for a total of 32,587 smokers) were meta-analyzed by population and results were compared across all three populations. We demonstrate association between smoking quantity and markers in the chromosome 15q25 region across all three populations, and narrow the region of association. Of the variants tested, only rs16969968 is associated with smoking (P < 0.01) in each of these three populations (odds ratio [OR] = 1.33, 95% CI = 1.25-1.42, P = 1.1 × 10(-17) in meta-analysis across all population samples). Additional variants displayed a consistent signal in both European ancestry and Asian datasets, but not in African Americans. The observed consistent association of rs16969968 with heavy smoking across multiple populations, combined with its known biological significance, suggests rs16969968 is most likely a functional variant that alters risk for heavy smoking. We interpret additional association results that differ across populations as providing evidence for additional functional variants, but we are unable to further localize the source of this association. Using the cross-population study paradigm provides valuable insights to narrow regions of interest and inform future biological experiments.