SERPINA1 PiZ and PiS Heterozygotes and Lung Function Decline in the SAPALDIA Cohort

Abstract

Background

Severe alpha1-antitrypsin (AAT) deficiency is a strong risk factor for COPD. But the impact of gene variants resulting in mild or intermediate AAT deficiency on the longitudinal course of respiratory health remains controversial. There is indication from experimental studies that pro-inflammatory agents like cigarette smoke can interact with these variants and thus increase the risk of adverse respiratory health effects. Therefore, we tested the effect of the presence of a protease inhibitor (Pi) S or Z allele (PiMS and PiMZ) on the change in lung function in different inflammation-exposed subgroups of a large, population-based cohort study.

Methodology and Principal Findings

The SAPALDIA population includes over 4600 subjects from whom SERPINA1 genotypes for S and Z alleles, spirometry and respiratory symptoms at baseline and after 11 years follow-up, as well as proxies for inflammatory conditions, such as detailed smoking history, obesity and high sensitivity C-reactive protein (hs-CRP), were available. All analyses were performed by applying multivariate regression models. There was no overall unfavourable effect of PiMS or PiMZ genotype on lung function change. We found indication that PiZ heterozygosity interacted with inflammatory stimuli leading to an accelerated decline in measures in use as indices for assessing mild airway obstruction. Obese individuals with genotype PiMM had an average annual decline in the forced mid expiratory flow (ΔFEF25-75%) of 58.4 ml whereas in obese individuals with PiMZ it amounted to 92.2 ml (p = 0.03). Corresponding numbers for persistent smokers differed even more strongly (66.8 ml (PiMM) vs. 108.2 ml (PiMZ), p = 0.005). Equivalent, but less strong associations were observed for the change in the FEV1/FVC ratio.

Conclusions

We suggest that, in addition to the well established impact of the rare PiZZ genotype, one Z allele may be sufficient to accelerate lung function decline in population subgroups characterized by elevated levels of low grade inflammation.

About the Authors:

Gian-Andri Thun

Affiliations Swiss Tropical and Public Health Institute, Basel, Switzerland, University of Basel, Basel, Switzerland

Ilaria Ferrarotti

Affiliation: Center for Diagnosis of Inherited Alpha1-antitrypsin Deficiency, Institute for Respiratory Disease, IRCCS San Matteo Hospital Foundation, University of Pavia, Pavia, Italy

Medea Imboden

Affiliations Swiss Tropical and Public Health Institute, Basel, Switzerland, University of Basel, Basel, Switzerland

Thierry Rochat

Affiliation: Division of Pulmonary Medicine, University Hospital of Geneva, Geneva, Switzerland

Margaret Gerbase

Affiliation: Division of Pulmonary Medicine, University Hospital of Geneva, Geneva, Switzerland

Florian Kronenberg

Affiliation: Division of Genetic Epidemiology, Department of Medical Genetics, Molecular and Clinical Pharmacology, Innsbruck Medical University, Innsbruck, Austria

Pierre-Olivier Bridevaux

Affiliation: Division of Pulmonary Medicine, University Hospital of Geneva, Geneva, Switzerland

Elisabeth Zemp

Affiliations Swiss Tropical and Public Health Institute, Basel, Switzerland, University of Basel, Basel, Switzerland

Michele Zorzetto

Affiliation: Center for Diagnosis of Inherited Alpha1-antitrypsin Deficiency, Institute for Respiratory Disease, IRCCS San Matteo Hospital Foundation, University of Pavia, Pavia, Italy

Stefania Ottaviani

Affiliation: Center for Diagnosis of Inherited Alpha1-antitrypsin Deficiency, Institute for Respiratory Disease, IRCCS San Matteo Hospital Foundation, University of Pavia, Pavia, Italy

Erich W. Russi

Affiliation: Pulmonary Division, University Hospital of Zurich, Zurich, Switzerland

Maurizio Luisetti

Affiliation: Center for Diagnosis of Inherited Alpha1-antitrypsin Deficiency, Institute for Respiratory Disease, IRCCS San Matteo Hospital Foundation, University of Pavia, Pavia, Italy

Nicole M. Probst-Hensch

* E-mail: [email protected]

Affiliations Swiss Tropical and Public Health Institute, Basel, Switzerland, University of Basel, Basel, Switzerland

Introduction

Reduced lung function measures are the most common diagnostic parameter to detect airway obstruction which is the main characteristic of asthma and chronic obstructive pulmonary disease (COPD). The only functionally characterized gene variants associated with airway obstruction and increased loss of lung function are infrequent polymorphisms in the SERPINA1 gene causing deficiency of alpha1-antitrypsin (AAT) [1]. This antiprotease inhibits neutrophil elastase, an enzyme that degrades pulmonary elastic fibers. Homozygosity for the protease inhibitor deficiency variant Z (PiZZ, also referred to as severe AAT deficiency) and compound heterozygosity for both deficiency variants S and Z (PiSZ) are widely accepted risk factors for airway obstruction and accelerated lung function decline, particularly among smokers [2]. But since these allele combinations all have frequencies below 0.1% in the general European population [3], they only account for 2-5% of all COPD cases. The more prevalent heterozygous genotypes PiMS and PiMZ (M stands for the wildtype allele) reduce the AAT blood levels only slightly [4] and are therefore referred to as mild (for PiMS) and intermediate (for PiMZ) AAT deficiencies.

While PiMS is generally believed not to be associated with low lung function or a higher risk of COPD [5], the evidence for PiMZ remains unclear even in the light of a meta-analysis [6]. The few population-based longitudinal studies have not shown adverse health effects, but they varied with regard to the phenotype studied and the inclusion of gene-environment interactions [7], [8], [9]. An investigation restricted to smokers showed that PiMZ was overrepresented in the group with rapid FEV1 (forced expiratory volume in one second) decline, suggesting that susceptibility may be refined to population subgroups with elevated inflammatory and proteolytic stress in the lungs [10]. These processes may locally increase cleavage, as well as oxidant-induced inactivation [11] and polymerization [12] of AAT, leading to a further reduction of this enzyme in PiMZ carriers. Apart from inhalant triggers like smoking, systemic inflammation may also compromise pulmonary health [13], [14] and therefore particularly affect individuals with reduced AAT levels.

We hypothesized therefore that the intrapulmonary anti-proteolytic capacity in people with mild or intermediate AAT deficiency may not be sufficient to counterbalance an excess of inflammatory triggers targeting the respiratory system. We used the SAPALDIA cohort (Swiss Cohort Study on Air Pollution and Lung and Heart Diseases in Adults) to test such possibly unfavourable effects of PiS or PiZ heterozygosity on the longitudinal course of lung function over 11 years of follow-up in the general population. The large and well characterized study population allowed us to particularly study subgroups exposed to elevated local airway or systemic inflammatory conditions, such as active and passive smokers, and people suffering from obesity.

Results

Study Sample

A comparison of characteristics between SAPALDIA follow-up participants included and excluded in this study is provided in Table 1. The study population consisted of a healthier sample with a higher percentage of never-smokers and fewer obese individuals. A comparison between the three genotype classes showed highly significant differences in AAT serum levels (Table 1). Compared to PiM homozygotes, unadjusted AAT serum concentrations were more than 16% and 38% reduced in PiS and PiZ heterozygotes, respectively. There were no differences in circulating high-sensitivity C-reactive protein (hs-CRP), the main marker for systemic inflammation, between the genotype groups. Furthermore, the genotype classes differed slightly in the unadjusted forced expiratory flow over the middle half of FVC (FEF25-75%) at baseline and in unadjusted declines of FEF25-75%.

[Figure omitted. See PDF.]

Table 1. Characteristics of SAPALDIA follow-up participants at baseline (1991; upper part), follow-up (2002; middle part), and between the two examinations (lower part).

https://doi.org/10.1371/journal.pone.0042728.t001

Adjusted Spirometric Decline Rates according to Genotype

In adjusted models, neither PiMS nor PiMZ subjects exhibited statistically significant steeper annual declines than PiMM individuals in any measure of lung function (all p≥0.07, Table 2 and Table 3). Hypothesizing that PiMS and PiMZ carriers can only compensate their reduced anti-proteolytic capacity in pulmonary tissue if no excess of pulmonary or systemic inflammation is present, we tested if stratification by smoking or obesity status may alter these associations. While we could not find any significant association between the presence of a S or Z allele and ΔFEV1 or ΔFVC (forced vital capacity) irrespective of the smoking or obesity category (Table 2), smokers and obese individuals with PiMZ genotype showed elevated declines in FEF25-75% (Table 3). In ever smokers, PiMZ carriers lost on average additional 17.4 ml per year compared to PiMM in FEF25-75% (p = 0.05), and this difference became more pronounced in people who smoked at baseline and follow-up (41.4 ml in persistent smokers, p = 0.005). A similar pattern could be observed in obese participants (ΔFEF25-75% = 58.4 ml (PiMM) vs. 92.2 ml (PiMZ), p = 0.03). There was no such effect in never-smokers exposed to environmental tobacco smoke. Values for Δ (FEV1/FVC) consistently confirmed this trend, but associations did not reach statistical significance. For ΔFEF25-75%, the presence of a Z allele interacted statistically significant with smoking (p_interaction = 0.04 with smoking status (persistent vs. never) and p_interaction = 0.002 with packyears between the two surveys), but not with obesity status (p = 0.14). Statistically significant modification of the Z allele effect on Δ (FEV1/FVC) could be observed for packyears between the two surveys and obesity status (p_interaction = 0.04 and 0.08, respectively). As we had previously found sex differences in the association between circulating AAT concentrations and lung function [15], we further evaluated a possible effect modification by sex. However, sex did not significantly modify the allele effect on any lung function measure (all p_interaction ≥ 0.15).

[Figure omitted. See PDF.]

Table 2. Adjusted mean values in ΔFEV1 and ΔFVC over 11 years of follow-up comparing different SERPINA1 genotypes.

https://doi.org/10.1371/journal.pone.0042728.t002

[Figure omitted. See PDF.]

Table 3. Adjusted mean values in Δ (FEV1/FVC) and ΔFEF25-75% over 11 years of follow-up comparing different SERPINA1 genotypes.

https://doi.org/10.1371/journal.pone.0042728.t003

High-sensitivity C-reactive Protein (hs-CRP) as a Proxy for Low Systemic Inflammation

In order to strengthen the hypothesis that low grade systemic inflammation drove the observed Z allele effect on pulmonary function, we investigated the associations in participants with elevated hs-CRP levels at follow-up (defined as the upper tertile of the study population, i.e.≥1.8 mg/l). Consistent with our findings for smokers and obese people, declines in PiZ heterozygotes were significantly enhanced for FEF25-75% (99.3 ml (PiMZ) vs. 71.2 ml (PiMM), p = 0.006) and suggestively for FEV1/FVC (5.46‰ (PiMZ) vs. 4.00‰ (PiMM), p = 0.07) (Table 3), but not for FEV1 or FVC. Interestingly, also PiS heterozygotes in the upper tertile of hs-CRP values showed statistically significant larger declines in FEF25-75% (p = 0.003) and FEV1/FVC (p = 0.04) than PiMM subjects.

Incidence of Airway Obstruction and Respiratory Symptoms

As a next step we tried to expand our analyses on the change of spirometry measures with clinically relevant outcomes like the incidence of airway obstruction or the development of respiratory symptoms. 103 PiMZ carriers did not show airway obstruction at baseline which was defined by FEV1/FVC < 0.7. In 17 of them a smaller ratio than 0.7 was calculated at follow-up, which classified them as incident cases. Comparing with PiMM individuals this resulted in an adjusted odds ratio of 1.11 (95% CI = 0.64 to 1.93, p = 0.70; Table 4). Corresponding investigations with individuals exposed to elevated inflammatory conditions did not show an increased obstruction risk for Z heterozygotes, but were based on very few numbers of cases. As far as the development of respiratory symptoms like regular cough, phlegm or shortness of breath during sleep was concerned, there seemed to be a trend that PiMZ individuals were more susceptible than PiMM carriers, in particular in subgroups exposed to pro-inflammatory conditions (Table 4). But the limited statistical power did unfortunately not allow a comprehensive evaluation of our hypothesis.

[Figure omitted. See PDF.]

Table 4. Adjusted odds ratios for developing airway obstruction (FEV1/FVC<0.7) and respiratory symptoms over 11 years of follow-up comparing different SERPINA1 genotypes.

https://doi.org/10.1371/journal.pone.0042728.t004

Validity

As chronic asthma is suggested to be associated with accelerated loss of pulmonary function [16] and in the absence of post-bronchodilation spirometry, we conducted a sensitivity analysis by including only subjects who did not report physician-diagnosed asthma. Furthermore, level of lung function change in adults is also determined by lung function growth during early adulthood, and to account for that, we performed another sensitivity analysis by only including participants older than 30 years of age at baseline. Neither of the two restrictions did essentially alter the results (Table S1).

Replacing the change of FEF25-75% by the change of the ratio FEF25-75%/FVC led to very similar conclusions (Table S2). Finally, in order to detect potential participation bias, we weighted each observation inverse to the probability of being included in the study sample. None of the results of the regression analyses did materially change (Table S3).

Discussion

In the present study, neither PiS nor PiZ heterozygosity influenced longitudinal lung function measured by ΔFEV1 or ΔFVC, independent of smoking or obesity status. However, PiMZ genotype was associated with an accelerated FEF25-75% decline in smoking and obesity subgroups from the general population. Results for participants in the upper tertile of hs-CRP values strengthened the notion that PiMZ carriers might be more susceptible to systemic pro-inflammatory conditions with respect to lung function parameters indicating narrowing of small airways.

The Role of Smoking and Pulmonary Oxidative Stress

It is well established that smokers suffering from severe AAT deficiency, a condition accompanied by only 15% of normal AAT blood concentrations, are particularly vulnerable to developing early onset COPD [2]. There is less evidence that intermediate or even mild AAT deficiency modify the effect of smoking or other oxidative inhalants on lung function. In a small study of 56-year-old men, a higher mean annual decrease in FEV1 in smoking PiMZ individuals was reported as compared with non-smoking PiMZ or smoking PiMM individuals [17], but larger studies could not find such an interaction [7], [9]. Results of studies investigating the impact of environmental and occupational exposure are heterogeneous [18]. Passive smoking has been shown in school children to be associated with cross-sectional lower lung function only in children having low levels of AAT in the blood, and particularly in measures of mid- to end-expiratory flow rates [19]. A recent longitudinal study found accelerated lung function declines in New York City firefighters of PiMZ genotype compared to PiMM in the years after World Trade Center collapse accompanied by massive air pollution, and no such difference could be observed prior to September 11, 2001 [20].

Support for a gene-environment interaction comes also from genome-wide as well as from experimental studies. While genome-wide association studies on COPD and lung function have not found the SERPINA1 locus among the top hit signals [21], [22], [23], it was the most strongly associated candidate gene in ever smokers in a comprehensive evaluation of potential lung function associated genes in more than 20,000 individuals from the general population [24]. Oxidation of the Z form by cigarette smoke induced its polymerization in lung tissue of transgenic PiZZ mice [12]. Such polymers could also be detected on the lung surface in human PiZ homozygotes, originating either from lung epithelial cells or macrophages and independent of the main hepatocytic secretion [25]. In addition to being an ineffective elastase inhibitor, these polymers attracted a higher number of neutrophils, further shifting the protease-antiprotease equilibrium towards proteolysis. As subjects with PiMZ genotype also exhibit elevated levels of intrapulmonary neutrophilic inflammation [26], they are potentially more susceptible to oxidative inhalants. The observed dose-response trend in the interaction between smoking and the presence of a Z allele is in agreement with this biological concept.

The Role of Obesity and Systemic Low Grade Inflammation

The observed accelerated loss of pulmonary function in obese subjects carrying the PiMZ genotype is novel. The often observed inverse association between obesity and lung function is in part explained by a mechanical effect of obesity on lung volume and airway caliber. Yet more limited evidence suggests an additional effect of obesity on peripheral airway obstruction that may be related to systemic low grade inflammation [27]. Adipose tissue from subjects who are overweight or obese produces pro-inflammatory adipokines that spill over to the blood stream. Elevated circulating concentrations of inflammation markers such as hs-CRP or interleukin-6 (IL-6) are both, higher in obese persons [28], [29] and associated with accelerated lung function decline [30], [31]. Also in the SAPALDIA cohort, we have previously reported that accelerated lung function decline and obesity are both associated with increased hs-CRP, particularly in women [32]. Mendelian randomization approaches to examine causality between increased hs-CRP and lower respiratory function led to conflicting results [30], [33]. However, these studies did not assess the effect of CRP gene variants on lung function separately in groups exposed to pro-inflammatory agents such as smokers.

This Study in Context

A longitudinal study in smokers reported an overrepresentation of PiMZ participants in rapid lung function decliners [10]. However, population-based cohort studies did not find consistent PiMZ effects on lung function decline, neither generally, nor in smoking strata [7], [8], [9]. The only marginally statistical significant result was found in non-smoking PiMZ carriers who showed a steeper unadjusted lung function decline than non-smoking PiMM subjects [7]. Our study was generally in good agreement with these results as we could not detect any statistical significant effect of SERPINA1 genotypes on ΔFEV1, which was the primary focus of the previously mentioned studies. Compared to those existing publications we included in our analysis a wider range of spirometry measures and found associations with the decline of FEF25-75% as well as weaker, but relatively consistent associations with the decline of the FEV1/FVC ratio. Both these two measures are in use to assess early airway obstruction and may probably best reflect the volume of the small airways [34], [35]. There are hardly any large studies which used these measures in connection with SERPINA1 genotypes so far, apart from a recent study of two large populations that found PiMZ genotypes associated with lower FEV1/FVC ratio and with more severe emphysema on chest computer tomography scan, but not with COPD status [36]. Flow related spirometric characteristics such as FEF25-75% may be decreased in the presence of airway abnormalities including inflammation or alterations in elastic recoil, two important correlates of AAT deficiency [37]. For example, interactions between glutathione S-transferase (GST) deficiency genotypes and passive smoking were strongest for mid expiratory flow measures in children [38]. However, these measures have often been criticized for being more variable and therefore less reliable than FEV1 [39]. We observed a correlation coefficient of 0.82 between baseline and follow-up FEF25-75% in SAPALDIA which is smaller than the one for FEV1 (0.92) and FVC (0.91), but larger than the one for the more commonly used FEV1/FVC ratio (0.74). Since we consistently found main and interacting effects of air pollution strongest for this mid flow parameter [40], [41], [42], it seems unlikely that the results of the present study are driven by measurement error. Moreover, FEF25-75% was found to have a high heritability in families with severe COPD [37].

Strengths and Weaknesses

The strength of this study is its large sample size, its detailed characterization of subjects, and its stringent quality control of spirometry [43]. The credibility of our results is supported by the fact that the reduction of AAT serum levels in PiMS and PiMZ compared to PiMM subjects was similar to that described by others [7]. Compared to the hitherto existing publications, we carefully excluded carriers of additional, rare mutations influencing AAT serum levels in order to diminish misclassification of wildtype alleles.

Our study has some limitations. First, a possible selection of healthy individuals may limit the generalizability of the results. However, giving more weight to underrepresented groups within the study sample did not alter the results. Moreover, if persons with low levels of lung function were preferentially lost among PiMZ carriers, stated effects may be an underestimation of the true effect. Second, PiMZ individuals showed slightly higher baseline FEF25-75% and FEV1/FVC values which can be partially explained by the younger age and the reduced number of smokers in this group, but which may question the clinical relevance of the accelerated decline in these measures. Yet, the combination of a higher level of cross-sectional lung function and a steeper lung function decline after exposure to inflammatory agents parallels observations in New York City firefighters before and after the September 11 attacks [20]. Furthermore, spirometric measurements were carried out without bronchodilator [44], which hinders a clinically acceptable definition of airway obstruction. Unlike some comparable studies [9], [10], we performed spirometry at only two time points. This makes our change values susceptible to imprecision, but we do not expect substantial measurement error for several reasons. The same spirometers and stringent quality criteria at baseline and follow-up were applied, correlation coefficients between the measurements at the two time points were high, and the direction of most genotype effects were consistent for all lung function outcomes. In addition, since it would be unlikely that any measurement error is associated with the SERPINA1 genotype, misclassification would be non-differential, which indicates an underestimation rather than an overestimation of the real effect. Another limitation is the low statistical power for analyzing obese or persistent smoking PiMZ carriers despite the large cohort size. Therefore, our findings must be interpreted with caution. If a correction for multiple testing was applied, none of the results would remain statistically significant (since all observed p-values > 0.001). For the same reason, we could neither distinguish emerging from persistent obesity, nor could we form strata with enhanced hs-CRP levels of clinical relevance (i.e.>10 mg/l). Finally, the clinical meaning of the observed excess decline in PiZ heterozygotes could not be reliably estimated, as the numbers of incident cases with airway obstruction or respiratory symptoms were too small to investigate differences between SERPINA1 genotype classes with respect to inflammatory conditions. Nevertheless, we found indication that those subgroups are more likely to develop respiratory symptoms and we know from literature that mid expiratory flow rates have been shown to be a powerful predictor of mortality from COPD, independent of FEV1 [45].

Conclusion

We confirm in this population-based study that neither PiMS, nor PiMZ carriers have a substantial impact on longitudinal lung function. There is indication however, that the presence of one Z allele may be sufficient to accelerate loss of small airway volume and incidence of respiratory symptoms in defined population subgroups which are exposed to pro-inflammatory agents and conditions. This is a potentially relevant observation as the prevalence of PiMZ genotype is quite common in Western Europe [3] and was 2.4% in this study sample. In order to estimate the public health relevance of our findings, future studies must associate subgroups of PiMZ individuals to post-bronchodilator lung function and clinically relevant outcomes like respiratory symptoms, emphysema by chest computer tomography and hospitalization for COPD.

Materials and Methods

Ethics Statement

SAPALDIA was approved by the Swiss Academy of Medical Sciences, the supraregional ethics committee for clinical research (UREK, Project Approval Number 123/00) and the Cantonal Ethics Committees for each of the eight examination areas (Ethics commissions of the cantons Aargau, Basel, Geneva, Grisons, Ticino, Valais, Vaud and Zurich). Participants were required to give written consent before any part of the health examination was conducted either globally (for all health examinations) or separately for each investigation.

Study Population

In 1991, a random sample of 9651 adults, aged 18-60 years, from eight areas in Switzerland underwent a detailed health examination including a questionnaire about respiratory health, occupational and lifestyle exposures [46]. Participants were predominantly of European-Caucasian ethnicity and represented urban and rural areas. Eleven years later, 8047 persons were reassessed [47]. 6058 follow-up subjects provided blood samples and consented to DNA analysis. 5274 of these subjects underwent spirometry testing at baseline and follow-up. Not included in this analysis were participants with missing smoking history or body mass index (BMI) data (n = 525), subjects without valid hs-CRP (n = 18), and subjects for whom genotyping of the S or Z allele either failed (n = 6) or resulted in PiS homozygosity (n = 10), PiSZ compound heterozygosity (n = 10) or PiZ homozygosity (n = 1). Other SERPINA1 rare mutations which lower AAT blood levels were detected according to a procedure described elsewhere [48] in additional 29 samples which were also excluded. Our study sample included thus 4675 subjects.

Measurements

Spirometry was assessed according to American Thoracic Society criteria using the same spirometers in 1991 and 2002 (Sensormedics model 2200, USA) and by applying stringent quality control criteria [43]. The forced expiratory manoeuvre was obtained without bronchodilators. FEV1 and FVC had to originate from the same manoeuvre in order to provide a valid FEV1/FVC ratio. Information about the smoking history was collected by questionnaire. Passive smoking was positive if never smoking subjects gave an affirmative answer at baseline or follow-up to the question if they were exposed to environmental tobacco smoke in the 12 months prior to the examination on most days or nights. Height and weight were measured and BMI was calculated as weight divided by squared height. Incident cases of airflow obstruction were defined as persons with a FEV1/FVC ratio ≥ 0.7 at baseline and < 0.7 at follow-up and were compared to individuals without obstruction at both examinations. Incident cases of respiratory symptoms were defined as people with self-reported regular cough, phlegm or shortness of breath at follow-up, but not at baseline. They were compared with individuals without any of these symptoms at baseline and follow-up. Cough or phlegm had to be present during the day or at night on most days for as much as 3 months per year and shortness of breath had to occur during sleep in the past 12 months before the examination. Subjects who declared an asthma diagnosis by a physician at baseline or follow-up were defined as asthmatics.

Serum Analysis

AAT and hs-CRP concentrations were determined from blood serum aliquots by latex-enhanced immunoturbidimetric assays (Roche diagnostics, Germany). Lower detection thresholds for the AAT and CRP assays were 0.21 g/l and 0.1 mg/l, respectively.

Genotyping

Genotyping of SERPINA1 PiS (rs17580) and PiZ (rs28929474) polymorphisms was carried out using 5′ nuclease fluorescent real-time PCR (TaqMan Probes technology) on LightCycler480 (Roche) as described before [4]. Probes and primers are given in Table S4. Genotype distributions for PiS and PiZ were both in Hardy Weinberg equilibrium (p = 0.93, N = 6050, and p = 0.99, N = 6051, respectively).

Statistical Analysis

Statistical tests to evaluate differences in the characteristics among the different groups of genotype carriers encompassed Pearson's χ² for testing equal proportions, analysis of variance (ANOVA) for testing equal means of normally distributed continuous data, and Kruskal-Wallis for testing equal distributions of continuous data which were not normally distributed. Main effects of SERPINA1 alleles on lung function decline were assessed using multiple unconditional linear regression models adjusted for sex, age, recruiting area, smoking history (packyears at baseline and between baseline and follow-up), height, baseline BMI and BMI change between baseline and follow-up. Age and packyears between baseline and follow-up were modeled with linear and squared terms to better fit to spirometry data. Interactions between genotypes and other covariates were tested by integrating multiplicative terms in the regression models. Two-sided p-values of <0.05 (and of <0.10 for interactions) were considered as statistically significant. We performed 56 different linear regression tests (4 respiratory outcomes * 2 genotype comparisons * 7 categories). Bonferroni correction would thus lower the significance threshold to p = 0.05/56 = 0.001. However, since all analyses were hypothesis-driven and most tests not independent of each other, we decided to give the results uncorrected for multiple testing.

Logistic regression models were used to compare the odds of developing airflow obstruction or respiratory symptoms between baseline and follow-up among the SERPINA1 genotype classes. The models were adjusted for the same covariates mentioned above. All statistical analyses were performed with STATA, release 10.1 IC (STATA corporation, USA).

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Sensitivity analyses for adjusted mean values in Δ(FEV1/FVC) and ΔFEF25-75% over 11 years of follow-up comparing different SERPINA1 genotypes.

https://doi.org/10.1371/journal.pone.0042728.s001

(PDF)

Table S2.

Adjusted mean values in Δ (FEF25-75%/FVC) over 11 years of follow-up comparing different SERPINA1 genotypes.

https://doi.org/10.1371/journal.pone.0042728.s002

(PDF)

Table S3.

Adjusted mean values in lung function change over 11 years of follow-up comparing different SERPINA1 genotypes in unweighted and weighted models.

https://doi.org/10.1371/journal.pone.0042728.s003

(PDF)

Table S4.

Primers and probes for genotyping the SERPINA1 PiS and PiZ polymorphisms (rs17580 and rs28929474) using 5′ nuclease fluorescent real-time PCR (TaqMan Probes technology) on LightCycler480 (Roche).

https://doi.org/10.1371/journal.pone.0042728.s004

(PDF)

Acknowledgments

The SAPALDIA Team:

Study directorate: T Rochat (p), JM Gaspoz (c), N Künzli (e/exp), LJS Liu (exp), NM Probst Hensch (e/g), C Schindler (s).

Scientific team: JC Barthélémy (c), W Berger (g), R Bettschart (p), A Bircher (a), G Bolognini (p), O Brändli (p), C Brombach (n), M Brutsche (p), L Burdet (p), M Frey (p), U Frey (pd), MW Gerbase (p), D Gold (e/c/p), E de Groot (c), W Karrer (p), R Keller (p), B Knöpfli (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), JM Tschopp (p), A von Eckardstein (cc), E Zemp Stutz (e).

Scientific team at coordinating centers: M Adam (e/g), E Boes (g), PO Bridevaux (p), D Carballo (c), E Corradi (e), I Curjuric (e), J Dratva (e), A Di Pasquale (s), L Grize (s), D Keidel (s), S Kriemler (pa), A Kumar (g), M Imboden (g), N Maire (s), A Mehta (e), F Meier (e), H Phuleria (exp), E Schaffner (s), GA Thun (g) A Ineichen (exp), M Ragettli (exp), M Ritter (exp), T Schikowski (e), G Stern (pd), M Tarantino (s), M Tsai (exp), M Wanner (pa)

(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.

Administrative staff: C Gabriel, R Gutknecht.

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. Local fieldworkers : Aarau: S Brun, G Giger, M Sperisen, M Stahel. Basel: C Bürli, C Dahler, N Oertli, I Harreh, F Karrer, G Novicic, N Wyttenbacher. Davos: A Saner, P Senn, R Winzeler, Geneva: F Bonfils, B Blicharz, C Landolt, J Rochat. Lugano: S Boccia, E Gehrig, MT Mandia, G Solari, B Viscardi. Montana: AP Bieri, C Darioly, M Maire. Payerne: F Ding, P Danieli A Vonnez. Wald: D Bodmer, E Hochstrasser, R Kunz, C Meier, J Rakic, U Schafroth, A Walder.

Author Contributions

Conceived and designed the experiments: TR EWR ML NPH. Performed the experiments: GAT IF MI MZ SO. Analyzed the data: GAT IF MI TR NPH. Contributed reagents/materials/analysis tools: MG FK POB EZ. Wrote the paper: GAT NPH. Contributed to the interpretation of the data: IF MI TR NPH .

Citation: Thun G-A, 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(8): e42728. https://doi.org/10.1371/journal.pone.0042728

References

1. Hall IP, Lomas DA (2010) The genetics of obstructive lung disease: big is beautiful. Thorax 65: 760–761.

2. Janus ED, Phillips NT, Carrell RW (1985) Smoking, lung function, and alpha 1-antitrypsin deficiency. Lancet 1: 152–154.

3. Luisetti M, Seersholm N (2004) Alpha1-antitrypsin deficiency. 1: epidemiology of alpha1-antitrypsin deficiency. Thorax 59: 164–169.

4. Ferrarotti I, Thun GA, Zorzetto M, Ottaviani S, Imboden M, et al.. Serum levels and genotype distribution of alpha1-antitrypsin in the general population. Thorax(2012).

5. Dahl M, Hersh CP, Ly NP, Berkey CS, Silverman EK, et al. (2005) The protease inhibitor PI*S allele and COPD: a meta-analysis. Eur Respir J 26: 67–76.

6. Hersh CP, Dahl M, Ly NP, Berkey CS, Nordestgaard BG, et al. (2004) Chronic obstructive pulmonary disease in alpha1-antitrypsin PI MZ heterozygotes: a meta-analysis. Thorax 59: 843–849.

7. Dahl M, Tybjaerg-Hansen A, Lange P, Vestbo J, Nordestgaard BG (2002) Change in lung function and morbidity from chronic obstructive pulmonary disease in alpha1-antitrypsin MZ heterozygotes: A longitudinal study of the general population. Ann Intern Med 136: 270–279.

8. Wadsworth ME, Vinall LE, Jones AL, Hardy RJ, Whitehouse DB, et al. (2004) Alpha1-antitrypsin as a risk for infant and adult respiratory outcomes in a national birth cohort. Am J Respir Cell Mol Biol 31: 559–564.

9. Silva GE, Sherrill DL, Guerra S, Barbee RA (2003) A longitudinal study of alpha1-antitrypsin phenotypes and decline in FEV1 in a community population. Chest 123: 1435–1440.

10. Sandford AJ, Chagani T, Weir TD, Connett JE, Anthonisen NR, et al. (2001) Susceptibility genes for rapid decline of lung function in the lung health study. Am J Respir Crit Care Med 163: 469–473.

11. Johnson D, Travis J (1979) The oxidative inactivation of human alpha-1-proteinase inhibitor. Further evidence for methionine at the reactive center. J Biol Chem 254: 4022–4026.

12. Alam S, Li Z, Janciauskiene S, Mahadeva R (2011) Oxidation of Z alpha1-antitrypsin by cigarette smoke induces polymerization: a novel mechanism of early-onset emphysema. Am J Respir Cell Mol Biol 45: 261–269.

13. Gan WQ, Man SF, Senthilselvan A, Sin DD (2004) Association between chronic obstructive pulmonary disease and systemic inflammation: a systematic review and a meta-analysis. Thorax 59: 574–580.

14. McKeever T, Saha S, Fogarty AW (2011) The association between systemic inflammatory cellular levels and lung function: a population-based study. PLoS One 6: e21593.

15. Senn O, Russi EW, Schindler C, Imboden M, von Eckardstein A, et al. (2008) Circulating alpha1-antitrypsin in the general population: determinants and association with lung function. Respir Res 9: 35.

16. Vonk JM, Jongepier H, Panhuysen CI, Schouten JP, Bleecker ER, et al. (2003) Risk factors associated with the presence of irreversible airflow limitation and reduced transfer coefficient in patients with asthma after 26 years of follow up. Thorax 58: 322–327.

17. Eriksson S, Lindell SE, Wiberg R (1985) Effects of smoking and intermediate alpha 1-antitrypsin deficiency (PiMZ) on lung function. Eur J Respir Dis 67: 279–285.

18. Senn O, Russi EW, Imboden M, Probst-Hensch NM (2005) alpha1-Antitrypsin deficiency and lung disease: risk modification by occupational and environmental inhalants. Eur Respir J 26: 909–917.

19. von Ehrenstein OS, von Mutius E, Maier E, Hirsch T, Carr D, et al. (2002) Lung function of school children with low levels of alpha1-antitrypsin and tobacco smoke exposure. Eur Respir J 19: 1099–1106.

20. Banauch GI, Brantly M, Izbicki G, Hall C, Shanske A, et al. (2010) Accelerated spirometric decline in New York City firefighters with alpha-antitrypsin deficiency. Chest 138: 1116–1124.

21. Cho MH, Castaldi PJ, Wan ES, Siedlinski M, Hersh CP, et al. (2012) A genome-wide association study of COPD identifies a susceptibility locus on chromosome 19q13. Hum Mol Genet 21: 947–957.

22. Soler Artigas M, Loth DW, Wain LV, Gharib SA, Obeidat M, et al. (2011) Genome-wide association and large-scale follow up identifies 16 new loci influencing lung function. Nat Genet 43: 1082–1090.

23. Imboden M, Bouzigon E, Curjuric I, Ramasamy A, Kumar A, et al. (2012) Genome-wide association study of lung function decline in adults with and without asthma. J Allergy Clin Immunol 129: 1218–1228.

24. Obeidat M, Wain LV, Shrine N, Kalsheker N, Soler Artigas M, et al. (2011) A comprehensive evaluation of potential lung function associated genes in the SpiroMeta general population sample. PLoS One 6: e19382.

25. Mulgrew AT, Taggart CC, Lawless MW, Greene CM, Brantly ML, et al. (2004) Z alpha1-antitrypsin polymerizes in the lung and acts as a neutrophil chemoattractant. Chest 125: 1952–1957.

26. Malerba M, Ricciardolo F, Radaeli A, Torregiani C, Ceriani L, et al. (2006) Neutrophilic inflammation and IL-8 levels in induced sputum of alpha-1-antitrypsin PiMZ subjects. Thorax 61: 129–133.

27. Salome CM, King GG, Berend N (2010) Physiology of obesity and effects on lung function. J Appl Physiol 108: 206–211.

28. Roytblat L, Rachinsky M, Fisher A, Greemberg L, Shapira Y, et al. (2000) Raised interleukin-6 levels in obese patients. Obes Res 8: 673–675.

29. Bullo M, Garcia-Lorda P, Megias I, Salas-Salvado J (2003) Systemic inflammation, adipose tissue tumor necrosis factor, and leptin expression. Obes Res 11: 525–531.

30. Sunyer J, Pistelli R, Plana E, Andreani M, Baldari F, et al. (2008) Systemic inflammation, genetic susceptibility and lung function. Eur Respir J 32: 92–97.

31. Margretardottir OB, Thorleifsson SJ, Gudmundsson G, Olafsson I, Benediktsdottir B, et al. (2009) Hypertension, systemic inflammation and body weight in relation to lung function impairment-an epidemiological study. COPD 6: 250–255.

32. Bridevaux PO, Gerbase MW, Schindler C, Dietrich DF, Curjuric I, et al. (2009) Sex-specific effect of body weight gain on systemic inflammation in subjects with COPD: results from the SAPALDIA cohort study 2. Eur Respir J 34: 332–339.

33. Dahl M, Vestbo J, Zacho J, Lange P, Tybjaerg-Hansen A, et al. (2011) C reactive protein and chronic obstructive pulmonary disease: a Mendelian randomisation approach. Thorax 66: 197–204.

34. Cochrane GM, Prieto F, Hickey B, Benatar SR, Clark TJ (1974) Early diagnosis of airways obstruction. Thorax 29: 389–393.

35. Simon MR, Chinchilli VM, Phillips BR, Sorkness CA, Lemanske RF Jr, et al.. Forced expiratory flow between 25% and 75% of vital capacity and FEV1/forced vital capacity ratio in relation to clinical and physiological parameters in asthmatic children with normal FEV1 values. J Allergy Clin Immunol (2010) 126: :527-534 e521-528.

36. Sorheim IC, Bakke P, Gulsvik A, Pillai SG, Johannessen A, et al. (2010) alpha-Antitrypsin protease inhibitor MZ heterozygosity is associated with airflow obstruction in two large cohorts. Chest 138: 1125–1132.

37. DeMeo DL, Carey VJ, Chapman HA, Reilly JJ, Ginns LC, et al. (2004) Familial aggregation of FEF(25-75) and FEF(25-75)/FVC in families with severe, early onset COPD. Thorax 59: 396–400.

38. Kabesch M, Hoefler C, Carr D, Leupold W, Weiland SK, et al. (2004) Glutathione S transferase deficiency and passive smoking increase childhood asthma. Thorax 59: 569–573.

39. Miller MR, Hankinson J, Brusasco V, Burgos F, Casaburi R, et al. (2005) Standardisation of spirometry. Eur Respir J 26: 319–338.

40. Downs SH, Schindler C, Liu LJ, Keidel D, Bayer-Oglesby L, et al. (2007) Reduced exposure to PM10 and attenuated age-related decline in lung function. N Engl J Med 357: 2338–2347.

41. Imboden M, Schwartz J, Schindler C, Curjuric I, Berger W, 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.

42. Curjuric I, Imboden M, Schindler C, Downs SH, Hersberger M, et al. (2010) HMOX1 and GST variants modify attenuation of FEF25-75% decline due to PM10 reduction. Eur Respir J 35: 505–514.

43. Kunzli N, Kuna-Dibbert B, Keidel D, Keller R, Brandli O, et al. (2005) Longitudinal validity of spirometers–a challenge in longitudinal studies. Swiss Med Wkly 135: 503–508.

44. Probst-Hensch NM, Curjuric I, Pierre-Olivier B, Ackermann-Liebrich U, Bettschart RW, et al. (2010) Longitudinal change of prebronchodilator spirometric obstruction and health outcomes: results from the SAPALDIA cohort. Thorax 65: 150–156.

45. Thomason MJ, Strachan DP (2000) Which spirometric indices best predict subsequent death from chronic obstructive pulmonary disease? Thorax 55: 785–788.

46. Martin BW, Ackermann-Liebrich U, Leuenberger P, Kunzli N, Stutz EZ, 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.

47. Ackermann-Liebrich U, Kuna-Dibbert B, Probst-Hensch NM, Schindler C, Felber Dietrich D, 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.

48. Zorzetto M, Russi E, Senn O, Imboden M, Ferrarotti I, et al. (2008) SERPINA1 gene variants in individuals from the general population with reduced alpha1-antitrypsin concentrations. Clin Chem 54: 1331–1338.