Due to recent changes in operational resources, EHP and JHP are currently unable to consider new manuscript submissions for peer review or publication.

Skip to main content
Skip to main content
Open access
Research
15 August 2014

Modification of the Association between PM10 and Lung Function Decline by Cadherin 13 Polymorphisms in the SAPALDIA Cohort: A Genome-Wide Interaction Analysis

Publication: Environmental Health Perspectives
Volume 123, Issue 1
Pages 72 - 79

Abstract

Background: Both air pollution and genetic variation have been shown to affect lung function. Their interaction has not been studied on a genome-wide scale to date.
Objectives: We aimed to identify, in an agnostic fashion, genes that modify the association between long-term air pollution exposure and annual lung function decline in an adult population-based sample.
Methods: A two-stage genome-wide interaction study was performed. The discovery (n = 763) and replication (n = 3,896) samples were derived from the multi-center SAPALDIA cohort (Swiss Cohort Study on Air Pollution and Lung Disease in Adults). Annual rate of decline in the forced mid-expiratory flow (FEF25–75%) was the main end point. Multivariate linear regression analyses were used to identify potential multiplicative interactions between genotypes and 11-year cumulative PM10 exposure.
Results: We identified a cluster of variants intronic to the CDH13 gene as the only locus with genome-wide significant interactions. The strongest interaction was observed for rs2325934 (p = 8.8 × 10–10). Replication of the interaction between this CDH13 variant and cumulative PM10 exposure on annual decline in FEF25–75% was successful (p = 0.008). The interaction was not sensitive to adjustment for smoking or body weight.
Conclusions: CDH13 is functionally linked to the adipokine adiponectin, an inflammatory regulator. Future studies need to confirm the interaction and assess how the result relates to previously observed interactions between air pollution and obesity on respiratory function.
Citation: Imboden M, Kumar A, Curjuric I, Adam M, Thun GA, Haun M, Tsai MY, Pons M, Bettschart R, Turk A, Rochat T, Künzli N, Schindler C, Kronenberg F, Probst-Hensch NM. 2015. Modification of the association between PM10 and lung function decline by cadherin 13 polymorphisms in the SAPALDIA cohort: a genome-wide interaction analysis. Environ Health Perspect 123:72–79; https://doi.org/10.1289/ehp.1307398

Introduction

Lung function is a complex phenotype influenced by lifestyle, environmental, and genetic factors. Inverse associations between chronic exposure to air pollution, such as particulate matter (PM), and respiratory function level as well as its decline have been reported in independent settings (Downs et al. 2007; Katanoda et al. 2011; Romieu et al. 2009; Schikowski et al. 2010). Air pollutants are thought to mediate their acute and chronic effects through an increase in oxidative stress, inflammation, and cytotoxicity (Andreau et al. 2012; Huang et al. 2012). However, mechanisms and differences in susceptibility remain poorly characterized (Brook et al. 2010). Only few candidate gene–air pollution interaction studies have been published. These reports also point to the oxidative and inflammatory effects of air pollution in mediating adverse respiratory health effects (Breton et al. 2011; Curjuric et al. 2012; Imboden et al. 2009; Melén et al. 2008; Romieu et al. 2006; Yang et al. 2005).
Genome-wide association studies (GWAS) on lung function were mostly cross-sectional in nature (Artigas et al. 2011; Hancock et al. 2010, 2012; Obeidat et al. 2011; Repapi et al. 2010; Wilk et al. 2009), and more recent reports have shown that the overlap in genetic determinants of the level of lung function and its decline is minimal (Hansel et al. 2013; Imboden et al. 2012). None of these lung function GWAS studies has considered ambient air pollution.
In the present study, we used a genome-wide interaction study (GWIS) approach to uncover novel genetic loci modifying the association between particulate matter exposure and 11-year lung function decline. We applied a two-stage approach with a discovery sample (n = 763) and a replication sample (n = 3,896). Both are subpopulations of the SAPALDIA cohort study (Swiss Cohort Study on Air Pollution And Lung Diseases In Adults) (Martin et al. 1997). This multi-center population-based cohort was specifically designed to investigate long-term effects of air pollution on respiratory health. We a priori chose forced mid-expiratory flow (FEF25–75%) as the dependant lung function phenotype because it was the outcome most strongly associated with ambient particulate matter air pollution exposure in SAPALDIA (Downs et al. 2007). As a proxy for long-term exposure to complex air pollution mixtures, we chose personal estimates of 11-year cumulative exposure to home outdoor PM mass with ≤ 10 μm in aerodynamic diameter (PM10) (Curjuric et al. 2012; Liu et al. 2007).

Methods

SAPALDIA cohort study. SAPALDIA was initiated in 1991. Participants, 18–60 years of age, were randomly selected from the population registries of eight geographic Swiss regions, with varying degrees of urbanization and different environmental and cultural characteristics. Participants of the baseline examination (n = 9,651) were invited in 2002 (n = 8,047) for a second examination. Ethical approval was obtained from the Swiss Academy of Medical Sciences and the Regional Ethics Committees; written informed consent was obtained from all participants before health examination and biological sample collection at each survey. Study design and data collected have been described elsewhere (Ackermann-Liebrich et al. 2005). Briefly, health examinations and standardized questionnaires focused on respiratory and cardiovascular health. Formal collection of fractioned blood and DNA samples was established at the follow-up survey.
Study population. Nonparticipation at follow-up (n = 1,604) and missing information on lung function phenotype data (n = 2,302), genotype data (n = 476), or covariates (n = 43) led to the exclusion of some SAPALDIA cohort participants from the present study. Participants reporting a history of asthma were excluded (n = 567) because of evidence of genetic heterogeneity of lung function decline in asthmatic and nonasthmatic subjects (Imboden et al. 2012). The final study population included participants with blood samples available for genetic testing and complete baseline and follow-up data on spirometry, smoking history, weight, weight change, height, PM10 exposure, and residential history, as well as statistical model covariates (n = 4,659) (Downs et al. 2007). The discovery sample with genome-wide data was a random sample of the nonasthmatic SAPALDIA study population (n = 763) (Moffatt et al. 2010). The replication sample consisted of the remainder of SAPALDIA participants with complete data (n = 3,896) and was subjected to targeted genotyping for replication of promising discovery interaction signals.
Phenotype and covariate assessment. For lung function assessment, identical spirometer devices (Sensormedics model 2200; Sensormedics, Yorba Linda, CA, USA) and protocols were used at both examinations (Ackermann-Liebrich et al. 2005). Comparability of devices was ascertained (Künzli et al. 2005). Each participant performed three to a maximum of eight forced expiratory lung function maneuvers to obtain a minimum of two acceptable forced expiratory flows, forced vital capacity (FVC), forced expiratory volume in the first second (FEV1) complying with American Thoracic Society (1995) criteria. Expiratory flow measures during the middle half of the FVC (FEF25–75%) were taken from the flow-volume curves with the highest sum of FVC and FEV1. Given evidence from the SAPALDIA cohort, we focused in the present study on the annual rate of decline in FEF25–75% as a sensitive marker of age-related decline because it was more strongly associated with the PM10 exposure than was FEV1 or FEV1/FVC decline (Curjuric et al. 2010; Downs et al. 2007; Imboden et al. 2009; Thun et al. 2012). Annual decline in FEF25–75% was calculated as the difference between follow-up and baseline measure in milliters per second, divided by length of follow-up in years. Accordingly, declines in FEV1, FEV1/FVC, and FEF25–75%/FVC were calculated for sensitivity analyses of the FEF25–75% GWIS top hits. Covariate information was assessed including a computer-assisted personal interview at baseline and follow-up examinations, including age, sex, current and past smoking status, and smoking history (number of cigarettes/day, years of smoking). The exposure to other inhaled pollutants such as environmental tobacco smoke or occupational exposure to dust and fumes, and respiratory symptoms were assessed with the same questions at both surveys. Participants who reported smoking < 20 packs of cigarettes and < 360 g of tobacco in their lifetime at both time points were defined as never-smokers. Cumulative cigarette exposure of participants was assessed by pack-years smoked before the first examination and pack-years smoked during follow-up. Height was measured (without shoes) at baseline and follow-up. Weight was self-reported at baseline and measured at follow-up (without shoes and coat). Weight change was calculated as weight at follow-up minus weight at baseline, with positive values reflecting weight gain during follow-up period.
Home outdoor PM10 exposure assessment. We used PM10 as the air pollution exposure measure. Air pollution exposure assessment, dispersion model validation, as well as attribution of individual 11-year cumulative PM10 exposure have been described elsewhere (Liu et al. 2007). Briefly, a hybrid exposure model incorporated geocoded data on seasonal, meteorological, and traffic, industrial, regional, and agricultural emission activities. Hourly concentrations of PM10 were calculated on a spatial resolution of 200 × 200 m grid cells over the follow-up period. Annual averages of the modeled PM10 concentrations were obtained for each grid cell. We estimated the cumulative PM10 exposure for study participants using their residential history, in geocoded data format, assigning annual PM10 exposure averages derived from the grid cells generated by the dispersion model, and adding up the averages over the 11-year follow-up period (Liu et al. 2007).
Genotyping, imputation, population stratification, and quality control. DNA extraction from EDTA-buffered whole blood has been previously described (Ackermann-Liebrich et al. 2005; Imboden et al. 2006). Genome-wide genotyping was obtained using the platform Illumina 610K quad Bead Chip. Discovery genotyping quality control, imputation, and correction for population stratification was applied as previously described (Moffatt et al. 2010). Briefly, genome-wide genotyping was centrally performed for the GABRIEL asthma Consortium at the Centre National de Génotypage (CNG, Evry, France). We satisfactorily genotyped 567,589 autosomal single nucleotide polymorphisms (SNPs) (mean call rate, 99.7%). We obtained 2,588,592 autosomal HapMap-based SNPs by cohort-specific imputation using the MACH v1.00 software and the HapMap2 release 22 CEU reference sample (Moffatt et al. 2010). Statistical power to detect gene–environment interaction is expected to be limited, so we excluded SNPs with minor allele frequencies < 5% to avoid inflation of false positive findings produced by rare genetic variants. Final number of SNPs used for interaction association testing was 2,198,793. To account for population stratification, we relied on inferred ancestry-informative principal components (Moffatt et al. 2010) that were previously carried out using EIGENSTRAT 2.0 software and the all ethnicity HapMap data, as well as additional European reference samples (Heath et al. 2008). Subjects of non-European descent were excluded based on the first and second principal components. Adjustment for population stratification in the linear regression analyses was done by incorporating the third and fourth principal component in the statistical model. Cryptic relatedness was detected based on identity-by-state analysis, and one participant per family cluster was retained in the study population.
Statistical analysis. Discovery sample. We performed agnostic GWIS analysis in the discovery sample using an additive genetic model, with a priori selection of potential confounders based on previous analyses of the association between air pollution and lung function decline (Curjuric et al. 2010; Downs et al. 2007; Imboden et al. 2009). We regressed each SNP with cumulative PM10 on FEF25–75% annual decline adjusting for study center, age, sex, height, never-smoking status, seasonal effects (sine and cosine function of day of examination), and population stratification factors. In addition, models were adjusted for weight at baseline, weight change during follow-up, and the multiplicative interaction between baseline weight and weight change, based on our recent analysis demonstrating an interaction between air pollution and obesity on lung function in the study population (Schikowski et al. 2013). Potential interaction effects between genotype and PM10 were captured by the inclusion of a multiplicative interaction term in the linear multivariate regression analyses. We used a joint test with two degrees of freedom to derive p-values for the joint effects of gene marginal and gene–environment interactions. This approach has been shown under a range of scenarios to have greater power for identifying novel genetic candidates than tests of the gene marginal effect or gene–environment interaction effects alone (Hancock et al. 2012; Kraft et al. 2007). We used the following terminology to report the results of the genetic effects related to the gene marginal (pmain), the gene–environment (pint), and the joint (pjoint) effects referring to their respective null hypothesis of gene marginal (βmain = 0), the gene–environment (βint = 0), and the joint (βmain = 0 and βint = 0) effects. We defined the genome-wide significance level at p < 5 × 10–8 using Bonferroni adjustment for one million independent tests. The lambda (λ) for the main GWIS—a metric for estimating genomic inflation of the observed associations—was calculated as the ratio of the observed versus expected median of the chi-square distribution with 2 degrees of freedom [ΣChi2(2df); medianobserved divided by 1.386 (medianexpected)]. GWIS sensitivity analyses were additionally performed on FEV1, FEV1/FVC, and FEF25–75%/FVC decline to determine the genome-wide ranking of the FEF25–75% top hits in the GWIS results for other lung function phenotypes.
Replication sample. The replication analyses was performed on two intronic CDH13 SNPs exhibiting the lowest or very low p-values in the discovery phase without being in high linkage disequilibrium with one other (top hit, rs232593, and rs17284098). Replication de novo genotyping of rs2325934 and rs17284098 was performed on a 7900HT Fast Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) by using 5´-nuclease allelic discrimination assays. A random sample of approximately 10% of all DNA samples was re-genotyped, and all genotypes were confirmed. The genotype call rate was > 99%. The same adjustments as for the discovery GWIS were used, except that adjustment for population stratification was not possible in the replication sample. We do not, however, expect associations to be confounded in the replication sample because adjustment for population stratification did not influence associations in the discovery sample (data not shown). Given the gene–environment interaction identified, we performed genotype stratified analyses in the combined sample (discovery and replication) as well as additional explorative analyses assessing the robustness of the observed PM10CDH13 interaction with a particular focus on smoking and on weight-related variables.
Post hoc analysis on CDH13. In a post hoc analysis, we first searched the dbGaP database (http://www.ncbi.nlm.nih.gov/projects/gapplusprev/sgap_plus.htm) for reported associations of CDH13 genetic variants using “CDH13” as the search term and looked up the GWIS result of these SNPs for interaction with PM10 on decline in FEF25–75%. Second, we used the imputed data obtained in the SAPALDIA discovery sample to construct haplotypes in a 200-kb chromosomal window centered on the GWIS top hit using the software Haploview (Barrett et al. 2005). Third, based on the strong functional link between CDH13 and adiponectin, we looked up the GWIS result of SNPs in the ADIPOQ gene, the adiponectin precursor protein, for interaction with PM10 on decline in FEF25–75% and made a regional association plot of the ADIPOQ locus using the software LocusZoom (Pruim et al. 2010).

Results

Baseline characteristics of the SAPALDIA cohort study participants included in the current GWIS analysis of the discovery and replication sample are presented in Table 1. We observed a highly comparable distribution of sex, age, baseline body mass index, weight change during follow-up, baseline lung function level, and average PM10 exposure at baseline and during follow-up, except for a small difference in proportion in smokers and smoking intensity between the two samples (Table 1).
Table 1 Baseline characteristics of the study population of the SAPALDIA cohort, discovery and replication sample (mean ± SD or %).
CharacteristicDiscoveryReplication
n7633,896
Female (%)51.151.6
Age (years)41.6 ± 11.041.1 ± 11.4
Body mass index (kg/m2)23.5 ± 3.523.7 ± 3.6
Height (cm)169.6 ± 9.0169.5 ± 8.7
Weight (kg)67.9 ± 12.668.4 ± 12.6
Weight changea (kg)5.4 ± 6.25.6 ± 6.0
Baseline lung function
FEF25–75% (L/sec)3.5 ± 1.23.5 ± 1.2
FEV1 (L)3.6 ± 0.83.6 ± 0.8
FEV1/FVC (%)79.6 ± 7.079.4 ± 7.2
FEF25–75%/FVC (%)78.1 ± 24.577.6 ± 24.7
Air pollution exposure
PM10 annual mean (μg/m3)27.4 ± 9.427.3 ± 9.7
PM10 cumulativea (μg/m3)246.5 ± 79.0245.2 ± 81.4
Smoking status
Never-smokersb (%)42.943.8
Pack-yearsc16.5 ± 17.317.0 ± 18.2
Pack-yearsa,c5.4 ± 6.56.3 ± 8.3
Abbreviations: FEF25–75%, forced mid-expiratory flow; FEF25–75%/FVC, ratio of forced mid-expiratory flow and of forced vital capacity; FEV1, forced expiratory volume in the first second; FEV1/FVC, ratio of forced expiratory volume in the first second and of forced vital capacity; FVC, forced vital capacity. aDuring 11-year follow-up. bNever-smokers defined as nonsmoker at baseline and at follow-up survey. cDiscovery sample: missing data on pack-years at baseline (n = 20, 2.6%) and during follow-up (n = 62, 7.3%); replication sample: missing data on pack-years at baseline (n = 67, 1.7%) and during follow-up (n = 326, 8.4%).
GWIS discovery results. We observed significant association signals with a group of 13 SNPs, interacting with cumulative PM10 exposure on annual decline in FEF25–75%, at a single locus on chromosome 16 (Figure 1A). The quantile–quantile plot of interaction p-values showed evidence for a higher number of significant signals than expected by chance (Figure 1B). Based on the lambda observed (λ = 1.0476), adequate genomic control of the genome-wide associations had been applied.
Figure 1 Genome-wide interactions between cumulative PM10 exposure on annual decline in FEF25–75% in the discovery sample (= 763) of the SAPALDIA cohort study. (A) Manhattan plot of the negative log of the p-values (pint) of 2,198,793 SNPs used for interaction association testing. CDH13 SNP cluster with interaction p-values reaching genome-wide significance are above the line and SNPs in the CDH13 locus are highlighted in green. (B) Quantile–quantile plot representing calculated p-values (pint) for each PM10 by SNP interaction tested plotted against the expected chi-square–distributed p-values. Deviation from the diagonal identity line points to the presence of potentially true associations.
The association signal at 16q23.3, located intronic to gene CDH13 (Table 2), had the strongest interaction (pint = 8.8 × 10–10) for rs2325934, an uncommon variant [MAF (frequency of the least common allele in the study sample), 9.6%]. There was evidence for additional potentially independent interaction signals in this locus, as variants with varying MAFs and differing linkage disequilibrium (LD) values also showed significant interactions with cumulative PM10 on annual FEF25–75% decline (e.g., rs11643197: MAF, 13.4%, pint = 6.87 × 10–8, LD with the CDH13 top hit rs2325934 r2 = 0.585 and D´ = 1) (Figure 2; see also Supplemental Material, Figure S1). p-Values for gene main effects (pmain), gene–environment interaction effects (pint), and the joint test (pjoint) are provided in Table 2 for intronic CHD13 SNPs, and in Supplemental Material, Table S1, for the top 1,000 SNPs associated with the decline in FEF25–75% in the discovery GWIS, ranked by gene–environment interaction effects (pint). The joint test of the SNP main effect and the interaction effect (pjoint, null hypothesis: βmain = 0 and βint = 0) did not identify additional genetic modifiers of the association between cumulative PM10 and lung function decline beyond those already identified based on pint for the gene–environment interaction, and the CDH13 variants ranked high according to gene main effect, interaction, and joint tests (maximum genome-wide rank = 22 for the joint test; Table 2).
Table 2 Discovery GWIS top hits with p-values for interaction (pint) < 10–7, clustered intronically to the CDH13 gene: adjusted interaction association with individualized cumulative PM10 exposure on annual decline in FEF25–75% in the SAPALDIA cohort.
dbSNP IDChromosomePositionMinor allele frequency (%)MainpmainInteractionpintJointpjointGenome-wide test rankStrong LD groupsa
MainInteractionJoint
rs2325934b16819000009.646.94 × 10–118.80 × 10–105.80 × 10–10111Reference (A)
rs17282232168190582411.125.26 × 10–97.28 × 10–83.64 × 10–8131713Reference (B)
rs1051458216819104328.663.71 × 10–93.93 × 10–82.72 × 10–8111211A
rs1051458016819108729.711.89 × 10–104.10 × 10–91.47 × 10–9474A
rs1696023416819135129.959.90 × 10–111.82 × 10–98.05 × 10–10232A
rs12325503168191724811.215.98 × 10–96.23 × 10–84.30 × 10–8141514B
rs10514578168191731211.206.01 × 10–96.07 × 10–84.34 × 10–8151415B
rs1721059916819185689.931.18 × 10–101.49 × 10–99.77 × 10–10323A
rs1051457516819313209.853.95 × 10–102.61 × 10–93.12 × 10–9545A
rs1721137116819330409.987.20 × 10–104.03 × 10–95.56 × 10–9666A
rs1424168168193560010.101.36 × 10–96.95 × 10–91.03 × 10–8797A
rs17211581168193724010.111.44 × 10–97.40 × 10–91.09 × 10–88109A
rs17284098b168194757612.871.56 × 10–83.98 × 10–81.02 × 10–7161316Reference (C)
rs17284265168194979212.145.13 × 10–91.57 × 10–83.55 × 10–8121112C
rs17284390168195478411.902.01 × 10–95.08 × 10–91.33 × 10–810810C
rs17212165168195568811.961.71 × 10–93.30 × 10–91.05 × 10–8958C
rs11643197168196479213.366.35 × 10–86.87 × 10–83.17 × 10–7171622C
Ranks shown in the table refer to the genome-wide ranking over 2,198,793 SNPs. The following terminology defined genetic effects referring to their respective null hypothesis: pmain: gene marginal (βmain = 0); pint: gene–environment (βint = 0); pjoint: joint effect (βmain = 0 and βint = 0).GWIS was adjusted for study center, age, sex, height, never-smoking status, weight at baseline, weight change during follow-up, interaction between baseline weight and weight change, seasonal effects of time point of baseline and follow-up examination date (sine and cosine function of day of examination) and population stratification components. Cohort participants with self-report of asthma history had been excluded from the analysis. Discovery sample size was n = 763. aLD (linkage disequilibrium): A indicates strong LD with rs2325934 (r2 > 0.85; D´ = 1), B indicates strong LD with rs17282232 (r2 > 0.85; D´ = 1), and C indicates strong LD with rs17284098 (r2 > 0.85; D´ = 1). The replication SNPs, rs2325934 and rs17284098, were in the moderate linkage disequilibrium with each other (r2 = 0.685, D´ = 1). bSNPs selected for replication; they were in the moderate linkage disequilibrium with respect to r2 (r2 = 0.685, D´ = 1).
Figure 2 Regional association plot showing the p-values of interaction between cumulative PM10 exposure and CDH13 SNPs on annual decline in FEF25–75% in the discovery sample (= 763) of the SAPALDIA cohort study. Shown is the regional association plot for the genome-wide significant GWIS association signal located in the CDH13 gene at 16q23.3. Negative log of the p-values are plotted on the y-axis. Genomic coordinates (Mb) of the plotted SNPs refer to genome build 36/hg18 and dbSNP128 and are given on the x-axis. Linkage disequilibrium information (r2 values) refers to HapMap Phase II data of Caucasian samples. Recombination rate shown over this chromosomal window indicates recombination sites as vertical lines. The plot was generated using LocusZoom (Pruim et al. 2010). Genes in the genomic vicinity are HSBP1, heat-shock factor-binding protein 1; MBTPS1, membrane-bound transcription factor protease, site 1; MIR3182, microRNA 3182; MLYCD, malonyl-CoA decarboxylase; NECAB2, N-terminal EF-hand calcium binding protein 2; OSGIN1, oxidative stress-induced growth inhibitor 1; SLC38A8, solute carrier family 38 (amino acid transporter), member 8.
Next, we performed sensitivity GWIS analyses for annual decline in FEV1, in FEV1/FVC, and in FEF25–75%/FVC to assess the ranking of the CDH13 variants (see Supplemental Material, Table S2). Briefly, although genome-wide significance was not reached, the CDH13 locus was the top-ranking locus in the GWIS for decline in ratios of both FEV1/FVC and FEF25–75%/FVC, but it was not strongly associated with annual decline in FEV1 (pint ≥ 0.001). For FEV1/FVC and FEF25–75%/FVC, the CDH13 SNP with the strongest interaction p-value was rs2325934 (for FEV1/FVC: pint = 1.99 × 10–6; for FEF25–75%/FVC: pint = 1.47 × 10–6).
GWIS replication results. We selected two genome-wide significant CDH13 SNPs for de novo genotyping in the replication study sample (n = 3,896). The rs2325934 variant was selected because it exhibited the lowest p-value of association (pint) in the discovery analysis. A second SNP, rs17284098 (MAF, 12.9%; discovery pint = 3.98 × 10–8) was chosen for its higher MAF compared with the top hit. Both replication SNPs were in moderate LD (r2 = 0.685 and D´ = 1). They replicated yielding p-values below the Bonferroni corrected significance level for two tests (p < 0.025; rs2325934: pint = 0.008; rs17284098: pint = 0.016; Table 3). Interactions between the CDH13 replication SNPs and PM10 were robust to adjustment for different covariates (Table 3). The observed association became slightly stronger despite diminished sample size (n = 3,504) when adjusted for history of smoking intensity, including pack-years at baseline and pack-years smoked during follow-up. However, it is not possible to determine whether the change was attributable to adjustment or to a difference in the sample. Omitting smoking adjustment resulted in weakened signal, as did omitting the interaction term between weight and weight change. In contrast, adjusting additionally for age-squared slightly strengthened the associations. The PM10 effect modification by the CDH13 SNPs remained significant even in minimally adjusted (age, sex, and study area) models (rs2325934: pint = 0.019; rs17284098: pint = 0.020; Table 3).
Table 3 Replication results of adjusteda interaction of CDH13 intronic SNPs (rs2325934 and rs17284098) with cumulative PM10 during 11-year follow-up on annual decline in FEF25–75%, the SAPALDIA cohort study.
dbSNP IDnCoefficientb (95% CI)SEpint
Adjusted modela
rs23259343,8790.0742 (0.0191, 0.1294)0.02810.008
rs172840983,8780.0632 (0.0117, 0.1147)0.02630.016
Adjusted + smoking historyc
rs23259343,5040.0766 (0.0198, 0.1335)0.0290.008
rs172840983,5040.0723 (0.0188, 0.1257)0.02730.008
Adjusted + age-squaredc
rs23259343,8790.0755 (0.0205, 0.1306)0.02810.007
rs172840983,8780.0645 (0.0130, 0.1160)0.02630.014
Adjusted: smoking statusd
rs23259343,8790.074 (0.0188, 0.1291)0.02810.009
rs172840983,8780.0632 (0.0118, 0.1147)0.02620.016
Adjusted: interaction between baseline weight and weight changed
rs23259343,8790.0742 (0.0191, 0.1293)0.02810.008
rs172840983,8780.0629 (0.0114, 0.1144)0.02630.017
Adjusted: weight change and interaction between baseline weight and weight changed
rs23259343,8790.0757 (0.0212, 0.1301)0.02780.006
rs172840983,8780.0626 (0.0116, 0.1136)0.0260.016
Adjusted: weight at baseline and interaction between baseline weight and weight changed
rs23259343,8790.0766 (0.0211, 0.1320)0.02830.007
rs172840983,8780.0657 (0.0141, 0.1173)0.02630.013
Minimal adjustmente
rs23259343,8790.0659 (0.0108, 0.1211)0.02820.019
rs172840983,8780.0606 (0.0095, 0.1117)0.02610.02
aSame adjustments applied as for discovery GWIS, including study center, age, sex, height, never-smoking status, weight at baseline, weight change during follow-up, interaction between baseline weight and weight change, seasonal effects of time point of baseline and follow-up examination date (sine and cosine function of day of examination). No adjustment for population stratification was available. bCoefficient refers to the additive SNP effect in annual change in FEF25–75% (mL/sec) per 1-μg/m3 change in PM10 exposure. cSame adjustments as for the discovery GWIS, adding indicated additional covariate(s) in the model. dSame adjustments as for the discovery GWIS, omitting indicated covariate(s) in the model. eBasic adjustment including only study center, age, and sex in the model.
In a genotype-stratified analysis, combining discovery and replication sample, for both CDH13 SNPs, the PM10–FEF25–75% association appeared to be restricted to participants who were homozygous for the major allele (see Supplemental Material, Table S3). For rs2325934, the major homozygous genotype strata (n = 3,750) was estimated to have an average annual change in FEF25–75% of –0.102 mL per increase of 1 μg/m3 PM10 [95% confidence interval (CI): –0.19, –0.01; p = 0.03] in contrast with the group carrying at least one minor allele (estimated average annual change of 0.074 mL; 95% CI: –0.16, 0.31; p = 0.53, n = 886).
Previously published CDH13 genome-wide association results and linkage disequilibrium in the 200-kb GWIS window. In GWAS, variants of the CDH13 gene were previously associated with a number of different phenotypes (see Supplemental Material, Table S4). Genetic variants in the 5´ end of the CDH13 gene have been repeatedly associated with circulating adiponectin levels (Chung et al. 2011; Dastani et al. 2012; Jee et al. 2010; Morisaki et al. 2012; Wu et al. 2010). We thus looked up the interactions of these previous CDH13 GWAS hits with PM10 on decline in FEF25–75%, but none of these SNPs ranked high in the current GWIS (pint ≥ 0.02; see Supplemental Material, Table S4). Pairwise LD in the CDH13 gene between the GWAS SNPs and the PM10 interacting SNPs was low (r2 ≤ 0.33).
The structure of linkage disequilibrium (see Supplemental Material, Figure S1) and haplotypes (see Supplemental Material, Figure S2) in the 200-kb window centered on rs2325934 were constructed using the imputed genotype data of the discovery sample. The LD pattern and derived haplotypes suggested that the GWIS top hit tagged specifically one 33 kb–long haplotype (block 9, stretching from rs2352934 to rs1426166). The second replication variant, rs17284098, was located in a different 11 kb–long haplotype (block 11, stretching from rs1424168 to rs17284098) in an intron downstream of the GWIS top hit.

Discussion

To our knowledge, this is the first report presenting a genome-wide interaction study aiming to identify novel genes modifying the association of PM on lung function decline. We identified a cluster of SNPs intronic to the gene CDH13 that modified the estimated effect of cumulative PM10 on the decline in FEF25–75% in our study population. We estimated that participants who were homozygous for the major allele of rs2325934 experienced an excess average decline of 11 mL/sec in FEF25–75% per 10-μg/m3 increase in cumulative PM10 exposure over 11 years. Interestingly, cadherin 13, the protein encoded by CDH13 is functionally linked to adiponectin, a predominantly anti-inflammatory adipokine.
Experimental animal studies have provided strong evidence that major cellular responses to PM exposure include oxidative stress (Manzo et al. 2012) and inflammation (Uski et al. 2012). Results from 125 subjects monitored before, during, and after the Beijing Olympics were consistent with oxidative and inflammatory effects of ambient air pollution in the respiratory tract (Huang et al. 2012). The relationship between PM exposure and systemic inflammation, as indicated by serum C-reactive protein (CRP), was the subject of a recent systematic review reporting more consistent results for a positive PM–CRP association in longitudinal studies of healthy subjects than in short-term studies or longitudinal studies of subjects with chronic inflammatory conditions (Li et al. 2012). The few candidate gene–air pollution interaction studies published to date also support the oxidative and inflammatory effects of air pollution in mediating adverse respiratory health effects (Breton et al. 2011; Curjuric et al. 2010; Imboden et al. 2009; Melén et al. 2008; Romieu et al. 2006; Yang et al. 2005).
In humans, CDH13 is expressed in various lung cell types, including bronchial epithelium and airway smooth muscle cells. The CDH13 gene, spanning 1.17 Mb, encodes 15 different transcripts with alternate exons that produce structural proteins, which are expressed in endothelia, epithelia (including bronchial epithelial cells), smooth muscle cells, and in nervous tissue. The major CDH13 transcript contains 14 exons (Figure 3), encoding an open reading frame for a 713-amino acid polypeptide, cadherin 13, also known as T-cadherin, H-cadherin, or vascular adiponectin receptor. Previous GWAS have indicated that CDH13 genetic variants may contribute to various phenotypes. The predicted molecular and cellular functions of cadherin 13 are congruent with some of the GWAS findings. The strongest and most consistent GWAS signals have been SNPs in the 5´ untranslated region or in intron 1 associated with circulating adiponectin levels, consistent with the molecular function of adiponectin binding (GO:0055100) (Chung et al. 2011; Dastani et al. 2012; Jee et al. 2010; Morisaki et al. 2012; Wu et al. 2010). CDH13 SNPs have also been reported to be associated with body height (Okada et al. 2010) and with respiratory function (http://www.ncbi.nlm.nih.gov/projects/gapplusprev/sgap_plus.htm).
Figure 3 Schematic representation of CDH13 gene and genetic variants identified by GWAS to be associated with various phenotypes. The CDH13 gene, spanning 1.17 Mb, at 82.6 Mb (build 36) on chromosome 16, encodes 15 different transcripts with alternate exons which produce structural proteins. The major CDH13 transcript contains 14 exons, encoding an open reading frame for a 713-amino acid polypeptide. In this schematic view of the CDH13 gene, we pinpoint SNPs identified in the dbGaP database (http://www.ncbi.nlm.nih.gov/projects/gapplusprev/sgap_plus.htm) for reported associations with various phenotypes using “CDH13” as the search term. For a list of the associated phenotypes, see Supplemental Material, Table S4. *SNPs, rs2325934 and rs17284098, identified in the present report to interact with PM10 exposure on decline in FEF25–75%.
Cadherin 13 is one type of adiponectin-binding protein (Hug et al. 2004), among others such as adiponectin receptors (AdipoR1, AdipoR2) (Yamauchi et al. 2003) or calreticulin (Takemura et al. 2007), and might exert its role in respiratory health through adiponectin. Adiponectin—a 244 amino acid protein resembling collagen VII, X, and complement factor C1—has been identified as a potent and pleiotropic regulator of inflammation (Ohashi et al. 2012). Experimental evidence in mice demonstrated that cadherin 13 was required to mediate the protective effect of adiponectin on allergen-induced airway inflammation (Williams et al. 2012). In human studies, CDH13 has been consistently identified by GWAS as a determinant of circulating adiponectin (Dastani et al. 2012), and serum adiponectin concentrations were positively associated with peak lung function in a prospective study of young healthy adults (Thyagarajan et al. 2010).
Adiponectin is secreted primarily by visceral adipocytes (Arita et al. 1999). Body composition and especially visceral adiposity have been associated with lower lung function and accelerated age-related decline (Rossi et al. 2011; Wehrmeister et al. 2012). We were among the first to report evidence of a modifying effect of obesity on the association of air pollution with lung function decline (Schikowski et al. 2013).
A strength of the current analysis is the detailed characterization of the cohort participants, as well as the cohort’s prospective design to specifically investigate longitudinally air pollution health effects. The fact that discovery and replication samples derive from the SAPALDIA cohort pool is a strength of this analysis. Both samples were recruited at the same time by the same field workers using the same standardized procedures. Nonetheless, larger studies in independent populations with different environmental and ethnic characteristics are needed to confirm the observed interactions between CDH13 genetic variants and PM10. Genome-wide interactions between environmental exposures and genetic variants on complex health phenotypes form an active field of investigation, and novel methodologies are being developed to address analytical challenges associated with this research (Ege et al. 2011; Gauderman et al. 2013; Hutter et al. 2013; Sohns et al. 2013). Recently, a genome-wide analysis investigating genetic modifiers of associations between occupational exposures and lung function combined the GWIS approach with an in silico pathway analysis that indicated the involvement of inflammatory pathways (Liao et al. 2013).
There are several limitations, in addition to low power, in the present study. First, less than half of all baseline cohort participants were included in the analysis, leaving room for potential bias. Second, given the pleiotropic health outcomes associated with CDH13 in previous GWAS, it is conceivable that modification of the air pollution–lung function association reflected an underlying susceptibility caused by health conditions associated with CDH13 genotypes, such as cardiovascular phenotypes, rather than a causal interaction between the genotypes and PM10. In the absence of measured adiponectin levels we cannot verify whether the apparent modifying effect of CDH13 is mediated through this adipokine. Genetic variants of the adiponectin precursor protein encoded by the gene ADIPOQ were by definition included in our GWIS analysis. We looked up the interactions between PM10 exposure and nine haplotype tagging variants of the ADIPOQ gene (Peters et al. 2013). None of these SNPs showed significant interactions (pint = 0.07 to pint = 0.95; for regional association plot of the ADIPOQ locus, see Supplemental Material, Figure S3). The CDH13 SNPs associated with adiponectin circulating level (rs3865188, rs4783244, rs12051272) (Chung et al. 2011; Dastani et al. 2012; Jee et al. 2010; Morisaki et al. 2012; Wu et al. 2010) were in very low LD (r2 < 0.1) and thus contained in a different haplotype block than the CDH13 SNPs interacting with PM10 (Figure 3). Nevertheless, rs3865188 previously associated with adiponectin (Jee et al. 2010; Wu et al. 2010) interacted with PM10 at a nominal p-value of 0.06 in our study.
Because we focused a priori on FEF25–75% decline for this GWIS, we did not evaluate interactions between PM10 and previously identified GWAS signals on associations with FEV1 or FVC (Artigas et al. 2011; Hancock et al. 2010, 2012; Obeidat et al. 2011; Repapi et al. 2010; Wilk et al. 2009). Another limitation is that although asthmatic subjects appear to be more vulnerable to effects of air pollution exposure (Trasande and Thurston 2005), we restricted the current analyses to nonasthmatics. This sample restriction was based on our previous finding of an extended heterogeneity in the GWAS-identified determinants of lung function decline of healthy individuals compared with asthmatics (Imboden et al. 2012). Finally, comparing genotype-stratified analyses in the SAPALDIA cohort, the size of effect modification by the CDH13 SNPs reported here on the PM10–FEF25–75% association is substantially smaller than the size of effect modification by SERPINA1 genotypes (underlying intermediate alpha 1 antitrypsin deficiency) on the association of FEF25–75% decline with occupational exposure to vapors, dusts, gases, and fumes (Mehta et al. 2012). This is not unexpected for genome-wide signals of unknown functional relevance. The limitation of genome-wide signals with regard to clinical utility is a well-known problem of GWAS.
In conclusion, the mechanistic link between adiponectin (its modulating action on inflammatory processes systemically and locally in the lung) and cadherin 13 (its sequestering action on circulating adiponectin levels) make our GWIS finding, CDH13, a biologically plausible candidate gene for modifying the air pollution exposure effect. Follow-up studies need to confirm the observed interaction with CDH13 SNPs and must assess whether the finding is related to recent evidence on the modifying effect of obesity on the association between PM10 and decline in lung function.

Acknowledgments

The SAPALDIA Team: Study directorate: N.M. Probst Hensch (e/g), T. Rochat (p), N. Künzli (e/exp), C. Schindler (s), J.M. Gaspoz (c). Scientific team: J.C. Barthélémy (c), W. Berger (g), R. Bettschart (p), A. Bircher (a), O. Brändli (p), C. Brombach (n), M. Brutsche (p), L. Burdet (p), D. Felber Dietrich (exp), M. Frey (p), U. Frey (pd), M.W. Gerbase (p), D. Gold (e/c/p), E. de Groot (c), W. Karrer (p), R. Keller (p), B. Martin (pa), D. Miedinger (o), U. Neu (exp), L. Nicod (p), M. Pons (p), F. Roche (c), T. Rothe (p), E. Russi (p), P. Schmid-Grendelmeyer (a), A. Schmidt-Trucksäss (pa), A. Turk (p), J. Schwartz (e), D. Stolz (p), P. Straehl (exp), J.M. Tschopp (p), A. von Eckardstein (cc), E. Zemp Stutz (e). Scientific team at coordinating centers: M. Adam (e/g), C. Autenrieth (pa), P.O. Bridevaux (p), D. Carballo (c), E. Corradi (exp), I. Curjuric (e), J. Dratva (e), A. Di Pasquale (s), E. Dupuis Lozeron (s), E. Fischer (e), M. Germond (s), L. Grize (s), D. Keidel (s), S. Kriemler (pa), A. Kumar (g), M. Imboden (g), N. Maire (s), A. Mehta (e), H. Phuleria (exp), E. Schaffner (s), G.A. Thun (g), A. Ineichen (exp), M. Ragettli (e), M. Ritter (exp), T. Schikowski (e), M. Tarantino (s), M. Tsai (exp). Abbreviations: (a) allergology, (c) cardiology, (cc) clinical chemistry, (e) epidemiology, (exp) exposure, (g) genetic and molecular biology, (m) meteorology, (n) nutrition, (o) occupational health, (p) pneumology, (pa) physical activity, (pd) pediatrics, (s) statistics.
The study could not have been done without the help of the study participants, technical and administrative support and the medical teams, and field workers at the local study sites.

Supplemental Material

(1.7 MB) PDF
Click here for additional data file.

References

Ackermann-Liebrich U, Kuna-Dibbert B, Probst-Hensch N, Schindler C, Felber Dietrich D, Zemp Stutz Eet al. 2005. Follow-up of the Swiss Cohort Study on Air Pollution and Lung Diseases in Adults (SAPALDIA 2) 1991–2003: methods and characterization of participants. Soz Praventivmed 50:245-263 https://pubmed.ncbi.nlm.nih.gov/16167509/.
American Thoracic Society. 1995. Standardization of Spirometry, 1994 update. Am J Respir Crit Care Med 152(3):1107-1136 https://pubmed.ncbi.nlm.nih.gov/7663792/.
Andreau K, Leroux M, Bouharrour A. 2012. Health and cellular impacts of air pollutants: from cytoprotection to cytotoxicity. Biochem Res Int 2012. 493894; https://doi.org/10.1155/2012/493894.
Arita Y, Kihara S, Ouchi N, Takahashi M, Maeda K, Miyagawa Jet al. 1999. Paradoxical decrease of an adipose-specific protein, adiponectin, in obesity. Biochem Biophys Res Commun 257:79-83 https://pubmed.ncbi.nlm.nih.gov/10092513/.
Artigas MS, Loth DW, Wain LV, Gharib SA, Obeidat M, Tang Wet al. 2011. Genome-wide association and large-scale follow up identifies 16 new loci influencing lung function. Nat Genet 43:1082-1090 https://pubmed.ncbi.nlm.nih.gov/21946350/.
Barrett JC, Fry B, Maller J, Daly MJ. 2005. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21:263-265 https://pubmed.ncbi.nlm.nih.gov/15297300/.
Breton CV, Salam MT, Vora H, Gauderman WJ, Gilliland FD. 2011. Genetic variation in the glutathione synthesis pathway, air pollution, and children’s lung function growth. Am J Respir Crit Care Med 183:243-248 https://pubmed.ncbi.nlm.nih.gov/20802163/.
Brook RD, Rajagopalan S, Pope CA, Brook JR, Bhatnagar A, Diez-Roux AVet al. 2010. Particulate matter air pollution and cardiovascular disease: an update to the scientific statement from the American Heart Association. Circulation 121:2331-2378 https://pubmed.ncbi.nlm.nih.gov/20458016/.
Chung CM, Lin TH, Chen JW, Leu HB, Yang HC, Ho HYet al. 2011. A genome-wide association study reveals a quantitative trait locus of adiponectin on CDH13 that predicts cardiometabolic outcomes. Diabetes 60:2417-2423 https://pubmed.ncbi.nlm.nih.gov/21771975/.
Curjuric I, Imboden M, Nadif R, Kumar A, Schindler C, Haun Met al. 2012. Different genes interact with particulate matter and tobacco smoke exposure in affecting lung function decline in the general population. PLoS One 7:e40175; https://doi.org/10.1371/journal.pone.0040175 https://pubmed.ncbi.nlm.nih.gov/22792237/.
Curjuric I, Imboden M, Schindler C, Downs SH, Hersberger M, Liu SLet al. 2010. HMOX1 and GST variants modify attenuation of FEF25–75% decline due to PM10 reduction. Eur Respir J 35:505-514 https://pubmed.ncbi.nlm.nih.gov/20190330/.
Dastani Z, Hivert MF, Timpson N, Perry JR, Yuan X, Scott RAet al. 2012. Novel loci for adiponectin levels and their influence on type 2 diabetes and metabolic traits: a multi-ethnic meta-analysis of 45,891 individuals. PLoS Genet 8:e1002607; https://doi.org/10.1371/journal.pgen.1002607 https://pubmed.ncbi.nlm.nih.gov/22479202/.
Downs SH, Schindler C, Liu LJ, Keidel D, Bayer-Oglesby L, Brutsche MHet al. 2007. Reduced exposure to PM10 and attenuated age-related decline in lung function. N Engl J Med 357:2338-2347 https://pubmed.ncbi.nlm.nih.gov/18057336/.
Ege MJ, Strachan DP, Cookson WO, Moffatt MF, Gut I, Lathrop M, et al. 2011. Gene-environment interaction for childhood asthma and exposure to farming in Central Europe. J Allergy Clin Immunol 127:138-144.e4; https://doi.org/10.1016/j.jaci.2010.09.041 https://pubmed.ncbi.nlm.nih.gov/21211648/.
Gauderman WJ, Zhang P, Morrison JL, Lewinger JP. 2013. Finding novel genes by testing G x E interactions in a genome-wide association study. Genet Epidemiol 37:603-613 https://pubmed.ncbi.nlm.nih.gov/23873611/.
Hancock DB, Artigas MS, Gharib SA, Henry A, Manichaikul A, Ramasamy Aet al. 2012. Genome-wide joint meta-analysis of SNP and SNP-by-smoking interaction identifies novel loci for pulmonary function. PLoS Genet 8:e1003098; https://doi.org/10.1371/journal.pgen.1003098 https://pubmed.ncbi.nlm.nih.gov/23284291/.
Hancock DB, Eijgelsheim M, Wilk JB, Gharib SA, Loehr LR, Marciante KDet al. 2010. Meta-analyses of genome-wide association studies identify multiple loci associated with pulmonary function. Nat Genet 42:45-52 https://pubmed.ncbi.nlm.nih.gov/20010835/.
Hansel NN, Ruczinski I, Rafaels N, Sin DD, Daley D, Malinina Aet al. 2013. Genome-wide study identifies two loci associated with lung function decline in mild to moderate COPD. Hum Genet 132:79-90 https://pubmed.ncbi.nlm.nih.gov/22986903/.
Heath SC, Gut IG, Brennan P, McKay JD, Bencko V, Fabianova Eet al. 2008. Investigation of the fine structure of European populations with applications to disease association studies. Eur J Hum Genet 16:1413-1429 https://pubmed.ncbi.nlm.nih.gov/19020537/.
Huang W, Wang G, Lu SE, Kipen H, Wang Y, Hu Met al. 2012. Inflammatory and oxidative stress responses of healthy young adults to changes in air quality during the Beijing Olympics. Am J Respir Crit Care Med 186:1150-1159 https://pubmed.ncbi.nlm.nih.gov/22936356/.
Hug C, Wang J, Ahmad NS, Bogan JS, Tsao TS, Lodish HF. 2004. T-cadherin is a receptor for hexameric and high-molecular-weight forms of Acrp30/adiponectin. Proc Natl Acad Sci USA 101:10308-10313 https://pubmed.ncbi.nlm.nih.gov/15210937/.
Hutter CM, Mechanic LE, Chatterjee N, Kraft P, Gillanders EM, Gene-Environment Think Tank NCI. 2013. Gene-environment interactions in cancer epidemiology: a National Cancer Institute Think Tank report. Genet Epidemiol 37:643-657 https://pubmed.ncbi.nlm.nih.gov/24123198/.
Imboden M, Bouzigon E, Curjuric I, Ramasamy A, Kumar A, Hancock DBet al. 2012. Genome-wide association study of lung function decline in adults with and without asthma. J Allergy Clin Immunol 129:1218-1228 https://pubmed.ncbi.nlm.nih.gov/22424883/.
Imboden M, Nieters A, Bircher AJ, Brutsche M, Becker N, Wjst Met al. 2006. Cytokine gene polymorphisms and atopic disease in two European cohorts. (ECRHS-Basel and SAPALDIA). Clin Mol Allergy 4:9; https://doi.org/10.1186/1476-7961-4-9 https://pubmed.ncbi.nlm.nih.gov/16759385/.
Imboden M, Schwartz J, Schindler C, Curjuric I, Berger W, Liu SLet al. 2009. Decreased PM10 exposure attenuates age-related lung function decline: genetic variants in p53, p21, and CCND1 modify this effect. Environ Health Perspect 117:1420-1427; https://doi.org/10.1289/ehp.0800430 https://pubmed.ncbi.nlm.nih.gov/19750108/.
Jee SH, Sull JW, Lee JE, Shin C, Park J, Kimm Het al. 2010. Adiponectin concentrations: a genome-wide association study. Am J Hum Genet 87:545-552 https://pubmed.ncbi.nlm.nih.gov/20887962/.
Katanoda K, Sobue T, Satoh H, Tajima K, Suzuki T, Nakatsuka Het al. 2011. An association between long-term exposure to ambient air pollution and mortality from lung cancer and respiratory diseases in Japan. J Epidemiol 21:132-143 https://pubmed.ncbi.nlm.nih.gov/21325732/.
Kraft P, Yen YC, Stram DO, Morrison J, Gauderman WJ. 2007. Exploiting gene-environment interaction to detect genetic associations. Hum Hered 63:111-119 https://pubmed.ncbi.nlm.nih.gov/17283440/.
Künzli N, Kuna-Dibbert B, Keidel D, Keller R, Brändli O, Schindler Cet al. 2005. Longitudinal validity of spirometers—a challenge in lung function follow-up studies. Swiss Med Wkly 135(33–34):503-508 https://pubmed.ncbi.nlm.nih.gov/16208589/.
Li Y, Rittenhouse-Olson K, Scheider WL, Mu L. 2012. Effect of particulate matter air pollution on C-reactive protein: a review of epidemiologic studies. Rev Environ Health 27:133-149 https://pubmed.ncbi.nlm.nih.gov/23023922/.
Liao SY, Lin X, Christiani DC. 2013. Gene-environment interaction effects on lung function- a genome-wide association study within the Framingham heart study. Environ Health 12:101; https://doi.org/10.1186/1476-069X-12-101 https://pubmed.ncbi.nlm.nih.gov/24289273/.
Liu LJ, Curjuric I, Keidel D, Heldstab J, Künzli N, Bayer-Oglesby Let al. 2007. Characterization of source-specific air pollution exposure for a large population-based Swiss cohort (SAPALDIA). Environ Health Perspect 115:1638-1645; https://doi.org/10.1289/ehp.10177 https://pubmed.ncbi.nlm.nih.gov/18007997/.
Manzo ND, LaGier AJ, Slade R, Ledbetter AD, Richards JH, Dye JA. 2012. Nitric oxide and superoxide mediate diesel particle effects in cytokine-treated mice and murine lung epithelial cells—implications for susceptibility to traffic-related air pollution. Part Fibre Toxicol 9:43; https://doi.org/10.1186/1743-8977-9-43 https://pubmed.ncbi.nlm.nih.gov/23151036/.
Martin BW, Ackermann-Liebrich U, Leuenberger P, Künzli N, Stutz EZ, Keller Ret al. 1997. SAPALDIA: methods and participation in the cross-sectional part of the Swiss Study on Air Pollution and Lung Diseases in Adults. Soz Praventivmed 42:67-84 https://pubmed.ncbi.nlm.nih.gov/9151378/.
Mehta AJ, Adam M, Schaffner E, Barthélémy JC, Carballo D, Gaspoz JMet al. 2012. Heart rate variability in association with frequent use of household sprays and scented products in SAPALDIA. Environ Health Perspect 120:958-964; https://doi.org/10.1289/ehp.1104567 https://pubmed.ncbi.nlm.nih.gov/22538298/.
Melén E, Nyberg F, Lindgren CM, Berglind N, Zucchelli M, Nordling Eet al. 2008. Interactions between glutathione S-transferase P1, tumor necrosis factor, and traffic-related air pollution for development of childhood allergic disease. Environ Health Perspect 116:1077-1084; https://doi.org/10.1289/ehp.11117 https://pubmed.ncbi.nlm.nih.gov/18709160/.
Moffatt MF, Gut IG, Demenais F, Strachan DP, Bouzigon E, Heath Set al. 2010. A large-scale, consortium-based genomewide association study of asthma. N Engl J Med 363:1211-1221 https://pubmed.ncbi.nlm.nih.gov/20860503/.
Morisaki H, Yamanaka I, Iwai N, Miyamoto Y, Kokubo Y, Okamura Tet al. 2012. CDH13 gene coding T-cadherin influences variations in plasma adiponectin levels in the Japanese population. Hum Mutat 33:402-410 https://pubmed.ncbi.nlm.nih.gov/22065538/.
Obeidat M, Wain LV, Shrine N, Kalsheker N, Soler Artigas M, Repapi Eet al. 2011. A comprehensive evaluation of potential lung function associated genes in the SpiroMeta general population sample. PLoS One 6:e19382; https://doi.org/10.1371/journal.pone.0019382 https://pubmed.ncbi.nlm.nih.gov/21625484/.
Ohashi K, Ouchi N, Matsuzawa Y. 2012. Anti-inflammatory and anti-atherogenic properties of adiponectin. Biochimie 94:2137-2142 https://pubmed.ncbi.nlm.nih.gov/22713764/.
Okada Y, Kamatani Y, Takahashi A, Matsuda K, Hosono N, Ohmiya Het al. 2010. A genome-wide association study in 19 633 Japanese subjects identified LHX3-QSOX2 and IGF1 as adult height loci. Hum Mol Genet 19:2303-2312 https://pubmed.ncbi.nlm.nih.gov/20189936/.
Peters KE, Beilby J, Cadby G, Warrington NM, Bruce DG, Davis WAet al. 2013. A comprehensive investigation of variants in genes encoding adiponectin (ADIPOQ) and its receptors (ADIPOR1/R2), and their association with serum adiponectin, type 2 diabetes, insulin resistance and the metabolic syndrome. BMC Med Genet 14:15; https://doi.org/10.1186/1471-2350-14-15 https://pubmed.ncbi.nlm.nih.gov/23351195/.
Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TPet al. 2010. LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics 26:2336-2337 https://pubmed.ncbi.nlm.nih.gov/20634204/.
Repapi E, Sayers I, Wain LV, Burton PR, Johnson T, Obeidat Met al. 2010. Genome-wide association study identifies five loci associated with lung function. Nat Genet 42:36-44 https://pubmed.ncbi.nlm.nih.gov/20010834/.
Romieu I, Ramirez-Aguilar M, Sienra-Monge JJ, Moreno-Macías H, del Rio-Navarro BE, David Get al. 2006. GSTM1 and GSTP1 and respiratory health in asthmatic children exposed to ozone. Eur Respir J 28:953-959 https://pubmed.ncbi.nlm.nih.gov/16870661/.
Romieu I, Riojas-Rodríguez H, Marrón-Mares AT, Schilmann A, Perez-Padilla R, Masera O. 2009. Improved biomass stove intervention in rural Mexico: impact on the respiratory health of women. Am J Respir Crit Care Med 180:649-656 https://pubmed.ncbi.nlm.nih.gov/19556519/.
Rossi AP, Watson NL, Newman AB, Harris TB, Kritchevsky SB, Bauer DCet al. 2011. Effects of body composition and adipose tissue distribution on respiratory function in elderly men and women: the health, aging, and body composition study. J Gerontol A Biol Sci Med Sci 66:801-808 https://pubmed.ncbi.nlm.nih.gov/21498841/.
Schikowski T, Ranft U, Sugiri D, Vierkötter A, Brüning T, Harth Vet al. 2010. Decline in air pollution and change in prevalence in respiratory symptoms and chronic obstructive pulmonary disease in elderly women. Respir Res 11:113; https://doi.org/10.1186/1465-9921-11-113 https://pubmed.ncbi.nlm.nih.gov/20727210/.
Schikowski T, Schaffner E, Meier F, Phuleria HC, Vierkötter A, Schindler Cet al. 2013. Improved air quality and attenuated lung function decline: modification by obesity in the SAPALDIA cohort. Environ Health Perspect 121:1034-1039; https://doi.org/10.1289/ehp.1206145 https://pubmed.ncbi.nlm.nih.gov/23820868/.
Sohns M, Viktorova E, Amos CI, Brennan P, Fehringer G, Gaborieau Vet al. 2013. Empirical hierarchical Bayes approach to gene-environment interactions: development and application to genome-wide association studies of lung cancer in TRICL. Genet Epidemiol 37:551-559 https://pubmed.ncbi.nlm.nih.gov/23893921/.
Takemura Y, Ouchi N, Shibata R, Aprahamian T, Kirber MT, Summer RSet al. 2007. Adiponectin modulates inflammatory reactions via calreticulin receptor-dependent clearance of early apoptotic bodies. J Clin Invest 117:375-386 https://pubmed.ncbi.nlm.nih.gov/17256056/.
Thun GA, Ferrarotti I, Imboden M, Rochat T, Gerbase M, Kronenberg Fet al. 2012. SERPINA1 PiZ and PiS heterozygotes and lung function decline in the SAPALDIA cohort. PLoS One 7:e42728; https://doi.org/10.1371/journal.pone.0042728 https://pubmed.ncbi.nlm.nih.gov/22912729/.
Thyagarajan B, Jacobs DR, Smith LJ, Kalhan R, Gross MD, Sood A. 2010. Serum adiponectin is positively associated with lung function in young adults, independent of obesity: the CARDIA study. Respir Res 11:176; https://doi.org/10.1186/1465-9921-11-176 https://pubmed.ncbi.nlm.nih.gov/21143922/.
Trasande L, Thurston GD. 2005. The role of air pollution in asthma and other pediatric morbidities. J Allergy Clin Immunol 115:689-699 https://pubmed.ncbi.nlm.nih.gov/15805986/.
Uski OJ, Happo MS, Jalava PI, Brunner T, Kelz J, Obernberger Iet al. 2012. Acute systemic and lung inflammation in C57Bl/6J mice after intratracheal aspiration of particulate matter from small-scale biomass combustion appliances based on old and modern technologies. Inhal Toxicol 24:952-965 https://pubmed.ncbi.nlm.nih.gov/23216156/.
Wehrmeister FC, Menezes AM, Muniz LC, Martínez-Mesa J, Domingues MR, Horta BL. 2012. Waist circumference and pulmonary function: a systematic review and meta-analysis. Syst Rev 1:55; https://doi.org/10.1186/2046-4053-1-55 https://pubmed.ncbi.nlm.nih.gov/23153289/.
Wilk JB, Chen TH, Gottlieb DJ, Walter RE, Nagle MW, Brandler BJet al. 2009. A genome-wide association study of pulmonary function measures in the Framingham Heart Study. PLoS Genet 5:e1000429; https://doi.org/10.1371/journal.pgen.1000429 https://pubmed.ncbi.nlm.nih.gov/19300500/.
Williams AS, Kasahara DI, Verbout NG, Fedulov AV, Zhu M, Si Het al. 2012. Role of the adiponectin binding protein, T-cadherin (Cdh13), in allergic airways responses in mice. PLoS One 7:e41088; https://doi.org/10.1371/journal.pone.0041088 https://pubmed.ncbi.nlm.nih.gov/22815927/.
Wu Y, Li Y, Lange EM, Croteau-Chonka DC, Kuzawa CW, McDade TWet al. 2010. Genome-wide association study for adiponectin levels in Filipino women identifies CDH13 and a novel uncommon haplotype at KNG1-ADIPOQ. Hum Mol Genet 19:4955-4964 https://pubmed.ncbi.nlm.nih.gov/20876611/.
Yamauchi T, Kamon J, Ito Y, Tsuchida A, Yokomizo T, Kita Set al. 2003. Cloning of adiponectin receptors that mediate antidiabetic metabolic effects. Nature 423:762-769 https://pubmed.ncbi.nlm.nih.gov/12802337/.
Yang IA, Holz O, Jörres RA, Magnussen H, Barton SJ, Rodríguez Set al. 2005. Association of tumor necrosis factor-α polymorphisms and ozone-induced change in lung function. Am J Respir Crit Care Med 171:171-176 https://pubmed.ncbi.nlm.nih.gov/15486341/.

Information & Authors

Information

Published In

Environmental Health Perspectives
Volume 123Issue 1January 2015
Pages: 72 - 79
PubMed: 25127211

History

Received: 20 July 2013
Accepted: 13 August 2014
Published online: 15 August 2014

Authors

Affiliations

Medea Imboden
Swiss Tropical and Public Health Institute, Basel, Switzerland
University of Basel, Basel, Switzerland
Ashish Kumar
Swiss Tropical and Public Health Institute, Basel, Switzerland
University of Basel, Basel, Switzerland
Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, United Kingdom
Ivan Curjuric
Swiss Tropical and Public Health Institute, Basel, Switzerland
University of Basel, Basel, Switzerland
Martin Adam
Swiss Tropical and Public Health Institute, Basel, Switzerland
University of Basel, Basel, Switzerland
Gian Andri Thun
Swiss Tropical and Public Health Institute, Basel, Switzerland
University of Basel, Basel, Switzerland
Margot Haun
Division of Genetic Epidemiology, Department of Medical Genetics, Molecular and Clinical Pharmacology, Innsbruck Medical University, Innsbruck, Austria
Ming-Yi Tsai
Swiss Tropical and Public Health Institute, Basel, Switzerland
University of Basel, Basel, Switzerland
Marco Pons
Division of Pulmonary Medicine, Regional Hospital of Lugano, Lugano, Switzerland
Robert Bettschart
Lungenzentrum, Hirslanden Klinik, Aarau, Switzerland
Alexander Turk
Zürcher Höhenklinik, Wald, Faltigberg-Wald, Switzerland
Thierry Rochat
Division of Pulmonary Medicine, University Hospitals of Geneva, Geneva, Switzerland
Nino Künzli
Swiss Tropical and Public Health Institute, Basel, Switzerland
University of Basel, Basel, Switzerland
Christian Schindler
Swiss Tropical and Public Health Institute, Basel, Switzerland
University of Basel, Basel, Switzerland
Florian Kronenberg
Division of Genetic Epidemiology, Department of Medical Genetics, Molecular and Clinical Pharmacology, Innsbruck Medical University, Innsbruck, Austria
Nicole M. Probst-Hensch
Swiss Tropical and Public Health Institute, Basel, Switzerland
University of Basel, Basel, Switzerland

Notes

Address correspondence to M. Imboden, SwissTPH, Socinstr. 59, 4002 Basel, Switzerland. Telephone: 41 61 284 82 40. E-mail: [email protected]

Competing Interests

The authors declare they have no actual or potential competing financial interests.

Funding Information

This research was supported by the Swiss National Science Foundation (grants 33CSCO-134276/1, 33CSCO-108796, 3247BO-104283, 3247BO-104288, 3247BO-104284, 3247-065896, 3100-059302, 3200-052720, 3200-042532, 4026-028099); the Federal Office for Forest, Environment and Landscape; the Federal Office of Public Health; the Federal Office of Roads and Transport; the canton’s government of Aargau, Basel-Stadt, Basel-Land, Geneva, Luzern, Ticino, Valais, and Zürich; the Swiss Lung League; the canton’s Lung League of Basel Stadt/Basel Landschaft, Geneva, Ticino, Valais, and Zurich; SUVA (Schweizerische Versicherungsanstalt); Freiwillige Akademische Gesellschaft; UBS Wealth Foundation; Talecris Biotherapeutics GmbH; Abbott Diagnostics; European Commission 018996 (GABRIEL); Wellcome Trust (WT 084703MA).

Metrics & Citations

Metrics

About Article Metrics


Citations

Download citation

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click DOWNLOAD.

Cited by

  • Air Pollution and Its Impact on Health and Performance in Football Players, Sports, 10.3390/sports13060170, 13, 6, (170), (2025).
  • Joint association of air pollution exposure and inflammation-related proteins in relation to infant lung function, International Journal of Hygiene and Environmental Health, 10.1016/j.ijheh.2023.114294, 255, (114294), (2024).
  • Body composition modify the association between ambient particulate matter and lung function among asthma patients, Environmental Science and Pollution Research, 10.1007/s11356-023-28597-4, 30, 37, (88060-88071), (2023).
  • Swiss Cohort & Biobank – The White Paper, Public Health Reviews, 10.3389/phrs.2022.1605660, 43, (2022).
  • Inflammation resolution in environmental pulmonary health and morbidity, Toxicology and Applied Pharmacology, 10.1016/j.taap.2022.116070, 449, (116070), (2022).
  • Association Between Metabolic and Hormonal Derangements and Professional Exposure to Urban Pollution in a High Intensity Traffic Area, Frontiers in Endocrinology, 10.3389/fendo.2020.00509, 11, (2020).
  • Antioxidant genes and susceptibility to air pollution for respiratory and cardiovascular health, Free Radical Biology and Medicine, 10.1016/j.freeradbiomed.2020.01.181, 151, (88-98), (2020).
  • A genome-wide by PM10 interaction study identifies novel loci for lung function near BICD1 and IL1RN-IL1F10 genes in Korean adults, Chemosphere, 10.1016/j.chemosphere.2019.125581, 245, (125581), (2020).
  • From Genetic Epidemiology to Exposome and Systems Epidemiology, Applied Genomics and Public Health, 10.1016/B978-0-12-813695-9.00002-9, (11-35), (2020).
  • Pathway analysis of a genome-wide gene by air pollution interaction study in asthmatic children, Journal of Exposure Science & Environmental Epidemiology, 10.1038/s41370-019-0136-3, 29, 4, (539-547), (2019).

View Options

View options

PDF

View PDF

Restore your content access

Enter your email address to restore your content access:

Note: This functionality works only for purchases done as a guest. If you already have an account, log in to access the content to which you are entitled.

Figures

Tables

Share

Share

Copy the content Link

Share on social media