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.

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.

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.

Genetic variants in ARID5B and CEBPE are childhood ALL susceptibility loci in Hispanics.

Recent genome-wide studies conducted in European Whites have identified novel susceptibility genes for childhood acute lymphoblastic leukemia (ALL). We sought to examine whether these loci are susceptibility genes among Hispanics, whose reported incidence of childhood ALL is the highest of all ethnic groups in California, and whether their effects differ between Hispanics and non-Hispanic Whites (NHWs). We genotyped 13 variants in these genes among 706 Hispanic (300 cases, 406 controls) and 594 NHW (225 cases, 369 controls) participants in a matched population-based case-control study in California. We found significant associations for the five studied ARID5B variants in both Hispanics (p values of 1.0 × 10(-9) to 0.004) and NHWs (p values of 2.2 × 10(-6) to 0.018). Risk estimates were in the same direction in both groups (ORs of 1.53-1.99 and 1.37-1.84, respectively) and strengthened when restricted to B-cell precursor high-hyperdiploid ALL (>50 chromosomes; ORs of 2.21-3.22 and 1.67-2.71, respectively). Similar results were observed for the single CEBPE variant. Hispanics and NHWs exhibited different susceptibility loci at CDKN2A. Although IKZF1 loci showed significant susceptibility effects among NHWs (p < 1 × 10(-5)), their effects among Hispanics were in the same direction but nonsignificant, despite similar minor allele frequencies. Future studies should examine whether the observed effects vary by environmental, immunological, or lifestyle factors.

Socioeconomic status and lung cancer: unraveling the contribution of genetic admixture.

OBJECTIVES: We examined the relationship between genetic ancestry, socioeconomic status (SES), and lung cancer among African Americans and Latinos. METHODS: We evaluated SES and genetic ancestry in a Northern California lung cancer case-control study (1998-2003) of African Americans and Latinos. Lung cancer case and control participants were frequency matched on age, gender, and race/ethnicity. We assessed case-control differences in individual admixture proportions using the 2-sample t test and analysis of covariance. Logistic regression models examined associations among genetic ancestry, socioeconomic characteristics, and lung cancer. RESULTS: Decreased Amerindian ancestry was associated with higher education among Latino control participants and greater African ancestry was associated with decreased education among African lung cancer case participants. Education was associated with lung cancer among both Latinos and African Americans, independent of smoking, ancestry, age, and gender. Genetic ancestry was not associated with lung cancer among African Americans. CONCLUSIONS: Findings suggest that socioeconomic factors may have a greater impact than genetic ancestry on lung cancer among African Americans. The genetic heterogeneity and recent dynamic migration and acculturation of Latinos complicate recruitment; thus, epidemiological analyses and findings should be interpreted cautiously.

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.

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.

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.