Intermittent exposure to whole cigarette smoke alters the differentiation of primary small airway epithelial cells in the air-liquid interface culture

Gindele, Julia A.; Kiechle, Tobias; Benediktus, Kerstin; Birk, Gerald; Brendel, Michael; Heinemann, Fabian; Wohnhaas, Christian T.; LeBlanc, Michelle; Zhang, Haijun; Strulovici-Barel, Yael; Crystal, Ronald G.; Thomas, Matthew J.; Stierstorfer, Birgit; Quast, Karsten; Schymeinsky, Jürgen

doi:10.1038/s41598-020-63345-5

Download PDF

Article
Open access
Published: 10 April 2020

Intermittent exposure to whole cigarette smoke alters the differentiation of primary small airway epithelial cells in the air-liquid interface culture

Julia A. Gindele^1,2,
Tobias Kiechle¹,
Kerstin Benediktus¹,
Gerald Birk³,
Michael Brendel⁴,
Fabian Heinemann³,
Christian T. Wohnhaas^5,6,
Michelle LeBlanc⁷,
Haijun Zhang⁷,
Yael Strulovici-Barel⁷,
Ronald G. Crystal⁷,
Matthew J. Thomas¹,
Birgit Stierstorfer³,
Karsten Quast⁵ &
…
Jürgen Schymeinsky^1,2

Scientific Reports volume 10, Article number: 6257 (2020) Cite this article

6773 Accesses
41 Citations
1 Altmetric
Metrics details

Subjects

Abstract

Cigarette smoke (CS) is the leading risk factor to develop COPD. Therefore, the pathologic effects of whole CS on the differentiation of primary small airway epithelial cells (SAEC) were investigated, using cells from three healthy donors and three COPD patients, cultured under ALI (air-liquid interface) conditions. The analysis of the epithelial physiology demonstrated that CS impaired barrier formation and reduced cilia beat activity. Although, COPD-derived ALI cultures preserved some features known from COPD patients, CS-induced effects were similarly pronounced in ALI cultures from patients compared to healthy controls. RNA sequencing analyses revealed the deregulation of marker genes for basal and secretory cells upon CS exposure. The comparison between gene signatures obtained from the in vitro model (CS vs. air) with a published data set from human epithelial brushes (smoker vs. non-smoker) revealed a high degree of similarity between deregulated genes and pathways induced by CS. Taken together, whole cigarette smoke alters the differentiation of small airway basal cells in vitro. The established model showed a good translatability to the situation in vivo. Thus, the model can help to identify and test novel therapeutic approaches to restore the impaired epithelial repair mechanisms in COPD, which is still a high medical need.

Smoking shifts human small airway epithelium club cells toward a lesser differentiated population

Article Open access 08 September 2021

Mahboubeh R. Rostami, Michelle G. LeBlanc, … Ronald G. Crystal

Electronic cigarette smoke reduces ribosomal protein gene expression to impair protein synthesis in primary human airway epithelial cells

Article Open access 01 September 2021

Hae-Ryung Park, Jose Vallarino, … Quan Lu

Altered generation of ciliated cells in chronic obstructive pulmonary disease

Article Open access 29 November 2019

Sophie Gohy, François M. Carlier, … Charles Pilette

Introduction

Chronic obstructive pulmonary disease (COPD) is the third leading cause of death worldwide and its prevalence continues to rise¹. The main risk factor to develop COPD is cigarette smoke^2,3. Smoking induces epithelial injury and this repeated injury of the epithelium triggers a pathophysiologic response, which leads to tissue remodeling of the airways that is characteristic for COPD^4,5. These changes of the small airway epithelium in COPD include: Goblet cell metaplasia^6,7,8, reduced cilia function^{9,10,11,12,13}, reduced club cell numbers^7,14,15, basal membrane thickening^16,17, epithelial barrier dysfunction^18,19,20 and squamous metaplasia^8,21,22,23. Furthermore, the epithelial defense mechanisms against inhaled particles and pathogens are impaired enabling sub-epithelial penetration of pathogens that increases the risk of COPD patients suffering from bacterial and viral infections and subsequent exacerbations^24,25,26.

To address cigarette smoke (CS)-induced damage on epithelial cells in vitro, previous studies used primary epithelial cells or cell lines that were mostly exposed to cigarette smoke extract (CSE) or to whole CS. These studies demonstrate smoke effects e.g. on epithelial barrier integrity, mucus production and cilia toxicity^{27,28,29,30,31,32,33,34}. The majority of these studies focus on the pathophysiology of large airways, i.e. bronchial or tracheal epithelial cells.

One hallmark of COPD is chronic bronchiolitis that mainly affects the small airways¹⁷. Therefore, a novel protocol was developed to culture human small airway epithelial cells (SAEC) at the air-liquid interface (ALI). To model CS induced alterations of the small airway epithelium, SAEC were repeatedly exposed to whole cigarette smoke during the differentiation process. In order to address the question if smoking and health status has an impact on the sensitivity to cigarette smoke exposure, cells from three healthy non-smokers and three COPD-smokers were compared.

To assess, how the in vitro model translates to the situation in vivo, the gene expression profiles of the ALI cultures were compared to the profiles of human small airway samples from healthy non-smokers and healthy smokers.

Results

Establishment of SAEC ALI culture

To further understand the molecular mechanisms underlying small airway disease and airway remodeling in COPD a new in vitro protocol was established for the cultivation of primary SAEC at ALI. First, the differentiation of SAEC under ALI conditions was characterized by analyzing the histology, the cellular composition (expression of cellular markers by RT-PCR technique), and the epithelial physiology (measuring transepithelial resistance and cilia beating). Removal of the apical medium (initiation of ALI culture) induced the differentiation of SAEC into a pseudostratified two-layered epithelium in 28 days (Fig. 1a). By means of immunohistochemistry, the main small airway cell types: KRT5⁺ basal cells, FOXJ1⁺ ciliated cells and SCGB1A1⁺ club cells were detected (Fig. 1a). MUC5AC⁺ goblet cells are normally absent in the human small airways³⁵, however an intermediate cell population of SCGB1A1⁺/MUC5AC⁺ double positive secretory cells and a few SCGB1A1⁻/MUC5AC⁺ goblet cells were detected in the ALI cultures (Fig. 1a). For the characterization of the epithelial physiology, transepithelial electrical resistance (TEER; Fig. 1b) was measured and cilia beating was investigated using a high-speed camera, enabling the quantification of the area of ciliated cells and the frequency of cilia beating (Fig. S1). At day 28 post airlift, the ALI cultures established from SAECs have a typical TEER between 400–600 Ω*cm² (Fig. 1b), a cilia beat frequency between 7–10 Hz and show an area of 5–20% covered by actively beating cilia.

A COPD-like phenotype is conserved in ALI cultures derived from donors with COPD

SAECs from three different COPD donors and three healthy controls were used to test whether pathologic changes in the epithelial histology are preserved in the differentiated ALI cultures. Of note, histological slices from ALI cultures obtained from COPD patients show a reduced number of ciliated cells and an increased number of secretory cells in comparison to control ALI cultures from healthy donors (Fig. 1a). To characterize the differential mRNA expression in detail, a Next Generation Sequencing (NGS) analysis was performed of ALI cultures obtained from healthy donors and COPD patients at day 28 upon air-lift. 248 transcripts were detected that were significantly changed in COPD-derived cells compared to healthy volunteer-derived cells. Hierarchical clustering of these transcripts reveals different mRNA signatures in healthy and COPD cultures (Fig. 1c). To assess the ability of the COPD ALI signature to distinguish between human small airway epithelial (SAE) samples of healthy non-smokers and COPD smokers as well as between healthy smokers and COPD smokers, the signature was used as feature vector in a random forest classification analysis³⁶. The power of the signature to separate the respective samples was tested by determining the ROC (receiver operating characteristic) curves (Fig. 1d). Classifying samples from healthy non-smokers and COPD smokers resulted in an AUC of 1.0, indicating a nearly perfect separation of SAE samples by utilizing the ALI culture derived signature. To ensure, that the high classification power is not driven solely by very few well separating transcripts, but rather by a large proportion of the signature, the best correlating transcripts were identified. Subsequently, ten transcripts were removed from the signature, which showed the highest correlation. After reiterating the classification, an AUC of 0.98 was obtained. This indicates a clear relationship between the transcriptional changes within the ALI culture and the SAE samples, which is not dependent upon the top ten classifying transcripts. The high overlap between the in vitro derived signature and the clinical situation allows a good prediction of the health status. In the task of separating smokers from COPD smokers an AUC of 0.8 was obtained. Upon removal of the ten highest correlating transcripts, this drops to 0.65. Even though the classification power is lower in this case, it still indicates an overlap of COPD driven transcriptional changes in the ALI samples with the changes observed in SAE from human subjects.

Cigarette smoke exposure impairs epithelial barrier integrity

To investigate the harmful effects of CS on SAEC differentiation and physiology the ALI cultures were intermittently exposed to well-defined doses of whole CS during the differentiation period of 28 days. Unlike commonly used exposure models that use cigarette smoke extract (CSE) in which many volatile components are lost during the extraction process, using whole CS exposes the cells to all volatile toxins contained within CS. Furthermore, this is likely to focus CS exposure to the apical cell layer and reduce exposure to the basal cell layer, in contrast to a stimulation with CSE added to the basal medium. It is assumed that this setting is an optimal framework to mimic the direct toxic effects of smoking on the small airway epithelium.

In order to assess effects of CS on epithelial barrier integrity, TEER was repeatedly measured in the ALI cultures during differentiation. Air control treated ALI cultures of healthy donors and COPD patients exhibited similar values at days 7, 14, 21 and 28 (Fig. 2a). However, the TEER values of CS exposed cultures were reduced compared to air controls (Fig. 2a). To quantify this effect the area under the curve (AUC) was calculated and standardized with respect to the time span. The adjusted AUC dropped in the healthy ALI cultures from 535 ± 91 Ω*cm² in air controls to 354 ± 36 Ω*cm² (p = 0.01) upon CS exposure, and in the COPD ALI cultures from 500 ± 143 Ω*cm² to 337 ± 50 Ω*cm² (p = 0,017), respectively (Fig. 2a). To investigate the underlying molecular mechanisms mRNA was extracted from the ALI cultures and the expression of junctional protein transcripts was analyzed. Some genes associated with cell adhesion, permeability and tight junction regulation were deregulated by CS. JAM3, CLDN11 and CLDN18 were down-regulated after 28 days of intermittent CS treatment. CLDN7 and CLDN10 were up-regulated in both healthy and COPD cultures (Fig. 2b). The upregulation of the strongest deregulated junctional protein on mRNA level, namely Claudin-10, was confirmed on protein level by means of immunohistochemistry (Fig. 2c). Of note, most of the prevalent classic claudins in the lung³⁷, i.e. claudin-1, -3, -4, and -5, were not significantly deregulated on mRNA level (data not shown).

Cigarette smoke exposure induces squamous differentiation

One hallmark of COPD is squamous metaplasia²¹. Therefore, the expression of different squamous cell marker genes was investigated during the differentiation of smoked and control ALI cultures (RNA-seq, Fig. 3c). During the differentiation of the ALI cultures, the difference of threshold cycle values (∆Ct: Ct target gene – Ct endogenous control) for Keratin 5 (KRT5) increased from the day of airlift until day 14, where it reached a plateau, indicating that KRT5 mRNA expression is downregulated during the normal differentiation process. This was true for ALI cultures from healthy and COPD SAECs. However, the exposure to CS led to a decrease of the ∆Ct values which reflects an increase of KRT5 mRNA expression in healthy as well as in COPD ALI cultures (Fig. 3a). To quantify this effect over the observed time the AUC was calculated and standardized with respect to the time span. The exposure of ALI cultures to CS led to a decreased AUC which was similar in healthy (from −2.3 ± 0.4 ∆Ct to −2.8 ± 0.3 ∆Ct, p = 0.0004) and COPD derived cultures (from −2.6 ± 0.1 ∆Ct to −3.0 ± 0.1 ∆Ct, p = 0.0004) (Fig. 3b). The ∆Ct values for the squamous differentiation marker KRT14 increased during the ALI differentiation indicating that KRT14 mRNA expression is downregulated during the normal differentiation process. However, the ∆Ct values were not down regulated during the differentiation of CS exposed cultures (Fig. 3b) indicating that CS-exposure is blocking the physiologic downregulation of KRT14 expression. Accordingly, the AUC was downregulated upon CS exposure in ALI cultures from COPD patients (from 5.0 ± 0.3 ∆Ct to 3.3 ± 0.3 ∆Ct, p < 0.0001) and healthy controls (from 5.5 ± 0.8 ∆Ct to 3.1 ± 0.9 ∆Ct, p < 0.0001) (Fig. 3b). RNA sequencing analysis revealed deregulation of several cell marker genes, i.e. genes associated with basal cells, goblet and club cells as well as ciliated cells (Fig. 3c). KRT5 and KRT14 up-regulation was confirmed and the expression of other genes associated with squamous differentiation, namely Keratin 6 A (KRT6A), Keratin 13 (KRT13), and Stratifin (SFN) were found to be also up-regulated upon CS exposure, indicated by a positive fold-change (Fig. 3c). Besides evaluation on RNA level, KRT5 expression was assessed on protein level, qualitatively and representatively (Fig. 3d). Deregulation of secretory and ciliated markers were further evaluated in the following sections.

Cigarette smoke exposure stimulates the differentiation of MUC5AC ⁺ cells and reduces SCGB1A1 ⁺ cells

To investigate the effect of CS on basal cell differentiation into secretory cells, the club cell marker SCGB1A1 and the goblet cell marker MUC5AC were analyzed on mRNA level at day 7, 14, 21 and 28 after air-lift. The ΔCt values decreased after the day of air-lift, indicating an induction of MUC5AC and SCGB1A1 mRNA expression, which goes along with the normal differentiation process. The CS exposed cultures exhibit lower ΔCt values for MUC5AC, indicating an increased MUC5AC mRNA expression. At the same time, ΔCt values for SCGB1A1 were increased during differentiation, indicating a reduced SCGB1A1 mRNA expression (Fig. 4a). The adjusted AUC was calculated to compare the differentiation into secretory cells of smoke exposed cultures and air controls. The AUC for MUC5AC was down-regulated upon CS exposure in ALI cultures from COPD patients (from −1.1 ± 0.9 ∆Ct to –1.5 ± 0.5 ∆Ct, p = 0.3252) and healthy controls (from −1.3 ± 0.6 ∆Ct to −2.0 ± 0.5 ∆Ct, p = 0.2080). This indicates an approximate increase of MUC5AC mRNA expression of 1.3 fold in the COPD cultures and 1.6 fold in the healthy ALI cultures. The AUC for SCGB1A1 was up-regulated in CS-exposed ALI cultures from COPD patients (from −10.3 ± 0.03 ∆Ct to −8.9 ± 0.8 ∆Ct, p = 0.0289) and healthy controls (from −9.9 ± 0.1 ∆Ct to −9.3 ± 0.1 ∆Ct, p = 0.1317). This indicates an approximate reduction of SCGB1A1 mRNA expression of 2.6 fold in the COPD cultures and 1.6 fold in the healthy ALI cultures. Furthermore, a RNA sequencing analysis indicated an increase of goblet cell associated genes (MUC5AC, MUC1) and a reduction of club cell genes (SCGB1A1, SCGB3A1, CYP2F1, TMEM45A) upon CS treatment (Fig. 3c).

To investigate if the changes on mRNA level led to altered protein expression and cell composition MUC5AC and SCGB1A1 protein expression were analyzed by means of FACS analysis. The analysis revealed for the majority of donors (five out of six) a decrease in SCGB1A1⁺ single positive cells after 28 days with intermittent CS treatment (Fig. 4b,c). The effect on MUC5AC⁺ single positive and MUC5AC⁺/SCGB1A1⁺ double positive cells was not significant. Still half of the donors showed a clear increased expression of MUC5AC (however, two donors showed a decreased number of MUC5AC⁺/SCGB1A1⁻ cells). Next, it was assessed whether CS exposure induces histological differences of the ALI cultures by performing immunohistochemical stainings of sections. Unexpectedly, CS exposed cultures as well as air controls exhibited a differentiated pseudostratified epithelium without significant differences in the epithelial thickness. However, the analysis confirmed the previous findings, that the number of SCGB1A1⁺ cells decreased in the smoke exposed ALI cultures compared to air controls, whereas the number of MUC5AC⁺ cells increased, respectively (Fig. 4d).

Cigarette smoke exposure impairs the development of ciliated cells

The development of ciliated cells was monitored by measuring the area covered by beating cilia and their corresponding beat frequency at days 7, 14, 21 and 28. Moreover, the expression of the transcription factor FOXJ1 and the presence of acetylated tubulin was measured, since proper cilia function, changes in transcript levels and changes in protein levels may not necessarily correspond to each other. COPD-derived cultures exhibited a delayed and reduced differentiation of ciliated cells compared to healthy cultures, indicated by a reduced area covered by actively moving cilia, and reduced acetylated tubulin expression (data not shown). These parameters were further reduced by intermittent CS treatment (Fig. 5a–c). Interestingly, cigarette smoke exposure impaired the development of ciliated cells in healthy controls in almost the same manner. To quantify the reduction in ciliary beating, the AUC was calculated and standardized with respect to the time span. Both, the area covered by actively beating cilia and the cilia beat frequency were reduced by intermittent CS treatment (Fig. 5b). FOXJ1 mRNA increased during the normal differentiation process, indicated by a decrease in the ΔCt values (Fig. 5c). This induction was delayed in CS exposed cultures. To quantify the CS-induced effect on FOXJ1 mRNA expression the AUC was calculated and standardized with respect to the time span (Fig. 5c). Furthermore, FOXJ1 and DNAI1 expression was assessed by RNA sequencing on day 28. The analysis revealed a non-significant deregulation of the ciliated cell markers on RNA level upon 28 days of intermittent CS exposure (Fig. 3c). To investigate if CS exposure not only affects the cilia function and mRNA expression but also the cell composition acetylated tubulin was stained as a cilia marker in the histological slices. Air control cultures exhibited a well-differentiated ciliated epithelium with many acetylated tubulin-positive cilia. CS-exposed ALI cultures revealed less cilia (Fig. 5d).

CS exposure lead to the deregulation of a similar subset of genes in healthy and COPD derived ALI cultures

To elucidate the transcriptional response of the small airway epithelium to CS a NGS analysis was performed. To address the question whether COPD derived cultures are differently affected by smoking compared to healthy cultures the transcriptional response upon CS exposure was compared by performing a differential expression analysis of the NGS data. CS modulated the expression of 2442 genes in healthy cultures and 3249 genes in COPD cultures. Notably, 1945 genes were regulated in healthy as well as in COPD derived cultures. That all but one of the 1945 overlapping transcripts were regulated the same way indicates that the main stimulus of altered gene expression in the study was CS and not the disease status of the donors, i.e. COPD vs. healthy.

To strengthen the translational link of the in vitro system with the situation in vivo, the CS-induced transcriptional changes of the healthy ALI cultures were compared with a published data set of epithelial brushes from healthy smokers and non-smokers (GSE11784)³⁸ to assess the overlap of deregulated genes. Therefore, the normalized gene expression data (range from 0 to 1) from SAE and ALI cultures were combined and filtered for those transcripts that were significantly deregulated in both comparisons: a) smoked ALI cultures versus air controls and b) SAE from healthy smokers versus SAE from non-smokers. Subsequently an unsupervised, hierarchical clustering was performed using correlation as distance measure on the combined and filtered data-set (Fig. 6a). This approach allowed the separation of different clusters according to their similarities. Two distinct clusters representing the in vitro samples from CS and air-treated ALI cultures were separated. Of note, most of the SAE samples from smokers clustered with the smoked ALI cultures whereas the SAE from non-smokers clustered with the ALI air-controls. This unbiased approach (which is only focusing on transcripts that were significantly deregulated but not filtering for common directionality of deregulation) indicates that CS alters overlapping molecular pathways in SAEC ALI cultures and in SAE from smokers.

To quantify to what extent the smoked ALI culture resembles the SAE of smokers a correlation analysis was performed (Fig. 6b). More than 86% of the genes that were significantly deregulated in the CS-exposed ALI cultures and the SAE of healthy smokers showed a consistent direction of deregulation (Fig. 6c). Furthermore, a pathway analysis was performed to elucidate the common mechanisms affected by CS-treated ALI in vitro and SAE from smokers. In order to allow for identifying pathways that were altered in both data sets, the significantly deregulated transcripts that exhibit the same direction of change were combined. With this set of transcripts, a core analysis was performed utilizing Ingenuity Pathways Analysis (IPA). A large number of top scoring canonical pathways were categorized to “xenobiotic metabolism” (supplementary Table S3). Therefore, all transcripts that were reported to be associated to xenobiotic metabolism were extracted and used to generate a combined network representing the in vivo and in vitro data (Fig. 7). Among the top deregulated canonical pathways associated to the network and consistent with the in vivo SAE data were Nuclear factor erythroid 2–related factor 2 (NRF2)-mediated oxidative stress, xenobiotic metabolism and aryl hydrocarbon (AHR) signaling (Figs. 7 and 8).

Discussion

COPD is a disease that strongly affects the small airways^17,39. To investigate the underlying pathophysiological and molecular mechanisms in vitro, a new protocol for SAEC ALI cultures and novel methods for their analysis were established. In contrast to multi-layered bronchial epithelial cell cultures, the SAEC ALI culture exhibit a two-layered pseudostratified epithelium closely resembling the human small airway in vivo⁴⁰. This model allows culturing SAEC from healthy donors and COPD patients with the latter displaying some conserved pathological characteristics, such as reduced numbers of ciliated cells, increased numbers of secretory cells and a distinct gene signature. The analysis of the gene signature of the ALI cultures derived from healthy or COPD donors with in vivo samples from human SAE revealed a high classification power and therefore a high predictability of the health status.

A pathologic hallmark of COPD is a reduced airway epithelial barrier function leading to an increased susceptibility to infections^{18,19,20,24,25,26}. In this model, CS exposure during SAEC differentiation impaired epithelial barrier formation characterized by a decreased TEER. This is consistent with previous studies investigating the effect of cigarette smoke extract on HBEC ALI cultures²⁷. The observed functional decline of epithelial barrier function may be linked to the deregulation of several genes that are associated with the maintenance and control of cell-cell adhesion, permeability and tight junctions. An important molecular role in these aspects play the Claudin family members, which are tight junction membrane proteins that are conceptually divided into barrier-forming (e.g. Claudin-11, -18) and pore-forming (e.g. Claudin-7, -10) proteins^41,42,43. Notably, Claudin-11 and -18, which are known to seal paracellular permeability⁴¹, were down-regulated in the CS exposed cultures. Claudin-18 deficiency was also described to be associated with airway epithelial barrier dysfunction and asthma⁴⁴. Furthermore, Claudin-10, which has been unequivocally defined as pore-forming claudin^41,45 was up-regulated. Claudin-10 and Claudin-7 were also shown to be up-regulated in smokers and COPD smokers, respectively⁴⁶. Apart from claudins, the junctional adhesion molecule C (JAM3), a component of desmosomes, which mediates cell-cell adhesion⁴⁷ was described to be down-regulated in COPD patients. JAM3 was also decreased in the ALI cultures after 28 days of CS exposure. Taken together, different genes associated with the regulation of epithelial barrier function were deregulated in the CS-treated ALI cultures that may be associated with the pathogenesis in COPD.

Basal cell hyperplasia and squamous metaplasia play an essential role in the emergence of airway obstruction in COPD¹⁷ and increase the risk to develop squamous cell carcinoma^48,49,50. Previous studies have found increased squamous metaplasia in the bronchial epithelium of smokers with COPD compared to healthy smokers²¹. In this study, an increase in several genes associated with squamous metaplasia, namely KRT5, KRT6A, KRT13, KRT14 and SFN was observed indicating that CS plays a role in squamous differentiation. These in vitro findings are in accordance with histological (known for a long time⁵¹) and expression analyses (more recent ones⁵²) of COPD patients indicating that basal cells are the first cell population affected in smokers and leading to an altered differentiation pattern and architecture of the small airway epithelium. Taken together, the model can be used to investigate the molecular mechanisms underlying the early effects of smoking on basal cell proliferation.

Smoking induces changes in secretory cell populations leading to hyperplasia of mucous-producing goblet cells and a reduction of club cells, which are associated with epithelial defense mechanisms^7,53. These findings resemble the goblet cell hyperplasia in COPD patients, which was shown to be associated with current smoking^6,54. Recent findings indicate a shift in the differentiation pattern of the small airways of smokers and COPD patients towards a more proximal airway epithelium^35,55; an effect recently described as smoking-dependent distal-to-proximal repatterning⁵⁶. Goblet cells for example were found in the small airways of smokers⁶, where they are normally absent. This aspect of the disease was found in the smoked SAEC ALI cultures. CS exposure led to a reduction of club cell markers (SCGB1A1, CYP2F1) and a downward trend of MUC5AC⁻/SCGB1A1⁺ cells in the ALI cultures. This may reflect the pathophysiological decline in club cells observed in COPD. The reduced club cell numbers in the airways of smokers and COPD patients, respectively^7,15, are associated with reduced uteroglobin (SCGB1A1) serum levels, which were correlated to a decline in FEV1 (Forced Expiratory Volume in 1 Second) and COPD disease progression^14,52. In contrast to our findings, the exposure to CSE did not show a reduction in club cell numbers²⁷. Therefore, we think that the apical exposure of the SAEC ALI culture with whole CS leads to a better representation of the pathophysiological processes induced in human smokers. The model may help to dissect the molecular mechanisms, leading to the impaired differentiation of basal cells in COPD patients^{14,53,57,58,59}.

Coordinated and efficient ciliary beating is essential to clear the airways from inhaled particles and pathogens and is therefore an important defense mechanism⁶⁰. COPD patients often suffer from severe exacerbations caused by bacterial or viral infections, due to insufficient clearance^26,61. We show that CS exposure during SAEC differentiation impairs the development and function of ciliated cells. These findings are consistent with studies, investigating the effect of cigarette smoke extract on human bronchial epithelial cell differentiation²⁷. To evaluate to which extent the CS ALI in vitro model reflects smoking induced effects on the lung epithelium of small airways in living humans, a transcriptomic comparison analysis of SAEC ALI cultures and human SAE brushings was conducted. Of note, SAEC ALI cultures stimulated with CS exhibit a gene expression pattern overlapping to epithelial samples from healthy smokers. When combining both datasets and performing a hierarchical clustering on the overlapping de-regulated transcripts, it would be expected to see similar patterns but not to find the in vivo data and the in vitro data forming separate clusters. Remarkably, CS stimulated SAEC ALI cultures clustered with samples from smokers, whereas air treated SAEC ALI cultures clustered with healthy control samples in an unexpected fashion. By this means, we did not include mechanisms that were only affected in vitro and not in vivo or vice versa, but we focused on processes that were modulated in both situations. Our results highlight the presence of shared processes that show a corresponding de-regulation, i.e. a similar up or downregulation in vitro and in vivo. Thus, we think that the CS treated ALI culture from SAEC is a valuable in vitro method to investigate smoke induced effects on epithelial differentiation. However, our data does not allow distinguishing CS-induced effects in ALI cultures from healthy donors compared to COPD patients, since the CS is a very dominant driver of differential changes. The dissection of potential differences in smoked healthy versus smoked COPD ALI cultures would require a much higher donor number. Here we see a clear limitation of our complex experimental in vitro setup; nevertheless, ALI cultures from COPD patients may be extremely helpful to test whether potential therapeutic concepts work also in a diseased setting.

Among the top deregulated canonical pathways were Nuclear factor erythroid 2–related factor 2 (NRF2)-mediated oxidative stress, xenobiotic metabolism and aryl hydrocarbon receptor (AHR) signaling^62,63,64,65. It is known that CS induces reactive oxygen species, which activates the NRF2 pathway^63,66. Moreover, CS comprises a number of toxic chemicals including polycyclic aromatic hydrocarbons (PAH)⁶⁷. A prominent member of PAH is benzo[a]pyrene, which activates the AHR pathway, inducing the expression of several metabolic enzymes, including CYP1A1, CYP1A2 and CYP1B1^62,68. The model reflects these CS induced pathways in close resemblance to the in vivo situation and offers novel opportunities to study new therapeutic targets addressing oxidative stress responses and xenobiotic metabolism.

Inconsistencies with the in vivo data were strongly associated with lipid metabolism (e.g. synthesis of cholesterol).

Our data demonstrate that CS specifically alters the small airway epithelium by affecting basal cell differentiation. Thereby, important pathways are consistently modulated compared to the small airways of human smokers. Taken together, critical pathophysiologic and molecular effects could be mimicked in vitro that were described in small airways of smokers, which might lead to novel essential therapeutic approaches to restore the impaired epithelial repair mechanisms in a smoke-induced chronic disease like COPD.

Methods

SAEC ALI culture

Human Small Airway Epithelial Cells (SAEC, CC-2547, Lonza, Basel, Switzerland) of three healthy donors and three COPD patients (see Supplementary Table S1) were cultured as described in Supplementary Information. As stated by Lonza, the cells were isolated from donated human tissue after obtaining permission for their use in research applications by informed consent or legal authorization. Established ethical practices of the donation and transplantation organizations in the US (American Association of Tissue Banks, Association of Organ Procurement Organizations, Eye Bank Association of America) are followed at Lonza. All methods were carried out in accordance with relevant guidelines and regulations. Human SAEC were exclusively used for in vitro experiments. SAEC were grown in PneumaCult™-Ex Plus Medium (Stemcell Technologies, Vancouver, Canada) on rat tail collagen type 1 (Corning Life Sciences B.V., Amsterdam, the Netherlands) coated transwells (#3460, Corning Life Sciences B.V., Amsterdam, the Netherlands) until confluent. Cells were further cultivated in air-liquid interface using Pneumacult-ALI-S Medium (Stemcell Technologies, Vancouver, Canada) for four weeks until fully differentiated.

Cigarette smoke exposure

SAECs were exposed to whole cigarette smoke of K3R4F reference cigarettes (University of Kentucky, Lexington, KY, USA) using an automated cigarette smoking machine (In-Expose smoking robot, Scireq Montreal, QC, Canada) and the P.R.I.T. (Professional In-Vitro Technologies) ExpoCube in the experimental setup developed by Fraunhofer Institute for Toxicology and Experimental Medicine. Four cigarettes were smoked in parallel in compliance to ISO 3308, drawing every 15 seconds a puff from the sequentially smoked cigarettes (9 puffs per cigarette). The smoke was diluted with ambient air (dilution rate: 0.5 l/min). The cells were exposed to CS three times a week during the differentiation phase (28 days) beginning on day 0 (= day of air-lift). Cells were cultured on 12-well plates with the upper four wells exposed to smoke and the lower four wells exposed to air⁶⁹. All wells were uniformly exposed to CS due to computational fluid dynamics-optimized geometry of the flow paths. This was previously tested and described in detail by Ritter et al., Fraunhofer Institute for Toxicology and Experimental Medicine⁷⁰. Readouts were performed 24 hours post last exposure.

Cytotoxicity of smoke-exposed cells was assessed on day 28, based on the measurement of lactate dehydrogenase (LDH) activity released from the cytosol of damaged cells into the supernatant (compare Supplementary Methods). There was no significant difference in LDH release between smoke-exposed and air-exposed cells (Supplementary Fig. S4).

Transepithelial Electrical Resistance (TEER)

Epithelial barrier integrity was monitored by measuring Transepithelial Electrical Resistance (TEER) using an EVOM2 epithelial volt-ohmmeter (World precision Instruments, Sarasota, FL, USA) according to the manufacturer’s instructions. To perform the TEER measurement, phosphate buffered saline (PBS, 300 µl) was added to the apical compartment of the insert. The long stem of the electrode was inserted through the gap of the transwell to be in contact with the basolateral medium and the short stem of the electrode was placed above the apical surface to be in contact with the apical PBS. The measured resistance was multiplied by the surface area of the epithelium (1.12 cm²) to obtain TEER (Ω*cm²).

Histology and Immunohistochemistry

Cells, growing on transwell membranes (pore size of 0.4 µm), were fixed in 4% paraformaldehyde over night and embedded in paraffin. After sectioning (section slice thickness approximately 3–4 µm), slides were deparaffinized and hematoxylin and eosin (H&E) stainings were performed according to standard procedures. For immunohistochemical (IHC) investigation, paraffin-embedded tissue slides were subjected to heat induced antigen retrieval in citrate buffer (pH 6.0, H.I.E.R, BioLegend San Diego, CA, USA) for 30 minutes. Antibodies used for immunofluorescence or chromogenic antigen detection are listed in Supplementary Information, Table S2. After staining, tissues slides were mounted either with ProLong™ Diamond Antifade Mountant (Thermo Fisher Scientific, Waltham, MA, USA) or with Entellan® new (Merck KGaA, Darmstadt, Germany). For immunofluorescence samples, stained slides were stored at 4 °C in the dark and images were taken using the Zeiss Laser Scanning Microscope 710 or an AxioImager M2 (both Carl Zeiss, Oberkochen, Germany) in case of chromogenic IHC. Stainings for KRT5, MUC5AC, SCGB1A1 and acetylated Tubulin were used to assess protein expression qualitatively and representatively. Claudin-10 protein expression was semi-quantitatively analyzed using digital image analysis (compare section “Image Analysis”).

Image Analysis

ALI cultures with Claudin-10 IHC staining (ab52234, Rabbit polyclonal to Claudin-10, Abcam, Cambridge, MA, USA) were scanned with an Axio Scan.Z1 scanner (Carl Zeiss Microscopy GmbH, Jena, Germany) using a 20x objective (0.22 µm/px) in bright field illumination. Semi-quantitative image analysis was performed with the Halo 3.0 software using the CytoNuclear 1.6 module (Indicalabs, Corales, NM, USA). In preparation for the automated analysis, regions with damaged ALIs or insufficiently stained cells were manually excluded. For automated analysis, a random forest tissue classifier was trained to recognize the cellular layer of the ALIs to limit the following analysis to this region and exclude the supporting material and background. Subsequently, the CytoNuclear module was applied within this region. As first step, the channels of nuclear (hematoxylin) and IHC stain (3’Diaminobenzidine (DAB) signal of Claudin-10) were unmixed using color deconvolution. Cell segmentation used the nuclear stain and empirical parameters on size and shape of the cells. Positive cells were identified using a manually optimized threshold value on the DAB channel, within the segmented region of a cell. The generated readout per ALI culture was the number of positive cells n per analyzed area of the cell layer A, i.e. n/A.

Cilia beat measurement of SAEC ALI cultures

Ciliary function was analyzed using video-microscopy as described in Supplementary Information. Briefly, cilia beat frequency was calculated from image stacks of 2D + time. For each pixel in the field of view, the signal in time domain is searched for a meaningful frequency. The results were mapped back to image space forming an image, where each pixel represents the beat frequency of the cilia covering that area.

RT-qPCR

RT-qPCRs were conducted using the High-Capacity cDNA Reverse Transcription Kit and TaqMan Assay-on-Demand kits from Applied Biosystems (Darmstadt, Germany) for KRT5 (Hs00361185_m1), MUC5AC (Hs01365616_m1), SCGB1A1 (Hs00171092_m1), FOXJ1 (HS00230964_m1), KRT14 (Hs00265033_m1) and POLR2A (Hs00172187_m1). The relative expression of the RNAs was determined by the comparative Ct method. Expression values were normalized to the control RNA (POLR2A).

Next generation sequencing and detection of differentially expressed genes

RNA extraction, Illumina library preparation and sequencing was performed as previously described⁷¹. Comparative analysis was done using the limma R-package⁷² by applying a the following linear model:

anno.sub < - data.frame(group = pre.design$group,donorID = pre.design$donorID)

dmat < - model.matrix(~0+group+donorID,data = anno.sub)

Benjamini-Hochberg correction was used to adjust for multiple testing. Expression data have been deposited to Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/; accession no. GSE135188). The previously published data set GSE11784³⁸ was used for comparative analysis.

Hierarchical clustering

Hierarchical clustering was performed with TIBCO Spotfire Analyst 7.11.1 LTS HF-013. As clustering method UPGM with correlation as distance measure was used for both columns and rows, respectively. Expression values were normalized to a per transcript range between 0 and 1 prior to clustering.

Classification

The Random Forest classifier offered by WEKA Workbench 3.8.3 (Weka, University of Waikato, New Zealand) with parameters -P 100 -I 100 -num-slots 1 -K 10 -M 1.0 -V 0.001 -S 1 was used for classifying samples and generating ROC curves and confusion matrices. Classification was performed on expression values normalized to a per transcript range between 0 and 1. The WEKA CorrelationAttributeEval attribute evaluator (Weka, University of Waikato, New Zealand) with standard parameters and WEKA Ranker as search method with standard parameters was used for ranking features according to their classification power.

Network generation

Ingenuity Pathway Analysis (IPA; Qiagen) version 48207413 was used to perform all network related analyses. Molecules associated with “Xenobiotic Metabolism Signaling” were compiled, utilizing the “Pathways and Tox Lists” related search tool, and gathered in one pathway view. Subsequently, molecules were overlaid with ALI (GSE135188) and SAE (GSE11784) expression data representing significantly deregulated transcripts in one or the other experiment. Molecules, which did not show deregulation in either experiment, were removed. The “Connect” tool was used to introduce relationships as reported within the IPA knowledge base among the remaining transcripts. Orphans and smaller graphs (less than four nodes) were removed in order to focus on the largest interconnected graph. To focus on specific transcripts and increase readability, a number of hubs were removed from the graph. Finally, either ALI or SAE expression values were laid over the resulting network and transcripts representing selected biological processes according to the IPA knowledge base were highlighted.

Flow cytometry

Single cell suspensions of SAEC ALI cultures (see Supplementary Information) were fixed in 4% paraformaldehyde for 20 minutes at room temperature. Cells were washed twice using Perm/Wash Buffer (BD Biosciences, San Jose, CA, USA) and blocked to prevent non-specific binding using Human BD Fc Block (BD Biosciences, San Jose, CA, USA). Cells were stained in Perm/Wash Buffer for 1 hour at room temperature. Anti-human antibodies used for staining are listed in the supplement. Flow cytometry was performed using a BD LSRFortessa X-20 cytometer equipped with DIVA-software (BD Biosciences, San Jose, CA, USA).

Statistics

Statistical testing used for RNA sequencing analysis is described in section “Next generation sequencing and detection of differentially expressed genes”. All the other analyses were performed as followed. Instead of modelling the longitudinal data, the analyses were performed with adjusted areas under the curve (AUC) or the logarithmized adjusted AUC. The adjusted AUC was obtained by dividing the AUC by the length of the considered time interval. Each of the endpoints was analyzed separately. The data was analyzed with a repeated measurement model, as both treatments (CS/Air) were applied to cells of each donor. Compound symmetry was assumed as underlying covariance structure in order to account for the repeated measurement within a donor. The fixed factors status of the cells (HC/COPD), treatment (CS/Air) and the interaction of status of the cells and treatment were included in the model. The degrees of freedom were calculated by dividing the residual degrees of freedom into between-subject and within-subject portions. If a fixed effect changed within any subject, the within-subject degrees of freedom were assigned to the considered fixed effect (cell status, treatment, interaction effect). Otherwise, the between-subject degrees of freedom were used. The applied repeated measurement model assumes that the errors follow a normal distribution.

Using linear combinations of the parameter estimates obtained from the repeated measurement model, it was tested whether treatment with CS compared to Air had an impact on the considered endpoint. This analysis was performed separately for each status of cells (HC/COPD). Depending on the expected change of the endpoint under the treatment with smoke, the following hypotheses were tested:

${ {\mathcal H} }_{0,{\rm{HC}},{\rm{Lower}}}:{\mu }_{{\rm{HC}},{\rm{CS}}}-{\mu }_{{\rm{HC}},{\rm{Air}}}\ge 0$ vs. ${ {\mathcal H} }_{1,{\rm{Lower}}}:{\mu }_{{\rm{HC}},{\rm{CS}}}-{\mu }_{{\rm{HC}},{\rm{Air}}} < 0$ and ${ {\mathcal H} }_{0,{\rm{COPD}},{\rm{Lower}}}:{\mu }_{{\rm{COPD}},{\rm{CS}}}-$ ${\mu }_{{\rm{COPD}},{\rm{Air}}}\ge 0$ vs. ${{\mathcal{H}}}_{1,{\rm{L}}{\rm{o}}{\rm{w}}{\rm{e}}{\rm{r}}}:{\mu }_{{\rm{C}}{\rm{O}}{\rm{P}}{\rm{D}},{\rm{C}}{\rm{S}}}-{\mu }_{{\rm{C}}{\rm{O}}{\rm{P}}{\rm{D}},{\rm{A}}{\rm{i}}{\rm{r}}} < 0$, where μ_i,j denotes the mean value of the investigated endpoint under treatment j for the status of cells i.
For the following endpoints the above stated hypotheses were tested as a decrease was assumed when cells were treated with smoke: Logarithmized AUC of Area covered by active cilia, AUC of Cilia beat frequency, logarithmized AUC of TEER, AUC of KRT5 mRNA (Ct), AUC of KRT14 mRNA (Ct), AUC of MUC5AC mRNA (Ct).
${ {\mathcal H} }_{0,{\rm{HC}},{\rm{Upper}}}:{\mu }_{{\rm{HC}},{\rm{CS}}}-{\mu }_{{\rm{HC}},{\rm{Air}}}\le 0$ vs. ${ {\mathcal H} }_{1,{\rm{HC}},{\rm{Upper}}}:{\mu }_{{\rm{HC}},{\rm{CS}}}-{\mu }_{{\rm{HC}},{\rm{Air}}} > 0$ and ${ {\mathcal H} }_{0,{\rm{COPD}},{\rm{Upper}}}:{\mu }_{{\rm{COPD}},{\rm{CS}}}-$ ${\mu }_{{\rm{COPD}},{\rm{Air}}}\le 0$ vs. ${ {\mathcal H} }_{1,{\rm{COPD}},{\rm{Upper}}}:{\mu }_{{\rm{COPD}},{\rm{CS}}}-{\mu }_{{\rm{COPD}},{\rm{Air}}} > 0.$

For the following endpoints the above stated hypotheses were tested as an increase was assumed when cells were treated with smoke: Claudin-10 protein expression, AUC of SCGB1A1 mRNA (Ct), AUC of FOXJ1 mRNA (Ct).

All tests are t-test like tests based on the linear combinations of the parameter estimates from the repeated measurement model and were performed post-hoc. The p-values were not adjusted for multiplicity induced by the number of tested endpoints and hypotheses and the selection of the finally applied statistical model.

The statistical evaluation was prepared using the software package SAS Version 9.4 (SAS Institute Inc., Cary, North Carolina, USA).

Data availability

All relevant data are within the paper and its Supplemental Information file. Sequencing data that support the findings of this study have been deposited in Gene Expression Omnibus (GEO; https://www.ncbi.nlm.nih.gov/geo/) with the accession code GSE135188.

References

Lozano, R. et al. Global and regional mortality from 235 causes of death for 20 age groups in 1990 and 2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet (London, England) 380, 2095–2128, https://doi.org/10.1016/s0140-6736(12)61728-0 (2012).
Article Google Scholar
Yoshida, T. & Tuder, R. M. Pathobiology of cigarette smoke-induced chronic obstructive pulmonary disease. Physiological reviews 87, 1047–1082, https://doi.org/10.1152/physrev.00048.2006 (2007).
Article CAS PubMed Google Scholar
Laniado-Laborin, R. Smoking and chronic obstructive pulmonary disease (COPD). Parallel epidemics of the 21 century. International journal of environmental research and public health 6, 209–224, https://doi.org/10.3390/ijerph6010209 (2009).
Article CAS PubMed PubMed Central Google Scholar
Puchelle, E., Zahm, J. M., Tournier, J. M. & Coraux, C. Airway epithelial repair, regeneration, and remodeling after injury in chronic obstructive pulmonary disease. Proceedings of the American Thoracic Society 3, 726–733, https://doi.org/10.1513/pats.200605-126SF (2006).
Article CAS PubMed Google Scholar
Perotin, J. M. et al. Delay of airway epithelial wound repair in COPD is associated with airflow obstruction severity. Respiratory research 15, 151, https://doi.org/10.1186/s12931-014-0151-9 (2014).
Article PubMed PubMed Central Google Scholar
Saetta, M. et al. Goblet cell hyperplasia and epithelial inflammation in peripheral airways of smokers with both symptoms of chronic bronchitis and chronic airflow limitation. American journal of respiratory and critical care medicine 161, 1016–1021, https://doi.org/10.1164/ajrccm.161.3.9907080 (2000).
Article CAS PubMed Google Scholar
Lumsden, A. B., McLean, A. & Lamb, D. Goblet and Clara cells of human distal airways: evidence for smoking induced changes in their numbers. Thorax 39, 844–849 (1984).
Article CAS Google Scholar
Polosukhin, V. V. et al. Secretory IgA Deficiency in Individual Small Airways Is Associated with Persistent Inflammation and Remodeling. American journal of respiratory and critical care medicine 195, 1010–1021, https://doi.org/10.1164/rccm.201604-0759OC (2017).
Article CAS PubMed PubMed Central Google Scholar
Yaghi, A. & Dolovich, M. B. Airway Epithelial Cell Cilia and Obstructive Lung Disease. Cells 5, https://doi.org/10.3390/cells5040040 (2016).
Rennard, S. I. Pathogenesis of COPD. Clinical cornerstone 5, 11–16 (2003).
Article Google Scholar
Hessel, J. et al. Intraflagellar transport gene expression associated with short cilia in smoking and COPD. PLoS One 9, e85453, https://doi.org/10.1371/journal.pone.0085453 (2014).
Article ADS CAS PubMed PubMed Central Google Scholar
Lam, H. C. et al. Histone deacetylase 6-mediated selective autophagy regulates COPD-associated cilia dysfunction. The Journal of clinical investigation 123, 5212–5230, https://doi.org/10.1172/jci69636 (2013).
Article CAS PubMed PubMed Central Google Scholar
Yaghi, A., Zaman, A., Cox, G. & Dolovich, M. B. Ciliary beating is depressed in nasal cilia from chronic obstructive pulmonary disease subjects. Respiratory medicine 106, 1139–1147, https://doi.org/10.1016/j.rmed.2012.04.001 (2012).
Article PubMed Google Scholar
Park, H. Y. et al. Club cell protein 16 and disease progression in chronic obstructive pulmonary disease. American journal of respiratory and critical care medicine 188, 1413–1419, https://doi.org/10.1164/rccm.201305-0892OC (2013).
Article CAS PubMed PubMed Central Google Scholar
Lomas, D. A. et al. Evaluation of serum CC-16 as a biomarker for COPD in the ECLIPSE cohort. Thorax 63, 1058–1063, https://doi.org/10.1136/thx.2008.102574 (2008).
Article CAS PubMed Google Scholar
Gohy, S. T. et al. Imprinting of the COPD airway epithelium for dedifferentiation and mesenchymal transition. The European respiratory journal 45, 1258–1272, https://doi.org/10.1183/09031936.00135814 (2015).
Article PubMed Google Scholar
Hogg, J. C., Macklem, P. T. & Thurlbeck, W. M. Site and nature of airway obstruction in chronic obstructive lung disease. The New England journal of medicine 278, 1355–1360, https://doi.org/10.1056/nejm196806202782501 (1968).
Article CAS PubMed Google Scholar
Aghapour, M., Raee, P., Moghaddam, S. J., Hiemstra, P. S. & Heijink, I. H. Airway Epithelial Barrier Dysfunction in Chronic Obstructive Pulmonary Disease: Role of Cigarette Smoke Exposure. American journal of respiratory cell and molecular biology 58, 157–169, https://doi.org/10.1165/rcmb.2017-0200TR (2018).
Article CAS PubMed Google Scholar
Heijink, I. H., Noordhoek, J. A., Timens, W., van Oosterhout, A. J. & Postma, D. S. Abnormalities in airway epithelial junction formation in chronic obstructive pulmonary disease. American journal of respiratory and critical care medicine 189, 1439–1442, https://doi.org/10.1164/rccm.201311-1982LE (2014).
Article PubMed Google Scholar
Aghapour, M., Raee, P., Moghaddam, S. J., Hiemstra, P. S. & Heijink, I. H. Airway Epithelial Barrier Dysfunction in COPD: Role of Cigarette Smoke Exposure. American journal of respiratory cell and molecular biology, https://doi.org/10.1165/rcmb.2017-0200TR (2017).
Rigden, H. M. et al. Squamous Metaplasia Is Increased in the Bronchial Epithelium of Smokers with Chronic Obstructive Pulmonary Disease. PLoS One 11, e0156009, https://doi.org/10.1371/journal.pone.0156009 (2016).
Article CAS PubMed PubMed Central Google Scholar
Araya, J. et al. Squamous metaplasia amplifies pathologic epithelial-mesenchymal interactions in COPD patients. The Journal of clinical investigation 117, 3551–3562, https://doi.org/10.1172/jci32526 (2007).
Article CAS PubMed PubMed Central Google Scholar
Polosukhin, V. V. et al. Bronchial secretory immunoglobulin a deficiency correlates with airway inflammation and progression of chronic obstructive pulmonary disease. American journal of respiratory and critical care medicine 184, 317–327, https://doi.org/10.1164/rccm.201010-1629OC (2011).
Article CAS PubMed PubMed Central Google Scholar
Yi, G. et al. A large lung gene expression study identifying IL1B as a novel player in airway inflammation in COPD airway epithelial cells. Inflammation research: official journal of the European Histamine Research Society… [et al.], https://doi.org/10.1007/s00011-018-1145-8 (2018).
Jafarinejad, H. et al. Worldwide prevalence of viral infection in AECOPD patients: A meta-analysis. Microbial pathogenesis 113, 190–196, https://doi.org/10.1016/j.micpath.2017.10.021 (2017).
Article PubMed Google Scholar
Leung, J. M. et al. The role of acute and chronic respiratory colonization and infections in the pathogenesis of COPD. Respirology (Carlton, Vic.) 22, 634–650, https://doi.org/10.1111/resp.13032 (2017).
Article Google Scholar
Schamberger, A. C., Staab-Weijnitz, C. A., Mise-Racek, N. & Eickelberg, O. Cigarette smoke alters primary human bronchial epithelial cell differentiation at the air-liquid interface. Scientific reports 5, 8163, https://doi.org/10.1038/srep08163 (2015).
Article ADS CAS PubMed PubMed Central Google Scholar
Maunders, H., Patwardhan, S., Phillips, J., Clack, A. & Richter, A. Human bronchial epithelial cell transcriptome: gene expression changes following acute exposure to whole cigarette smoke in vitro. American journal of physiology. Lung cellular and molecular physiology 292, L1248–1256, https://doi.org/10.1152/ajplung.00290.2006 (2007).
Article CAS PubMed Google Scholar
Heijink, I. H., Brandenburg, S. M., Postma, D. S. & van Oosterhout, A. J. Cigarette smoke impairs airway epithelial barrier function and cell-cell contact recovery. The European respiratory journal 39, 419–428, https://doi.org/10.1183/09031936.00193810 (2012).
Article CAS PubMed Google Scholar
Kuehn, D. et al. Impact assessment of repeated exposure of organotypic 3D bronchial and nasal tissue culture models to whole cigarette smoke. Journal of visualized experiments: JoVE, https://doi.org/10.3791/52325 (2015).
Aufderheide, M., Scheffler, S., Ito, S., Ishikawa, S. & Emura, M. Ciliatoxicity in human primary bronchiolar epithelial cells after repeated exposure at the air-liquid interface with native mainstream smoke of K3R4F cigarettes with and without charcoal filter. Experimental and toxicologic pathology: official journal of the Gesellschaft fur Toxikologische Pathologie 67, 407–411, https://doi.org/10.1016/j.etp.2015.04.006 (2015).
Article CAS Google Scholar
Tatsuta, M. et al. Effects of cigarette smoke on barrier function and tight junction proteins in the bronchial epithelium: protective role of cathelicidin LL-37. Respiratory research 20, 251, https://doi.org/10.1186/s12931-019-1226-4 (2019).
Article CAS PubMed PubMed Central Google Scholar
Yu, Q., Chen, X., Fang, X., Chen, Q. & Hu, C. Caveolin-1 aggravates cigarette smoke extract-induced MUC5AC secretion in human airway epithelial cells. International journal of molecular medicine 35, 1435–1442, https://doi.org/10.3892/ijmm.2015.2133 (2015).
Article CAS PubMed Google Scholar
Brekman, A., Walters, M. S., Tilley, A. E. & Crystal, R. G. FOXJ1 prevents cilia growth inhibition by cigarette smoke in human airway epithelium in vitro. American journal of respiratory cell and molecular biology 51, 688–700, https://doi.org/10.1165/rcmb.2013-0363OC (2014).
Article CAS PubMed PubMed Central Google Scholar
Burgel, P. R. et al. Update on the roles of distal airways in COPD. European respiratory review: an official journal of the European Respiratory Society 20, 7–22, https://doi.org/10.1183/09059180.10010610 (2011).
Article Google Scholar
Breiman, L. Random Forests. Machine Learning 45, 5–32, https://doi.org/10.1023/A:1010933404324 (2001).
Article MATH Google Scholar
Schlingmann, B., Molina, S. A. & Koval, M. Claudins: Gatekeepers of lung epithelial function. Seminars in cell & developmental biology 42, 47–57, https://doi.org/10.1016/j.semcdb.2015.04.009 (2015).
Article CAS Google Scholar
Tilley, A. E. et al. Biologic phenotyping of the human small airway epithelial response to cigarette smoking. PLoS One 6, e22798, https://doi.org/10.1371/journal.pone.0022798 (2011).
Article ADS CAS PubMed PubMed Central Google Scholar
Higham, A., Quinn, A. M., Cancado, J. E. D. & Singh, D. The pathology of small airways disease in COPD: historical aspects and future directions. Respiratory research 20, 49, https://doi.org/10.1186/s12931-019-1017-y (2019).
Article PubMed PubMed Central Google Scholar
Min, K. A., Rosania, G. R. & Shin, M. C. Human Airway Primary Epithelial Cells Show Distinct Architectures on Membrane Supports Under Different Culture Conditions. Cell biochemistry and biophysics 74, 191–203, https://doi.org/10.1007/s12013-016-0719-8 (2016).
Article CAS PubMed PubMed Central Google Scholar
Koval, M. Claudin heterogeneity and control of lung tight junctions. Annual review of physiology 75, 551–567, https://doi.org/10.1146/annurev-physiol-030212-183809 (2013).
Article CAS PubMed Google Scholar
Soini, Y. Claudins in lung diseases. Respiratory research 12, 70, https://doi.org/10.1186/1465-9921-12-70 (2011).
Article CAS PubMed PubMed Central Google Scholar
Gunzel, D. & Yu, A. S. Claudins and the modulation of tight junction permeability. Physiological reviews 93, 525–569, https://doi.org/10.1152/physrev.00019.2012 (2013).
Article CAS PubMed PubMed Central Google Scholar
Sweerus, K. et al. Claudin-18 deficiency is associated with airway epithelial barrier dysfunction and asthma. The Journal of allergy and clinical immunology 139, 72–81.e71, https://doi.org/10.1016/j.jaci.2016.02.035 (2017).
Article CAS PubMed Google Scholar
Gunzel, D. et al. Claudin-10 exists in six alternatively spliced isoforms that exhibit distinct localization and function. Journal of cell science 122, 1507–1517, https://doi.org/10.1242/jcs.040113 (2009).
Article CAS PubMed Google Scholar
Shaykhiev, R. et al. Cigarette smoking reprograms apical junctional complex molecular architecture in the human airway epithelium in vivo. Cellular and molecular life sciences: CMLS 68, 877–892, https://doi.org/10.1007/s00018-010-0500-x (2011).
Article CAS PubMed Google Scholar
Zen, K. et al. JAM-C is a component of desmosomes and a ligand for CD11b/CD18-mediated neutrophil transepithelial migration. Molecular biology of the cell 15, 3926–3937, https://doi.org/10.1091/mbc.e04-04-0317 (2004).
Article CAS PubMed PubMed Central Google Scholar
Barnes, P. J. & Adcock, I. M. Chronic obstructive pulmonary disease and lung cancer: a lethal association. American journal of respiratory and critical care medicine 184, 866–867, https://doi.org/10.1164/rccm.201108-1436ED (2011).
Article PubMed Google Scholar
Papi, A. et al. COPD increases the risk of squamous histological subtype in smokers who develop non-small cell lung carcinoma. Thorax 59, 679–681, https://doi.org/10.1136/thx.2003.018291 (2004).
Article CAS PubMed PubMed Central Google Scholar
de Torres, J. P. et al. Lung cancer in patients with chronic obstructive pulmonary disease–incidence and predicting factors. American journal of respiratory and critical care medicine 184, 913–919, https://doi.org/10.1164/rccm.201103-0430OC (2011).
Article PubMed Google Scholar
Auerbach, O., Stout, A. P., Hammond, E. C. & Garfinkel, L. Changes in bronchial epithelium in relation to cigarette smoking and in relation to lung cancer. The New England journal of medicine 265, 253–267, https://doi.org/10.1056/nejm196108102650601 (1961).
Article CAS PubMed Google Scholar
Ghosh, M. et al. Exhaustion of Airway Basal Progenitor Cells in Early and Established COPD. American journal of respiratory and critical care medicine, https://doi.org/10.1164/rccm.201704-0667OC (2017).
Rokicki, W., Rokicki, M., Wojtacha, J. & Dzeljijli, A. The role and importance of club cells (Clara cells) in the pathogenesis of some respiratory diseases. Kardiochirurgia i torakochirurgia polska = Polish journal of cardio-thoracic surgery 13, 26–30, https://doi.org/10.5114/kitp.2016.58961 (2016).
Article PubMed PubMed Central Google Scholar
Kim, V. et al. Chronic bronchitis and current smoking are associated with more goblet cells in moderate to severe COPD and smokers without airflow obstruction. PLoS One 10, e0116108, https://doi.org/10.1371/journal.pone.0116108 (2015).
Article CAS PubMed PubMed Central Google Scholar
Cosio, M. G., Hale, K. A. & Niewoehner, D. E. Morphologic and morphometric effects of prolonged cigarette smoking on the small airways. The American review of respiratory disease 122, 265–221, https://doi.org/10.1164/arrd.1980.122.2.265 (1980).
Article CAS PubMed Google Scholar
Yang, J. et al. Smoking-Dependent Distal-to-Proximal Repatterning of the Adult Human Small Airway Epithelium. American journal of respiratory and critical care medicine 196, 340–352, https://doi.org/10.1164/rccm.201608-1672OC (2017).
Article CAS PubMed PubMed Central Google Scholar
Shijubo, N. et al. Serum and BAL Clara cell 10 kDa protein (CC10) levels and CC10-positive bronchiolar cells are decreased in smokers. The European respiratory journal 10, 1108–1114 (1997).
Article CAS Google Scholar
Laucho-Contreras, M. E. et al. Protective role for club cell secretory protein-16 (CC16) in the development of COPD. The European respiratory journal 45, 1544–1556, https://doi.org/10.1183/09031936.00134214 (2015).
Article PubMed PubMed Central Google Scholar
Knabe, L., Fort, A., Chanez, P. & Bourdin, A. Club cells and CC16: another “smoking gun”? (With potential bullets against COPD). The European respiratory journal 45, 1519–1520, https://doi.org/10.1183/09031936.00010515 (2015).
Article PubMed Google Scholar
Wanner, A., Salathe, M. & O’Riordan, T. G. Mucociliary clearance in the airways. American journal of respiratory and critical care medicine 154, 1868–1902, https://doi.org/10.1164/ajrccm.154.6.8970383 (1996).
Article CAS PubMed Google Scholar
Randell, S. H. & Boucher, R. C. Effective mucus clearance is essential for respiratory health. American journal of respiratory cell and molecular biology 35, 20–28, https://doi.org/10.1165/rcmb.2006-0082SF (2006).
Article CAS PubMed PubMed Central Google Scholar
Nebert, D. W. et al. Role of the aromatic hydrocarbon receptor and [Ah] gene battery in the oxidative stress response, cell cycle control, and apoptosis. Biochemical pharmacology 59, 65–85 (2000).
Article CAS Google Scholar
Huang, M. F., Lin, W. L. & Ma, Y. C. A study of reactive oxygen species in mainstream of cigarette. Indoor air 15, 135–140, https://doi.org/10.1111/j.1600-0668.2005.00330.x (2005).
Article CAS PubMed Google Scholar
Guerrina, N., Traboulsi, H., Eidelman, D. H. & Baglole, C. J. The Aryl Hydrocarbon Receptor and the Maintenance of Lung Health. International journal of molecular sciences 19, https://doi.org/10.3390/ijms19123882 (2018).
Yun, J. H. et al. Transcriptomic Analysis of Lung Tissue from Cigarette Smoke-Induced Emphysema Murine Models and Human Chronic Obstructive Pulmonary Disease Show Shared and Distinct Pathways. American journal of respiratory cell and molecular biology 57, 47–58, https://doi.org/10.1165/rcmb.2016-0328OC (2017).
Article CAS PubMed PubMed Central Google Scholar
Sajadimajd, S. & Khazaei, M. Oxidative Stress and Cancer: The Role of Nrf2. Current cancer drug targets 18, 538–557, https://doi.org/10.2174/1568009617666171002144228 (2018).
Article CAS PubMed Google Scholar
Stedman, R. L. The chemical composition of tobacco and tobacco smoke. Chemical reviews 68, 153–207 (1968).
Article CAS Google Scholar
Miller, K. P. & Ramos, K. S. Impact of cellular metabolism on the biological effects of benzo[a]pyrene and related hydrocarbons. Drug metabolism reviews 33, 1–35, https://doi.org/10.1081/dmr-100000138 (2001).
Article CAS PubMed Google Scholar
Ritter, D. & Knebel, J. Investigations of the Biological Effects of Airborne and Inhalable Substances by Cell-Based In Vitro Methods: Fundamental Improvements to the ALI Concept. Advances in Toxicology 2014, 11, https://doi.org/10.1155/2014/185201 (2014).
Article Google Scholar
Ritter, D. et al. Development and Evaluation of an In Vitro Test System for Toxicity Screening of Aerosols Released from Consumer Products and First Application to Aerosols from a Hair Straightening Process. Applied In Vitro Toxicology 4, 180–192, https://doi.org/10.1089/aivt.2017.0036 (2018).
Article CAS Google Scholar
Sollner, J. F. et al. An RNA-Seq atlas of gene expression in mouse and rat normal tissues. Scientific data 4, 170185, https://doi.org/10.1038/sdata.2017.185 (2017).
Article CAS PubMed PubMed Central Google Scholar
Ritchie, M. E. et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic acids research 43, e47, https://doi.org/10.1093/nar/gkv007 (2015).
Article CAS PubMed PubMed Central Google Scholar

Download references

Acknowledgements

The authors thank Nadine Rehm for excellent technical assistance and David J. Lamb for comments that greatly improved the manuscript.

Author information

Authors and Affiliations

Immunology & Respiratory Diseases Research, Boehringer Ingelheim Pharma GmbH & Co. KG, Biberach an der Riß, Germany
Julia A. Gindele, Tobias Kiechle, Kerstin Benediktus, Matthew J. Thomas & Jürgen Schymeinsky
Department of General Physiology, University of Ulm, Ulm, Germany
Julia A. Gindele & Jürgen Schymeinsky
Drug Discovery Sciences, Boehringer Ingelheim Pharma GmbH & Co. KG, Biberach an der Riß, Germany
Gerald Birk, Fabian Heinemann & Birgit Stierstorfer
Non-Clinical Statistics Biberach, Boehringer Ingelheim Pharma GmbH & Co. KG, Biberach an der Riß, Germany
Michael Brendel
Global Computational Biology and Digital Sciences, Boehringer Ingelheim Pharma GmbH & Co. KG, Biberach an der Riß, Germany
Christian T. Wohnhaas & Karsten Quast
Department of Biology, University of Konstanz, Konstanz, Germany
Christian T. Wohnhaas
Department of Genetic Medicine, Weill Cornell Medical College, New York, USA
Michelle LeBlanc, Haijun Zhang, Yael Strulovici-Barel & Ronald G. Crystal

Authors

Julia A. Gindele
View author publications
You can also search for this author in PubMed Google Scholar
Tobias Kiechle
View author publications
You can also search for this author in PubMed Google Scholar
Kerstin Benediktus
View author publications
You can also search for this author in PubMed Google Scholar
Gerald Birk
View author publications
You can also search for this author in PubMed Google Scholar
Michael Brendel
View author publications
You can also search for this author in PubMed Google Scholar
Fabian Heinemann
View author publications
You can also search for this author in PubMed Google Scholar
Christian T. Wohnhaas
View author publications
You can also search for this author in PubMed Google Scholar
Michelle LeBlanc
View author publications
You can also search for this author in PubMed Google Scholar
Haijun Zhang
View author publications
You can also search for this author in PubMed Google Scholar
Yael Strulovici-Barel
View author publications
You can also search for this author in PubMed Google Scholar
Ronald G. Crystal
View author publications
You can also search for this author in PubMed Google Scholar
Matthew J. Thomas
View author publications
You can also search for this author in PubMed Google Scholar
Birgit Stierstorfer
View author publications
You can also search for this author in PubMed Google Scholar
Karsten Quast
View author publications
You can also search for this author in PubMed Google Scholar
Jürgen Schymeinsky
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

J.A.G. and J.S. designed the study. J.A.G., T.K., K.B., G.B. and C.T.W. acquired the data. J.A.G., T.K., G.B., M.B., F.H. and K.Q. analyzed the data. J.A.G., B.S., J.S. and K.Q. discussed the results. J.A.G., J.S. and K.Q. wrote the manuscript. R.G.C., M.J.T., B.S., M.L., H.Z., Y.S. and all authors consulted on the final version.

Corresponding author

Correspondence to Jürgen Schymeinsky.

Ethics declarations

Competing interests

T.K., K.B., G.B., M.B., F.H., M.J.T., B.S., K.Q. and J.S. are employees of Boehringer Ingelheim Pharma GmbH & Co. K.G. J.A.G. and C.T.W. receive grant support from Boehringer Ingelheim Pharma GmbH & Co. K.G. The funder provided support in the form of salaries for authors [T.K., K.B., G.B., M.B., F.H., M.J.T., B.S., K.Q. and J.S.] and research materials, but did not have any additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Supplementary Information.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Gindele, J.A., Kiechle, T., Benediktus, K. et al. Intermittent exposure to whole cigarette smoke alters the differentiation of primary small airway epithelial cells in the air-liquid interface culture. Sci Rep 10, 6257 (2020). https://doi.org/10.1038/s41598-020-63345-5

Download citation

Received: 29 August 2019
Accepted: 30 March 2020
Published: 10 April 2020
DOI: https://doi.org/10.1038/s41598-020-63345-5

This article is cited by

Alterations in the molecular control of mitochondrial turnover in COPD lung and airway epithelial cells
- Christy B. M. Tulen
- Cheryl van de Wetering
- Alexander H. V. Remels
Scientific Reports (2024)
Azithromycin ameliorated cigarette smoke-induced airway epithelial barrier dysfunction by activating Nrf2/GCL/GSH signaling pathway
- Yun Song
- Wenhuan Fu
- Bin Wang
Respiratory Research (2023)
Effect of sub-chronic exposure to cigarette smoke, electronic cigarette and waterpipe on human lung epithelial barrier function
- Baishakhi Ghosh
- Hermes Reyes-Caballero
- Venkataramana K. Sidhaye
BMC Pulmonary Medicine (2020)
Development of a miniaturized 96-Transwell air–liquid interface human small airway epithelial model
- Teresa Bluhmki
- Sarah Bitzer
- Ralf Heilker
Scientific Reports (2020)

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.