Abstract
E-cadherins form intercellular junctions that maintain epithelial integrity. Epithelial integrity is impaired in asthma and can be restored by inhaled corticosteroids (ICSs). Our aim was to investigate the association of CDH1 gene polymorphisms (single-nucleotide polymorphisms (SNPs)) with airway remodelling, inflammation and forced expiratory volume in 1 s (FEV1) decline in asthma patients and assess whether ICSs modulate these effects.
Bronchial biopsies of 138 asthmatics were available (population 1). Associations of 17 haplotype-tagging SNPs with epithelial E-cadherin expression, biopsy parameters and FEV1/vital capacity (VC) ratio were tested. FEV1 and VC data were collected in 281 asthmatics with 30-yr follow-up (population 2). Linear mixed-effect models were used to assess associations of SNPs with FEV1 decline.
Seven out of the 17 SNPs were associated with airway remodelling, three with CD8+ T-cell counts, two with eosinophil counts and seven with FEV1 decline. All associations occurred only in patients using ICS. In general, alleles associated with less remodelling correlated with less FEV1 decline and higher FEV1/VC. Decreased epithelial E-cadherin expression was associated with five SNPs in non-ICS users.
In conclusion, our data show that CDH1 polymorphisms are associated with epithelial E-cadherin expression and suggest that epithelial adhesion is an important contributor to airway remodelling and lung function in asthma. These effects are modified by the use of inhaled corticosteroids.
Asthma is a chronic inflammatory disease of the airways characterised by bronchial hyperresponsiveness (BHR), and pathological changes in the epithelium and submucosal area of the bronchi. Histological studies indicate that the epithelium is interrupted and fragile in bronchial biopsies from asthmatics [1, 2]. Both genetic and environmental factors as well as their interactions are involved in asthma development and its course during lifetime. Recent advances in the understanding of the pathophysiology of asthma have pointed towards a prominent role of the airway epithelium in asthma development and severity [3, 4].
Airway epithelium is the natural barrier between the environment and underlying tissue. Its integrity is important for protecting the airways against noxious inhalants, such as environmental tobacco smoke, particles and biological agents. Regulation of cell–cell junction stability and dynamics is crucial to maintaining tissue integrity and allowing tissue remodelling throughout development [5, 6]. E-cadherin, the major cadherin expressed in epithelial cells, is an adhesion molecule that is able to establish and stabilise cellular junctions between adjacent cells in the presence of Ca2+ [7] and is coded by the CDH1 gene on chromosome 16q22.1. Hirosako et al. [8] suggested that human bronchial intraepithelial lymphocytes have roles distinct from subsets of other lymphocytes, and that CD8+ cells and CD103+ (ligand of E-cadherin) cells have potentially important functions in the bronchial epithelium.
Longitudinal studies have demonstrated that some asthma patients develop irreversible airway obstruction and progressive loss of lung function [9]. Inhaled corticosteroids (ICSs) have been widely used for suppressing inflammation and asthma symptoms. Early use of ICSs has been found to be associated with less accelerated forced expiratory volume in 1 s (FEV1) decline [10], and with improvement of BHR and airway remodelling [11]. Furthermore, corticosteroids are thought to modulate epithelial repair [12] and there is evidence that they upregulate epithelial E-cadherin expression [13].
We hypothesised that loss of epithelial integrity resulting from loss of epithelial E-cadherin expression or disruption of E-cadherin-mediated cell–cell contacts predisposes airways to abnormal responses to inhaled substances. Subsequent remodelling processes and airway inflammation may lead to airway narrowing and accelerated lung function decline.
Therefore, we investigated the associations between CDH1 polymorphisms and airway remodelling (basement membrane (BM) thickness and subepithelial vasculature), inflammation and annual FEV1 decline in asthma patients and assessed whether ICS use modulates these associations.
MATERIALS AND METHODS
Populations
Population 1
Population 1 consisted of 138 asthmatics who had participated in cohort studies conducted by our research group in previous years. They were re-examined for lung function and BHR, and underwent bronchoscopy with biopsy collection [14, 15]. The main exclusion criteria were: FEV1 <1.2 L; bronchiectasis; upper respiratory tract infection (e.g. colds); and/or use of antibiotics or oral corticosteroids within the 2 months before inclusion. DNA samples and biopsy data were available in 137 subjects.
Population 2
Population 2 consisted of 281 asthma patients who had been referred to the asthma clinic at Beatrixoord Hospital in Haren, the Netherlands, over the period 1966–1975. Patients who, at their first visit, were <45 yrs of age, had ≥20% fall in FEV1 during a histamine challenge test (provocative concentration causing a 20% fall in FEV1 ≤32 mg·mL−1) and were symptomatic according to current American Thoracic Society criteria [16] were re-examined during the period 1991–1999 [9]. At both visits, patients filled in a questionnaire on respiratory symptoms and underwent lung function, histamine challenge and skin-prick tests. The clinical assessment was performed as previously described [10]. Before testing, participants had to be in a stable condition without any exacerbation in the previous 6 weeks. After their first visit, subjects had annual routine check-ups for their asthma. Their medical records provided information on lung function and corticosteroid use during these check-ups. DNA samples have been collected from 253 subjects. 32 subjects in population 1 were derived from population 2. All participants originated from the northern region of the Netherlands. The study was approved by the medical ethics committee of the University Medical Center Groningen (Groningen, the Netherlands), and all participants gave signed, informed consent.
Biopsy collection and processing
Bronchoscopy, collection and processing of the bronchial biopsies were performed as previously described [14]. Immunohistochemistry was performed on 3-μm, formalin-fixed, paraffin-embedded tissue specimens that were deparaffinised with xylene, dehydrated in ethanol and after antigen retrieval, incubated with the primary antibodies: anti-E-cadherin antibody (#610181; BD Bioscience, Breda, the Netherlands) for epithelial E-cadherin expression; anti-eosinophilic peroxidase antibody for eosinophil detection (laboratories of N.A. Lee and J.J. Lee, Mayo Clinic, Scottsdale, AZ, USA); anti-CD8 (DAKO, Heverlee, the Netherlands) antibody for T-cells; and anti-CD31 monoclonal antibody for vessel endothelial cells. BM thickness was calculated by dividing the BM surface area by BM length (μm). The number of positively stained inflammatory cells was counted in a total area of 0.1 mm2 in the submucosa, 100 μm beneath the BM. The number of CD31+ vessels in the submucosal area was counted in the whole section (excluding epithelium, muscle and mucous gland areas); therefore, we measured the number of vessels per area (0.1 mm2). Finally, E-cadherin expression was determined as the percentage of BM covered with E-cadherin-positive intact epithelium (fig. 1). Further details on the immunochemistry and quantification procedures used are presented in the online supplementary material.
Immunostaining of airway wall biopsies of asthma patients. Representative images of a) strong and b) weak CD31 staining of vessels, and c) strong, d) moderate and e) weak epithelial E-cadherin staining. Magnification: objective 20×. Scale bars = 150 μm.
Genotyping of single-nucleotide polymorphisms
We genotyped 17 haplotype-tagging single-nucleotide polymorphisms (SNPs) in CDH1 according to HapMap CEU (Utah residents with Northern and Western European ancestry) genotype data (r2<0.8; minor allele frequency >10%). 12 of the SNPs were genotyped by KBioscience Ltd (Hoddesdon, UK) using their competitive allele-specific PCR system (KASPar), and five SNPs were derived from a genome-wide association study on asthma using the 370-kb Illumina chip (Illumina, Eindhoven, the Netherlands) (see figure 2 for the genotyped SNPs).
Linkage disequilibrium plot and correlation coefficients (100r2) of 17 genotyped CDH1 gene single-nucleotide polymorphisms (SNPs) in population 2 (n=281). The location of SNPs is given for the HapMap (www.hapmap.org).
Statistics
Normalisation of the distribution of the variables was performed by natural logarithm transformation where necessary. Multiple linear regression analysis adjusted for sex, age and smoking was used to assess the effect of the SNPs on airway remodelling parameters (BM thickness and submucosal vessel numbers) and inflammatory cells (eosinophils and CD8+T-cells) in population 1. Additionally, we used FEV1/vital capacity (VC) ratio post-bronchodilation (BD) adjusted for sex, age, height and smoking as a marker of remodelling in both populations [17]. Since E-cadherin expression is known to be upregulated by ICSs [13, 18, 19], we tested for interaction between ICS and genotypes by introducing interaction terms into the models. General (heterozygotes versus wild types and homozygous mutants versus wild types) genetic models were used. When the number of subjects with the homozygote mutant genotype was less than five, only dominant (heterozygotes and homozygous mutants pooled versus wild types) models were used. Association between epithelial E-cadherin expression (percentage of BM covered with intact epithelium positive for E-cadherin) and genotypes was assessed with nonparametric tests using dominant genetic models. The association between denuded epithelium (percentage BM covered with denuded epithelium) and epithelial E-cadherin expression was tested with multiple linear regression adjusted for sex and age. In population 2, FEV1 decline was analysed with linear mixed-effect (LME) models as described previously [10]. Only FEV1 measurements after the age of 30 yrs were included, since this is the age where the lung function is in a plateau phase or decline begins in normal subjects [20]. LME models were adjusted for sex, age, first FEV1 after age 30 yrs, height, pack-years of smoking and use of oral corticosteroids. Inclusion criteria required subjects having FEV1 measurements before and after the introduction of ICS, with a minimum of three FEV1 measurements over a period of ≥2 yrs. Subjects who never used ICSs were excluded. FEV1 levels and annual FEV1 declines were estimated for the two periods (without and with ICS use) by introducing the following variables: time; genotype; ICS use; and their interactions. LME models were conducted in S-plus 7.0 (Insightful Corp., Seattle, WA, USA). All other analyses were conducted with SPSS (version 16; SPSS Inc., Chicago, IL, USA). Two-tailed p-values of <0.05 were considered statistically significant.
RESULTS
Clinical characteristics and genotype distribution
Clinical characteristics of both populations are presented in table 1. All subjects were atopic and hyperresponsive at first clinical assessment. Population 1 had a higher mean lung function and 57 (41.3%) subjects used ICS, while 146 (52%) subjects were on ICS treatment in population 2. For the analysis on FEV1 decline 125 subjects met the inclusion criteria. Table 2 presents details on corticosteroid use in populations 1 and 2. Genotype distributions in the two asthmatic populations are presented in table E1 in the online supplementary material. All SNPs were in Hardy-Weinberg equilibrium (p>0.05). The linkage disequilibrium (LD) plot of the SNPs in population 2 is shown in figure 2.
Epithelial E-cadherin expression
Five CDH1 SNPs were significantly associated with epithelial E-cadherin expression (rs8056633, rs16958383, rs2276330, rs3785078 and rs7203904). The minor alleles of these SNPs were significantly associated with lower expression of E-cadherin in individuals without ICS and generally with similar expression as the wild types in those with ICS use (table 3 and fig. 3; table E2 in the online supplementary material). There were no significant differences in epithelial E-cadherin expression between never- and ever- (ex- and current pooled) smokers or among the three smoking categories (data not shown). The association between denuded epithelium with epithelial E-cadherin expression was very close to significance (b=-0.003; p=0.05; natural logarithm scale).
CDH1 rs16958383 associations with epithelial E-cadherin expression, airway remodelling, inflammation and lung function decline in adult asthma patients. CDH1 rs16958383 represents single-nucleotide polymorphisms associated with airway remodelling (rs7188750, rs8056633, rs16958383, rs2276330, rs3785878, rs7203904 and rs17690554). General linear models were used to estimate the means of each dependent variable adjusted for sex, age and smoking. a) Epithelial E-cadherin expression, expressed as % basement membrane (BM) covered with E-cadherin-positive intact epithelium. •: subjects; —: median. b) CD31+ vessels per 0.1 mm2 submucosal area. c) BM thickness. d) Post-bronchodilation (BD) forced expiratory volume in 1 s (FEV1)/vital capacity (VC) ratio. e) Annual change in FEV1. f) CD8+ T-cells per 0.1 mm2 submucosal area. g) Eosinophils per 0.1 mm2 submucosal area. •: mean; whiskers: 95% CI. ICS: inhaled corticosteroid. Bold indicates significance.
Airway remodelling: BM thickness and subepithelial vasculature
In population 1, seven out of 17 SNPs were associated with BM thickness and subepithelial vasculature in the presence of ICS (rs7188750, rs8056633, rs16958383, rs2276630, rs3785078, rs7203904 and rs17690554). The minor alleles of these polymorphisms were all associated with thinner BM and/or a lower number of vessels in the submucosal area (table 3 and fig. 3). Without ICS use, these effects of the minor alleles were less prominent or even reversed, but not significantly so. Interactions between these SNPs and ICSs with respect to airway remodelling were significant for rs2276330 and rs3785078 (tables E3.1 and E3.2 in the online supplementary material).
Airway inflammation: submucosal CD8+ T-cell and eosinophil counts
With ICS use, rs11075699, rs2276630 and rs1125557 were significantly associated with CD8+ T-cell counts in the submucosal area. CD8+ T-cell counts were higher in subjects heterozygous for rs11075699 and in carriers of the minor C-allele of rs2276330, while subjects with a homozygous mutant genotype for rs1125557 had lower CD8+T-cell counts compared with wild-types. Minor alleles of rs2902185 and rs10431924 were associated with higher eosinophil counts. These associations were absent in asthmatics without ICS. Interactions between SNPs and ICSs on inflammatory cell counts (CD8+ T-cells and eosinophils) were statistically significant for three out of five SNPs (rs1125557, rs10431924 and rs2276330) and for rs7188750, respectively (fig. 3; tables E4.1 and E4.2 in the online supplementary material).
FEV1 decline and FEV1/VC
Seven out of 17 SNPs were significantly associated with FEV1 decline in the presence of ICSs. The minor alleles for rs8056633, rs16958383, rs7203904 and rs17690554 were associated with less FEV1 decline. For example, during ICS use, FEV1 decline was 20.4 mL·yr−1 for subjects with the wild-type genotype of rs16958383, while it was 36.6 mL·yr−1 less in homozygous carriers of the minor allele (p=0.004; fig. 3). The minor alleles for rs1125557, rs7199991 and rs7186053 were associated with accelerated FEV1 decline. A significant interaction of ICS with CDH1 SNPs on FEV1 decline was present only for rs7199991 and rs3785078. There were no associations with lung function decline in the absence of ICS use (table 3; table E5 in the online supplementary material).
In population 2, SNPs that were associated with less FEV1 decline were also associated with higher post-BD FEV1/VC with ICS use. Additionally, rs4783573, rs7188750 and rs3785078 were associated with higher post-BD FEV1/VC (table 3 and fig. 3). Interactions of these SNPs with ICSs on post-BD FEV1/VC were significant. In population 1, these SNPs had the same direction of association with post-BD FEV1/VC, but only the interaction with rs7203904 was significant (tables E6.1 and E6.2 in the online supplementary material).
DISCUSSION
Our results indicate that E-cadherin (CDH1) gene polymorphisms are associated with airway remodelling, inflammation and lung function decline in individuals with asthma. In summary, with ICS use, seven out of the 17 SNPs were associated with airway remodelling, three with CD8+ T-cells counts and two with eosinophil counts in the submucosa, and seven with FEV1 decline. Consistently, alleles associated with less airway remodelling correlated with less FEV1 decline and higher post-BD FEV1/VC. While all associations occurred during ICS use, associations with epithelial E-cadherin expression were significant for five SNPs in the absence of ICSs. These results may indicate that ICSs may influence the way these polymorphisms express their effect and are consistent with in vitro results showing that administration of glucocorticosteroids can upregulate cadherin expression, improve cell–cell contact and strengthen the epithelial barrier [19, 21, 22].
Altered expression of epithelial junction proteins on bronchial epithelium may contribute to a defective epithelial barrier [1, 23] with subsequent functional and clinical manifestations. In our longitudinal study (population 2), the minor alleles of three SNPs (rs1125557, rs7199991 and rs7186053) were associated with accelerated FEV1 decline during ICS use. rs7199991 is in complete LD (r2 = 1 according to www.hapmap.org, version 4.0) with a promoter-region SNP (rs16260), suggesting that this may have functional significance. To date, E-cadherin gene polymorphisms have not been studied in asthma, but recent findings have shown associations of CDH1 SNPs with Crohn's disease [24], another impaired epithelial barrier disorder, which may share pathways with asthma [25, 26]. Moreover, in genome-wide association studies, CDH1 has recently been associated with susceptibility to ulcerative colitis [27] and colorectal cancer [28]. The new associations suggest that changes in the integrity of the intestinal epithelial barrier may contribute to the pathogenesis of these diseases; hence, CDH1 polymorphisms are important for understanding the pathogenesis of impaired epithelial barrier diseases.
Additionally, our study showed that the minor alleles of four other SNPs (rs8056633, rs16958383, rs7203904 and rs17690554) were associated with less FEV1 decline. The latter four alleles were also associated with less airway remodelling in biopsies of asthmatics using ICSs in population 1. The fact that we observed a significant decrease in BM thickness in the presence of ICSs is consistent with the findings of Broekema et al. [15], who previously showed that symptomatic asthmatics using ICS in our population 1 had significantly thinner BMs than asymptomatic asthmatics who did not use ICSs (5.7 and 6.5 μm, respectively; p<0.05). However, our results highlight that this decrease in BM thickness is probably not present in all asthmatics using ICSs, but only in those with a specific CDH1 genotype.
Rasmussen et al. [17] used post-BD FEV1/VC as an indirect measure of airway remodelling, and our results point in the same direction, i.e. alleles associated with less airway remodelling (rs7188750, rs8056633, rs16958383, rs3785078, rs7203904 and rs17690554) were also associated with higher post-BD FEV1/VC and less FEV1 decline in population 2. In population 1, only rs7203904 showed a significant association with FEV1/VC, which is probably due to lower power: lower number of ICS users (57 (41%) versus 146 (52%) subjects) and less post-BD FEV1/VC variance (sd 6.5 versus 13.5) compared with population 2. In contrast, Kosciuch et al. [29] concluded that there is no relationship between BM thickness and lung function tests, such as FEV1, FEV1/VC, residual volume % predicted and total lung capacity % predicted, but they included only asthmatic patients who never used ICS or who were withdrawn from ICSs for ≥3 months. Our study suggests that discrepant findings may, at least partially, be due to a modulating effect of both ICSs and CDH1 SNPs on airway remodelling and lung function. As previous studies have shown, treatment with ICSs may reduce or control the intensity of airway remodelling, improve FEV1 and reduce loss of lung function [30–32]. It has yet to be elucidated how ICSs prevent these adverse effects. Based on our findings and on evidence that corticosteroids upregulate cadherin expression [18, 33], we speculate that ICS interaction with CDH1 SNPs might be one possible explanation.
The fact that we do not observe significant differences in epithelial E-cadherin expression between the mutant and wild-type alleles during ICS use could be attributed to the fact that the mutant alleles interact with ICSs and promote E-cadherin upregulation to a similar level as the wild-types. However, in the absence of ICS the mutant alleles are not capable of producing the same amount of protein compared to wild types. The nearly significant, inverse association we found between damaged/denuded epithelium and epithelial E-cadherin expression indicates that E-cadherin is involved in epithelial integrity, and it is probable that epithelial damage is increased in subjects with lower E-cadherin expression.
In our study, we also observed that specific CDH1 polymorphisms, in the presence of ICSs, tended to be associated with recruitment of CD8+ T-cells in the submucosa. Previously, Sont et al. [34] found that despite ICS treatment, CD8+ T-cell infiltration might persist in asthma patients and that CD8+, but not CD4+, cells are associated with BHR [34] and lung function decline [35], and recently, Hirosako et al. [8] showed that CD8 and CD103 (ligand for E-cadherin) are highly expressed in asthmatic bronchial intraepithelial lymphocytes. In that study, the percentage of CD8+ cells was higher than the percentage of CD4+ cells in intraepithelial lymphocytes in asthma. The expression of CD103 was significantly higher in CD8+ cells compared with CD4+ cells, suggesting that the interaction between E-cadherin (CD103) and CD8+ cells is the reason for the higher percentage of CD8+ cells in asthmatic intraepithelial lymphocytes [8]. Furthermore, Cepek et al. [36] demonstrated that T-cells express a member of the cadherin superfamily and this may contribute to T-cell-mediated immune surveillance via heterotypic adhesions with mucosal epithelial cells. Combining these findings, we could speculate that there might be a link between specific CDH1 polymorphisms and expression or binding of E-cadherin on T-cells, a finding that needs further investigation.
One could suggest that the lack of multiple testing correction is responsible for the current results. We decided not to apply a sequential (classical) Bonferroni correction for a number of reasons. First, our choice for the current study was explicitly driven by previous observations suggesting E-cadherin expression is related to epithelial integrity, airway remodelling and disease progression. Secondly, the independent variables in our analyses (e.g. airway remodelling and FEV1 decline) are mutually related, indicating that a rigid statistical procedure, such as Bonferroni correction for multiple testing, would not do justice to their biologically linked nature. Finally, although adjustment for multiple testing will decrease the chance of a type I error, it will also increase the chance of a type II error, so that a true association is not found. This is especially possible in a relatively small study such as ours. Thus, we followed the advice given by Perneger [37]: “Simply describing what was done and why, and discussing the possible interpretations of each result, should enable the reader to reach a reasonable conclusion without the help of Bonferroni adjustments”.
In our study, the associations between the SNPs and the biopsy parameters in population 1 were confirmed in population 2, where the same risk alleles were associated with FEV1 decline. This is called a loose replication [38]. The observed interaction of CDH1 SNPs with ICS use makes it difficult to find another study that is suitable for replication of our results. First, no other study on adult asthma patients exists with longitudinal data on lung function both with and without ICS use during such a long follow-up period. Secondly, bronchial biopsy studies in asthmatics usually investigate fewer patients than our study, which decreases the power to find associations.
In conclusion, our data show that CDH1 SNPs are associated with epithelial E-cadherin expression and suggest that epithelial adhesion is an important contributor to airway remodelling and lung function in asthma. These effects are modified by the use of ICSs.
Acknowledgments
We thank F. Volbeda and A. Dijkstra (Dept of Pulmonology, University Medical Center Groningen (UMCG), Groningen, the Netherlands) for collecting the biopsies, M. Broekema and M. Lodewijk (Dept of Pathology and Medical Biology, UMCG) for processing and staining of the biopsies, and M. Olivier (Dept of Genetics, UMCG) for her help with genotyping.
Footnotes
This article has supplementary material available from www.erj.ersjournals.com
Support Statement
This study was supported, in part, by a grant from the Netherlands Asthma Foundation (AF 3.2.07.15).
Statement of Interest
Statements of interest for D.S. Postma and G.H. Koppelman can be found at www.erj.ersjournals.com/site/misc/statements.xhtml
- Received December 16, 2010.
- Accepted April 6, 2011.
- ©ERS 2011