Association of transforming growth factor-β1 gene variants with risk of coal workers' pneumoconiosis ☆

Hiyng Qin, Zhifng Song, Meilin Wng, Xiomin Ji, Aiping Li, Ye Yng, Linlin Shen,Shsh Wng, Chunhui Ni*, Jinwei Zhou*

aDepartment of Occupational Medicine and Environmental Health, School of Public Health, Nanjing Medical University, Nanjing 210029, Jiangsu Province, ChinabDepartment of Occupational Health, General Hospital of Xuzhou Mining Business Group Co., Ltd. Xuzhou 221006,Jiangsu Province, China

INTRODUCTION

Pneumoconiosis is a global occupational health problem, especially in China with an annual incidence of nearly 80% of all occupational diseases. Inhaling respirable coal dust and silica particulates frequently develop coal workers' pneumoconiosis (CWP), a dustassociated pneumoconiosis characterized by chronic pulmonary inflammation and fibrotic nodular lesions that usually lead to progressive fibrosis[1]. Epidemiological and animal studies have found that silicosis may develop even when exposure was terminated prior to its initial development[2]. Although individuals may be exposed to similar levels of coal dust,only some of them develop lung fibrosis, suggesting that genetic predisposition may influence individual susceptibility to the development of lung fibrosis[3,4].Therefore, understanding genetic variability, and the interaction between genetic variation and dust exposure may aid in the identification of high-risk individuals and prevent them from developing CWP.

Transforming growth factor-β1 (TGF-β1) is an important multifunctional cytokine that modulates myriad cellular and tissue processes, including cell growth,differentiation, apoptosis, and inflammation and is involved in the pathogenesis of the lung fibrosis[5,6]. The profibrotic effects of TGF-β1 are numerous, including induction of myofibroblasts, increase of matrix synthesis, and inhibition of collagen breakdown. During lung fibrosis, TGF-β1 contributes to the influx and activation of inflammatory cells, transdifferentiation of epithelium to mesenchyme, and influx of fibroblasts and their subsequent elaboration of extracellular matrix[7,8].

Recently, studies also showed the contribution of TGF-β1 genetic variation to individual fibrosis susceptibility. Several polymorphic variants in the TGF-β1 gene were examined in pulmonary diseases,such as chronic obstructive pulmonary disease (COPD)[9], asthma[10], idiopathic pulmonary fibrosis (IPF)[11],and cystic fibrosis (CF)[12]. Further, most investigators reported a positive correlation between TGF-β1 gene variants and risk of development of these diseases.

Since many single nucleotide polymorphisms(SNPs) are in linkage disequilibrium with other nearby SNPs in the genome, it is feasible to select a small number of SNPs as representative genetic markers,also termed tagging single nucleotide polymorphisms(tSNPs), to capture the common variations in the gene. The use of tSNPs might improve the effectiveness of studies on association, because tSNPs can provide information about nearby SNPs that are not genotyped[13,14]. However, there are currently only a few reports on the role of TGF-β1 in CWP[15,16]. In the present study, five tSNPs were selected from data for Chinese subjects in the HapMap (http：//www.hapmap.org/) to evaluate the association between genetic variants in TGF-β1 and risk of CWP development in a Chinese population.

MATERIALS AND METHODS

Study population

This study consisted of 508 CWP patients and 526 controls, as described previously[17]. Briefly, all subjects were Han Chinese, underground coal miners recruited from coal mines of Xuzhou Mining Business Group Co., Ltd. between January 2006 and December 2008. The high kilovolt chest X-ray and physical examinations were performed for reconfirming the diagnoses based on the China National Diagnostic Criteria for Pneumoconiosis (GBZ 70-2002), which is the same as the 1980 International Labor Office Classification of Pneumoconioses in the judgment of opacity profusion. The controls were matched with CWP cases for age, dust exposure period and job type.Blood sample of 5 ml was obtained from all studied subjects, and used for routine lab tests, including hapetic function tests. Isolated leucocytes were used for DNA extraction. The research protocol was approved by the Institutional Review Board of Nanjing Medical University.

tSNPs identification and genotyping

TGF-β1 is located on chromosome 19q13.1-q13.3 in humans and consists of 7 exons (Fig. 1). To examine the gene extensively, this TGF-β1 gene was searched by including its flanking regions of 1,000-bp both upstream and downstream (from 23103114 to 23139491)and a total of 53 SNPs located, including 20 common SNPs [i.e., minor allele frequency (MAF)≥0.05]among a Chinese population included in the HapMap database (HapMap Data Rel 21a/PhaseⅡ, Jan07, on NCBI B35 assembly, dbSNP b125). For the genotyping, a set of tSNPs was selected in the TGFβ1 gene with the following criterion： A minimal set of haplotypes that ensured an R2h of at least 0.8 to cover all possible haplotypes that had a frequency of at least 5% as evaluated by the tSNPs program[13]. The reported SNPs in coding region (i.e., rs20577 Thr33Ile in exon1, rs6557634 His141Ar g in exon3, rs20575 Thr209Arg in exon4, rs20576 Glu228Ala in exon5)were not chosen in this study because of their low MAF (＜ 0.05) in Asian populations. As a result, five tSNPs, which could accurately predict the common (＞0.05) haplotypes with a minimum R2h of 0.817, were selected. The selected SNP IDs, locations and allele frequencies are summarized in Table 1.

The selected SNPs were genotyped in all 1,034 subjects by the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method.The tSNPs information, primers and restriction enzymes of polymorphisms are all listed in Table 2. The polymorphism analysis was performed independently by two individuals in a blinded fashion. About 1% ofPCR products were randomly selected and confirmed by sequencing (data not shown), and more than 10%of the samples were randomly selected for repeated genotyping. The results were 100% concordant.

Fig. 1 Location of five tSNPs in the TGF-β1 gene.

Table 1 Primary information on five tSNPs of the TGF-β1 gene among the CWP cases and controls

Table 2 The tSNP information, primer sequences, annealing temperatures and restriction enzymes for genotyping of TGF-β1 polymorphisms

Interviewer-administered questionnaire

Information for each person was obtained from the questionnaire, including personal characteristics (age,height, weight, education and family history), occupational history (work type, dust exposure period and personnel protection), occupational diseases (CWP stage, the date of first diagnoses, progress and hearing loss) and life-style (alcohol intake and tobacco smoking). Trained interviewers administered the questionnaires according to standard operating procedures prepared in advance.

Statistical analysis

Chi-square test was used to evaluate the differences in the frequency distributions of select demographic variables, dust exposure period, smoking status, alleles and genotypes of the TGF-β1 polymorphisms between the CWP cases and controls. Student's t-test was used to compare the differences in means of age and dust exposure period between CWP cases and controls.Hardy-Weinberg equilibrium of the controls' geno-type distributions was used to test by a goodness-offit Chi-square test. Univariate and multivariate logistic regression analysis was performed to obtain crude and adjusted ORs for risk of CWP and their 95% CIs.

RESULTS

Characteristics of the study subjects

The selected characteristics of the CWP cases and controls included in the analyses are shown in Table 3. There were no significant differences in the distribution of work types and in the means of exposure period between CWP cases and controls. The average age of CWP cases was older than that of controls (P ＜0.001). The smoking prevalence of CWP was similar to the controls (P = 0.123); however, the pack-years smoked was significantly less than that of controls since the rate of smoking cessation was significantly higher in CWP cases (P ＜ 0.001).

Table 3 Distribution of selected variables between the CWP cases and controls

Individual SNP association analysis and risk of CWP

The SNP IDs, locations and allele frequencies are shown in Table 1. SNP1 (rs1800469), SNP2(rs1800470) and SNP5(rs11466345) of TGF-β1 genotype distributions in the control population deviated significantly from those expected for a population in Hardy-Weinberg equilibrium, while the frequencies of SNP1 and SNP5 genotype in CWP cases were in agreement with Hardy-Weinberg equilibrium (P＞ 0.05). The genotype distributions of the TGF-β1 polymorphisms in cases and controls are summarized in Table 4. The single SNP analysis indicated that the genotype frequencies of SNP2 rs1800470 and SNP5 rs11466345 were significantly different between cases and controls (P = 0.003 and P = 0.017, respectively).Multivariate logistic regression analysis revealed that the SNP2 (rs1800470) CC genotype was associated with decreased risk of CWP (OR = 0.50, 95% CI =0.32-0.78) and the SNP5 (rs11466345) GG genotype was associated with an increased risk of CWP (OR =2.24, 95% CI = 1.37-3.68) . No significant association was found between the other SNPs and CWP risk.

Stratification analysis between the genotypes of SNPs in TGF-β1 and risk of CWP

As shown in Table 5, the SNP2 (rs1800470) TT genotype was significantly decreased with risk of CWP in subgroups of those never smoke (OR = 0.40, 95%CI= 0.24-0.66), cases with stageⅡ (OR = 0.41, 95%CI= 0.22-0.76) and exposure period (＜ 28 y： OR = 0.54,95%CI = 0.31-0.95; ≥ 28 y： OR=0.52, 95%CI = 0.32-0.96). However, the SNP5 (rs11466345) GG genotype was associated with an increased risk of CWP in both smoking (OR = 0.40, 95%CI = 0.24-0.66) and nonsmoking groups, shorter and longer exposure groups,while the risk of CWP was only increased in patients with stage Ⅰ and Ⅱ (Table 6 ).

DISCUSSION

In this case-control study, five tSNPs in the TGF-β1 gene were investigated with respect to an association with the risk of CWP occurrence in Chinese subjects.It was found that the SNP2 (rs1800470) CC genotype decreased the risk of CWP compared to the TT genotype and the further stratification analysis showed that those with the TT genotype had a decreased risk of CWP in subgroups of never smoking, exposure period(＜ 28 y or≥ 28 y) and stageⅡ. However, the SNP5(rs11466345) GG genotype was associated with an increased risk of CWP compared to the AA genotype;the stratification analysis showed similar trends where GG genotype increased risks of CWP in all subgroups except stage Ⅲ.

TGF-β1 was found to play a key role in the development of lung fibrosis, and increased expression of TGF-β1 occurred in lung tissue in patients with lung fibrosis and animal models of pulmonary fibrosis[8,18-20]. Yao et al[15]noted the potential role of TGF-β1 genetic variants in the pathogenesis of CWP by demonstrating that the TGF-β1 gene -509 site(rs1800469) polymorphism influenced the concentra-tion of TGF-β1 in serum, but was not different among stages and exposure period of CWP. In contrast, Wu et al[21]did not find an association between the TGF-β1 gene polymorphisms at positions -509 (rs1800469),+869 (rs1800470), and +915 (rs1800471) with silicosis risk in Chinese iron miners.

Table 5 Stratification analysis between the genotypes of SNP2 rs1800470 polymorphism and CWP risk

In the present study, we found that SNP2(rs1800470, T ＞ C) and SNP5 (rs11466345, A ＞ G)were significantly associated with CWP risk. SNP2(rs1800470) is an exon polymorphism whose T ＞ Ctransition is at nucleotide +869, which leads to the Leu10Pro amino acid substitution. Previous studies reported that TGF-β1 gene variants were important genetic modifiers of lung disease progression since the genotypes of +869T/C (rs1800470) can influence TGF-β1 plasma levels[22,23]. The SNP5 is an intronic polymorphism whose functional consequences are less intuitive, but intronic polymorphisms have been reported to be associated with a variety of chronic diseases including breast cancer[24], essential hypertension[25]and typeⅡ diabetes[26]. By analyzing the putative transcription factors for TGF-β1 from the Genomatix program (http：//www.genomatix.de), the SNP5 A locus is found to be centered in a 21 bp that is nearly identical to the Tax/CREB complex binding site, and the SNP5 G locus is in 25bp nearly identical to the arylhydrocarbon-receptor nuclear translocator (ARNT). SNP5 (rs11466345) variants may influence the affinities of Tax/CREB complex and/or ARNT with the binding sites. Although exact molecular mechanisms underlying variation in SNP5 and the ability to increase risk of CWP are unknown, it is possible that this SNP may exert an effect on gene expression and be in linkage disequilibrium with other functional variants. However, further investigation is required to confirm these hypotheses.

Table 6 Stratification analysis between the genotypes of SNP5 rs11466345 polymorphism and CWP risk

Studies on the associations between the TGF-β1 gene polymorphisms and susceptibility of lung fibrosis have been reported[11,12,27], but with few involved CWP. Yucesoy et al[16]found no association in coal miners between TGF-β1 genetic variations of promoter region (rs1800469) in patients with susceptibility to progressive massive fibrosis. In the present casecontrol study, 5 tSNPs of TGF-β1 were examined,and data showed that although the result of SNPs in promoter region was similar to the study reported by Yucesoy et al[16], the variants in exon 1 (rs1800470)and intron 5 (rs11466345) were associated with susceptibility to CWP development.

In this study, cases and controls matched very well in work type and exposure period, which were the most important matching indicators. The smoking frequency of cases was not significantly different from controls, but over 50% smoking cases stopped smoking, much higher than that of controls (13.86%), due to the breathing problems. The age was not matched very well and this needs to be improved in later studies. The frequencies of SNP3 and SNP4 genotypes in controls and the frequencies of SNP1 and SNP5 genotypes in cases were in agreement with Hardy-Weinberg equilibrium, which indicated that our controls may differ from the general population since the controls were from underground miners who had the same dust exposure as the cases, but without lung fibrosis. Therefore, a control group from the local general population is needed in a future study.

In conclusion, our study demonstrated that some representative genetic variants in TGF-β1 may exert a role in the risk of CWP development. Studies with ethnically diverse populations and functional evaluation are warranted to confirm our findings.

[1] Huang X, Finkelman RB. Understanding the chemical properties of macerals and minerals in coal and its potential application for occupational lung disease prevention. J Toxicol Environ Health B Crit Rev 2008;11：45-67.

[2] Porter DW, Hubbs AF, Mercer R, Robinson VA, Ramsey D, McLaurin J, et al. Progression of lung inflammation and damage in rats after cessation of silica inhalation.Toxicol Sci 2004;79：370-80.

[3] Hu W, Zhang Q, Su WC, Feng Z, Rom W, Chen LC,et al. Gene expression of primary human bronchial epithelial cells in response to coal dusts with different prevalence of coal workers' pneumoconiosis. J Toxicol Environ Health A 2003;66：1249-65.

[4] Danila E, Sileikiene V, Nargela R, Zurauskas E,Loskutoviene G. Different course of silicosis in four brothers of one family. Int J Occup Med Environ Health 2009;22：51-7.

[5] Bartram U, Speer CP. The role of transforming growth factor beta in lung development and disease. Chest 2004;125：754-65.

[6] Sheppard D. Transforming growth factor beta： a central modulator of pulmonary and airway inflammation and fibrosis. Proc Am Thorac Soc 2006;3：413-7.

[7] Ihn H. Pathogenesis of fibrosis： role of TGF-beta and CTGF. Curr Opin Rheumatol 2002;14：681-5.

[8] Tarantal AF, Chen H, Shi TT, Lu CH, Fang AB, Buckley S,et al. Overexpression of TGF-{beta}1 in foetal monkey lung results in Prenatal Pulmonary Fibrosis. Eur Respir J 2010. Mar.

[9] Ito M, Hanaoka M, Droma Y, Hatayama O, Sato E,Katsuyama Y, et al. The association of transforming growth factor beta 1 gene polymorphisms with the emphysema phenotype of COPD in Japanese. Intern Med 2008;47：1387-94.

[10] Silverman ES, Palmer LJ, Subramaniam V, Hallock A,Mathew S, Vallone J, et al. Transforming growth factorbeta1 promoter polymorphism C-509T is associated with asthma. Am J Respir Crit Care Med 2004;169：214-9.

[11] Xaubet A, Marin-Arguedas A, Lario S, Ancochea J,Morell F, Ruiz-Manzano J, et al. Transforming growth factor-beta1 gene polymorphisms are associated with disease progression in idiopathic pulmonary fibrosis. Am J Respir Crit Care Med 2003;168：431-5.

[12] Brazova J, Sismova K, Vavrova V, Bartosova J,Macek M, Jr., Lauschman H, et al. Polymorphisms of TGF-beta1 in cystic fibrosis patients. Clin Immunol 2006;121：350-7.

[13] Carlson CS, Eberle MA, Rieder MJ, Yi Q, Kruglyak L,Nickerson DA. Selecting a maximally informative set of single-nucleotide polymorphisms for association analyses using linkage disequilibrium. Am J Hum Genet 2004;74：106-20.

[14] Long JR, Zhao LJ, Liu PY, Lu Y, Dvornyk V, Shen H,et al. Patterns of linkage disequilibrium and haplotype distribution in disease candidate genes. BMC Genet 2004; 5：11.

[15] Yao W, Wang ZM, Wang MZ, Hao CF, Wang N. The relationship between the concentration of TGF-beta1 in serum and its gene polymorphisms in CWP. Sichuan Da Xue Xue Bao Yi Xue Ban (in Chinese) 2005;36：827-89.

[16] Yucesoy B, Johnson VJ, Kissling GE, Fluharty K,Kashon ML, Slaven J, et al. Genetic susceptibility to progressive massive fibrosis in coal miners. Eur Respir J 2008; 31：1177-82.

[17] Wang M, Ye Y, Qian H, Song Z, Jia X, Zhang Z, et al.Common genetic variants in pre-microRNAs are associated with risk of coal workers' pneumoconiosis. J Hum Genet 2010;55：13-7.

[18] Khalil N, O'Connor RN, Unruh HW, Warren PW,Flanders KC, Kemp A, et al. Increased production and immunohistochemical localization of transforming growth factor-beta in idiopathic pulmonary fibrosis. Am J Respir Cell Mol Biol 1991;5：155-62.

[19] Santana A, Saxena B, Noble NA, Gold LI, Marshall BC. Increased expression of transforming growth factor beta isoforms (beta 1, beta 2, beta 3) in bleomycininduced pulmonary fibrosis. Am J Respir Cell Mol Biol 1995;13：34-44.

[20] Williams AO, Knapton AD, Ifon ET, Saffiotti U. Transforming growth factor beta expression and transformation of rat lung epithelial cells by crystalline silica(quartz). Int J Cancer 1996;65：639-49.

[21] Wu F, Qu Y, Tang Y, Cao D, Sun P, Xia Z. Lack of association between cytokine gene polymorphisms and silicosis and pulmonary tuberculosis in Chinese iron miners. J Occup Health 2008;50：445-54.

[22] Corvol H, Boelle PY, Brouard J, Knauer N, Chadelat K,Henrion-Caude A, et al. Genetic variations in inflammatory mediators influence lung disease progression in cystic fibrosis. Pediatr Pulmonol 2008;43：1224-32.

[23] Ogawa E, Ruan J, Connett JE, Anthonisen NR, Pare PD,Sandford AJ. Transforming growth factor-beta1 polymorphisms, airway responsiveness and lung function decline in smokers. Respir Med 2007;101：938-43.

[24] Freedman ML, Penney KL, Stram DO, Le Marchand L,Hirschhorn JN, Kolonel LN, et al. Common variation in BRCA2 and breast cancer risk： a haplotype-based analysis in the Multiethnic Cohort. Hum Mol Genet 2004;13：2431-41.

[25] Sano M, Kuroi N, Nakayama T, Sato N, Izumi Y, Soma M, et al. Association study of calcitonin-receptor-like receptor gene in essential hypertension. Am J Hypertens 2005; 18：403-8.

[26] Lehman DM, Fu DJ, Freeman AB, Hunt KJ, Leach RJ,Johnson-Pais T, et al. A single nucleotide polymorphism in MGEA5 encoding O-GlcNAc-selective N-acetyl-beta-D glucosaminidase is associated with type 2 diabetes in Mexican Americans. Diabetes 2005;54：1214-21.

[27] Grutters JC, du Bois RM. Genetics of fibrosing lung diseases. Eur Respir J 2005;25：915-27.