The role of IREB2 and transforming growth factor beta-1 genetic variants in COPD: a replication case-control study
- Sally L Chappell1,
- Leslie Daly2,
- Juzer Lotya2,
- Aiman Alsaegh1,
- Tamar Guetta-Baranes1,
- Josep Roca3,
- Roberto Rabinovich7,
- Kevin Morgan1,
- Ann B Millar4,
- Seamas C Donnelly5,
- Vera Keatings6,
- William MacNee7,
- Jan Stolk8,
- Pieter S Hiemstra8,
- Massimo Miniati9,
- Simonetta Monti9,
- Clare M O'Connor†5 and
- Noor Kalsheker†1Email author
© Chappell et al; licensee BioMed Central Ltd. 2011
Received: 5 August 2010
Accepted: 14 February 2011
Published: 14 February 2011
Genetic factors are known to contribute to COPD susceptibility and these factors are not fully understood. Conflicting results have been reported for many genetic studies of candidate genes based on their role in the disease. Genome-wide association studies in combination with expression profiling have identified a number of new candidates including IREB2. A meta-analysis has implicated transforming growth factor beta-1 (TGFbeta1) as a contributor to disease susceptibility.
We have examined previously reported associations in both genes in a collection of 1017 white COPD patients and 912 non-diseased smoking controls. Genotype information was obtained for seven SNPs in the IREB2 gene, and for four SNPs in the TGFbeta1 gene. Allele and genotype frequencies were compared between COPD cases and controls, and odds ratios were calculated. The analysis was adjusted for age, sex, smoking and centre, including interactions of age, sex and smoking with centre.
Our data replicate the association of IREB2 SNPs in association with COPD for SNP rs2568494, rs2656069 and rs12593229 with respective adjusted p-values of 0.0018, 0.0039 and 0.0053. No significant associations were identified for TGFbeta1.
These studies have therefore confirmed that the IREB2 locus is a contributor to COPD susceptibility and suggests a new pathway in COPD pathogenesis invoking iron homeostasis.
Chronic obstructive pulmonary disease (COPD) is a leading cause of morbidity and mortality, and is predicted to become the 4th leading cause of death by the year 2030 . Whilst smoking is a significant environmental cause of COPD, not all smokers will develop disease. It is well recognised that COPD has a genetic component as well as environmental, and that this may account for these differences in susceptibility.
Many candidate gene studies have been carried out over the past few years, with varying degrees of reproducibility. Conflicting results may be due to population differences, spurious results caused by small sample sizes and the subsequent low power of the study to detect true associations or to variation in the phenotype. Meta-analysis can be used to pool results from genetic studies and give an overall conclusion, and this approach has been used to summarise the results for several COPD candidates . This study revealed significant associations for three polymorphisms in transforming growth factor beta 1 (TGFB1), although it was acknowledged that this was a limited analysis based on a small number of studies.
Other strategies for identifying the genetic basis of complex human disease include integrative approaches which combine gene expression data with association studies, or unbiased genome-wide approaches. Both of these strategies have also been employed within the study of COPD.
A recent genome-wide association study highlighted a region on chromosome 15q25 which shows strong association to COPD . This region contains several genes, including members of the nicotinic acetylcholine receptors (CHRNA5, CHRNA3 and CHRNB4) and the iron regulatory protein 2 (IREB2). Although this region has also been linked to lung cancer and nicotine addiction as well as COPD [4, 5], the strong levels of linkage disequilibrium make it difficult to refine the location of the true association signal and the functional variants. However, the expression of IREB2 has also been shown to be altered in lung tissue from COPD patients compared to controls , adding to the evidence for IREB2 acting as a COPD susceptibility gene.
The aim of the current study was to investigate polymorphisms in the TGFB1 and IREB2 genes in our existing large collection of COPD cases and controls.
COPD cases and control subjects were recruited at six European centres. Numbers from each centre were: Barcelona - 70 controls and 138 cases; Bristol - 152 controls and 129 cases; Dublin - 195 controls and 196 cases; Edinburgh - 81 controls and 168 cases; Leiden - 216 controls and 188 cases and Pisa - 198 controls and 198 cases. Approval for the study was obtained from the appropriate committees at each recruitment centre. Informed consent was obtained from all subjects. Criteria for patient recruitment were a firm clinical diagnosis of stable COPD; airflow limitation as indicated by FEV1 ≤70% normal predicted values and FEV1/FVC <70%; no significant reversibility on bronchodilation and a smoking history of ≥20 pack years. Patients were excluded from the study if they had an established diagnosis of asthma, lung cancer, a history of atopy, known AAT deficiency or a serum AAT level of less than 1.0 g/L. They were also excluded if they had had an acute exacerbation in the 4 weeks preceding assessment for the study. Disease severity was classified according to GOLD as follows: GOLD Stage II (moderate disease): FEV1/FVC ratio of <0.70 and FEV1 predicted ≥50% and <80%; GOLD Stage III (severe disease): FEV1/FVC ratio of <0.70 and FEV1 predicted ≥30% and <50%; and GOLD Stage IV (very severe disease): FEV1/FVC ratio of <0.70 and FEV1 predicted <30% .
Control subjects were recruited at each centre to match COPD patients for age, gender and smoking history. Exclusion criteria were as described for cases and also included a family history of COPD. Only individuals with no evidence of airflow obstruction (FEV1 and FVC ≥80% and FEV1/FVC >70%) were included in the control group. Only white Caucasians were recruited for cases and controls. Complete matching between cases and controls was not achieved, but this was taken into account during the analysis.
Genotyping of Study Samples
Four SNPs were chosen for TGFB1, which were the same as the ones included in the recent meta-analysis . These were rs2241712, rs1800469, rs1800470 (formerly known as rs1982073) and rs6957. Seven SNPs were included for IREB2, as identified by a recent publication . These were rs2568494, rs2656069, rs1964678, rs12593229, rs10851906, rs965604 and rs13180. Genotyping was carried out commercially by K-Bioscience (Hertfordshire, UK). As a quality control measure, approximately 5% of samples were genotyped in duplicate to check for concordance. In addition, a selection of samples were also genotyped using restriction enzyme digestion, allele-specific PCR or direct sequencing to confirm the genotyping results from K-Bioscience.
Statistical Analysis of Genetic variation
Each of the SNPs in the IREB2 and TGF1 genes was analyzed for Hardy-Weinberg equilibrium (HWE) using SAS/Genetics PROC ALLELE software . HWE analysis was performed on all controls and by each centre. The analysis of allele and genotype frequencies in cases and controls was performed using the same program. The p-values and odds ratios for genotype and allele frequencies are obtained by using SAS/Genetics PROC LOGISTIC software. Odds ratios for individual SNP allele distributions are relative to the common allele and for genotypes they are relative to the major homozygous genotype. The adjusted p-values and odds ratios with 95% confidence intervals for genotype and allele frequencies adjust for age, sex, smoking and centre, including interactions of age, sex and smoking with centre. This adjustment was done to eliminate residual confounding due to age, sex, smoking and centre using logistic regression because matching was not completely achieved on recruitment. Further, to assess the sensitivity of the analysis, we examined the adjusted odds ratios, for individual SNP allele distributions relative to the common allele, for each centre separately adjusting for age, sex and smoking. Pairwise linkage disequilibrium coefficients were calculated in controls using Haploview .
Because there are significant differences in allele and genotype frequencies between cases and controls for the IREB2 gene, using FAMHAP18 [10, 11] we examined all possible two SNP and three SNP haplotypes for SNPs 1, 2 and 4 to explore if the haplotypes present significantly stronger associations than single SNPs. We tested the association of major and minor alleles with disease severity, comparing each of the GOLD Stages II, III and IV with Controls for TGF1 gene (SNPs 2, 3 and 4) and for IREB2 gene (SNPs 1, 2 and 4). The chi-square test of association of Major and Minor Alleles with severity was done using SAS PROC FREQ software. An apparent trend in the IREB2 gene for SNP allele frequency with increasing severity was tested using SAS PROG GLM software with CONTRAST statement. Using PROC GLM in SAS, the quantitative trait associations between SNPs and the phenotype FEV1 were tested by multivariate regression, adjusting for age, sex, smoking and centre.
Characteristics of Controls and COPD subjects
60.8 ± 8.9
65.9 ± 8.2
Smoking Pack Years
38.6 ± 17.3
48.8 ± 22.9
Predicted FEV 1 (%)
95.3 ± 10.9
43.1 ± 15.1
FEV 1 /FVC (%)
77.8 ± 4.9
47.5 ± 12.1
Characteristics of COPD patients according to GOLD classification of disease severity
Very Severe IV
66.0 ± 7.9
66.6 ± 8.4
64.3 ± 8.0
49.5 ± 21.3
48.4 ± 22.4
48.6 ± 26.0
Predicted FEV 1 (%)
60.3 ± 5.9
39.9 ± 5.8
23.5 ± 4.3
FEV 1 /FVC (%)
57.1 ± 7.8
45.9 ± 9.9
36.2 ± 9.4
Characteristics of the SNPs
Location in Gene
Exon 1 (Pro10Leu)
Exon 21 (Ala872Ala)
TGFΒ1 and IREB2 Allele frequencies and odds ratios in Controls and Cases
(95% C. I.)
1.02 (0.87, 1.18)
0.99 (0.86, 1.15)
0.90 (0.74, 1.08)
1.30 (1.12, 1.50)
0.77 (0.64, 0.92)
0.79 (0.68, 0.91)
IREB2 Genotype frequencies and odds ratios in Controls and Cases.
Minor Allele frequencies by case, control and severity status for IREB2 SNPs 1, 2 and 4
Though there were significant differences between centres for the allele distribution of some of the SNPs, our results are unlikely to be due to population stratification since cases and controls were matched within centres. We also examined the genotype odds ratios for all six SNPs in each centre (Additional file 1: Table S3). The within-centre odds ratios for the three TGFB1 SNPs were all non-significant, showed no particular pattern and fluctuated around the overall (all-centre) odds ratios. The within-centre analysis for the IREB2 SNPs showed that, though non-significant due to the smaller sample size, the direction of the relationship was the same in each centre as the all-centre results. Thus each centre is showing the same IREB2 SNP relationships with COPD.
This study has replicated the associations seen for polymorphisms in IREB2 which have been previously reported, but failed to reproduce the associations seen in a recent meta-analysis for TGFB1. Recent studies using cluster analysis suggest that SNP rs1800470 in TGFB1 is associated with the emphysema-predominant phenotype . There is a plausible biological explanation for this based on the critical role that TGFB1 plays in reducing matrix metalloproteinase-12 activity, a key mediator of emphysematous change. This is further supported by recent replicated studies showing that MMP-12 variants are associated with a protective effect on lung function and COPD [14, 15]. The supposition is that low levels of MMP-12 reduce the likelihood of emphysematous change. Over-expression of TGFB1 would therefore predictably have the same protective effect and reduced expression would be a predisposing factor. This is also supported by work in a murine model of COPD, where activation of the TGF-Β signalling pathway was observed alongside a partial rescue of the emphysema phenotype . We did not have phenotypic data to assess emphysema scores and so were unable to test this proposition. It is also possible that other genetic variants in the TGFB1 region which were not investigated as part of the current study may be involved in COPD susceptibility.
We have confirmed previous observations which suggest that IREB2 variants may play a role in COPD. This appears to be independent of an effect on lung function, as we failed to detect any association between FEV1 and any of the SNPs included in this study. This is in agreement with the genome-wide association study which looked for variants associated with lung function, which also fails to find evidence for a role of IREB2 . The IREB2 gene is located within the region identified by the COPD GWAS, and the high levels of linkage disequilibrium in this area make it difficult to identify the specific functional variant or gene which is underlying this association. Only SNP 1 (rs2568494) is in relatively high levels of linkage disequilibrium with SNPs identified by the genome-wide study of COPD (HapMap CEU data: r2 = 0.790 with rs8034191 and r2 = 0.692 with rs1051730), although the other IREB2 SNPs are also in LD with variants in other genes in the region, including CHRNA5. Identification of the genes and variants contributing to the association with COPD will require targeted resequencing of the region and further functional work. In support of IREB2 involvement, a previous publication has shown increased levels of IREB2 mRNA in COPD patients versus controls . This is an intriguing observation as this gene codes for an iron binding protein. Iron homeostasis and free iron concentration are likely to be important mediators of oxidative stress and iron could therefore contribute to local damage by this mechanism. IREB2 protein is expressed in the lungs and cigarette smoking has been associated with higher levels of iron in the lung . IREB2 knock-out mouse models have a predisposition to developing neurodegenerative disease due to aberrant cellular iron homeostasis , though the lungs of these animals were not examined in any detail. A potential mechanism for distorted ion homeostasis relates to membrane serine proteases (matriptases) which regulate a number of biological effects. Genome-wide association studies have shown that mutations in matriptase result in refractory iron deficiency anaemia and the proteolytic activity of matriptase is critical in regulating iron in this form of anaemia as loss of function mutations result in the disease . The catalytic domain of matriptase is inhibited by alpha1-antitrypsin  and it is therefore conceivable that two major pathways associated with progressive lung damage, namely proteolysis and oxidative stress, are linked through iron homeostasis. This could also explain the tendency for greater oxidative stress in alpha1-antitrypsin deficiency as more active matriptase may result in excessive free iron.
In summary, this study failed to replicate previous reports of associations between SNPs in TGFB1 and susceptibility to COPD, despite being adequately powered. This may be due to the possibility that the association is actually with specific emphysema phenotypes rather than overall susceptibility to disease, which could not be assessed in this patient group. The replication of the association between variants of IREB2 and COPD provides further evidence to support the role of this genomic region in COPD pathogenesis, and the role of iron regulation deserves further investigation.
We thank Ann Hann (Bristol), Breda Callaghan and Gemma Hogan (Dublin), Joyce Barr (Edinburgh), Clara Kolster-Bijdevaate (Leiden) for assisting with recruitment; G. Catapano E. Fornai and C. Carli for clinical and technical assistance in Pisa. This work was supported by European Union 5th Framework Programme Grant QLG1-CT-200-01012 (COPD GENE SCAN project) and is protected as intellectual property. SCD is supported by Science Foundation Ireland (SFI); J. Lotya was supported by a grant from the Health Research Board (Ireland).
- Mathers CD, Loncar D: Projections of global mortality and burden of disease from 2002 to 2030. PLoS medicine. 2006, 3 (11): e442-10.1371/journal.pmed.0030442.View ArticlePubMedPubMed CentralGoogle Scholar
- Smolonska J, Wijmenga C, Postma DS, Boezen HM: Meta-analyses on suspected chronic obstructive pulmonary disease genes: a summary of 20 years' research. American journal of respiratory and critical care medicine. 2009, 180 (7): 618-631. 10.1164/rccm.200905-0722OC.View ArticlePubMedGoogle Scholar
- Pillai SG, Ge D, Zhu G, Kong X, Shianna KV, Need AC, Feng S, Hersh CP, Bakke P, Gulsvik A, et al: A genome-wide association study in chronic obstructive pulmonary disease (COPD): identification of two major susceptibility loci. PLoS genetics. 2009, 5 (3): e1000421-10.1371/journal.pgen.1000421.View ArticlePubMedPubMed CentralGoogle Scholar
- Amos CI, Wu X, Broderick P, Gorlov IP, Gu J, Eisen T, Dong Q, Zhang Q, Gu X, Vijayakrishnan J, et al: Genome-wide association scan of tag SNPs identifies a susceptibility locus for lung cancer at 15q25.1. Nature genetics. 2008, 40 (5): 616-622. 10.1038/ng.109.View ArticlePubMedPubMed CentralGoogle Scholar
- Spitz MR, Amos CI, Dong Q, Lin J, Wu X: The CHRNA5-A3 region on chromosome 15q24-25.1 is a risk factor both for nicotine dependence and for lung cancer. Journal of the National Cancer Institute. 2008, 100 (21): 1552-1556. 10.1093/jnci/djn363.View ArticlePubMedPubMed CentralGoogle Scholar
- DeMeo DL, Mariani T, Bhattacharya S, Srisuma S, Lange C, Litonjua A, Bueno R, Pillai SG, Lomas DA, Sparrow D, et al: Integration of genomic and genetic approaches implicates IREB2 as a COPD susceptibility gene. American journal of human genetics. 2009, 85 (4): 493-502. 10.1016/j.ajhg.2009.09.004.View ArticlePubMedPubMed CentralGoogle Scholar
- Rabe KF, Hurd S, Anzueto A, Barnes PJ, Buist SA, Calverley P, Fukuchi Y, Jenkins C, Rodriguez-Roisin R, van Weel C, et al: Global strategy for the diagnosis, management, and prevention of chronic obstructive pulmonary disease: GOLD executive summary. American journal of respiratory and critical care medicine. 2007, 176 (6): 532-555. 10.1164/rccm.200703-456SO.View ArticlePubMedGoogle Scholar
- SAS I: SAS software. 2002, SAS Institute Inc, 9.1.3Google Scholar
- Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics (Oxford, England). 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.View ArticleGoogle Scholar
- Becker T, Knapp M: A powerful strategy to account for multiple testing in the context of haplotype analysis. American journal of human genetics. 2004, 75 (4): 561-570. 10.1086/424390.View ArticlePubMedPubMed CentralGoogle Scholar
- Herold C, Becker T: Genetic association analysis with FAMHAP: a major program update. Bioinformatics (Oxford, England). 2009, 25 (1): 134-136. 10.1093/bioinformatics/btn581.View ArticleGoogle Scholar
- Sidak Z: Rectangular Confidence Regions for the Means of Multivariate Normal Distributions. Journal of the American Statistical Association. 1967, 62: 626-633. 10.2307/2283989.Google Scholar
- Cho MH, Washko GR, Hoffmann TJ, Criner GJ, Hoffman EA, Martinez FJ, Laird N, Reilly JJ, Silverman EK: Cluster analysis in severe emphysema subjects using phenotype and genotype data: an exploratory investigation. Respiratory research. 2010, 11: 30-10.1186/1465-9921-11-30.View ArticlePubMedPubMed CentralGoogle Scholar
- Hunninghake GM, Cho MH, Tesfaigzi Y, Soto-Quiros ME, Avila L, Lasky-Su J, Stidley C, Melen E, Soderhall C, Hallberg J, et al: MMP12, lung function, and COPD in high-risk populations. The New England journal of medicine. 2009, 361 (27): 2599-2608. 10.1056/NEJMoa0904006.View ArticlePubMedPubMed CentralGoogle Scholar
- Haq I, Chappell S, Johnson SR, Lotya J, Daly L, Morgan K, Guetta-Baranes T, Roca J, Rabinovich R, Millar AB, et al: Association of MMP-2 polymorphisms with severe and very severe COPD: a case control study of MMPs-1, 9 and 12 in a European population. BMC medical genetics. 2010, 11: 7-10.1186/1471-2350-11-7.View ArticlePubMedPubMed CentralGoogle Scholar
- Wempe F, De-Zolt S, Koli K, Bangsow T, Parajuli N, Dumitrascu R, Sterner-Kock A, Weissmann N, Keski-Oja J, von Melchner H: Inactivation of sestrin 2 induces TGF-beta signaling and partially rescues pulmonary emphysema in a mouse model of COPD. Disease models & mechanisms. 3 (3-4): 246-253.Google Scholar
- Repapi E, Sayers I, Wain LV, Burton PR, Johnson T, Obeidat M, Zhao JH, Ramasamy A, Zhai G, Vitart V, et al: Genome-wide association study identifies five loci associated with lung function. Nature genetics. 2010, 42 (1): 36-44. 10.1038/ng.501.View ArticlePubMedGoogle Scholar
- Ghio AJ, Hilborn ED, Stonehuerner JG, Dailey LA, Carter JD, Richards JH, Crissman KM, Foronjy RF, Uyeminami DL, Pinkerton KE: Particulate matter in cigarette smoke alters iron homeostasis to produce a biological effect. American journal of respiratory and critical care medicine. 2008, 178 (11): 1130-1138. 10.1164/rccm.200802-334OC.View ArticlePubMedGoogle Scholar
- LaVaute T, Smith S, Cooperman S, Iwai K, Land W, Meyron-Holtz E, Drake SK, Miller G, Abu-Asab M, Tsokos M, et al: Targeted deletion of the gene encoding iron regulatory protein-2 causes misregulation of iron metabolism and neurodegenerative disease in mice. Nature genetics. 2001, 27 (2): 209-214. 10.1038/84859.View ArticlePubMedGoogle Scholar
- De Falco L, Totaro F, Nai A, Pagani A, Girelli D, Silvestri L, Piscopo C, Campostrini N, Dufour C, Al Manjomi F, et al: Novel TMPRSS6 mutations associated with iron-refractory iron deficiency anemia (IRIDA). Human mutation. 2010, 31 (5): E1390-1405. 10.1002/humu.21243.View ArticlePubMedGoogle Scholar
- Janciauskiene S, Nita I, Subramaniyam D, Li Q, Lancaster JR, Matalon S: Alpha1-antitrypsin inhibits the activity of the matriptase catalytic domain in vitro. American journal of respiratory cell and molecular biology. 2008, 39 (6): 631-637. 10.1165/rcmb.2008-0015RC.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/12/24/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.