Peroxisome proliferator-activated receptor gamma gene variants modify human airway and systemic responses to indoor dibutyl phthalate exposure
Respiratory Research volume 23, Article number: 248 (2022)
Single nucleotide polymorphisms (SNPs) of peroxisome proliferator-activated receptor gamma (PPAR-γ; gene: PPARG) and oxidative stress genes are associated with asthma risk. However, whether such variants modulate responses to dibutyl phthalate (DBP), a common plasticizer associated with increased asthma development, remains unknown. The purpose of this study is to investigate how SNPs in PPARG and oxidative stress genes, as represented by two separate genetic risk scores, modify the impact of DBP exposure on lung function and the airway and systemic response after an inhaled allergen challenge.
We conducted a double-blinded human crossover study with sixteen allergen-sensitized participants exposed for three hours to DBP and control air on distinct occasions separated by a 4-week washout. Each exposure was followed by an allergen inhalation challenge; subsequently, lung function was measured, and blood and bronchoalveolar lavage (BAL) were collected and analyzed for cell counts and allergen-specific immunoglobulin E (IgE). Genetic risk scores for PPAR-γ (P-GRS; weighted sum of PPARG SNPs rs10865710, rs709158, and rs3856806) and oxidative stress (OS-GRS; unweighted sum of 16 SNPs across multiple genes) were developed, and their ability to modify DBP effects were assessed using linear mixed-effects models.
P-GRS and OS-GRS modified DBP effects on allergen-specific IgE in blood at 20 h (interaction effect [95% CI]: 1.43 [1.13 to 1.80], p = 0.005) and 3 h (0.99 [0.98 to 1], p = 0.03), respectively. P-GRS also modified DBP effects on Th2 cells in blood at 3 h (− 25.2 [− 47.7 to − 2.70], p = 0.03) and 20 h (− 39.1 [− 57.9 to − 20.3], p = 0.0005), and Th2 cells in BAL at 24 h (− 4.99 [− 8.97 to − 1.01], p = 0.02). An increasing P-GRS associated with reduced DBP effect on Th2 cells. Neither GRS significantly modified DBP effects on lung function parameters.
PPAR-γ variants modulated several airway and systemic immune responses to the ubiquitous chemical plasticizer DBP. Our results suggest that PPAR-γ variants may play a greater role than those in oxidative stress-related genes in airway allergic responses to DBP.
Trial registration: This study reports results from The Phthalate-Allergen Immune Response Study that was registered on ClinicalTrials.gov with identification NCT02688478.
Peroxisome proliferator-activated receptors (PPAR) are ligand-activated transcription factors expressed in the airway epithelium, airway smooth muscle cells, alveolar macrophages, T lymphocytes, and eosinophils [1,2,3]. Notably, PPAR-γ (NR1C3) regulates systemic and airway allergic immune responses relevant in asthma [4, 5].
Genetic variants or single nucleotide polymorphisms (SNPs) in PPARG and oxidative stress-related genes have been associated with increased risk of asthma development and exacerbations [6,7,8,9,10,11], highlighting gene-environment interactions in asthma pathophysiology. In one study of Koreans, the PPARG polymorphism rs3856806 was associated with asthma development . However, Palmer et al. reported that rare alleles of PPARG (rs1805192 and rs3856806) were associated with a reduced number of asthma exacerbations in a Caucasian population . In a case–control study by Li et al., asthma risk was influenced by the rs1805192 and rs10865710 PPARG polymorphisms in a Chinese population . Variants of the oxidative stress-related genes are associated with phthalate urine concentrations and modified the association of phthalate exposure with asthma in children . In a group of Korean elders, variants of certain oxidative stress-related genes influenced the association between urinary phthalate levels and decreasing lung function .
Phthalates are suspected to interact with PPAR-γ receptors . Phthalates are synthetic diesters of phthalic acid (1,2-benzenedicarboxylic acid) used to enhance properties of plastics such as flexibility, transparency, and durability . Humans are exposed to these environmental contaminants through a wide variety of common consumer products, such as household furnishings, food packaging, and cosmetics [15, 16]. Phthalate exposure has been associated with allergic diseases in epidemiological studies . In animal models, inhalation toxicity and oral intake of phthalates have demonstrated increased airway inflammation following exposures . In vitro studies to date indicate that exposure to environmentally relevant concentrations of certain phthalates can increase the release of inflammatory mediators from macrophages and lymphocytes, decrease phagocytic activity, and modulate differentiation of dendritic cells [17, 18].
We have previously published results from this randomized, crossover human exposure study which demonstrated that acute indoor exposure to dibutyl phthalate (DBP) enhanced the allergen-associated drop in forced expiratory volume in 1 s (FEV1) . DBP and allergen co-exposure significantly increased the percentage of M2 macrophages and macrophage activation markers in BAL, but had modest effects on airway immune mediators .
This novel study explores whether the effects of DBP on lung function, immune cells, and allergen-specific IgE are associated with SNPs related to PPARG (P-GRS) or oxidative stress-related genes (OS-GRS), as represented by genetic risk scores. We hypothesized that both the P-GRS and OS-GRS modulate the immune response to DBP, with higher scores associated with larger allergen-driven immune responses.
The study was registered on ClinicalTrials.gov with identification NCT02688478 and was approved by the research ethics board at the University of British Columbia (H14-01119) and the Norwegian Regional Committees for Medical and Health Research Ethics (2014/1217).
Participants and study design
Sixteen allergen-sensitized, non-smoking participants between 19–49 years of age were enrolled in this double-blinded, order-randomized, crossover study . The individuals were either healthy or had mild asthma (FEV1 greater than 70%); and were allergen-sensitized to grass, birch, or house dust mite. Airway hyperresponsiveness (AHR) status was categorized as either hyperresponsive (provocative concentration of methacholine resulting in 20% drop in FEV1 [PC20] ≤ 16 mg/ml) or normally responsive (PC20 > 16 mg/ml). Participants with asthma taking oral or inhaled corticosteroids were excluded. All participants discontinued their bronchodilator, antihistamine, and nonsteroidal anti-inflammatory medications 7 days prior to each exposure and during exposure days, except for salbutamol (200 µg was given routinely after methacholine challenge [by protocol, 20 h after exposure]).
Participants were randomized to the order of 3-h indoor air exposures to DBP (nominal concentration of 150 µg/m3) and control air (CA) at the Air Pollution Exposure Laboratory, University of British Columbia. Exposures were separated by at least 4 weeks for washout. Full details of the exposure methods were previously reported .
Outcome measures and analysis
The study design and data collection time are outlined in Fig. 1. Data were collected after each exposure (DBP and CA) in each crossover study arm. Blood samples were collected at − 4 h (baseline), 3 and 20 h after the allergen challenge post-DBP or CA exposure, while BAL was collected at 24 h after the allergen challenge post-DBP or CA exposure via bronchoscopy. These samples were analyzed for allergen-specific IgE, and blood and airway immune cells. The allergen-specific IgE level in blood and BAL was determined by ImmunoCAP Specific IgE analysis on a Phadia 250 system (Thermo Fisher Scientific). Immunophenotyping was performed by flow cytometry with analysis completed in FCS Express (v6.04.0034; De Novo Software). Further details are provided in the Additional file 1. The DBP effects on immune cells, which reflect primary endpoints rather than gene-exposure analyses, were previously reported .
FEV1 was measured at 3 and 20 h after exposure; while methacholine challenge and fractional exhaled nitric oxide (FeNO) were measured at 20 h after exposure. Delta values were calculated for each FEV1 timepoint and FeNO by subtracting the baseline values measured at -4 h. Allergen area under the curve (AUC) was determined as the percent fall in FEV1 spanning baseline to 3 h after inhaled allergen challenge.
Genotyping was performed on the DNA extracted from peripheral blood mononuclear cells (PBMC). Details of the genotyping are presented in the Additional file 1. The DNA sequences used to perform targeted pyrosequencing are listed in Additional file 1: Table S1.
Genetic risk scores
The GRS was developed from 3 SNPs of PPARG (P-GRS) that have been associated with an increased asthma risk and have minor allele frequency greater than 1% [6, 8]. The details, chromosome location, and risk alleles of the SNPs are illustrated in Table 1. The P-GRS was weighted as the sum of the effect estimates of the 3 SNPs: rs10865710, rs709158, and rs3856806 . The genotypes and corresponding weights for each SNP of the PPAR-γ GRS (P-GRS) are displayed in Additional file 1: Tables S1, S2. Genotyping was also performed for a 4th SNP of PPARG, rs1805192. This SNP did not vary among the participants, as shown in Additional file 1: Table S3. This SNP was therefore not included in the GRS.
Our analyses also included an oxidative stress genetic risk score (OS-GRS) based on 16 allelic variants that have been suggested to modulate the response to air pollution and tobacco smoke via oxidative stress pathway [20,21,22,23,24,25,26,27], Additional file 1: Table S4. Selected genes were derived from a number of individual studies with heterogeneous endpoints, thus, no adequate external weights of the relative contribution of each variant were available for the development of weighted risk scores in OS-GRS .
Statistical tests were performed using linear mixed-effects (LME) models (nlme package version 3.1-151) in R, version 4.0.3 (R Foundation for Statistical Computing) implemented in R Studio, version 1.3.1093 (RStudio). Confidence intervals (CI) were computed using the gmodels package (version 2.18.1) in RStudio. Interaction plots for study endpoints were generated using R packages: LME4 (version 1.1-27) and interplot (version 0.2.3). All LME analyses for analytes in blood were performed for delta values, calculated for 3 h and 20 h timepoints relative to the baseline (− 4 h), for each individual and each exposure condition. When necessary, data were log10-transformed to normalize the data distribution. Any effect estimates and associated confidence intervals are based on log10-transformed data.
In our primary LME model for allergen-specific IgE level, exposure (DBP relative to CA) was the fixed effect, and participant ID was set as a random effect. To assess whether either of the two GRS modified the effect of DBP on the study outcomes, a model with the interaction term exposure-by-GRS (exposure*GRS) as a fixed effect was evaluated. Participant ID was set as a random effect in these models. P-values < 0.05 were considered statistically significant.
Ten females and 6 males enrolled in the PAIR study were included in these analyses. Eleven out of the 16 participants underwent research bronchoscopies. Baseline participant profiles at screening, FEV1 measurements, and P-GRS are shown in Table 2. Eight participants were categorized as hyperresponsive AHR status based on baseline methacholine PC20. The mean P-GRS among participants was − 0.12 ± 0.15, with a more negative number indicating a lower risk score and a more positive number (maximum is 0) indicating a higher risk score, i.e. greater risk of having asthma.
The 3-h exposure to DBP did not significantly change allergen-specific IgE levels compared to CA in blood (3 h and 20 h) or BAL (24 h) (Fig. 2A–C). However, there was a significant exposure-by-P-GRS interaction effect at 20 h (interaction effect [95% CI]: 1.43 [1.13 to 1.80], p = 0.005). The positive slope of the interaction line indicates that increasing P-GRS was associated with increasing DBP effect on allergen-specific IgE in blood (Fig. 2E). Participants with a higher P-GRS experienced a larger DBP effect on IgE at 20 h compared to baseline. There was no significant interaction of P-GRS at 3 h and 24 h in blood and BAL, respectively (Fig. 2D, F).
DBP exposure did not significantly affect Th2 lymphocytes percentages compared to CA in the blood (3 h and 20 h) or BAL (24 h) (Fig. 3A–C). The DBP effect on % Th2 lymphocytes was modified by P-GRS, such that those participants with higher P-GRS were associated with decreasing DBP effect on % Th2 (Fig. 3D–F). This was observed for % Th2 in blood at 3 and 20 h (interaction effect [95% CI]: − 25.2 [− 47.7 to − 2.70], p = 0.03; and − 39.1 [− 57.9 to − 20.3], p = 0.0005), and in BAL at 24 h (− 4.99 [− 8.97 to − 1.01], p = 0.02). Participants with a higher P-GRS experienced a smaller DBP effect on Th2 lymphocytes at 3 h, 20 h, and 24 h compared to baseline. There was however no significant exposure-by-P-GRS interaction for other immune cells in blood or BAL, as shown in Additional file 1: Table S5.
Lung function outcomes included allergen AUC (from baseline to 3 h post allergen challenge), FEV1 and FeNO (∆20 h), and methacholine challenge (20 h). The DBP effect on these outcomes was not significantly modified by P-GRS (Fig. 4).
The unweighted OS-GRS composed of 16 genes for each participant are listed in Additional file 1: Table S6. The mean OS-GRS was 9.7 ± 2.6, with a more positive number indicating a higher risk score (theoretical maximum is 32). There was a significant exposure-by-OS-GRS interaction (interaction effect [95% CI]: 0.99 [0.98 to 1], p = 0.03) such that increasing OS-GRS was associated with decreasing DBP effect on blood allergen-specific IgE at 3 h (Additional file 1: Fig. S1A). There was no significant interaction of OS-GRS for DBP effect on allergen-specific IgE and % Th2 lymphocytes at other time points (Additional file 1: Fig. S1B–F). OS-GRS did not modify the DBP effect on lung function parameters (Additional file 1: Fig. S2).
This study demonstrated that a genetic risk score based on three relevant PPARG polymorphisms may modify the DBP effect on two important mediators of allergic responses, allergen-specific IgE levels and the proportion of Th2 cells amongst lymphocytes. None of these endpoints were affected by the DBP exposure alone. In contrast to P-GRS, the impact of OS-GRS was limited. Thus, our findings suggest that PPAR-γ may play a role in the airway allergic effects in response to DBP, while oxidative stress does not appear to be substantially involved.
IgE plays an important role in the acute airway allergic response to allergens, by binding to mast cells and basophils to release histamine, neutral proteases, and chemotactic factors . IgE and Th2 lymphocytes are relevant in the DBP-related airway inflammation, and several phthalates have been shown to increase the release of Th2-associated mediators and modulates Th1/Th2 balance . We observed that individuals with a higher P-GRS augmented the DBP effect on blood allergen-specific IgE at 20 h. The interaction of P-GRS with DBP-induced effects on IgE level was opposite to that of Th2 lymphocytes. Moreover, the DBP effect on B cells, which are the main cells that produce IgE, were not significantly modified by P-GRS. In a prior study where a bronchoscopic segmental allergen challenge was performed in atopic asthmatic subjects, there was a significant increase in allergen-specific IgE in the bronchoalveolar lavage fluid 24 h after the challenge . In contrast, the observed rises in specific IgE in the serum was much smaller, suggesting a predominantly local allergen response that is independent of systemic alterations . Furthermore, the production of IgE after an inhaled allergen challenge appears to be mediated through local expansion of IgE+ B cells in the airways rather than in the blood, with the latter being a more delayed response . Overall, the interaction of P-GRS with DBP-related Th2 lymphocyte response in the blood and airways was more uniform and substantial compared to the observation regarding blood allergen-specific IgE alone.
Given these observations, the overall trend is a reduction in DBP-induced effects for increasing P-GRS, or higher asthma risk . A possible interpretation is that individuals with a low P-GRS are more susceptible to DBP-induced effects in the airways. This is somewhat analogous with the main result in the previous PAIR study which demonstrated that airway responsiveness increased by 48% after DBP exposure specifically in participants without baseline hyperresponsiveness, exacerbating the mean response to levels typical of those with airway hyperresponsiveness .
In our report of primary endpoints in this study, we showed that DBP exposure augmented the allergen-induced recruitment of total macrophages in the airways and there was skewing towards M2 macrophage phenotype in the BAL . However, we did not find an exposure-by-P-GRS interaction on the DBP effect for macrophages in the current analysis, suggesting that augmented macrophages recruitment in the lung may not be impacted by the genetic variation in PPAR-γ. Furthermore, Maestre-Batlle et al. reported that exposure to DBP resulted in significant augmentation of allergen-induced AHR as measured by greater decrease methacholine PC20 and FEV1 change in participants without AHR . But, the present study did not find a significant exposure-by-P-GRS interaction on lung function. Impact on the lung function likely reflects cumulative immune responses in the airways which modulate airway smooth muscle function rather than reflecting modulation of Th2 cells or IgE production. Our data suggests that genetic variation in PPAR-γ only impacted IgE and Th2 cells.
Asthma is a complex disease with numerous genetic determinants. The weighted GRS is an intuitive approach for combining genetic data into a single measure of susceptibility to a disease or physiological response . The weighted score was preferred over a simple-count genetic score equal to the sum of its risk alleles because the SNPs used in our study were derived from a prior observational study with reported odds ratio and effect estimates, allowing us to approximate the relative weights of each SNP . However, we acknowledge that our estimate of genetic risk is still crude as it contains only three susceptibility polymorphisms that are associated with asthma. In our study, two of the three SNPs, rs10865710 and rs709158 are intron variants of the PPAR-γ gene [33, 34]. We found that polymorphism in rs10865710 alone was sufficient to show effect modification on Th2 lymphocytes and allergen-specific IgE levels in blood and BAL (Additional file 1: Fig. S1). We also note a bimodal distribution of P-GRS. This is due to the weighing of different SNPs based on asthma odds ratio reported in previous research . Specifically, rs10865710 G/G genotype conferred the greatest reduction in risk (odd ratio (OR) = 0.53) for asthma compared to rs709158 GG (OR = 0.91) and rs3856806 (OR = 0.97). Since those individuals with rs10865710 G/G genotype (n = 4) had the smallest OR and consequently, a much lower P-GRS than the other genotypes, this resulted in a bimodal distribution of P-GRS. It is unknown if the selected SNPs in this study, though associated with an important clinical outcome, are involved with the interaction of phthalates and the ligand-binding sites on PPAR-γ receptors.
Although previous studies suggest that even a single SNP can mediate a profound effect on the regulation of phenotype , the optimal number of SNPs to utilize in creating a weighted GRS is still under debate. The strategy used to construct a GRS, whether based on the protein biological function or structure, is also contested. In our case, another approach to create the GRS would be to select SNPs specifically involved in the receptor ligand-binding domain through which phthalates are thought to interact. Thus, our study only suggests a plausible, albeit indirect relationship between these polymorphisms and PPAR-γ function. Further studies with SNPs relevant to the ligand-binding domain should be performed to assess whether PPAR-γ receptors play a role in the molecular cascade initiated by DBP in human airway cells. Larger studies are required for validating the clinical relevance of PPAR-γ polymorphisms.
The current findings of effect modification of PPAR-γ polymorphisms on DBP induced effects should be used to focus future experimental efforts. The modest participant sample size and the weighted P-GRS may have limited variations in the GRS among the participants as evident in the P-GRS distribution. Moreover, the BAL and blood samples collected in this study only captured cellular and clinical outcomes at defined time points within 24 h post-allergen and DBP exposure, and the long-term risk of allergic and inflammatory airway responses is unknown. Finally, OS-GRS had a limited impact on DBP-induced effects, which may reflect differences in study design and endpoints, and phthalate type and concentration compared to other studies.
Several genetic variations in PPARG, as represented by a genetic risk score, modified the impact of dibutyl phthalate exposure on allergen-specific IgE and Th2 lymphocytes after allergen challenge in a controlled human exposure study. Our study focused specifically on DBP given its real-world propensity for exposure by inhalation, and our findings may motivate further study for potential application to personal and public health practice and policy.
Availability of data and materials
All relevant data generated or analysed during this study are included in this published article and its supplementary information file. Additional data may be available from the corresponding author on reasonable request.
Allergen area under the curve
- FEV1 :
Forced expiratory volume in 1 s
Fractional exhaled nitric oxide
Oxidative stress genetic risk score
Peripheral blood mononuclear cells
Peroxisome proliferator-activated receptor gamma
PPAR-γ genetic risk score
Single nucleotide polymorphism
Kocbach Bølling A, Holme JA, Bornehag CG, Nygaard UC, Bertelsen RJ, Nånberg E, et al. Pulmonary phthalate exposure and asthma—is PPAR a plausible mechanistic link? EXCLI J. 2013;12:733–59.
Rehan VK, Asotra K, Torday JS. The effects of smoking on the developing lung: insights from a biologic model for lung development, homeostasis, and repair. Lung. 2009;187(5):281–9.
Becker J, Delayre-Orthez C, Frossard N, Pons F. Regulation of inflammation by PPARs: a future approach to treat lung inflammatory diseases? Fundam Clin Pharmacol. 2006;20(5):429–47.
Benayoun L, Letuve S, Druilhe A, Boczkowski J, Dombret MC, Mechighel P, et al. Regulation of peroxisome proliferator-activated receptor gamma expression in human asthmatic airways: relationship with proliferation, apoptosis, and airway remodeling. Am J Respir Crit Care Med. 2001;164(8 Pt 1):1487–94.
Yessoufou A, Wahli W. Multifaceted roles of peroxisome proliferator-activated receptors (PPARs) at the cellular and whole organism levels. Swiss Med Wkly. 2010;140: w13071.
Li W, Dai W, Sun J, Zhang W, Jiang Y, Ma C, et al. Association of peroxisome proliferator-activated receptor-gamma gene polymorphisms and gene-gene interaction with asthma risk in a Chinese adults population. Int J Clin Exp Med. 2015;8(10):19346–52.
Oh SH, Park SM, Lee YH, Cha JY, Lee JY, Shin EK, et al. Association of peroxisome proliferator-activated receptor-gamma gene polymorphisms with the development of asthma. Respir Med. 2009;103(7):1020–4.
Palmer CNA, Doney ASF, Ismail T, Lee SP, Murrie I, Macgregor DF, et al. PPARG locus haplotype variation and exacerbations in asthma. Clin Pharmacol Ther. 2007;81(5):713–8.
You H, Chen S, Mao L, Li B, Yuan Y, Li R, et al. The adjuvant effect induced by di-(2-ethylhexyl) phthalate (DEHP) is mediated through oxidative stress in a mouse model of asthma. Food Chem Toxicol Int J Publ Br Ind Biol Res Assoc. 2014;71:272–81.
Zhou S, Han M, Ren Y, Yang X, Duan L, Zeng Y, et al. Dibutyl phthalate aggravated asthma-like symptoms through oxidative stress and increasing calcitonin gene-related peptide release. Ecotoxicol Environ Saf. 2020;199: 110740.
Franken C, Lambrechts N, Govarts E, Koppen G, Den Hond E, Ooms D, et al. Phthalate-induced oxidative stress and association with asthma-related airway inflammation in adolescents. Int J Hyg Environ Health. 2017;220(2 Pt B):468–77.
Wang IJ, Karmaus W. Oxidative stress-related genetic variants may modify associations of phthalate exposures with asthma. Int J Environ Res Public Health. 2017;14(2):162.
Park HY, Kim JH, Lim YH, Bae S, Hong YC. Influence of genetic polymorphisms on the association between phthalate exposure and pulmonary function in the elderly. Environ Res. 2013;122:18–24.
North ML, Takaro TK, Diamond ML, Ellis AK. Effects of phthalates on the development and expression of allergic disease and asthma. Ann Allergy Asthma Immunol. 2014;112(6):496–502.
Schettler T. Human exposure to phthalates via consumer products. Int J Androl. 2006;29(1):134–9.
Bølling AK, Sripada K, Becher R, Bekö G. Phthalate exposure and allergic diseases: review of epidemiological and experimental evidence. Environ Int. 2020;139: 105706.
Kuo CH, Hsieh CC, Kuo HF, Huang MY, Yang SN, Chen LC, et al. Phthalates suppress type I interferon in human plasmacytoid dendritic cells via epigenetic regulation. Allergy. 2013;68(7):870–9.
Lee MH, Park J, Chung SW, Kang BY, Kim SH, Kim TS. Enhancement of interleukin-4 production in activated CD4+ T cells by diphthalate plasticizers via increased NF-AT binding activity. Int Arch Allergy Immunol. 2004;134(3):213–22.
Maestre-Batlle D, Huff RD, Schwartz C, Alexis NE, Tebbutt SJ, Turvey S, et al. Dibutyl phthalate augments allergen-induced lung function decline and alters human airway immunology. A randomized crossover study. Am J Respir Crit Care Med. 2020;202(5):672–80.
Wooding DJ, Ryu MH, Hüls A, Lee AD, Lin DTS, Rider CF, et al. Particle depletion does not remediate acute effects of traffic-related air pollution and allergen. A randomized, double-blind crossover study. Am J Respir Crit Care Med. 2019;200(5):565–74.
Bowatte G, Erbas B, Lodge CJ, Knibbs LD, Gurrin LC, Marks GB, et al. Traffic-related air pollution exposure over a 5-year period is associated with increased risk of asthma and poor lung function in middle age. Eur Respir J. 2017;50(4):1602357.
Wenzlaff AS, Cote ML, Bock CH, Land SJ, Schwartz AG. GSTM1, GSTT1 and GSTP1 polymorphisms, environmental tobacco smoke exposure and risk of lung cancer among never smokers: a population-based study. Carcinogenesis. 2005;26(2):395–401.
Masuko H, Sakamoto T, Kaneko Y, Iijima H, Naito T, Noguchi E, et al. An interaction between Nrf2 polymorphisms and smoking status affects annual decline in FEV1: a longitudinal retrospective cohort study. BMC Med Genet. 2011;20(12):97.
Aoki T, Hirota T, Tamari M, Ichikawa K, Takeda K, Arinami T, et al. An association between asthma and TNF-308G/A polymorphism: meta-analysis. J Hum Genet. 2006;51(8):677–85.
Mordukhovich I, Wilker E, Suh H, Wright R, Sparrow D, Vokonas PS, et al. Black carbon exposure, oxidative stress genes, and blood pressure in a repeated-measures study. Environ Health Perspect. 2009;117(11):1767–72.
Canova C, Dunster C, Kelly FJ, Minelli C, Shah PL, Caneja C, et al. PM10-induced hospital admissions for asthma and chronic obstructive pulmonary disease: the modifying effect of individual characteristics. Epidemiol Camb Mass. 2012;23(4):607–15.
MacIntyre EA, Brauer M, Melén E, Bauer CP, Bauer M, Berdel D, et al. GSTP1 and TNF Gene variants and associations between air pollution and incident childhood asthma: the traffic, asthma and genetics (TAG) study. Environ Health Perspect. 2014;122(4):418–24.
Hüls A, Krämer U, Carlsten C, Schikowski T, Ickstadt K, Schwender H. Comparison of weighting approaches for genetic risk scores in gene-environment interaction studies. BMC Genet. 2017;18(1):115.
Gauvreau GM, El-Gammal AI, O’Byrne PM. Allergen-induced airway responses. Eur Respir J. 2015;46(3):819–31.
Wilson DR, Merrett TG, Varga EM, Smurthwaite L, Gould HJ, Kemp M, et al. Increases in allergen-specific IgE in BAL after segmental allergen challenge in atopic asthmatics. Am J Respir Crit Care Med. 2002;165(1):22–6.
Oliveria JP, Salter BM, MacLean J, Kotwal S, Smith A, Harris JM, et al. Increased IgE+ B cells in sputum, but not blood, bone marrow, or tonsils, after inhaled allergen challenge in subjects with asthma. Am J Respir Crit Care Med. 2017;196(1):107–9.
Igo RP, Kinzy TG, Bailey JNC. Genetic risk scores. Curr Protoc Hum Genet. 2019;104(1): e95.
Li T, Xu X, Xu Y, Jin P, Chen J, Shi Y, et al. PPARG polymorphisms are associated with unexplained mild vision loss in patients with type 2 diabetes mellitus. J Ophthalmol. 2019;12(2019):1–7.
Lu H, Wen D, Sun J, Zeng L, Du J, Du D, et al. Enhancer polymorphism rs10865710 associated with traumatic sepsis is a regulator of PPARG gene expression. Crit Care Lond Engl. 2019;23(1):430.
The authors declare that they have no financial or non-financial disclosures. The authors thank Agnes Yuen and Ryan Huff (Department of Medicine, University of British Columbia, Vancouver, British Columbia, Canada) for their assistance and support of this study.
This study was supported by AllerGen NCE, Legacy for Airway Health, the Norwegian Research Council (research grant number 228129) and The Norwegian Asthma and Allergy Association NAAF.
Ethics approval and consent to participate
This study was approved by the research ethics board at the University of British Columbia (H14-01119) and the Norwegian Regional Committees for Medical and Health Research Ethics (2014/1217). All participants recruited to The Phthalate-Allergen Immune Response Study had consented to participate.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
DNA sequences for PPARG used to perform targeted pyrosequencing. Table S2. The genotypes and GRS are illustrated for rs709158, rs3856806, and rs10865710. Table S3. The weights of each SNP in the construction of P-GRS. Table S4. Oxidative stress-related genes included in the associated OS-GRS and definition of their risk alleles. Table S5. Exposure-by-P-GRS interaction on DBP effects on (A) blood immune cells at 3h, (B) blood immune cells at 20h, and (C) BAL immune cells at 24h. Table S6. Participant characteristics and oxidative stress genetic risk score (OS-GRS). Figure S1. Effect of OS-GRS interaction on DBP effect on IgE and Th2 lymphocytes. Blood allergen-specific IgE at 3 and 20 h (A-B), and BAL allergen-specific IgE at 24 h (C); blood % Th2 lymphocytes at 3 and 20 h (D-E), and BAL % Th2 at 24 h (F) after the allergen challenge post-DBP. The histogram underlying the plot illustrates the distribution of the participants’ OS-GRS. Blood and BAL samples were collected for 16 and 11 participants, respectively. The histogram underlying each plot illustrates the distribution of the participants’ OS-GRS. Effect values represent the DBP-attributable difference between the measurement at the given time point and the baseline. The solid line represents the interaction effect, with the p value corresponding to the interaction significance. The shaded region represents 95% confidence intervals. Figure S2. Exposure-by-OS-GRS interaction on DBP effect on lung function outcomes. Lung function outcomes included: (A) FEV1 (% predicted), ∆20 h; (B) Methacholine challenge (LogPC20), 20 h; (C) FeNO, ∆20 h; (D) Allergen AUC, 3 h. The listed time points refer to the time after the allergen challenge post-DBP or CA exposure. FEV1 = forced expiratory volume in one second; PC20 = concentration leading to 20% fall in FEV1; FeNO = fractional exhaled nitric oxide; AUC = area under curve, measured as percent decline in FEV1 spanning baseline to 3 h after allergen challenge. The histogram underlying each plot illustrates the distribution of the participants’ OS-GRS. Effect values represent the DBP-attributable difference between the measurement at the given time point and the baseline. The solid line represents the interaction effect, with the p value corresponding to the interaction significance. The shaded region represents 95% confidence intervals. Figure S3. PPAR-γ single nucleotide polymorphism rs10865710. Effect modification by genotype G/C or C/C on (A) Blood IgE at 3 and 20 hours post-exposure to CA with allergen challenge (Air+Ag) or DBP with allergen challenge (DBP+Ag), (B) Blood Th2 lymphocytes at 3 and 20 hours post-exposure to CA (Air) or DBP, and (C) BAL Th2 lymphocytes at 24 hours post-exposure to CA (Air) or DBP.
About this article
Cite this article
Leung, C., Ryu, M.H., Bølling, A.K. et al. Peroxisome proliferator-activated receptor gamma gene variants modify human airway and systemic responses to indoor dibutyl phthalate exposure. Respir Res 23, 248 (2022). https://doi.org/10.1186/s12931-022-02174-8