Abstract
The aim of this study was to determine the frequencies of single-nucleotide polymorphisms (SNPs) in the Smad3 gene in the Chinese population and their possible association with chronic obstructive pulmonary disease (COPD). The frequency distribution of nine SNPs in the introns of the Smad3 gene was determined in both the COPD patients and control groups by the TaqMan polymerase chain reaction method using a minor groove binder probe. The genotype distribution of the rs28683050 polymorphism in control subjects was as follows: CC, 67.6%; CT, 27.0% and TT, 5.4%; and in COPD patients was as follows: CC, 47.0%; CT, 38.8% and TT, 14.2%. Therefore, a significant difference in allelic and genotypic frequencies between the COPD patients and control subjects is observed. Their genotypic distribution was not significantly different from that predicted by Hardy–Weinberg's equilibrium. The frequency of the TT genotype in the COPD patients was significantly higher than that in controls (14.2% versus 5.4%, odds ratio [OR] = 3.762, 95% confidence interval: 1.650–8.581, P = 0.002), and the frequency of the T allele in the COPD patients was significantly higher than that in controls (81.1% versus 66.4%, OR = 1.774, 95% confidence interval: 1.354–2.324, P = 0.001). These data show that the rs28683050 polymorphism of the Smad3 gene may be associated with COPD.
Introduction
Cigarette smoking is the major risk factor for chronic obstructive pulmonary disease (COPD); however, only approximately 15% of smokers develop clinically relevant airflow obstruction. 1 Multiple studies in diverse populations have shown a large genetic contribution to the variability in pulmonary function and for the familial aggregation of COPD patients. 2,3 As expected, segregation analysis suggests that multiple genes may be involved. At present, however, only a single gene, α1-antitrypsin, a potent inhibitor of inflammatory cell protease in the lung, has been unequivocally implicated in the development of COPD. 4 Recently, the associations between COPD and polymorphisms in several other genes have been studied, which include α1-antichymotrypsin, 5 microsomal epoxide hydrolase, 6 vitamin D-binding protein, 7 tumor necrosis factor-α, 8 SERPINE2 9 and a-nicotinic acetylcholine receptor (CHRNA 3/5). 10,11
Transforming growth factor-β (TGF-β) plays an important role in the pathogenesis of COPD. Signals from the activated TGF-β receptor complex are transduced to the nucleus of airway cells by SMAD proteins, which represent a family of transcription factors that have recently been implicated to play a major role as intracellular mediators of inflammation. The Smad family consists of the receptor-regulated Smads, a common pathway Smad, and inhibitory Smads. Receptor-regulated Smads are phosphorylated by the TGF-β type I receptor. They include Smad2 and Smad3, which are recognized by TGF-β and activin receptors, and Smads 1, 5, 8 and 9, which are recognized by bone morphogenetic protein receptors. Smad4 is a common pathway of Smad, which is also defined as co-operating Smad and is not phosphorylated by the TGF-β type I receptor. Inhibitory Smads (anti-Smads) include Smad6 and Smad7, which downregulate TGF-β signaling. To date, the Smads are the only TGF-β receptor substrates with a demonstrated ability to propagate signals, and with regard to the growing number of investigations of Smad-mediated effects in the airways, Smads may prove to be an important target for future development of new therapeutic strategies for asthma and COPD.
The human Smad3 gene has been assigned to chromosome 15q21–22 and spans over 5.7 kb with nine exons and eight introns. 12–14 So far, data on the association between DNA polymorphisms in the Smad3 gene and COPD have not been reported. The current study was to screen for DNA sequence variants in the Smad3 gene in Chinese Hans – which account for 95% of the Chinese population – and to determine the association between the polymorphisms and COPD (Table 1).
Nine SNPs in the Smad3 gene in the study
SNP, single-nucleotide polymorphism
To examine the interaction between a Smad3 gene polymorphism and COPD, we carried out a case-control study.
Materials and methods
Subjects
Two hundred and nineteen patients with COPD were recruited from West China Hospital of Sichuan University. The definition of COPD was consistent with that in the American Thoracic Society (ATS) consensus statement. The patients had a history of chronic or recurrent productive cough for >2 y and decreased maximum expiratory flow that had been slowly progressive and irreversible. The presence of other lung or cardiac diseases as the cause of patient symptoms was excluded by clinical and radiographic examinations. The criteria of enrollment were as follows: (1) individuals with a forced expiratory volume in 1 s (FEV1) <70% of predicted, an FEV1/forced vital capacity (FVC) ratio of 70% and an increase in FEV1 of <12% 15 min after the inhalation of 400 μg Fenoterol HBr MDI (Berotec; Boehringer Ingelheim, Ridgefield, CT, USA); and (2) patient consent to participate in the study. Most of the patients were receiving oral methylxanthine, inhaled anticholinergic agents and inhaled β2-agonists as needed. The patient's name, age, sex, family history, smoking habits, the number of cigarettes smoked, the duration of diseases and chest radiographic findings were recorded. Pulmonary function testing (CHESTAC-33-8800, Tokyo, Japan) was performed according to the ATS performance requirements. One hundred and forty-eight unrelated, age-matched healthy subjects, who had no known medical illness or family disorders and were taking no medications, acted as control subjects. This study was approved by the ethics committee of West China Hospital, Sichuan University, and signed informed consent forms were obtained from all subjects.
Candidate single-nucleotide polymorphism selection and polymorphism genotyping
The single-nucleotide polymorphisms (SNPs) of Smad3 used for the test are shown in Figure 1.Their chromosomal location are illustrated in Figure 1. SNPs of Smad3 were extracted from the database of SNP (

A schematic overview of the nine investigated single-nucleotide polymorphisms in the Smad3 gene
Genotyping method
Genomic DNA was extracted from the buffy coat using a QIAamp DNA Blood Kit (Qiagen Inc, Stanford, CA, USA). Analysis of genetics polymorphisms SNP genotyping was performed by TaqMan allelic discrimination assay. 15 Reagents were purchased from Applied Biosystems (Foster City, CA, USA). TaqMan probes were designed and synthesized by Applied Biosystems, and distinguished the SNPs at the end of a polymerase chain reaction (PCR). One allelic probe was labeled with the fluorescent FAM dye and the other with the fluorescent VIC dye. PCR was performed by TaqMan Universal Master Mix without UNG (Applied Biosystems) with PCR primers at concentrations of 200 nmol/L. Reactions were performed in 384-well formats in a total reaction volume of 3 μL using 3.0 ng of genomic DNA.
The PCR amplification protocol for the TaqMan assays included denaturation at 95°C for 10 min, followed by 40 cycles at 92°C for 15 s, 60°C for 1 min and 72°C for 45 s, followed by elongation at 72°C for 5 min. The TaqMan assays were then read on a 7900HT Fast Real-Time PCR System and alleles were called using the SDS software (Applied Biosystems).
Statistical analysis
Hardy–Weinberg analyses were performed using the Hardy–Weinberg equilibrium (HWE) program of the LINKUTIL software package. 16 Data analyses were performed with the Statistical Package for the Social Science 13.0 (SPSS Inc., Chicago, IL, USA), and the significance level for statistical tests was taken to be 0.05. The two-sided Student's t-test was used for checking significant differences in clinical data between the COPD patients and the control subjects. Differences between the patients with COPD and the controls with respect to the allele frequencies and genotype distributions were analyzed by the χ 2 test or the Fisher's exact test when necessary. Haplotype frequencies for pairs of alleles, as well as χ 2 values for allele associations, were estimated by SHEsis software. 17 Linkage disequilibrium (LD) coefficients D′ = D/D max were calculated by SHEsis software, 17 and D′ > 0.70 was taken to represent a significant difference.
Results
General characteristics
The study population consisted of 219 patients with COPD and 148 control subjects; its characteristics are described in Table 2. All the individuals included in the study were from southwestern China. The COPD cases and control subjects did not significantly differ in sex, age and smoking history characteristics. The parameters used for FEV1 and FEV1/FVC were significantly decreased in the COPD case subjects compared with the controls (P < 0.01).
Description of study population
FEV1, forced expiratory volume in 1 s; FVC, forced vital capacity
Distribution of the nine SNPs of the Smad3 gene in COPD patients and controls
To determine the prevalence of the nine SNPs of the Smad3 gene, we screened them in all 219 COPD patients and in the 148 controls. Genotype and allele frequencies by case-control status are given in Table 3. There was no deviation from HWE in controls and COPD patients (see Table 3).
Frequency distributions of the Smad3 gene in patients with COPD versus controls
COPD, chronic obstructive pulmonary disease; HWE, Hardy–Weinberg equilibrium; ND, not determined
χ 2 1 is the HWE χ 2 of the controls, χ 2 2 is the HWE χ 2 of the COPD patients
In our study, four known SNPs in the Smad3 gene, including rs28719801, rs28627002, rs28669671 and rs7166015, were not found in the COPD patients and control subjects, indicating that they are rare in the Chinese Han population.
Allelic or genotypic frequency distributions of rs3825977, rs2278546, rs2289259 and rs10152593 polymorphism in Smad3 gene between the COPD patients and controls are shown in Table 3. No significant difference of allelic or genotypic frequencies of rs3825977, rs2278546, rs2289259 and rs10152593 polymorphism in the Smad3 gene between the COPD patients and control subjects were observed.
Rs28683050 has a significant difference in allelic or genotypic frequencies between the COPD patients and control subjects.
The genotype frequency of CC, CT and TT was 67.6%, 27.0% and 5.4% in controls, respectively, and 7.0%, 38.8% and 14.2% in COPD patients, respectively (χ 2 = 57.38, P = 0.01). The allele frequency of T and C was 81.1% and 18.9% in controls and 66.4% and 33.6% in COPD patients. (χ 2 = 18.93, P < 0.05).
The frequency of the TT genotype in the COPD patients was significantly higher than that in controls (14.2% versus 5.4%, odds ratio [OR] = 3.762, 95% confidence interval [CI]: 1.650–8.581, P = 0.002). The frequency of the T allele in the COPD patients was significantly higher than that in controls (81.1% versus 66.4%, OR = 1.774, 95% CI: 1.354–2.324, P = 0.001).
Relationship between cigarette smoking and the distribution of the Smad3 gene polymorphisms
The allele and genotype distributions were compared in smokers (COPD patients and controls) to determine whether the prevalence of different alleles or genotypes was associated with smoking. Results showed that no significant difference of allelic or genotypic frequencies of the nine SNPs polymorphism in the Smad3 gene between the COPD patients and control subjects were observed (Table 4).
Frequency distributions of the Smad3 gene in smoked COPD versus smoked controls
COPD, chronic obstructive pulmonary disease; ND, not determined
Dominant and recessive genetic models
We applied both dominant and recessive genetic models in these statistical analyses. Our analyses of Fisher's exact test showed that all SNPs in this study in the Smad3 gene were not significantly associated with COPD in the recessive model and the dominant model (Table 5).
Fisher's exact test in each SNP examined in this study
SNP, single-nucleotide polymorphism; COPD, chronic obstructive pulmonary disease; OR, odds ratio; CI, confidence interval; ND, not determined
LD in nine SNPs of the Smad3 gene
Haplotype analysis, testing associations using several polymorphisms, sometimes demonstrates genetic influences that are not detected by the analysis of single polymorphisms. The relation in the nine SNPs of the Smad3 gene and their effects on COPD were analyzed. The extent of D in pairwise combinations of alleles in different variation sites was estimated by means of maximum likelihood from the frequency of diploid genotypes in the COPD and control groups. Haplotype frequencies and the coefficient of LD (D′) are given in Table 6 and Figure 1. There was no strong disequilibrium in these SNPs of the Smad3 gene.
Disequilibrium statistics in nine SNPs of the Smad3 gene
SNP, single-nucleotide polymorphism
We used estimated haplotype frequencies higher than 0.01 to construct phased multilocus genotypes for the Smad3 gene (Table 7). In Table 7 we found significant difference in some haplotypes of the Smad3 gene, e.g. ccacccaac (P < 0.05, OR = 0.249 [95% CI = 0.132–0.469]), ccacccagt (P < 0.05, OR = 2.768 [95% CI = 1.234–6.209]), ccactcaat (P < 0.05, OR = 0.331 [95% CI = 0.131–0.838]), ccatccagc (P < 0.05, OR = 0.419 [95% CI = 0.213–0.826]) and ccattcaat (P < 0.05, OR = 0.417 [95%CI = 0.229–0.759]). These results suggest that these haplotypes may be associated with COPD.
Estimate of pairwise haplotype frequencies in nine SNPs of the Smad3 gene
SNP, single-nucleotide polymorphism; CI, confidence interval
The SNP order of constructing haplotype was as follows: rs28683050, rs28669671 rs7166015, rs3825977, rs2289259, rs28627002, rs28719801, rs2278546, rs10152593
Discussion
In COPD, a marked fibrosis of the airways is observed, which is partly mediated by TGF-β1 18 that has been proposed to be a potent fibrogenic factor in numerous diseases. 19 The effects of the extensively studied TGF-β1 include fibroblast proliferation and an increased production and deposition of extracellular matrix proteins. Enzymes involved in extracellular matrix degradation and destruction can also be stimulated by TGF-β1. 19–21 While much interest has focused on the TGF-β1-induced effects in COPD, little attention has been paid to the intracellular signal transduction pathways. The pathways of the TGF-β family include the family of transcription factors termed SMADs that are directly activated by the receptors. The most predominant member of the TGF superfamily in COPD, TGF-β1, mediates its effects by the regulatory SMAD2 and SMAD3 proteins, which are phosphorylated by the type I receptor after binding of TGF-β1. The phosphorylated SMADs bind the common SMAD4 and translocate into the nucleus. In contrast, the inhibitory SMAD6 and SMAD7 22,23 function as a negative feedback loop of the TGF-β-signal by either competing for SMAD4 or terminating the signal and therefore may represent important regulators of TGF-β1-induced effects.
The data above confirmed the central role of TGF-β and the Smad pathway for proceeding to progressive fibrosis. In these animals, the administration of interleukin-1β caused extensive inflammation and tissue damage in both the Smad3 null as well as wild-type mice. If anything, the Smad3 null mouse showed an even greater degree of inflammation. However, after 20 d, only the wild-type animal had progressed to fibrosis, whereas the Smad3 null mouse had no indication of fibrogenesis and was not markedly different from animals receiving only control vector. 24 The Smad3 null mouse has been used to demonstrate that, in the absence of Smad3, the ALK5 receptor cannot transmit messages through to the nucleus and thus cannot upregulate matrix gene expression in either bleomycin- or TGF-β gene-mediated fibrotic stimulation. 25,26 Taken together, these data indicate that only when there is an intact TGF-βRI signaling mechanism and intact Smad3 signaling can the progressive nature of fibrosis proceed.
In large genome-wide association studies on COPD, Dr Silverman and Spitz identified two SNPs at the CHRNA 3/5 locus, identified earlier as a risk factor for both lung cancer and nicotine dependence, to be associated with COPD. Smad3 plays an important role in the pathogenesis of COPD, and the CHRNA 3/5 gene and the Smad3 gene are both located in chromosome 15q. 10,11 Taking these backgrounds into account, we hypothesize that SNPs in the Smad3 may render susceptibility to COPD, and thereby to contribute in predicting its occurrence during the course of the disease. Although introns were originally believed to be non-functional because they do not code for proteins, it has been suggested that some of these sequences do indeed have relevance. Introns have also been implicated in regulating gene expression and DNA–protein interactions; mutations in intron sequences may affect these functions. 27 Otsuka et al. 28 performed the tag SNP analysis using 111 atopic dermatitis (AD) families (384 members) in Japan. Thirty SNPs to cover the Smad3 gene region were chosen for genotyping using an Illumina GoldenGate assay. One intron SNP (rs4147358) was associated with AD in the pedigree disequilibrium test (PDT) and transmission disequilibrium test (TDT) analysis (significance level of P < 0.05). The rs4147358 was further genotyped in independent case-control samples and it was found that rs4147358 was associated with AD; the association was strengthened when the AD samples were limited to those of patients with IgE ≥250 IU/mL. In another case-control study, three intron SNPs of the Smad3 gene (rs2289790, rs2289791 and rs3825977) were found to be in association with diffuse cutaneous scleroderma in a UK Caucasian population. 29 Just like COPD, the fibrosis induced by transforming growth factor ss-Smads signaling also plays a major role for pathogenesis of AD and diffuse cutaneous scleroderma. So in this study, we have investigated nine intron SNPs of the Samd3 gene in COPD patients and controls. And we only found rs28683050, as a known SNP in intron 1 of the Smad3 gene, has a significant difference in allelic or genotypic frequencies between the COPD patients and control subjects. Its frequencies of the TT genotype and T allele in the COPD patients were significantly higher than that in controls. Rs2289259 was not found to be in association with COPD, which was also not associated with AD in the PDT and TDT analysis. 28 Contrary to the study by Pushpakom et al., 29 rs3825977 was not associated with COPD in our study, and the different ethnic backgrounds may be one of the reasons.
In summary, we have performed a comprehensive investigation of Smad3 gene polymorphisms. The SNP showing significant differences may be useful for predicting COPD susceptibility or possible preventive intervention. No previous research reported the SNPs in the current study have a function, and the SNPs are not in spliced junction regions. The SMAD3 SNP (rs28683050) found to be associated with COPD in this study has no known function, and it is likely that the SNPs are in LD with other functional SNPs that determine disease susceptibility. In the current study, no extron SNP of SMAD3 was included. Considering they may also render susceptibility to COPD, further studies need to be conducted to clarify the association between extron SNPs of SMAD3 and COPD and the functional characterization of these SNPs.
Footnotes
ACKNOWLEDGEMENTS
This study was supported by Grant 30971327 from the National Natural Science Foundation of China; and Grants 00-722 and 06-834 from the China Medical Board of New York to Dr FQW.
