Association Between Environmental Factors and Asthma Using Mendelian Randomization: Increased Effect of Body Mass Index on Adult-Onset Moderate-to-Severe Asthma Subtypes

Although asthma is one of the most common chronic diseases throughout all age groups, its etiology remains unknown, primarily due to its heterogeneous characteristics. We examined the causal effects of various environmental factors on asthma using Mendelian randomization and determined whether the susceptibility to asthma due to the causal effect of a risk factor differs between asthma subtypes, based on age of onset, severity of asthma, and sex. We performed Mendelian randomization analyses (inverse variance weighted, weighted median, and generalized summary-data-based Mendelian randomization) using UK Biobank data to estimate the causal effects of 69 environmental factors on asthma. Additional sensitivity analyses (MR-Egger regression, Cochran’s Q test, clumping, and reverse Mendelian randomization) were performed to ensure minimal or no pleiotropy. For confirmation, two-sample setting analyses were replicated using BMI SNPs that had been reported by a meta-genome-wide association study in Japanese and European (GIANT) populations and a genome-wide association study in control individuals from the UK Biobank. We found that BMI causally affects the development of asthma and that the adult-onset moderate-to-severe asthma subtype is the most susceptible to causal inference by BMI. Further, it is likely that the female subtype is more susceptible to BMI than males among adult asthma cases. Our findings provide evidence that obesity is a considerable risk factor in asthma patients, particularly in adult-onset moderate-to-severe asthma cases, and that weight loss is beneficial for reducing the burden of asthma.

Introduction

Asthma is one of the most common chronic diseases, affecting children and adults; yet, much remains to be determined with regard to its etiology (Subbarao et al., 2009). Approximately 300 million persons have been afflicted by asthma worldwide (Dharmage et al., 2019), and it is estimated that 250,000 people die prematurely each year due to asthma (Bousquet et al., 2010). If current trends continue, roughly 100 million more persons will develop asthma by 2025 (Dharmage et al., 2019). However, asthma has considerable heterogeneity, including its pathophysiological mechanisms, environmental exposure, comorbidities, underlying disease severity, age of onset, medical accessibility, psychological factors, and medical responsiveness (Bousquet et al., 2010). For this reason, the identification of causal risk factors in a specific asthma case has been a challenging task.

Asthma is characterized by variable narrowing of the airway due to the interaction of airway inflammation and remodeling. Several studies have classified asthma by various clinical criteria, such as age of onset (early onset, late onset) (Ferreira et al., 2019), disease severity (mild, moderate-to-severe) (Shrine et al., 2019), and allergy (atopic, non-atopic) (Zhu et al., 2020). In addition, a gender-specific asthma cluster analysis has reported heterogeneous characteristics of asthma between sexes (Hsiao et al., 2019). T-helper type 2 (Th2) inflammation is regarded as the central molecular mechanism of asthma, and based on Th2 profiles, asthma is classified as Th2-high and Th2-low (Woodruff et al., 2009; Sterk and Lutter, 2014).

The severe form of asthma is aggravated by marked thickening of the airway walls and widespread inflammation (Busse et al., 2000). These changes can reduce lung function and cause the airway to be narrowed rapidly due to smooth muscle contraction (Busse et al., 2000). Thus, individuals with severe asthma present with symptoms that are distinct from those in mild asthma patients, including debilitating lung function, frequent exacerbation of asthma, and increased hospital admissions despite significant use of medicines (corticosteroids) (Shrine et al., 2019). Classification of asthma by phenotype provides a foundation on which to understand disease causality and develop management approaches that improve the control of asthma while avoiding adverse effects and decreasing the risk of serious outcomes (Moore et al., 2010).

Many environmental factors, such as BMI (Taylor et al., 2008), smoking (Hedman et al., 2011), physical inactivity (Cordova et al., 2017; Vlashki et al., 2018), and dietary habits (Protudjer et al., 2012), have been epidemiologically and clinically related to the occurrence or severity of asthma. But, observational studies encounter limitations in determining a causal link, due to their various potential confounding factors (Sun et al., 2020). Randomized controlled trials (RCTs), if appropriately conducted, are the standard for determining causal inferences in health sciences. However, RCTs are expensive and time-consuming, involving thousands of patients, some of whom might experience unwanted side effects of drugs (Roberts, 2018).

With some similarities to RCTs, Mendelian randomization (MR) is an alternative approach for establishing evidence of causal relationships for which RCTs are practically unavailable (Smith and Ebrahim, 2003). MR is a method from genetic epidemiology that uses randomly inherited genetic instruments [single-nucleotide polymorphisms (SNPs)] that are robustly associated with a risk factor as proxies for environmental exposure to assess causal inferences for the effects of exposure on an outcome (Smith and Ebrahim, 2003; Rosoff et al., 2019). The increasing availability of summary-level data from genome-wide association studies (GWASs) in the public domain allows MR to make inferences on causality by integrating summary-level GWAS data from various studies (Bowden et al., 2016).

The aim of our study was to estimate causal inferences of environmental factors on asthma using MR and to examine the differences in the susceptibility to asthma with regard to the causal effect of an environmental factor between asthma subtypes, classified by the age of onset, disease severity, and sex.

Materials and Methods

Study Population and Design

The UK Biobank is a population-based cohort that recruited over 487,409 individuals aged 40–69 years from 2006–2010 (Collins, 2012). For quality control of the samples, we used the following filter parameters of the Neale lab¹ : PCA calculation filter for selection of unrelated samples; sex chromosome filter for removal of aneuploidy; filter of principal components (PCs) for European sample selection to determine British ancestry; and filters for selection of self-reported ‘white British,’ ‘Irish,’ and ‘White.’ The UK Biobank has been granted ethical approval to collect data on participants by the North West Multicentre Research Ethics Committee, the National Information Governance Board for Health & Social Care, and the Community Health Index Advisory Group.

Asthma cases (n = 35,926) were determined as those that had been diagnosed with asthma by a doctor and had checked for age of onset. Participants who had been diagnosed with chronic obstructive pulmonary disease (COPD) were excluded. Controls (n = 227,924) were defined those who had not been diagnosed with asthma, rhinitis, eczema, allergy, or emphysema/chronic bronchitis. In addition, those who had diagnostic records of hay fever, allergic rhinitis, emphysema, chronic bronchitis, or COPD and those who had J40-47 records in the ICD 10 codes (the 10th revision of the International Statistical Classification of Diseases and Related Health Problems) were excluded (Supplementary Material and Supplementary Figure 1).

For the study of specific asthma subtypes, cases were divided into 4 groups by the age of onset and disease severity: child-onset mild (CM, n = 9,758), child-onset moderate-to-severe (CM-S, n = 1,875), adult-onset mild (AM, n = 19,415), and adult-onset moderate-to-severe (AM-S, n = 4,878). Individuals with an age of onset before 19 years were defined as child-onset asthma cases (Ferreira et al., 2019), and those with an age of onset after 20 years were considered adult-onset asthma. Moderate-to-severe asthma cases were selected from individuals for whom, in addition to the conditions above, medication information was available and the British Thoracic Society (BTS) stage 3-5 criteria were met, as described (Shrine et al., 2019). Asthma cases that did not satisfy these conditions were classified as mild asthma.

To determine whether gender differences had causal effects, adult-onset asthma cases and controls were separated by gender: female adult-onset mild (FAM, n = 12,208), female adult-onset moderate-to-severe (FAM-S, n = 3075), male adult-onset mild (MAM, n = 7207), male adult-onset moderate-to-severe (MAM-S, n = 1803), female control (FC, n = 119,515), and male control (MC, n = 108,409).

Genetic Instrumental Variants for Environmental Factors and Asthma

We studied 93 environmental factors in 18 categories that were associated with asthma (Supplementary Table 1, references attached). Association test was performed between asthma and each of the 93 factors by logistic linear regression in R, with adjustments for age, sex, and each environmental factor. As a result, 69 environmental factors were associated with asthma [P < 5.38E-04 (0.05/93)] (Supplementary Table 2). For each of the 69 environmental traits, we extracted SNPs that were significantly associated with each factor (e.g., SNP-environmental factor) (P < 5E-08) using genome-wide summary statistics, provided by the Neale lab UK Biobank GWAS summary data (Supplementary Table 3)². The SNPs-environmental factor were subject to clumping (r² > 0.05, 1-Mbp boundary distance) using FUMA (Watanabe et al., 2017) to ensure the independence of environmental factor-associated loci. As an additional quality control of SNPs for Mendelian randomization, we removed strand-ambiguous SNPs (e.g., A/T and C/G, MAF > 0.42) (Tang et al., 2020) and SNPs in the MHC region (chromosome 6:25-34M) due to their strong pleotropic effects (Supplementary Figure 2 and Supplementary Table 4; Zhu et al., 2020).

For asthma, we performed a genome-wide association analysis using UK Biobank data. Genotyping imputation was performed using the UK10K Project and 1000 Genome Project Phase 3 reference panels (UK10K Consortium, Walter et al., 2015; Genomes Project et al., 2015). General quality control procedures for exclusion (P for Hardy-Weinberg equilibrium test <1E-06, missing genotype call rate >0.05, minor allele frequency <0.01) were applied to 7,402,791 SNPs. In total, 5,664,578 SNPs were retained for further analysis. A GWAS for asthma was performed using PLINK 1.9, with adjustments for age, sex, genetic array, and 10PCs. A list of independent asthma-associated loci (e.g., asthma SNPs) were determined by clumping (P < 5E-08, r² > 0.05, 1-Mbp boundary distance), and SNPs that were strand-ambiguous and in the MHC region were excluded. The resulting quantile-quantile (QQ) plot and Manhattan plot are shown in Supplementary Figure 3. Thus, 142 SNPs were selected for genetic instruments of asthma (Supplementary Table 5).

For the two-sample MR setting of BMI → asthma, we extracted summary association statistics for the 158 genome-wide significant SNPs (P < 5E-08) that were associated with BMI in a trans-ethnic meta-GWAS of 173,430 Japanese subjects (the BioBank Japan project, the Japan Public Health Center-based Prospective Study, and the Tohoku Medical Megabank Project) and 339,224 Europeans (the GIANT consortium) (total N_max = 480,438) (Akiyama et al., 2017). Of 158 BMI SNPs, we removed SNPs with palindromes (e.g., A/T and C/G) (Tang et al., 2020) and SNPs in the MHC region (chromosome 6:25-34M) in the MR analyses (Zhu et al., 2020), retaining 149 SNPs.

To avoid the biases of one-sample settings, such as reverse causality and overfitting, we performed a GWAS for BMI only in controls (n = 227,924) from the UK Biobank data. The resulting genome-wide significant BMI SNPs (P < 5E-08) were subject to clumping using FUMA (Watanabe et al., 2017). The 170 independent BMI SNPs were further subject to the removal of SNPs with strand ambiguity, SNPs in the MHC region, and SNPs that were associated with asthma, leaving 159 SNPs (Supplementary Table 10). The resulting quantile-quantile (QQ) and Manhattan plots are shown in Supplementary Figure 4.

Mendelian Randomization

To assess the causal relationship between environmental factors and asthma, we applied 3 methods: inverse variance weighted (IVW) random effects model (Burgess et al., 2013), weighted median regression (Bowden et al., 2016; Censin et al., 2017), and generalized summary-data-based Mendelian randomization (GSMR) (Zhu et al., 2018). The causal effect estimate by IVW is liable to be biased if any SNP exhibits horizontal pleiotropy. As a complementary method to reduce heterogeneity, we performed GSMR, in which genetic variants were pruned at a high threshold of r² < 0.05 (Pasman et al., 2019) and filtered for pleiotropic effects on exposure and outcome [Heterogeneity In Dependent Instrument (HEIDI) filtering] (Zhu et al., 2018). The weighted median method provides an unbiased estimate of the causal effect even when up to 50% of the information comes from invalid genetic variants (Bowden et al., 2016; Censin et al., 2017).

To ensure minimal or no pleiotropy in our results, we performed additional sensitivity analyses. First, we estimated the intercept by MR-Egger regression, with an intercept that differs significantly from 0 (P < 0.05) as an indication of residual heterogeneity due to directional pleiotropy (Bowden et al., 2015). Next, we evaluated the residual heterogeneity using Cochran’s Q statistic, with significant heterogeneity (P < 0.05) due to horizontal pleiotropy. Then, we removed SNPs with any evidence of pleiotropy by clumping both environmental factor and asthma SNPs (r² > 0.05, 1-Mbp boundary distance) and excluded SNPs that were potentially associated with asthma (P < 0.05/number of environmental factor SNPs) for forward MR and SNPs that were linked to environmental factors (P < 0.05/number of asthma SNPs) for reverse MR and then repeated the MR analyses. Finally, only unidirectional causal effects were determined by performing the reverse MR of asthma → environmental factor (P < 5E-08).

The estimates from the IVW, weighted median, and GSMR were defined as causal effects only when meeting significance after Bonferroni correction for multiple tests as the threshold for the true causal estimate (P < 0.05/number of environmental factors x number of asthma subtypes).

Statistical Analysis

For the genome-wide association analyses, we used a logistic regression model and assumed an additive genetic model of trait status with genotype dose, fitted using the R package and adjusted for covariates, including age, sex, environmental factor, and PC10s. t-tests were used to compare characteristics between control and asthma cases using R.

We used the TwoSampleMR package for performing Mendelian randomization methods, such as IVW, weighted median, and MR-Egger; the gsmr package for GSMR, the qqman package for drawing Manhattan plot and quantile-quantile (QQ) plots; and the ggplot2 package for drawing odds ratio plots, available in the R stats package, version 3.6.3³.

Results

Characteristics of Study Population for Asthma Case and Control Subjects

The basic characteristics of the 263,850 UK Biobank participants (35,926 asthma cases and 227,924 controls) in this study are described in Table 1. BMI, obesity, eosinophil parameters, and female frequency were significantly higher in asthma cases than controls, and there was no difference in smoking status between groups.

TABLE 1

Table 1. Characteristics of the asthma cases and controls from the UK Biobank.

Effects of 69 Environmental Factors on Asthma

Forward MR analyses (environmental factor → asthma) were conducted using IVW, weighted median, and GSMR for 69 factors (Figure 1 and Supplementary Table 4). Factors that were related to white blood cells and anthropometry (leukocyte, eosinophil parameters, BMI, and waist circumference) satisfied the threshold by Bonferroni correction for multiple tests (P < 7.25E–04, 0.05/69) in all 3 methods. The intercept P-value from the MR-Egger regression suggests that there was relatively balanced pleiotropy (P = 0.08 for leukocyte count, P = 0.27 for eosinophil count, P = 0.47 for eosinophil percentage, P = 0.88 for BMI, and P = 0.19 for waist circumference). However, the Cochran’s Q values indicated significant residual heterogeneity in all 4 analyses (Q = 394.9, P-value = 1.65E–30 for leukocyte count; Q = 1095.7, P-value = 3.22E–149 for eosinophil count; Q = 1181.2, P-value = 6.41E–169 for eosinophil percentage; Q = 383.4, P-value = 8.25E–14 for BMI; and Q = 309.6, P-value = 2.46E–10 for waist circumference), which must be improved to determine the true causal effects.

FIGURE 1

Figure 1. MR estimates for associations between 5 environment factors and asthma. The odds ratio is plotted as dots with the 95% confidence interval. White dots are the IVW method, and black dots are the GSMR method. The P value indicates significance when P < 7.25E-04 (0.05/69) by Bonferroni correction. The Q value by Cochran’s Q statistic indicates residual heterogeneity in the MR estimates.

For the sensitivity analyses, we clumped all 873 environmental factor SNPs (123 for leukocyte count, 138 for eosinophil count, 131 for eosinophil percentage, 200 for BMI, 170 for waist circumference, and 111 for asthma; r² < 0.05, 1-Mbp boundary distance; Supplementary Table 15). In addition, we removed environmental factor SNPs that were associated with asthma [P < 5.05E–04 (0.05/99) for leukocyte count SNPs, P < 5.32E–04 (0.05/94) for eosinophil count SNPs, P < 6.17E–04 (0.05/81) for eosinophil percentage SNPs, P < 2.99E–04 (0.05/167) for BMI SNPs, and P < 5.21E–04 (0.05/96) for waist circumference SNPs] and asthma SNPs that were related to environmental factors [P < 5.38E–04 (0.05/93) for asthma SNPs]. RE–analysis of MR using SNPs that were retained after pruning showed that the overall heterogeneity improved, as indicated by the lower Cochran’s Q values (Figure 2A and Supplementary Table 6; Q = 221.7, P-value = 1.03E–12 for leukocyte count; Q = 173.1, P-value = 3.61E–09 for eosinophil count; Q = 98.0, P-value = 8.01E–04 for eosinophil percentage; Q = 254.7, P-value = 4.48E–06 for BMI; and Q = 167.0, P-value = 5.33E–06 for waist circumference). These additional analyses suggest that eosinophil count, eosinophil percentage, and BMI have causal inferences on asthma, whereas the effects of leukocyte count and waist circumference are not significant, because they failed to reach the threshold for significance [P < 1.00E–02 (0.05/5)].

FIGURE 2

Figure 2. MR estimates for associations between 5 environmental factors and asthma after sensitivity analysis. (A) Forward MR: environmental factors → asthma. (B) Reverse MR: asthma → environmental factors. The odds ratio is plotted as dots with the 95% confidence interval. White dots are the IVW method, and black dots are the GSMR method. The P value indicates significance when P < 1.00E-02 (0.05/5) by Bonferroni correction. The Q value by Cochran’s Q statistic indicates residual heterogeneity in the MR estimates.

To determine unidirectional causal effects (environmental factor → asthma), we conducted reverse MR analyses (asthma → environmental factor) using asthma SNPs (Figure 2B and Supplementary Table 6). We found that asthma had causal effects on leukocyte count, eosinophil count, and eosinophil percentage, whereas asthma did not causally affect BMI or waist circumference. The reverse MR results indicate that only BMI causally increases the risk of asthma.

Effects of BMI on the Susceptibility of Asthma Subtypes

Asthma is a highly heterogeneous disease that can be classified by various clinical criteria (Haldar et al., 2008; Zhu et al., 2020). Recent studies indicate that specific asthma subtypes are related to metabolic traits, such as obesity (Jeong et al., 2017). In our study, asthma cases were divided into 4 subsets, based on 2 criteria, combining age of onset and severity: child-onset mild (CM), child-onset moderate-to-severe (CM-S), adult-onset mild (AM), and adult-onset moderate-to-severe (AM-S) (Table 2). The child- and adult-onset moderate-to-severe subgroups had significantly higher BMI, obesity, and eosinophil parameters than their respective mild subgroups (Table 2).

TABLE 2

Table 2. Characteristics of the asthma subtypes and controls from the UK Biobank.

Mendelian randomization (BMI → asthma) was performed using 4 asthma subtypes. As a result, the causal impact of BMI on asthma was stronger in moderate-to-severe versus mild cases before sensitivity analysis, with the adult-onset moderate-to-severe subtype the most susceptible to BMI ([CM versus CM-S] IVW OR/CI95 = 1.02/0.99–1.05 versus 1.10/1.04–1.17, weighted median = 1.02/0.98–1.06 vs. 1.09/1.00–1.20, GSMR = 1.09/0.97–1.23 vs. 1.62/1.24–2.12; [AM vs. AM-S] IVW = 1.06/1.03–1.08 vs. 1.12/1.07–1.16, weighted median = 1.02/0.99–1.05 vs. 1.11/1.05–1.17, GSMR = 1.271.17–1.39 vs. 1.67/1.41–1.97; Supplementary Tables 7, 15). Additional sensitivity analyses improved the heterogeneity, as evidenced by the Cochran’s Q values, establishing BMI as having its strongest effect on asthma in the adult-onset moderate-to-severe subgroup ([CM vs. CM-S] IVW OR/CI95 = 1.02/0.99–1.05 vs. 1.12/1.05–1.20, weighted median = 1.01/0.97–1.06 vs. 1.10/1.00–1.22, GSMR = 1.09/0.96–1.25 vs. 1.74/1.30–2.32; [AM vs. AM-S] IVW = 1.06/1.03–1.08 vs. 1.10/1.06–1.15, weighted median = 1.02/0.99–1.06 vs. 1.09/1.03–1.16, GSMR = 1.25/1.14–1.38 vs. 1.51/1.26–1.82; Figure 3 and Supplementary Table 7). The child-onset moderate-to-severe and adult-onset mild subgroups were marginally susceptible to asthma, based on the IVW and GSMR analyses, but not weighted median analysis, suggesting the weak causality of BMI. These findings suggest that the causality of BMI in asthma increases as the severity of asthma rises in child- and adult-onset asthma cases and that the causal effect of BMI is the strongest in the adult-onset moderate-to-severe subgroup.

FIGURE 3

Figure 3. MR estimates of associations between BMI and asthma subtypes. The odds ratio is plotted as dots with the 95% confidence interval. White dots are the IVW method and black dots are the GSMR method. The P value indicates significance when P < 1.25E-02 (0.05/4) by a Bonferroni correction. The Q value by Cochran’s Q statistic indicates residual heterogeneity in the MR estimates.

Gender-Specific Effect of BMI in Adult Onset Asthma Cases

Observational studies have suggested that BMI causes the development of adult-onset asthma and that the effect of BMI on asthma is greater in female versus male adults (Chen et al., 2002; Guerra et al., 2002; Beuther and Sutherland, 2007), although the gender-specific effect of BMI on asthma is unknown in children (Castro-Rodriguez et al., 2001; Gold et al., 2003). To study the sex-specific causality of BMI, we divided adult-onset asthma cases into 4 subgroups and controls into 2 subgroups by gender: male adult-onset mild (MAM, n = 7207), male adult-onset moderate-to-severe (MAM-S, n = 1803), female adult-onset mild (FAM, n = 12,208), female adult-onset moderate-to-severe (FAM-S, n = 3075), male control (MC, n = 108,409), and female control (FC, n = 119,515). FAM-S and MAM-S patients had higher BMI, obesity, and eosinophil parameters compared with FAM and MAM cases, respectively (Table 3).

TABLE 3

Table 3. Characteristics of the adult-onset asthma cases and controls from the UK Biobank.

Prior to the sensitivity analysis, the results in the adult subgroups showed that female adult-onset cases were more susceptible to the causal effect of BMI than male cases ([MAM versus FAM] IVW OR = 1.05 versus 1.06, weighted median OR = 1.05 versus 1.04, GSMR OR = 1.26 versus 1.35; [MAM-S versus FAM-S] IVW OR = 1.08 versus 1.14, weighted median OR = 1.11 versus 1.09, GSMR OR = 1.43 versus 1.83; Supplementary Tables 8, 15). All 3 MR results in the female moderate-to-severe subgroup support the robust inference of BMI on asthma, with the heterogeneity by Cochran’s test implying little pleiotropy. In the sensitivity analyses using clumping of SNPs, the susceptibility to asthma in the female moderate-to-severe cases was marginally significant, because the weighted median regression failed to satisfy the threshold by Bonferroni correction (P < 1.25E-02, 0.05/4) (Figure 4 and Supplementary Table 8).

FIGURE 4

Figure 4. MR estimates of associations between BMI and adult-onset asthma subtypes. The odds ratio is plotted as dots with the 95% confidence interval. White dots are the IVW method, and black dots are the GSMR method. The P value indicates significance when P < 1.25E-02 (0.05/4) by Bonferroni correction. The Q value by Cochran’s Q statistic indicates residual heterogeneity in the MR estimates.

Effect of BMI on Asthma and Its Subtypes in Two-Sample Mendelian Randomization

To confirm the causal inference of BMI on asthma subtypes without the bias that often occurs in a one-sample setting, we performed two-sample MR using BMI SNPs from another data resource. A total of 149 BMI SNPs were selected from the 158 genome-wide significant SNPs (P < 5E-08) that were associated with BMI in a trans-ethnic meta-GWAS of Japanese individuals (the BioBank Japan consortium) and Europeans (the GIANT consortium) (total N_max = 480,438) (Akiyama et al., 2017). Then, we excluded 2 SNPs that were potentially associated with asthma [P < 3.36E-4 (0.05/149)] to mitigate heterogeneity as a sensitivity analysis.

The MR (BMI → asthma) results of the two-sample setting using 147 BMI SNPs were consistent with the one-sample MR result, implicating BMI as a causal risk factor for asthma [IVW OR = 1.18 (1.07–1.30), P = 5.97E–04; weighted median OR = 1.14 (1.01–1.29), P = 3.68E–02; GSMR OR = 1.17 (1.08–1.28), P = 2.15E–04; Figure 5A and Supplementary Table 9]. Further, in the two-sample MR of the 4 asthma subtypes, the most susceptible was adult-onset moderate-to-severe asthma [IVW OR = 1.35 (1.11–1.64), P = 2.81E–03; weighted median OR = 1.52 (1.12–2.06), P = 6.76E–03; GSMR OR = 1.40 (1.13–1.73), P = 2.24E–03; Figure 5B and Supplementary Table 9], replicating the one-sample MR findings.

FIGURE 5

Figure 5. Two-sample MR estimates of BMI → asthma (A,C) and BMI → 4 asthma subtypes (B,D). (A,B) BMI SNPs were obtained from a meta-GWAS of Japanese (JBB) and European individuals (GIANT). (C,D) BMI SNPs were obtained from a GWAS that was performed in only controls (no asthma) from the UK Biobank. The IVW and GSMR methods are indicated by white and black dots, respectively. The Q value by Cochran’s Q statistic indicates residual heterogeneity in the MR estimates.

To avoid reverse causality or overfitting bias in the one-sample setting, we performed an additional MR analysis using 170 BMI SNPs that were acquired from a genome-wide association analysis with BMI alone in controls (no asthma cases) from the UK Biobank. SNPs were further subjected to clumping by FUMA and exclusion of SNPs that were strand-ambiguous and in the MHC region and SNPs that were potentially associated with asthma (P < 3.05E-4 (0.05/164)) as a sensitivity analysis. Ultimately, 159 BMI SNPs were used for the MR (BMI → asthma). The results confirmed BMI as a causal factor for asthma [IVW OR = 1.16 (1.07–1.26), P = 2.89E–04; weighted median OR = 1.18 (1.07–1.32), P = 1.71E–03; GSMR OR = 1.18 (1.11–1.26), P = 4.37E–08; Figure 5C and Supplementary Table 11]. Additional MR analyses of the 4 asthma subtypes confirmed that the adult-onset moderate-to-severe cases were the most susceptible to the development of asthma due to the causal effect of BMI [IVW OR = 1.37 (1.16–1.62), P = 2.60E–04; weighted median OR = 1.60 (1.25–2.06), P = 2.41E–04; GSMR OR = 1.36 (1.17–1.58), P = 8.84E–05; Figure 5D and Supplementary Table 11].

To determine the effect of childhood BMI on asthma in our setting, we extracted 25 childhood BMI SNPs with associated statistical values from a previous study (Vogelezang et al., 2020). These BMI SNPs were subject to the removal of SNPs with strand ambiguity (e.g., A/T and C/G), SNPs in the MHC region (chromosome 6:25-34M), SNPs with no proxy (r² > 0.8) in UKB, and SNPs that were associated with asthma (P < 2.17E–03, 0.05/23) leaving 21 SNPs. Two-sample MR was performed using these 21 childhood BMI SNPs, and the results are described in Supplementary Table 12. The estimates of childhood BMI in total asthma showed significant effects by IVW (OR/CI95 = 1.11/1.03-1.19, P-value = 3.43E–03) and GSMR (OR/CI95 = 1.12/1.04–1.20, P-value = 3.43E-03) but not by the weighted median method (OR/CI95 = 1.09/0.99–1.21, P-value = 0.07). However, no significant causal effect was found between childhood BMI and the 4 asthma subtypes by IVW, GSMR, and weighted median, consistent with the previous report (Au Yeung et al., 2021). Based on the previous study (Au Yeung et al., 2021) and our result, childhood BMI has a weak causal impact on asthma but not on any specific subtype.

Discussion

In this study, we used MR to examine the causal relationship between 69 environmental factors and asthma and noted the following: BMI is a causal risk factor for asthma without reverse causation; the effect of BMI on asthma is strongest in the adult-onset moderate-to-severe asthma subgroup; and finally, female subtypes are more prone to asthma due to increased BMI than male subtypes in adults.

Epidemiological and genetic studies that used MR have suggested that BMI is causal factor in asthma (Skaaby et al., 2018; Xu et al., 2019; Sun et al., 2020). In addition, BMI is a risk factor for late-onset asthma (onset age >16) and atopic asthma (Zhu et al., 2020). Our study confirms that BMI is a risk factor for asthma and demonstrates that the causal effect of BMI increases significantly in individuals with child- and adult-onset asthma, exacerbating the asthma. Further, the MR result on the stronger effect of BMI in female subgroups is consistent with observational studies (Chen et al., 2002; Guerra et al., 2002; Beuther and Sutherland, 2007).

In the Epidemiology and Natural History of Asthma: Outcomes and Treatment Regimens (TENOR) study on severe asthma, approximately 57% of individuals with severe asthma were obese, implying a high prevalence of obesity in severe asthma cases compared with an obesity rate of 35% in non-asthma adults in the general United States population (Schatz et al., 2014). In our study, the obesity rate in the adult moderate-to-severe subgroup was significantly higher than in the adult mild subgroup by 7%. Another study found that lung function improved after weight loss in obese patients with asthma, suggesting that greater obesity is related to the severity of asthma (Hakala et al., 2000). There are several potential mechanisms by which BMI is linked to asthma. Obesity has been related to multiple traits of asthma, including eosinophil levels (Kim et al., 2014), lung function (Salome et al., 2010), and allergy (Luo et al., 2013). Studies have suggested that adipokines, such as leptin and adiponectin, are associated with the development and severity of asthma and mediate the exacerbation of asthma through the regulation of eosinophil survival and trafficking (Kim et al., 2014; Zhang et al., 2017; Zheng et al., 2018).

We initially aimed to identify environmental factors that cause the development of asthma. However, of 69 factors, only BMI was identified as a causal influence in asthma. We assume that our study was limited in obtaining the appropriate instruments for certain phenotypes. Environmental data from self-reported questionnaires (e.g., dietary intake, neuroticism, alcohol, smoking, sociodemographic factors, and physical activities) are prone to responder bias (Rask-Andersen et al., 2017). Thus, the presence of an interviewer is recommended to reduce the likelihood of responder bias when obtaining self-reported questionnaire data.

Further, certain anthropometric factors (e.g., fat and non-fat mass) were measured with a less accurate method (bioelectric impedance using a Tanita BC418MA body composition analyzer; UK Biobank) than such techniques as dual-energy x-ray absorptiometry (Speed et al., 2019). Our primary results show that fat mass and percentage, but not fat free mass, are causal factors in asthma by IVW and GSMR (Supplementary Table 4). However, the weighted median method does not support the causal effect of fat mass and percentage. A valid association analysis using precise values can improve the statistics of MR analyses.

Our MR analyses of 69 environmental factors → asthma used only the UK Biobank data as a resource. The one-sample MR setting has several benefits: MR and epidemiological findings can be compared in the same individuals, the validity of causal inferences is unaffected by differences in population characteristics when using 2 samples, and harmonization of genetic variants across datasets is not required (Burgess et al., 2019). However, a limitation of the one-sample analysis is that if the links between genetic instrument and exposure are weak, the causal estimation might suffer from reflection of the confounded associations between exposure and outcome and inflation of false positive (type 1 error) rates (Burgess et al., 2019). Bias in a one-sample analysis with a binary disease outcome can be avoided, such that genetic associations with the exposure are estimated in the controls only; consequently, genetic associations with exposure and outcome will not be correlated (Gharahkhani et al., 2019; Hindy et al., 2019; Wade et al., 2019). We replicated the two-sample MR analyses of BMI → asthma using BMI SNPs from another source of GWAS summary statistics and BMI SNPs from a GWAS only in the control, confirming the one-sample analysis findings.

Although our study stratified asthma by age of onset and sex, concerns remain, because these stratifications were not applied when obtaining genome-wide association statistics on BMI. Genetic studies on BMI have suggested significantly positive genetic correlations between childhood and adult BMI (r_g = 0.76, P-value = 1.45E-112) (Vogelezang et al., 2020) and between male and female BMI (r_g = 0.879, P-value = 5.9E-4) (Yang et al., 2015). Thus, it is likely that depending on the instrument source, BMI might have disparate causal effects on asthma subtypes. Notably, a recent study that used the IVW method confirmed the causal impact of adult BMI on asthma, whereas the possible impact of childhood BMI on the risk of asthma was less clear, mediated predominantly by its relationship with adult BMI, implicating that children with high BMI can reduce their risk of asthma by becoming normal-weight adults. This study was limited, in that there were far fewer childhood BMI SNPs (N = 14 and 5) than adult BMI SNPs (N = 323 and 115), decreasing the power of the MR estimation (Au Yeung et al., 2021). Our two-sample MR using 25 childhood BMI SNPs (Vogelezang et al., 2020) supports a causal relationship with asthma. However, based on the previous study (Au Yeung et al., 2021) and our result, childhood BMI has a weak causal impact on asthma but not on any specific subtype.

There is much evidence that suggests gender-specific effects of BMI on asthma. Previous epidemiological reports have suggested that the incidence and symptoms of adult asthma are higher and more severe in women than in males (Chen et al., 2003; Zein and Erzurum, 2015). A recent study that performed sex-specific transcriptomics in 5 tissues from asthma patients also showed sexual dimorphism in asthma, including sex-specific dysregulation of genes and signaling pathways (Gautam et al., 2019). Moreover, the effect of BMI on asthma is greater in female than male adults (Chen et al., 2002; Guerra et al., 2002; Beuther and Sutherland, 2007). Consistent with these reports, the asthma cases and M-S asthma subgroups in our study included more female than male adults (Table 3). Further, FAM and FAM-S subtypes had significantly increased obesity compared with MAM and MAM-S subtypes (Table 3). Although we observed that female adult subtypes are more susceptible to the causal effect of BMI than the male groups, the effects did not meet our strict criteria. We speculate that the genetic correlations between male and female asthma subtypes are too high to render any distinctive causal patterns (Supplementary Table 14) compared with correlations between the child- and adult-onset mild and moderate-to-severe subtypes (Supplementary Table 13). For these reasons, it is unlikely that gender-specific BMI instruments have a causal effect on gender-specific asthma subtypes.

In conclusion, our data indicate that elevated BMI levels are causally related to the risk of adult-onset moderate-to-severe asthma. Thus, reducing body weight can help alleviate the susceptibility to moderate-to-severe asthma.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics Statement

The studies involving human participants were reviewed and approved by the North West Multicentre Research Ethics Committee, the National Information Governance Board for Health & Social Care, and the Community Health Index Advisory Group. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

T-WH, BO, and J-OK designed the study. T-WH analyzed the data and wrote the first draft of the manuscript. J-OK revised the manuscript. BO, JL, HK, and J-OK collected the data and provided technical support. All authors contributed to the interpretation of the results and critical revision of the manuscript for important intellectual content and approved the final version of the manuscript.

Funding

This work was supported by a grant from Kyung Hee University in 2019 (KHU-20191224) and the National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT) [The Bio & Medical Technology Development Program of the NRF (2019M3E5D3073365) and NRF-2019R1A2C1083980].

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This study was conducted using the resources of UK Biobank (application 56987).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.639905/full#supplementary-material

Abbreviations

AM, adult-onset mild; AM-S, adult-onset moderate-to-severe; BMI, body mass index; BTS, British Thoracic Society; CI, confidence interval; CM, child-onset mild; CM-S, child-onset moderate-to-severe; FC, female control; FAM, female adult-onset mild; FAM-S, female adult-onset moderate-to-severe; GIANT Consortium, The genetic investigation of anthropometric traits consortium; GSMR, generalized summary-data-based mendelian randomization; GWAS, genome-wide association study; HEIDI, heterogeneity in dependent instrument; IVW, inverse-variance weighted; MR, mendelian randomization; MC, male control; MAM, male adult-onset mild; MAM-S, male adult-onset moderate-to-severe; OR, odds ratio; SD, standard deviation; SNP, single-nucleotide polymorphism.

Footnotes

1. ^ https://github.com/Nealelab/UK_Biobank_GWAS
2. ^ http://www.nealelab.is/uk-biobank
3. ^ www.r-project.org

References

Akiyama, M., Okada, Y., Kanai, M., Takahashi, A., Momozawa, Y., Ikeda, M., et al. (2017). Genome-wide association study identifies 112 new loci for body mass index in the Japanese population. Nat. Genet. 49, 1458–1467. doi: 10.1038/ng.3951

PubMed Abstract | CrossRef Full Text | Google Scholar

Au Yeung, S. L., Li, A. M., and Schooling, C. M. (2021). A life course approach to elucidate the role of adiposity in asthma risk: evidence from a Mendelian randomisation study. J. Epidemiol. Commun. Health. 75, 277–281.

Google Scholar

Beuther, D. A., and Sutherland, E. R. (2007). Overweight, obesity, and incident asthma - A meta-analysis of prospective epidemiologic studies. Am. J. Resp. Crit. Care 175, 661–666. doi: 10.1164/rccm.200611-1717oc

PubMed Abstract | CrossRef Full Text | Google Scholar

Bousquet, J., Mantzouranis, E., Cruz, A. A., Ait-Khaled, N., Baena-Cagnani, C. E., Bleecker, E. R., et al. (2010). Uniform definition of asthma severity, control, and exacerbations: document presented for the World Health Organization Consultation on Severe Asthma. J. Allergy Clin. Immun. 126, 926–938. doi: 10.1016/j.jaci.2010.07.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowden, J., Smith, G. D., and Burgess, S. (2015). Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int. J. Epidemiol. 44, 512–525. doi: 10.1093/ije/dyv080

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowden, J., Smith, G. D., Haycock, P. C., and Burgess, S. (2016). Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genet. Epidemiol. 40, 304–314. doi: 10.1002/gepi.21965

PubMed Abstract | CrossRef Full Text | Google Scholar

Burgess, S., Butterworth, A., and Thompson, S. G. (2013). Mendelian randomization analysis with multiple genetic variants using summarized data. Genet. Epidemiol. 37, 658–665. doi: 10.1002/gepi.21758

PubMed Abstract | CrossRef Full Text | Google Scholar

Burgess, S., Davey Smith, G., Davies, N. M., Dudbridge, F., Gill, D., Glymour, M. M., et al. (2019). Guidelines for performing Mendelian randomization investigations. Wellcome Open Res. 4:186. doi: 10.12688/wellcomeopenres.15555.1

CrossRef Full Text | Google Scholar

Busse, W. W., Banks-Schlegel, S., and Wenzel, S. E. (2000). Pathophysiology of severe asthma. J. Allergy Clin. Immunol. 106, 1033–1042.

Google Scholar

Castro-Rodriguez, L. A., Holberg, C. J., Morgan, W. J., Wright, A. L., and Martinez, F. D. (2001). Increased incidence of asthmalike symptoms in girls who become overweight or obese during the school years. Am. J. Resp. Crit. Care 163, 1344–1349. doi: 10.1164/ajrccm.163.6.2006140

PubMed Abstract | CrossRef Full Text | Google Scholar

Censin, J. C., Nowak, C., Cooper, N., Bergsten, P., Todd, J. A., and Fall, T. (2017). Childhood adiposity and risk of type 1 diabetes: a mendelian randomization study. Plos Med. 14:e1002362. doi: 10.1371/journal.pmed.1002362

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Dales, R., Tang, M., and Krewski, D. (2002). Obesity may increase the incidence of asthma in women but not in men: longitudinal observations from the Canadian National Population Health Surveys. Am. J. Epidemiol. 155, 191–197. doi: 10.1093/aje/155.3.191

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Stewart, P., Johansen, H., McRae, L., and Taylor, G. (2003). Sex difference in hospitalization due to asthma in relation to age. J. Clin. Epidemiol. 56, 180–187. doi: 10.1016/s0895-4356(02)00593-0

CrossRef Full Text | Google Scholar

Collins, R. (2012). What makes UK biobank special? Lancet 379, 1173–1174. doi: 10.1016/s0140-6736(12)60404-8

CrossRef Full Text | Google Scholar

Cordova, L., Gibson, P., Gardiner, P., and McDonald, V. (2017). Physical inactivity and sedentary time in severe asthma: prevalence and associations. Eur. Respir J. 50:A775.

Google Scholar

Dharmage, S. C., Perret, J. L., and Custovic, A. (2019). Epidemiology of asthma in children and adults. Front. Pediatr. 7:246. doi: 10.3389/fped.2019.00246

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferreira, M. A. R., Mathur, R., Vonk, J. M., Szwajda, A., Brumpton, B., Granell, R., et al. (2019). Genetic architectures of childhood- and adult-onset asthma are partly distinct. Am. J. Hum. Genet. 104, 665–684. doi: 10.1016/j.ajhg.2019.02.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Gautam, Y., Afanador, Y., Abebe, T., Lopez, J. E., and Mersha, T. B. (2019). Genome-wide analysis revealed sex-specific gene expression in asthmatics. Hum. Mol. Genet. 28, 2600–2614. doi: 10.1093/hmg/ddz074

PubMed Abstract | CrossRef Full Text | Google Scholar

Genomes Project, C., Auton, A., Brooks, L. D., Durbin, R. M., Garrison, E. P., Kang, H. M., et al. (2015). A global reference for human genetic variation. Nature 526, 68–74. doi: 10.1038/nature15393

PubMed Abstract | CrossRef Full Text | Google Scholar

Gharahkhani, P., Ong, J. S., An, J. Y., Law, M. H., Whiteman, D. C., Neale, R. E., et al. (2019). Effect of increased body mass index on risk of diagnosis or death from cancer. Br. J. Cancer 120, 565–570. doi: 10.1038/s41416-019-0386-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Gold, D. R., Damokosh, A. I., Dockery, D. W., and Berkey, C. S. (2003). Body-mass index as a predictor of incident asthma in a prospective cohort of children. Pediatr. Pulm. 36, 514–521. doi: 10.1002/ppul.10376

PubMed Abstract | CrossRef Full Text | Google Scholar

Guerra, S., Sherrill, D. L., Bobadilla, A., Martinez, F. D., and Barbee, R. A. (2002). The relation of body mass index to asthma, chronic bronchitis, and emphysema. Chest 122, 1256–1263. doi: 10.1378/chest.122.4.1256

PubMed Abstract | CrossRef Full Text | Google Scholar

Hakala, K., Stenius-Aarniala, B., and Sovijarvi, A. (2000). Effects of weight loss on peak flow variability, airways obstruction, and lung volumes in obese patients with asthma. Chest 118, 1315–1321. doi: 10.1378/chest.118.5.1315

PubMed Abstract | CrossRef Full Text | Google Scholar

Haldar, P., Pavord, I. D., Shaw, D. E., Berry, M. A., Thomas, M., Brightling, C. E., et al. (2008). Cluster analysis and clinical asthma phenotypes. Am. J. Resp. Crit. Care 178, 218–224.

Google Scholar

Hedman, L., Bjerg, A., Sundberg, S., Forsberg, B., and Ronmark, E. (2011). Both environmental tobacco smoke and personal smoking is related to asthma and wheeze in teenagers. Thorax 66, 20–25. doi: 10.1136/thx.2010.143800

PubMed Abstract | CrossRef Full Text | Google Scholar

Hindy, G., Akesson, K. E., Melander, O., Aragam, K. G., Haas, M. E., Nilsson, P. M., et al. (2019). Cardiometabolic polygenic risk scores and osteoarthritis outcomes: a mendelian randomization study using data from the malmo diet and cancer study and the UK biobank. Arthritis Rheumatol. 71, 925–934. doi: 10.1002/art.40812

PubMed Abstract | CrossRef Full Text | Google Scholar

Hsiao, H. P., Lin, M. C., Wu, C. C., Wang, C. C., and Wang, T. N. (2019). Sex-specific asthma phenotypes, inflammatory patterns, and asthma control in a cluster analysis. J. Allergy Clin. Imm Pract. 7, 556.e15–567.e15.

Google Scholar

Jeong, A., Imboden, M., Hansen, S., Zemp, E., Bridevaux, P. O., Lovison, G., et al. (2017). Heterogeneity of obesity-asthma association disentangled by latent class analysis, the SAPALDIA cohort. Respir. Med. 125, 25–32. doi: 10.1016/j.rmed.2017.02.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, S. H., Sutherland, E. R., and Gelfand, E. W. (2014). Is there a link between obesity and asthma? Allergy Asthma Immun. 6, 189–195. doi: 10.4168/aair.2014.6.3.189

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, X., Xiang, J., Dong, X. H., Cai, F. W., Suo, J. N., Wang, Z. Q., et al. (2013). Association between obesity and atopic disorders in Chinese adults: an individually matched case-control study. BMC Public Health 13:12. doi: 10.1186/1471-2458-13-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Moore, W. C., Meyers, D. A., Wenzel, S. E., Teague, W. G., Li, H. S., Li, X. N., et al. (2010). Identification of Asthma phenotypes using cluster analysis in the severe asthma research program. Am. J. Resp. Crit. Care 181, 315–323.

Google Scholar

Pasman, J. A., Verweij, K. J. H., Gerring, Z., Stringer, S., Sanchez-Roige, S., Treur, J. L., et al. (2019). GWAS of lifetime cannabis use reveals new risk loci, genetic overlap with psychiatric traits, and a causal effect of schizophrenia liability. Nat. Neurosci. 22, 1196–1196. doi: 10.1038/s41593-019-0402-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Protudjer, J. L. P., Sevenhuysen, G. P., Ramsey, C. D., Kozyrskyj, A. L., and Becker, A. B. (2012). Low vegetable intake is associated with allergic asthma and moderate-to-severe airway hyperresponsiveness. Pediatr. Pulm. 47, 1159–1169. doi: 10.1002/ppul.22576

PubMed Abstract | CrossRef Full Text | Google Scholar

Rask-Andersen, M., Karlsson, T., Ek, W. E., and Johansson, A. (2017). Gene-environment interaction study for BMI reveals interactions between genetic factors and physical activity, alcohol consumption and socioeconomic status. PLoS Genet. 13:e1006977. doi: 10.1371/journal.pgen.1006977

PubMed Abstract | CrossRef Full Text | Google Scholar

Roberts, R. (2018). Mendelian randomization studies promise to shorten the journey to FDA approval. JACC 3, 690–703. doi: 10.1016/j.jacbts.2018.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosoff, D. B., Clarke, T. K., Adams, M. J., McIntosh, A. M., Davey Smith, G., Jung, J., et al. (2019). Educational attainment impacts drinking behaviors and risk for alcohol dependence: results from a two-sample Mendelian randomization study with ~780,000 participants. Mol. Psychiatry 26, 1119–1132. doi: 10.1038/s41380-019-0535-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Salome, C. M., King, G. G., and Berend, N. (2010). Physiology of obesity and effects on lung function. J. Appl. Physiol. 108, 206–211. doi: 10.1152/japplphysiol.00694.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Schatz, M., Hsu, J. W. Y., Zeiger, R. S., Chen, W. S., Dorenbaum, A., Chipps, B. E., et al. (2014). Phenotypes determined by cluster analysis in severe or difficult-to-treat asthma. J. Allergy Clin. Immun. 133, 1549–1556. doi: 10.1016/j.jaci.2013.10.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Shrine, N., Portelli, M. A., John, C., Artigas, M. S., Bennett, N., Hall, R., et al. (2019). Moderate-to-severe asthma in individuals of European ancestry: a genome-wide association study. Lancet Resp. Med. 7, 20–34.

Google Scholar

Skaaby, T., Taylor, A. E., Thuesen, B. H., Jacobsen, R. K., Friedrich, N., Mollehave, L. T., et al. (2018). Estimating the causal effect of body mass index on hay fever, asthma and lung function using Mendelian randomization. Allergy. 73, 153–164. doi: 10.1111/all.13242

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, G. D., and Ebrahim, S. (2003). ‘Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease? Int. J. Epidemiol. 32, 1–22. doi: 10.1093/ije/dyg070

PubMed Abstract | CrossRef Full Text | Google Scholar

Speed, M. S., Jefsen, O. H., Borglum, A. D., Speed, D., and Ostergaard, S. D. (2019). Investigating the association between body fat and depression via Mendelian randomization. Transl. Psychiatry 9:184.

Google Scholar

Sterk, P. J., and Lutter, R. (2014). Asthma phenotyping: TH2-high, TH2-low, and beyond. J. Allergy Clin. Immunol. 133, 395–396. doi: 10.1016/j.jaci.2013.10.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Subbarao, P., Mandhane, P. J., and Sears, M. R. (2009). Asthma: epidemiology, etiology and risk factors. Can. Med. Assoc. J. 181, E181–E190.

Google Scholar

Sun, Y. Q., Brumpton, B., Langhammer, A., Chen, Y., Kvaloy, K., and Mai, X. M. (2020). Adiposity and asthma in adults: a bidirectional Mendelian randomisation analysis of The HUNT Study. Thorax 75, 202–208. doi: 10.1136/thoraxjnl-2019-213678

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, B. W., Yuan, S., Xiong, Y., He, Q. Q., and Larsson, S. C. (2020). Major depressive disorder and cardiometabolic diseases: a bidirectional Mendelian randomisation study. Diabetologia 63, 1305–1311. doi: 10.1007/s00125-020-05131-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Taylor, B., Mannino, D., Brown, C., Crocker, D., Twum-Baah, N., and Holguin, F. (2008). Body mass index and asthma severity in the National Asthma Survey. Thorax 63, 14–20. doi: 10.1136/thx.2007.082784

PubMed Abstract | CrossRef Full Text | Google Scholar

UK10K Consortium, Walter, K., Min, J. L., Huang, J., Crooks, L., Memari, Y., et al. (2015). The UK10K project identifies rare variants in health and disease. Nature 526, 82–90. doi: 10.1038/nature14962

PubMed Abstract | CrossRef Full Text | Google Scholar

Vlashki, E., Cholakovska, V. C., Kimovska, M., Seckova, L., Ristevska, T., and Lawson, J. (2018). The association of physical activity and sedentary lifestyle with asthma in childhood. Eur. Respir. J. 52:A4603.

Google Scholar

Vogelezang, S., Bradfield, J. P., Ahluwalia, T. S., Curtin, J. A., Lakka, T. A., Grarup, N., et al. (2020). Novel loci for childhood body mass index and shared heritability with adult cardiometabolic traits. PLoS Genet. 16:e1008718. doi: 10.1371/journal.pgen.1008718

PubMed Abstract | CrossRef Full Text | Google Scholar

Wade, K. H., Carslake, D., Sattar, N., Smith, G. D., and Timpson, N. J. (2019). Obesity BMI and mortality in UK biobank: revised estimates using mendelian randomization. Obesity 27, 349–349. doi: 10.1002/oby.22397

PubMed Abstract | CrossRef Full Text | Google Scholar

Watanabe, K., Taskesen, E., van Bochoven, A., and Posthuma, D. (2017). Functional mapping and annotation of genetic associations with FUMA. Nat. Commun. 8:1826.

Google Scholar

Woodruff, P. G., Modrek, B., Choy, D. F., Jia, G. Q., Abbas, A. R., Ellwanger, A., et al. (2009). T-helper Type 2-driven Inflammation defines major subphenotypes of asthma. Am. J. Resp. Crit. Care 180, 388–395. doi: 10.1164/rccm.200903-0392oc

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, S. J., Gilliland, F. D., and Conti, D. V. (2019). Elucidation of causal direction between asthma and obesity: a bi-directional Mendelian randomization study. Int. J. Epidemiol. 48, 899–907. doi: 10.1093/ije/dyz070

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Bakshi, A., Zhu, Z., Hemani, G., Vinkhuyzen, A. A., Nolte, I. M., et al. (2015). Genome-wide genetic homogeneity between sexes and populations for human height and body mass index. Hum. Mol. Genet. 24, 7445–7449. doi: 10.1093/hmg/ddv443

PubMed Abstract | CrossRef Full Text | Google Scholar

Zein, J. G., and Erzurum, S. C. (2015). Asthma is different in women. Curr. Allergy Asthma Rep. 15:28.

Google Scholar

Zhang, L., Yin, Y., Zhang, H., Zhong, W., and Zhang, J. (2017). Association of asthma diagnosis with leptin and adiponectin: a systematic review and meta-analysis. J. Invest. Med. 65, 57–64. doi: 10.1136/jim-2016-000127

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, H., Wu, D., Wu, X., Zhang, X., Zhou, Q., Luo, Y., et al. (2018). Leptin promotes allergic airway inflammation through targeting the unfolded protein response pathway. Sci. Rep. 8:8905.

Google Scholar

Zhu, Z. H., Zheng, Z. L., Zhang, F. T., Wu, Y., Trzaskowski, M., Maier, R., et al. (2018). Causal associations between risk factors and common diseases inferred from GWAS summary data. Nat. Commun. 9:224.

Google Scholar

Zhu, Z. Z., Guo, Y. J., Shi, H., Liu, C. L., Panganiban, R. A., Chung, W., et al. (2020). Shared genetic and experimental links between obesity-related traits and asthma subtypes in UK Biobank. J Allergy Clin. Immun. 145, 537–549. doi: 10.1016/j.jaci.2019.09.035

PubMed Abstract | CrossRef Full Text | Google Scholar