Research Article Issue 1 | January 2015 | Vol. 123
Modification of the Association between PM10 and Lung Function Decline by Cadherin 13 Polymorphisms in the SAPALDIA Cohort: A Genome-Wide Interaction Analysis
Medea Imboden,1,2 Ashish Kumar,1,2,3 Ivan Curjuric,1,2 Martin Adam,1,2 Gian Andri Thun,1,2 Margot Haun,4 Ming-Yi Tsai,1,2 Marco Pons,5 Robert Bettschart,6 Alexander Turk,7 Thierry Rochat,8 Nino Künzli,1,2 Christian Schindler,1,2 Florian Kronenberg,4 and Nicole M. Probst-Hensch1,2
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; http://dx.doi.org/10.1289/ehp.1307398
Address correspondence to M. Imboden, SwissTPH, Socinstr. 59, 4002 Basel, Switzerland. Telephone: 41 61 284 82 40. E-mail: email@example.com
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.
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).
The authors declare they have no actual or potential competing financial interests.
Received: 20 July 2013
Accepted: 13 August 2014
Advance Publication: 15 August 2014
Final Publication: 1 January 2015
- Supplemental Material (1.7 MB) PDF
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).
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 PM10–CDH13 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).
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).
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 (n = 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.
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 (n = 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.
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.
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.
Ackermann-Liebrich U, Kuna-Dibbert B, Probst-Hensch N, Schindler C, Felber Dietrich D, Zemp Stutz E, et 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.
Andreau K, Leroux M, Bouharrour A. 2012. Health and cellular impacts of air pollutants: from cytoprotection to cytotoxicity. Biochem Res Int 2012:493894; doi: 10.1155/2012/493894.
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.
Brook RD, Rajagopalan S, Pope CA III, Brook JR, Bhatnagar A, Diez-Roux AV, et al. 2010. Particulate matter air pollution and cardiovascular disease: an update to the scientific statement from the American Heart Association. Circulation 121:2331–2378.
Chung CM, Lin TH, Chen JW, Leu HB, Yang HC, Ho HY, et al. 2011. A genome-wide association study reveals a quantitative trait locus of adiponectin on CDH13 that predicts cardiometabolic outcomes. Diabetes 60:2417–2423.
Curjuric I, Imboden M, Nadif R, Kumar A, Schindler C, Haun M, et 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; doi: 10.1371/journal.pone.0040175.
Dastani Z, Hivert MF, Timpson N, Perry JR, Yuan X, Scott RA, et 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; doi: 10.1371/journal.pgen.1002607.
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; doi: 10.1016/j.jaci.2010.09.041.
Hancock DB, Artigas MS, Gharib SA, Henry A, Manichaikul A, Ramasamy A, et al. 2012. Genome-wide joint meta-analysis of SNP and SNP-by-smoking interaction identifies novel loci for pulmonary function. PLoS Genet 8:e1003098; doi: 10.1371/journal.pgen.1003098.
Hancock DB, Eijgelsheim M, Wilk JB, Gharib SA, Loehr LR, Marciante KD, et al. 2010. Meta-analyses of genome-wide association studies identify multiple loci associated with pulmonary function. Nat Genet 42:45–52.
Heath SC, Gut IG, Brennan P, McKay JD, Bencko V, Fabianova E, et al. 2008. Investigation of the fine structure of European populations with applications to disease association studies. Eur J Hum Genet 16:1413–1429.
Huang W, Wang G, Lu SE, Kipen H, Wang Y, Hu M, et 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.
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.
Imboden M, Bouzigon E, Curjuric I, Ramasamy A, Kumar A, Hancock DB, et al. 2012. Genome-wide association study of lung function decline in adults with and without asthma. J Allergy Clin Immunol 129:1218–1228.
Imboden M, Nieters A, Bircher AJ, Brutsche M, Becker N, Wjst M, et al. 2006. Cytokine gene polymorphisms and atopic disease in two European cohorts. (ECRHS-Basel and SAPALDIA). Clin Mol Allergy 4:9; doi: 10.1186/1476-7961-4-9.
Imboden M, Schwartz J, Schindler C, Curjuric I, Berger W, Liu SL, et 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; doi: 10.1289/ehp.0800430.
Katanoda K, Sobue T, Satoh H, Tajima K, Suzuki T, Nakatsuka H, et 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.
Künzli N, Kuna-Dibbert B, Keidel D, Keller R, Brändli O, Schindler C, et al. 2005. Longitudinal validity of spirometers—a challenge in lung function follow-up studies. Swiss Med Wkly 135(33–34):503–508.
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; doi: 10.1186/1476-069X-12-101.
Liu LJ, Curjuric I, Keidel D, Heldstab J, Künzli N, Bayer-Oglesby L, et al. 2007. Characterization of source-specific air pollution exposure for a large population-based Swiss cohort (SAPALDIA). Environ Health Perspect 115:1638–1645; doi: 10.1289/ehp.10177.
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; doi: 10.1186/1743-8977-9-43.
Martin BW, Ackermann-Liebrich U, Leuenberger P, Künzli N, Stutz EZ, Keller R, et 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.
Mehta AJ, Adam M, Schaffner E, Barthélémy JC, Carballo D, Gaspoz JM, et al. 2012. Heart rate variability in association with frequent use of household sprays and scented products in SAPALDIA. Environ Health Perspect 120:958–964; doi: 10.1289/ehp.1104567.
Melén E, Nyberg F, Lindgren CM, Berglind N, Zucchelli M, Nordling E, et 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; doi: 10.1289/ehp.11117.
Morisaki H, Yamanaka I, Iwai N, Miyamoto Y, Kokubo Y, Okamura T, et al. 2012. CDH13 gene coding T-cadherin influences variations in plasma adiponectin levels in the Japanese population. Hum Mutat 33:402–410.
Obeidat M, Wain LV, Shrine N, Kalsheker N, Soler Artigas M, Repapi E, et al. 2011. A comprehensive evaluation of potential lung function associated genes in the SpiroMeta general population sample. PLoS One 6:e19382; doi: 10.1371/journal.pone.0019382.
Okada Y, Kamatani Y, Takahashi A, Matsuda K, Hosono N, Ohmiya H, et 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.
Peters KE, Beilby J, Cadby G, Warrington NM, Bruce DG, Davis WA, et 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; doi: 10.1186/1471-2350-14-15.
Romieu I, Ramirez-Aguilar M, Sienra-Monge JJ, Moreno-Macías H, del Rio-Navarro BE, David G, et al. 2006. GSTM1 and GSTP1 and respiratory health in asthmatic children exposed to ozone. Eur Respir J 28:953–959.
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.
Rossi AP, Watson NL, Newman AB, Harris TB, Kritchevsky SB, Bauer DC, et 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.
Schikowski T, Ranft U, Sugiri D, Vierkötter A, Brüning T, Harth V, et 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; doi: 10.1186/1465-9921-11-113.
Schikowski T, Schaffner E, Meier F, Phuleria HC, Vierkötter A, Schindler C, et al. 2013. Improved air quality and attenuated lung function decline: modification by obesity in the SAPALDIA cohort. Environ Health Perspect 121:1034–1039; doi: 10.1289/ehp.1206145.
Sohns M, Viktorova E, Amos CI, Brennan P, Fehringer G, Gaborieau V, et 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.
Takemura Y, Ouchi N, Shibata R, Aprahamian T, Kirber MT, Summer RS, et al. 2007. Adiponectin modulates inflammatory reactions via calreticulin receptor-dependent clearance of early apoptotic bodies. J Clin Invest 117:375–386.
Thun GA, Ferrarotti I, Imboden M, Rochat T, Gerbase M, Kronenberg F, et al. 2012. SERPINA1 PiZ and PiS heterozygotes and lung function decline in the SAPALDIA cohort. PLoS One 7:e42728; doi: 10.1371/journal.pone.0042728.
Thyagarajan B, Jacobs DR Jr, 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; doi: 10.1186/1465-9921-11-176.
Uski OJ, Happo MS, Jalava PI, Brunner T, Kelz J, Obernberger I, et 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.
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; doi: 10.1186/2046-4053-1-55.
Wilk JB, Chen TH, Gottlieb DJ, Walter RE, Nagle MW, Brandler BJ, et al. 2009. A genome-wide association study of pulmonary function measures in the Framingham Heart Study. PLoS Genet 5:e1000429; doi: 10.1371/journal.pgen.1000429.
Williams AS, Kasahara DI, Verbout NG, Fedulov AV, Zhu M, Si H, et al. 2012. Role of the adiponectin binding protein, T-cadherin (Cdh13), in allergic airways responses in mice. PLoS One 7:e41088; doi: 10.1371/journal.pone.0041088.
Wu Y, Li Y, Lange EM, Croteau-Chonka DC, Kuzawa CW, McDade TW, et 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.
Yang IA, Holz O, Jörres RA, Magnussen H, Barton SJ, Rodríguez S, et al. 2005. Association of tumor necrosis factor-α polymorphisms and ozone-induced change in lung function. Am J Respir Crit Care Med 171:171–176.
Come visit with EHP’s Science Editor, Jane Schroeder, at the 2016 ISEE Meeting in Rome, Italy, from 1–4 September 2016. This is a great opportunity to meet a member of our team, and to learn more about the journal. A number of our Associate Editors will be in attendance as well. Jane is also co-hosting a half-day writing and publishing workshop on 4 September, following the conference.
EHP is pleased to present the abstracts from the 28th annual meeting of the International Society for Environmental Epidemiology (ISEE), held in Rome, Italy, 1–4 September 2016, and hosted by the Department of Epidemiology Lazio Regional Health Service, ASL Roma 1, and the Italian Epidemiological Association. The focus of this year’s conference is current and future challenges in exposure assessment, study design, and data analyses.
Featured Children’s Health
Birgit Claus Henn, Adrienne S. Ettinger, Marianne R. Hopkins, Rebecca Jim, Chitra Amarasiriwardena, David C. Christiani, Brent A. Coull, David C. Bellinger, and Robert O. Wright
Diane Gilbert-Diamond, Jennifer A. Emond, Emily R. Baker, Susan A. Korrick, and Margaret R. Karagas
Benjamin B. Green, Margaret R. Karagas, Tracy Punshon, Brian P. Jackson, David J. Robbins, E. Andres Houseman, and Carmen J. Marsit