Skip to main content
  • Research article
  • Open access
  • Published:

Impact of atopy on risk of glioma: a Mendelian randomisation study

Abstract

Background

An inverse relationship between allergies with glioma risk has been reported in several but not all epidemiological observational studies. We performed an analysis of genetic variants associated with atopy to assess the relationship with glioma risk using Mendelian randomisation (MR), an approach unaffected by biases from temporal variability and reverse causation that might have affected earlier investigations.

Methods

Two-sample MR was undertaken using genome-wide association study data. We used single nucleotide polymorphisms (SNPs) associated with atopic dermatitis, asthma and hay fever, IgE levels, and self-reported allergy as instrumental variables. We calculated MR estimates for the odds ratio (OR) for each risk factor with glioma using SNP-glioma estimates from 12,488 cases and 18,169 controls, using inverse-variance weighting (IVW), maximum likelihood estimation (MLE), weighted median estimate (WME) and mode-based estimate (MBE) methods. Violation of MR assumptions due to directional pleiotropy were sought using MR-Egger regression and HEIDI-outlier analysis.

Results

Under IVW, MLE, WME and MBE methods, associations between glioma risk with asthma and hay fever, self-reported allergy and IgE levels were non-significant. An inverse relationship between atopic dermatitis and glioma risk was found by IVW (OR 0.96, 95% confidence interval (CI) 0.93–1.00, P = 0.041) and MLE (OR 0.96, 95% CI 0.94–0.99, P = 0.003), but not by WME (OR 0.96, 95% CI 0.91–1.01, P = 0.114) or MBE (OR 0.97, 95% CI 0.92–1.02, P = 0.194).

Conclusions

Our investigation does not provide strong evidence for relationship between atopy and the risk of developing glioma, but findings do not preclude a small effect in relation to atopic dermatitis. Our analysis also serves to illustrate the value of using several MR methods to derive robust conclusions.

Peer Review reports

Background

Although glioma accounts for approximately 80% of malignant primary brain tumours [1], to date, few aetiological risk factors are well established for the disease [2]. Over the past three decades the search for an immune-mediated risk factor that might influence risk has led to studies of a possible relationship between multiple allergic conditions and autoimmune disorders with glioma [3].

Several case-control studies have shown that self-reported allergic conditions may protect against glioma [4]. For example, in the International Adult Brain Tumour Study, based on 1178 glioma patients, an odds ratio (OR) of 0.59 was found for any self-reported allergy [5]. Other case-control studies have reported similar ORs, however, most have been reliant on substantial numbers of proxy informants (up to 44%) [4, 6] and have potential bias as a consequence of how controls were ascertained, thereby casting doubt on findings. In contrast to case-control studies, evidence for an association between glioma and allergy from cohort-based analyses has been less forthcoming [7], although such studies have been poorly powered to demonstrate a relationship.

Assaying IgE potentially reduces bias stemming from self-reporting despite levels not necessarily corresponding to specific allergies or equating to a single allergic response. Nevertheless, measurement of IgE has been explored by a number of researchers seeking to identify risk factors for glioma [8,9,10]. In a case-control study of 228 cases and 289 controls performed in 2004 [8], self-reported allergies and IgE levels were both inversely associated with glioma, but concordance between the two outcomes was poor. In a larger study of 535 cases and 532 controls [11], both self-reported allergies and IgE levels were inversely related to glioma risk; however, IgE levels in patients were affected by temozolomide treatment. A case-control study nested within the European Prospective Investigation into Cancer and Nutrition cohort based on prospectively collected serum IgE levels reported a non-significant OR of 0.73 [9]. A similar nested case-control study performed in the USA based on 181 cases reported a non-significant OR of 0.72 for high serum IgE [10].

Several mechanisms have been proposed to explain a possible association between atopic disease and glioma [12]. The findings could reflect a true causal effect of the heightened immune function reported for atopy on tumour development. Alternatively, the associations observed might be non-causal, arising as a consequence of methodological biases inherent in the study design. Imprecisely defined exposures, such as allergic disease, are likely to have affected the validity of the findings of both case-control and cohort studies. The heterogeneous description of allergy in studies and different levels of detail in self-reporting on individual allergies complicate the interpretation of results. Additional biases include possible selection bias in controls, recall bias from self-reported allergy assessment and reverse causation or confounding from unmeasured effects. Finally, the high frequency of exposure ascertainment by proxy for cases is also likely to have systematically biased findings.

Mendelian randomisation (MR) analysis can be used to minimise potential biases in conventional observational studies and to determine the causal association of an exposure with an outcome such as disease risk [13]. The causal association can also be manifested by common genetic and biological pathways that determine two sequentially developed phenotypes such as an atopic trait and glioma risk. Atopy has a strong heritable basis [14, 15] and, thus far, genome-wide association studies (GWAS) have identified over 50 loci associated with different atopy-related traits [16]. The alleles associated with atopy should be randomly assigned to offspring from parents during mitosis, a process analogous to the random assignment of subjects to an exposure of interest in randomised clinical trials. Thus, genetic scores summarising the effects of single nucleotide polymorphisms (SNPs) associated with atopy-related traits can serve as instrumental variables (IVs) in a MR analysis of atopy and glioma risk.

To examine the nature of the association between atopy and glioma, we implemented two-sample MR [17] to estimate associations between atopy-associated SNPs and glioma risk using summary data from the recent GWAS meta-analysis performed by the Glioma International Case-Control Consortium study [18].

Methods

Two-sample MR was undertaken using GWAS data. Ethical approval was not sought for this specific project because all data came from the summary statistics of published GWAS, and no individual-level data were used.

Glioma genotyping data

Glioma genotyping data were derived from the most recent meta-analysis of GWAS in glioma, which related > 10 million genetic variants (after imputation) to glioma, in 12,488 glioma patients and 18,169 controls from eight independent GWAS datasets of individuals of European descent [18] (Additional file 1: Table S1). Comprehensive details of the genotyping and quality control of the seven GWAS have been previously reported [18].

Genetic variant instruments for atopic traits

SNPs associated with each of the atopy-related traits investigated, namely atopic dermatitis (eczema), asthma and hay fever, IgE level, and self-reported allergy, by the NHGRI-EBI GWAS Catalog [19,20,21,22,23,24,25,26] at genome-wide significance (i.e. P ≤ 5.0 × 10− 8) in individuals with European ancestry were used as IVs. To avoid co-linearity between SNPs for each trait, we excluded SNPs that were correlated (i.e. r2 value of ≥ 0.001) within each trait, and only considered the SNPs with the strongest effect on the trait for use as IVs (Additional file 2: Table S2). For each SNP, we recovered the chromosome position, risk allele, association estimates (per-allele log-OR) and standard errors (Table 1). The allele that was associated with increased risk of the exposure was considered the effect allele. For IgE level, the allele associated with an increase in serum IgE was considered the effect allele. Allele frequencies for these SNPs were compared between the atopy-related trait and glioma datasets to ensure that the effect estimates were recorded with respect to the same allele. Gliomas are heterogeneous and different tumour subtypes, defined in part by malignancy grade (e.g. pilocytic astrocytoma World Health Organization (WHO) grade I, diffuse ‘low-grade’ glioma WHO grade II, anaplastic glioma WHO grade III and glioblastoma (GBM) WHO grade IV) can be distinguished [27]. For the sake of brevity we considered gliomas as being either GBM or non-GBM.

Table 1 Variant and effect allele with frequencies and magnitude of effect on each atopy-related trait and strength of association with glioma

Two-sample MR method

The association between each atopy-related trait and glioma was examined using MR on summary statistics using the inverse-variance weighting (IVW) method and maximum likelihood estimation (MLE) as per Burgess et al. [28]. The IVW ratio estimate \( \left(\widehat{\upbeta}\right) \) of all SNPs associated with each atopy-related trait on glioma risk was calculated as follows:

$$ \widehat{\upbeta}=\frac{\sum_k{X}_k{Y}_k{\sigma_Y}_k^{-2}}{\sum_k{X}_k^2{\sigma_Y}_k^{-2}} $$

Where X k corresponds to the association of SNP k (as log of the OR per risk allele) with the atopy-related trait Y k is the association between SNP k and glioma risk (as log OR) with standard error \( {\sigma}_{Y_k} \). The estimate for \( \left(\widehat{\upbeta}\right) \) represents the causal increase in the log odds of glioma for each trait. The standard error of the combined ratio estimate is given by:

$$ se\left(\widehat{\upbeta}\right)=\sqrt{\frac{1}{\sum_k{X}_k^2{\sigma_Y}_k^{-2}}} $$

For the MLE, a bivariate normal distribution for the genetic associations was assumed, and the R function optim was used to estimate β. \( se\left(\widehat{\upbeta}\right) \) was calculated using observed information. The correlation between the errors of Y k and X k was taken to be 0 as they were derived from independent studies.

A central tenet in MR is the absence of pleiotropy (i.e. a gene influencing multiple traits) between the SNPs influencing the exposure and outcome disease risk [13]. This would be revealed as deviation from a linear relationship between SNPs and their effect size for atopy and glioma risk. To examine for violation of the standard IV assumptions in our analysis we first performed MR-Egger regression, as well as HEIDI-outlier analysis, as per Zhu et al. [29], imposing the advocated threshold of P ≤ 0.01. Additionally, we derived weighted median estimates (WME) [30] and mode-based estimates (MBE) [31] to establish the robustness of findings.

Atopic dermatitis, asthma and hay fever, and self-reported allergy as well as all of the disease outcomes (all glioma, GBM and non-GBM glioma) are binary. The causal effect estimates therefore represent the odds for outcome disease risk per unit increase in the log OR of the exposure disease [32]. These ORs were converted to represent the OR for the outcome disease per doubling in odds of the exposure disease to aid interpretation [32].

For each statistical test we considered a global significance level of P < 0.05 as being satisfactory to derive conclusions. To assess the robustness of our conclusions, we initially imposed a conservative Bonferroni-corrected significance threshold of 0.0125 (i.e. 0.05/4 atopy-related traits). We considered a P value ≥ 0.05 as non-significant (i.e. no association), a P < 0.05 as evidence for a potential causal association, and a P < 0.0125 as significant evidence for an association. All statistical analyses were undertaken using R software (Version 3.1.2). The meta and gsmr packages were used to generate forest plots and perform HEIDI-outlier analysis [29].

The power of a MR investigation depends greatly on the proportion of variance in the risk factor that is explained by the IV. We estimated study power a priori using the methodology of Burgess et al. [33], making use of published estimates of the heritability of trait associated IV SNPs [34,35,36], as well as estimates found by direct calculation (Additional file 3: Table S3), and the reported effect of each trait on glioma risk reported in a meta-analysis of epidemiological studies [18]. Additional file 4: Table S4 shows the range of ORs for which we had less than 80% power to detect for each of the four atopy-related traits.

Simulation model

Through simulation we evaluated the suitability of using each employed MR method in a two-sample setting with binary-exposure and binary-outcome data. Let i index genetic variants, N be the total number of genetic variants, and j index individuals. Genetic variants g ij were generated independently by sampling from a Binomial(2,p j ) distribution with probability p j drawn from a Uniform(0.1,0.9) distribution, to mimic bi-allelic SNPs in Hardy–Weinberg equilibrium. Let w j correspond to the per-allele OR for the exposure disease, sampled from ORs reported for genome-wide significant SNPs reported in the GWAS Catalog [37], and v be the OR for the outcome disease per doubling in odds of the exposure disease. For each individual, exposure disease odds x j , outcome disease odds y j , exposure disease status a j , and outcome disease status b j were determined as follows:

$$ {x}_j={x}_0\prod \limits_{i=1}^N{w_i}^{g_{ij}} $$
$$ {y}_j={y}_0\times {2}^{{\mathsf{\log}}_2{x}_j\times {\mathsf{\log}}_2v} $$
$$ {a}_j\sim \mathsf{Binomial}\left(1,\frac{x_j}{1+{x}_j}\right) $$
$$ {b}_j\sim \mathsf{Binomial}\left(1,\frac{y_j}{1+{y}_j}\right) $$

Data for 1,000,000 individuals were simulated and partitioned at random to reflect the two-sample setting. Cases and controls for the exposure and outcome GWAS were sampled from each half of the dataset using the exposure and outcome disease statuses of each individual, and association statistics computed under an additive logistic regression model. To ensure the simulated data closely resembled the atopy-related trait and glioma data, the simulation analysis was repeated for each binary atopy-related trait using the same number of genetic variants as IVs and the same numbers of case and control individuals as used to estimate the atopy-related trait and glioma association statistics (Additional file 5: Table S5). Parameters x 0 = 0.0005 and y 0 = 0.01 were chosen to ensure the prevalence of the simulated exposure and outcome diseases were similar to that of the atopy-related traits and glioma, respectively (Additional file 5: Table S5). To determine the suitability of each MR method we considered two scenarios: (1) no causal relationship between exposure and outcome (v = 1.00) and (2) a causal relationship between exposure and outcome (v = 1.33). We performed 100 simulations for each scenario for each binary atopy-related trait.

Results

The atopic dermatitis risk SNP rs909341, which is highly correlated with the chromosome 20q13.33 glioma risk SNP rs2297440 (D’ = 0.89, r2 = 0.77), was strongly associated with risk of glioma (P = 2.10 × 10−34). Testing for pleiotropy using HEIDI-outlier analysis formally identified rs909341 as violating the assumption of the instrument on the outcome. Henceforth, we confined our analysis of the relationship between atopic dermatitis and glioma to a dataset excluding this SNP.

Figure 1 shows forest plots of ORs for glioma generated from the SNPs. There was minimal evidence of heterogeneity between variants for asthma and hay fever, atopic dermatitis, IgE levels and self-reported allergy (respective I2 and Phet values being 28% and 0.192, 8% and 0.377, 0% and 0.444, and 0% and 0.707). Including rs909341 in the analysis for atopic dermatitis, the I2 value was 90% and Phet < 10− 4 (Additional file 6: Figure S1), providing further evidence that inclusion of this SNP would invalidate the MR analysis.

Fig. 1
figure 1

Forest plot of Wald odds ratios (ORs) and 95% confidence intervals generated from single nucleotide polymorphisms (SNPs) associated with atopy-related traits. ORs for individual SNPs are listed according to magnitude of effect in the instrumental variable analysis and are presented with pooled effects using the inverse-variance weighting method. Squares represent the point estimate, and the bars are the 95% confidence intervals. a Asthma and hay fever, b atopic dermatitis, c IgE level, d self-reported allergy

The results of the IVW, MLE, WME, MBE and MR-Egger methods are summarised in Table 2. Using the IVW method to pool results from individual SNPs, no associations (i.e. P ≥ 0.05) were identified between genetically conferred risk of raised IgE level (OR 0.88, 95% CI 0.69–1.13, P = 0.319), asthma and hay fever (OR 0.96, 95% CI 0.90–1.03, P = 0.248), or self-reported allergy (OR 1.03, 95% CI 0.95–1.11, P = 0.534) with risk of all glioma. There was some support for an inverse relationship between atopic dermatitis and glioma risk (OR 0.96, 95% CI 0.93–1.00, P = 0.041), albeit not significant after adjustment for multiple testing.

Table 2 Inverse-variance weighting, maximum likelihood estimation, weighted median estimate, mode-based estimate and Mendelian randomisation-Egger test results for combined atopy-related instrumental variables

Using MLE, no associations were identified between asthma and hay fever (OR 0.96, 95% CI 0.93–1.00, P = 0.066), IgE levels (OR 0.88, 95% CI 0.74–1.05, P = 0.157) or self-reported allergy (OR 1.02, 95% CI 0.97–1.08, P = 0.429) with risk of all glioma. For atopic dermatitis, an OR of 0.96 (95% CI 0.94–0.99, P = 0.003) was shown, which remained significant after adjusting for multiple testing. Figure 2 shows relaxation of the assumption that the correlation between the errors in X k and Y k is zero for each of the atopy-related traits demonstrating the consistency of findings. Specifically, for a correlation in the range −0.15 to 0.15, the association between atopic dermatitis and glioma risk remained significant.

Fig. 2
figure 2

Plot of P value of maximum likelihood estimation associations with glioma against correlation between errors in X k and Y k . a Asthma and hay fever, b atopic dermatitis, c IgE level, d self-reported allergy

In contrast to findings from IVW and MLE, no significant support was provided by either the WME or MBE for an association between any of the atopy-related traits and glioma risk, including atopic dermatitis (WME: OR 0.96, 95% CI 0.91–1.01, P = 0.114; MBE: OR 0.97, 95% CI 0.92–1.02, P = 0.194; Table 2).

The respective effect estimated from MR-Egger regression (Fig. 3) were 0.97 for atopic dermatitis (95% CI 0.92–1.03; P = 0.375), 0.63 for IgE levels (95% CI 0.32–1.25; P = 0.184), 0.99 for asthma and hay fever (95% CI 0.72–1.36, P = 0.951) and 0.92 for self-reported allergy (95% CI 0.69–1.22; P = 0.540), with intercepts of −0.004 (95% CI −0.014 to 0.006, P = 0.396), 0.027 (95% CI 0.001 to 0.053, P = 0.042), −0.007 (95% CI −0.030 to 0.016, P = 0.542) and 0.017 (95% CI 0.003–0.031, P = 0.018). Collectively, these findings provide possible evidence of systematic bias in the IVW estimate for IgE level and self-reported allergy, which might have arisen through overall unbalanced horizontal pleiotropy. There was no such evidence for such pleiotropy in respect of atopic dermatitis.

Fig. 3
figure 3

Scatter plots of genetic associations with glioma against genetic associations with the exposure. a Asthma and hay fever, b atopic dermatitis, c IgE level, d self-reported allergy

We explored the possibility that a relationship between atopy and glioma might be subtype specific, considering GBM and non-GBM separately. Imposing a stronger significance threshold of P = 0.00625 (0.05/8, to correct for testing four traits over two outcomes), no histology-specific associations were shown by the IVW method between asthma and hay fever, IgE levels and self-reported allergy and glioma risk, with the respective ORs for the IVW method being 0.97, 0.92 and 1.04 for GBM tumours, and 0.96, 0.97 and 1.04 for non-GBM tumours (Additional file 7: Table S6). For atopic dermatitis, a significant OR of 0.94 (95% CI 0.90–0.98, P = 0.004) was shown for GBM but not for non-GBM (OR 0.98, 95% CI 0.93–1.03, P = 0.421). The association between atopic dermatitis and risk of GBM was also apparent in the MLE analysis, which provided an OR of 0.94 (95% CI 0.91–0.97, P = 2.17 × 10− 4). MR-Egger regression provided for an intercept of −0.007 (95% CI −0.019 to 0.005, P = 0.247). As with the analysis of all glioma, the association between atopic dermatitis and GBM was weaker under the WME (OR 0.96, 95% CI 0.91–1.02, P = 0.172) and MBE (OR 0.95, 95% CI 0.90–1.01, P = 0.096) frameworks.

Although previously implemented in other studies [32, 38], ratio estimators may not fully recapitulate an estimate of the causal OR in the case of binary exposures, such as atopic dermatitis, and binary outcomes such as glioma [39]. We therefore evaluated, through simulation, whether the IVW, MLE, WME, MBE and MR-Egger methods provide reliable estimates of causal ORs. When no causal relationship between exposure and outcome was simulated, each MR method provided accurate estimates of the null relationship (Additional file 5: Table S5). Conversely, when a causal relationship was simulated, the magnitudes of the relationship estimates were weakly inflated in some instances (Additional file 5: Table S5), indicating the importance of considering additional evidence when evaluating causal relationships between binary exposures and binary outcomes.

Discussion

To our knowledge, this is the first MR study evaluating a range of atopy-related traits with glioma risk. Overall, our results provide evidence for a causal protective effect of atopic dermatitis with GBM tumours, but do not provide evidence that asthma and hay fever, raised IgE levels, or self-reported allergy is protective against the risk of developing glioma.

Possible mechanisms explaining an observed inverse relation between the risk of atopic dermatitis and the risk of glioma have been suggested in previous papers [12], postulated to be the consequence of immune system hyperactivity. The question thus arises as to how such divergent findings for other atopic traits can be explained or reconciled, when they have been previously reported in high numbers.

A key assumption in MR is that the instrument affects glioma risk through its effect on a specific phenotype/exposure (i.e. atopic traits), and does not have a direct effect on glioma risk. We tested this assumption using MR-Egger regression and HEIDI-outlier analysis and found possible evidence of violation of this assumption for IgE and self-reported allergy. It is notable that self-reported allergy does not show an approximately quadratic response to correlation, in contrast to asthma and hay fever, atopic dermatitis and IgE level. This is likely to be a consequence of imprecise estimates of the association between SNPs and allergy, illustrating the inherent issue in attempting to make use of self-reported allergy data as an atopy-related trait.

The meta-analyses of published epidemiological observational studies has indeed provided strong evidence for an inverse relationship between atopy and glioma risk [40]. However, most of the support for such a relationship came from case-control studies [4]. A common limitation in retrospective studies of glioma has been the use of proxy respondents for patients with cognitive impairment, who may not remember past exposures accurately due to cognitive deficits [4]. Such issues are compounded by the fact that, across studies, multiple atopic traits have been assessed. The strength of support for a relationship seen across case-control studies contrasts markedly with the limited evidence for a relationship from prospective cohort-based analyses [7].

By inference, a relationship between long-term antihistamine use could theoretically provide supporting evidence, albeit indirect, that atopic-mediated mechanisms influence glioma risk. However, the impact of antihistamine use is difficult to disentangle from that of allergies, as these factors are highly correlated and few individuals without allergies use antihistamines regularly. Paradoxically, an increased risk for glioma associated with antihistamines, particularly among individuals with allergic conditions, has been found in some studies [41, 42].

Raised IgE levels and self-reported allergy suffer limitations as traits used to assess the effect of atopy on glioma risk as they are both variable over short time scales in their level of expression (in contrast to clinical diagnosis of atopic dermatitis). Further, allergies may develop later in life, and patients may not necessarily exhibit symptoms. This introduces the possibility of bias and error due to the time varying association of SNPs with the exposure. However, it has been suggested that seasonality does not have a significant effect [11].

An additional possible explanation for the lack of causal association between IgE levels and glioma risk seen in this study is that the causality is in fact reversed, which could result in epidemiological observational studies reporting inverse relationships [8, 9], but would not affect an MR analysis. Immunosuppression caused by glioblastoma is well documented [43, 44] and may lead to reduced expression of atopy. Furthermore, in addition to steroids, temozolomide therapy, routinely used to treat GBM nowadays, leads to reduced blood IgE levels [11].

Using data from large genetic consortia for multiple atopy-related traits and glioma risk has enabled us to more precisely test our study hypotheses than if we had used individual-level data from a smaller study. Through simulation scenarios, the IVW, MLE, WME, MBE and MR-Egger methods have been demonstrated to accurately estimate causal effects using summary-level data [28, 30, 31, 45]. However, using summary-level data instead of individual-level data limits the approaches that can be used to test the validity of genetic variants as IVs, as adjusting for measured covariates and assessing gene-environment interactions is generally not possible using summary-level data [46]. The first-stage F statistic was large (> 25 for all traits), and therefore weak instrument bias is unlikely.

Epidemiological observational studies have reported inverse relationships between atopy-related traits and glioma risk, with ORs in the range 0.43–0.96 for asthma [6, 47], 0.42–0.90 for atopic dermatitis [6, 47], 0.37–0.73 for IgE levels [8,9,10] and 0.47–0.69 for self-reported allergies [4, 5, 8]. Odds ratios for binary exposures estimated in this MR study represent the OR for the outcome disease per doubling in odds of the exposure disease, and the magnitudes of these causal effect estimates are therefore not directly comparable to those reported in observational studies.

Our MR analysis has several strengths. Firstly, by utilising the random allocation of genetic variants, we were able to overcome potential confounding and reverse causation that may bias estimates from observational studies. Secondly, given that a poor outcome from glioma is almost universal, it is unlikely that survival bias will have influenced study findings. Lastly, the findings from this study represent the association of a lifelong atopy with glioma in the general European population.

Nevertheless, our study does have limitations. Firstly, while it is entirely appropriate to implement different MR methods to assess the robustness of findings, they have a differing power to demonstrate associations, with the WME, MBE and MR-Egger methods having less power than IVW and MLE. Irrespective of such factors, our study only had 80% power to detect ORs of 1.16, 1.09, 1.16 and 1.22 for asthma and hay fever, atopic dermatitis, IgE level and self-reported allergy, respectively (Additional file 4: Table S4), due to the very low proportion of variability in the atopy-related traits explained by the SNPs used. Hence, we cannot exclude the possibility that these traits influence glioma risk, albeit modestly. To explore this possibility, will require additional IVs and larger sample sizes affording increased power. Furthermore, it is possible that an effect of atopy on glioma risk might be mediated through mechanisms associated with a trait that we have not captured by using MR to assess asthma and hay fever and self-reported allergy. Secondly, a weakness of the two-sample MR strategy is that it does not allow examination of non-linear relationships between exposures and outcomes. Finally, we have sought to examine whether bias could be introduced when considering a binary exposure for a binary outcome. Although in our simulation study we found no evidence of bias when estimating non-causal relationships, we did not extend our analysis to consider the potential impact of invalid SNPs.

Conclusions

In conclusion, our investigation does not provide strong evidence for a relationship between atopy-related diseases and risk of developing glioma, but findings do not preclude a small effect for atopic dermatitis. Our analysis also serves to illustrate the value of using several MR methods to derive robust conclusions.

Abbreviations

CI:

confidence interval

GBM:

glioblastoma

GWAS:

genome-wide association study

IV:

instrumental variable

IVW:

inverse-variance weighting

MBE:

mode-based estimate

MLE:

maximum likelihood estimation

MR:

Mendelian randomisation

OR:

odds ratio

SNP:

single nucleotide polymorphism

WHO:

World Health Organization

WME:

weighted median estimate

References

  1. Ostrom QT, Gittleman H, Xu J, Kromer C, Wolinsky Y, Kruchko C, Barnholtz-Sloan JS. CBTRUS Statistical Report: primary brain and other central nervous system tumors diagnosed in the United States in 2009-2013. Neuro Oncol. 2016;18(Suppl 5):v1–v75.

    Article  PubMed  Google Scholar 

  2. Ostrom QT, Bauchet L, Davis FG, Deltour I, Fisher JL, Langer CE, Pekmezci M, Schwartzbaum JA, Turner MC, Walsh KM, et al. The epidemiology of glioma in adults: a “state of the science” review. Neuro Oncol. 2014;16(7):896–913.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Wiemels JL, Wiencke JK, Sison JD, Miike R, McMillan A, Wrensch M. History of allergies among adults with glioma and controls. Int J Cancer. 2002;98(4):609–15.

    Article  CAS  PubMed  Google Scholar 

  4. Johansen C, Schüz J, Andreasen A-MS, Dalton SO. Study designs may influence results: the problems with questionnaire-based case–control studies on the epidemiology of glioma. Br J Cancer. 2017;116(7):841–8.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Schlehofer B, Blettner M, Preston-Martin S, Niehoff D, Wahrendorf J, Arslan A, Ahlbom A, Choi WN, Giles GG, Howe GR, et al. Role of medical history in brain tumour development. Results from the international adult brain tumour study. Int J Cancer. 1999;82(2):155–60.

    Article  CAS  PubMed  Google Scholar 

  6. Cicuttini FM, Hurley SF, Forbes A, Donnan GA, Salzberg M, Giles GG, McNeil JJ. Association of adult glioma with medical conditions, family and reproductive history. Int J Cancer. 1997;71(2):203–7.

    Article  CAS  PubMed  Google Scholar 

  7. Schwartzbaum J, Jonsson F, Ahlbom A, Preston-Martin S, Lonn S, Soderberg KC, Feychting M. Cohort studies of association between self-reported allergic conditions, immune-related diagnoses and glioma and meningioma risk. Int J Cancer. 2003;106(3):423–8.

    Article  CAS  PubMed  Google Scholar 

  8. Wiemels JL, Wiencke JK, Patoka J, Moghadassi M, Chew T, McMillan A, Miike R, Barger G, Wrensch M. Reduced immunoglobulin E and allergy among adults with glioma compared with controls. Cancer Res. 2004;64(22):8468–73.

    Article  CAS  PubMed  Google Scholar 

  9. Schlehofer B, Siegmund B, Linseisen J, Schuz J, Rohrmann S, Becker S, Michaud D, Melin B, Bas Bueno-de-Mesquita H, Peeters PH, et al. Primary brain tumours and specific serum immunoglobulin E: a case-control study nested in the European Prospective Investigation into Cancer and Nutrition cohort. Allergy. 2011;66(11):1434–41.

    Article  CAS  PubMed  Google Scholar 

  10. Calboli FC, Cox DG, Buring JE, Gaziano JM, Ma J, Stampfer M, Willett WC, Tworoger SS, Hunter DJ, Camargo CA Jr, et al. Prediagnostic plasma IgE levels and risk of adult glioma in four prospective cohort studies. J Natl Cancer Inst. 2011;103(21):1588–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Wiemels JL, Wilson D, Patel C, Patoka J, McCoy L, Rice T, Schwartzbaum J, Heimberger A, Sampson JH, Chang S, et al. IgE, allergy, and risk of glioma: update from the San Francisco Bay Area Adult Glioma Study in the temozolomide era. Int J Cancer. 2009;125(3):680–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Linos E, Raine T, Alonso A, Michaud D. Atopy and risk of brain tumors: a meta-analysis. J Natl Cancer Inst. 2007;99(20):1544–50.

    Article  PubMed  Google Scholar 

  13. Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet. 2014;23(R1):R89–98.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Meyers DA, Xu J, Postma DS, Levitt RC, Bleecker ER. Two locus segregation and linkage analysis for total serum IgE levels. Clin Exp Allergy. 1995;25:113–5.

    Article  PubMed  Google Scholar 

  15. Wilkinson J, Grimley S, Collins A, Simon Thomas N, Holgate ST, Morton N. Linkage of asthma to markers on chromosome 12 in a sample of 240 families using quantitative phenotype scores. Genomics. 1998;53(3):251–9.

    Article  CAS  PubMed  Google Scholar 

  16. Portelli MA, Hodge E, Sayers I. Genetic risk factors for the development of allergic disease identified by genome-wide association. Clin Exp Allergy. 2015;45(1):21–31.

    Article  CAS  PubMed  Google Scholar 

  17. Pierce BL, Burgess S. Efficient design for Mendelian randomization studies: subsample and 2-sample instrumental variable estimators. Am J Epidemiol. 2013;178(7):1177–84.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Melin BS, Barnholtz-Sloan JS, Wrensch MR, Johansen C, Il'yasova D, Kinnersley B, Ostrom QT, Labreche K, Chen Y, Armstrong G, et al. Genome-wide association study of glioma subtypes identifies specific differences in genetic susceptibility to glioblastoma and non-glioblastoma tumors. Nat Genet. 2017;49(5):789–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Granada M, Wilk JB, Tuzova M, Strachan DP, Weidinger S, Albrecht E, Gieger C, Heinrich J, Himes BE, Hunninghake GM, et al. A genome wide association study of plasma total IgE concentration in the Framingham Heart Study. J Allergy Clin Immunol. 2012;129(3):840–5.

    Article  CAS  PubMed  Google Scholar 

  20. Baurecht H, Hotze M, Brand S, Büning C, Cormican P, Corvin A, Ellinghaus D, Ellinghaus E, Esparza-Gordillo J, Fölster-Holst R, et al. Genome-wide comparative analysis of atopic dermatitis and psoriasis gives insight into opposing genetic mechanisms. Am J Hum Genet. 2015;96(1):104–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Paternoster L, Standl M, Waage J, Baurecht H, Hotze M, Strachan DP, Curtin JA, Bønnelykke K, Tian C, Takahashi A, et al. Multi-ethnic genome-wide association study of 21,000 cases and 95,000 controls identifies new risk loci for atopic dermatitis. Nat Genet. 2015;47(12):1449–56.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Schaarschmidt H, Ellinghaus D, Rodríguez E, Kretschmer A, Baurecht H, Lipinski S, Meyer-Hoffert U, Harder J, Lieb W, Novak N, et al. A genome-wide association study reveals 2 new susceptibility loci for atopic dermatitis. J Allergy Clin Immunol. 2015;136(3):802–6.

    Article  CAS  PubMed  Google Scholar 

  23. Weidinger S, Willis-Owen SAG, Kamatani Y, Baurecht H, Morar N, Liang L, Edser P, Street T, Rodriguez E, O'Regan GM, et al. A genome-wide association study of atopic dermatitis identifies loci with overlapping effects on asthma and psoriasis. Hum Mol Genet. 2013;22(23):4841–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Ferreira MA, Matheson MC, Tang CS, Granell R, Ang W, Hui J, Kiefer AK, Duffy DL, Baltic S, Danoy P, et al. Genome-wide association analysis identifies 11 risk variants associated with the asthma with hay fever phenotype. J Allergy Clin Immunol. 2014;133(6):1564–71.

    Article  CAS  PubMed  Google Scholar 

  25. Ramasamy A, Curjuric I, Coin LJ, Kumar A, McArdle WL, Imboden M, Leynaert B, Kogevinas M, Schmid-Grendelmeier P, Pekkanen J, et al. A genome-wide meta-analysis of genetic variants associated with allergic rhinitis and grass sensitization and their interaction with birth order. J Allergy Clin Immunol. 2011;128(5):996–1005.

    Article  CAS  PubMed  Google Scholar 

  26. Hinds DA, McMahon G, Kiefer AK, Do CB, Eriksson N, Evans DM, St Pourcain B, Ring SM, Mountain JL, Francke U, et al. A genome-wide association meta-analysis of self-reported allergy identifies shared and allergy-specific susceptibility loci. Nat Genet. 2013;45(8):907–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Louis DN, Perry A, Reifenberger G, von Deimling A, Figarella-Branger D, Cavenee WK, Ohgaki H, Wiestler OD, Kleihues P, Ellison DW. The 2016 World Health Organization Classification of Tumors of the Central Nervous System: a summary. Acta Neuropathol. 2016;131(6):803–20.

    Article  PubMed  Google Scholar 

  28. Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. 2013;37(7):658–65.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Zhu Z, Zheng Z, Zhang F, Wu Y, Trzaskowski M, Maier R, Robinson M, McGrath J, Visscher P, Wray N, et al. Causal associations between risk factors and common diseases inferred from GWAS summary data. Nat Commun. 2018;9(1):224.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40(4):304–14.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Hartwig FP, Smith GD, Bowden J. Robust inference in summary data Mendelian randomization via the zero model pleiotropy assumption. Int J Epidemiol. 2017;46(6):1985–98.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Gage SH, Jones HJ, Burgess S, Bowden J, Davey Smith G, Zammit S, Munafo MR. Assessing causality in associations between cannabis use and schizophrenia risk: a two-sample Mendelian randomization study. Psychol Med. 2017;47(5):971–80.

    Article  CAS  PubMed  Google Scholar 

  33. Burgess S. Sample size and power calculations in Mendelian randomization with a single instrumental variable and a binary outcome. Int J Epidemiol. 2014;43(3):922–9.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Paternoster L, Standl M, Johannes W, Baurecht H, Hotze M, Strachan DP, Curtin JA. Multi-ancestry genome-wide association study of 21,000 cases and 95,000 controls identifies new risk loci for atopic dermatitis. Nat Genet. 2015;47(12):1449–56.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Weidinger S, Gieger C, Rodriguez E, Baurecht H, Mempel M, Klopp N, Gohlke H, Wagenpfeil S, Ollert M, Ring J, et al. Genome-wide scan on total serum IgE levels identifies FCER1A as novel susceptibility locus. PLoS Genet. 2008;4(8):e1000166.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Ramasamy A, Kuokkanen M, Vedantam S, Gajdos ZK, Couto Alves A, Lyon HN, Ferreira MAR, Strachan DP, Zhao JH, Abramson MJ, et al. Genome-wide association studies of asthma in population-based cohorts confirm known and suggested loci and identify an additional association near HLA. PLoS One. 2012;7(9):e44008.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. MacArthur J, Bowler E, Cerezo M, Gil L, Hall P, Hastings E, Junkins H, McMahon A, Milano A, Morales J, et al. The new NHGRI-EBI Catalog of published genome-wide association studies (GWAS Catalog). Nucleic Acids Res. 2017;45(D1):D896–901.

    Article  CAS  PubMed  Google Scholar 

  38. Ahmad OS, Morris JA, Mujammami M, Forgetta V, Leong A, Li R, Turgeon M, Greenwood CM, Thanassoulis G, Meigs JB, et al. A Mendelian randomization study of the effect of type-2 diabetes on coronary heart disease. Nat Commun. 2015;6:7060.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Palmer TM, Sterne JA, Harbord RM, Lawlor DA, Sheehan NA, Meng S, Granell R, Smith GD, Didelez V. Instrumental variable estimation of causal risk ratios and causal odds ratios in Mendelian randomization analyses. Am J Epidemiol. 2011;173(12):1392–403.

    Article  PubMed  Google Scholar 

  40. Amirian ES, Zhou R, Wrensch MR, Olson SH, Scheurer ME, Il’yasova D, Lachance D, Armstrong GN, McCoy LS, et al. Approaching a scientific consensus on the association between allergies and glioma risk: a report from the Glioma International Case-Control Study. Cancer Epidemiol Biomarkers Prev. 2016;25(2):282–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Scheurer ME, El-Zein R, Thompson PA, Aldape KD, Levin VA, Gilbert MR, Weinberg JS, Bondy ML. Long-term anti-inflammatory and antihistamine medication use and adult glioma risk. Cancer Epidemiol Biomarkers Prev. 2008;17(5):1277–81.

    Article  CAS  PubMed  Google Scholar 

  42. Amirian ES, Marquez-Do D, Bondy ML, Scheurer ME. Antihistamine use and immunoglobulin E levels in glioma risk and prognosis. Cancer Epidemiol. 2013;37(6):908–12.

    Article  PubMed  Google Scholar 

  43. Razavi S-M, Lee KE, Jin BE, Aujla PS, Gholamin S, Li G. Immune evasion strategies of glioblastoma. Front Surg. 2016;3:11.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Gustafson MP, Lin Y, New KC, Bulur PA, O'Neill BP, Gastineau DA, Dietz AB. Systemic immune suppression in glioblastoma: the interplay between CD14(+)HLA-DR(lo/neg) monocytes, tumor factors, and dexamethasone. Neuro Oncol. 2010;12(7):631–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44(2):512–25.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Burgess S, Bowden J, Fall T, Ingelsson E, Thompson SG. Sensitivity analyses for robust causal inference from Mendelian randomization analyses with multiple genetic variants. Epidemiology. 2017;28(1):30–42.

    Article  PubMed  Google Scholar 

  47. Il'yasova D, McCarthy B, Marcello J, Schildkraut JM, Moorman PG, Krishnamachari B, Ali-Osman F, Bigner DD, Davis F. Association between glioma and history of allergies, asthma, and eczema: a case-control study with three groups of controls. Cancer Epidemiol Biomarkers Prev. 2009;18(4):1232–8.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

LD-H was supported by a Wellcome Trust Summer Student bursary. AS is supported by a Cancer Research UK clinical Fellowship. In the UK, funding was provided by Cancer Research UK (C1298/A8362) supported by the Bobby Moore Fund. The Glioma International Case-Control Consortium Study was supported by grants from the National Institutes of Health, Bethesda, Maryland (R01CA139020, R01CA52689, P50097257, P30CA125123). The UK Interphone Study was supported by the European Commission Fifth Framework Program “Quality of Life and Management of Living Resources” and the UK Mobile Telecommunications and Health Programme. The Mobile Manufacturers Forum and the GSM Association provided funding for the study through the scientifically independent International Union against Cancer (UICC).

Availability of data and materials

Genotype data from the Glioma International Case-Control Consortium Study GWAS are available from the database of Genotypes and Phenotypes (dbGaP) under accession phs001319.v1.p1. Additionally, genotypes from the GliomaScan GWAS can be accessed through dbGaP accession phs000652.v1.p1.

Author information

Authors and Affiliations

Authors

Contributions

RSH and AJC managed the project. LD-H, AJC, AS, PJL and RSH drafted the manuscript. LD-H and AJC performed statistical analyses. BK, KL, MJS and RSH acquired and analysed the UK data. MSi, PH, MMN and K-HJ acquired and analysed the German data. DIJ, QTO, JEE-P, GNA, EBC, DI, JS, JSB-S, SHO, JLB, RKL, CJ, RBJ, BSM, MRW, MLB and RSH acquired and analysed the Glioma International Case-Control Consortium Study data. SC and PR acquired and analysed the National Cancer Institute data. MSa acquired and analysed the French data. All authors reviewed the final manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Richard S. Houlston.

Ethics declarations

Ethics approval and consent to participate

Two-sample Mendelian randomisation was undertaken using GWAS data. Ethical approval was not sought for this specific project because all data came from the summary statistics of published GWAS, and no individual-level data were used.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figure S1. Forest plot of Wald odds ratios (ORs) and 95% confidence intervals generated from single nucleotide polymorphisms (SNPs) associated with atopic dermatitis, including rs909341. ORs for individual SNPs are listed according to magnitude of effect in the instrumental variable analysis and are presented with pooled effects using the inverse-variance weighting method. Squares represent the point estimate, and the bars are the 95% confidence intervals. (DOCX 89 kb)

Additional file 2:

Table S1. Summary of the eight glioma genome-wide association studies. (XLSX 29 kb)

Additional file 3:

Table S2. Table of single nucleotide polymorphisms (SNPs) reported in the NHGRI-EBI Genome-wide Association Studies Catalog for each trait, with correlations between SNPs. (XLSX 48 kb)

Additional file 4:

Table S3. Percentage of variance explained by the combined sets of single nucleotide polymorphisms used as instrumental variables. (XLSX 33 kb)

Additional file 5:

Table S4. Range of odds ratios for which study had < 80% power, for each atopy-related trait (P = 0.05, two-sided). (XLSX 9 kb)

Additional file 6:

Table S5. Simulation analyses. (XLSX 28 kb)

Additional file 7:

Table S6. Inverse-variance weighting, maximum likelihood estimation, weighted median estimate, mode-based estimate and Mendelian randomisation-Egger test results for combined atopy-related instrumental variables and glioma subtypes. (XLSX 39 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Disney-Hogg, L., Cornish, A.J., Sud, A. et al. Impact of atopy on risk of glioma: a Mendelian randomisation study. BMC Med 16, 42 (2018). https://doi.org/10.1186/s12916-018-1027-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-018-1027-5

Keywords