- Open Access
Differential gene expression of 3D primary human airway cultures exposed to cigarette smoke and electronic nicotine delivery system (ENDS) preparations
BMC Medical Genomics volume 15, Article number: 76 (2022)
Acute exposure to cigarette smoke alters gene expression in several biological pathways such as apoptosis, immune response, tumorigenesis and stress response, among others. However, the effects of electronic nicotine delivery systems (ENDS) on early changes in gene expression is relatively unknown. The objective of this study was to evaluate the early toxicogenomic changes using a fully-differentiated primary normal human bronchial epithelial (NHBE) culture model after an acute exposure to cigarette and ENDS preparations.
RNA sequencing and pathway enrichment analysis identified time and dose dependent changes in gene expression and several canonical pathways when exposed to cigarette preparations compared to vehicle control, including oxidative stress, xenobiotic metabolism, SPINK1 general cancer pathways and mucociliary clearance. No changes were observed with ENDS preparations containing up to 28 µg/mL nicotine. Full model hierarchical clustering revealed that ENDS preparations were similar to vehicle control.
This study revealed that while an acute exposure to cigarette preparations significantly and differentially regulated many genes and canonical pathways, ENDS preparations containing the same concentration of nicotine had very little effect on gene expression in fully-differentiated primary NHBE cultures.
Long-term cigarette smoking impacts lung biology, including induction of inflammation, altered microbial defense, and compromised immune surveillance. Several studies have demonstrated that cigarette smoke alters gene expression in different tissues including bronchial airway epithelia [1, 2], nasal epithelia , whole blood [4, 5] and peripheral blood mononuclear cells [6, 7]. Altered gene expression is associated with different biological processes such as xenobiotic metabolism, apoptosis, immune response, tumorigenesis and stress response, among others. Sustained perturbations in gene expression in key biological processes could drive/contribute to the development of smoking-related diseases such as lung cancer, cardiovascular disease and chronic obstructive pulmonary disease (COPD) [8,9,10,11,12,13].
Electronic nicotine delivery systems (ENDS) are battery operated devices designed to deliver nicotine without combustion . ENDS products in the current marketplace are a highly heterogeneous product category, with several types of devices, a wide range of formulations of e-liquids and variable use settings (voltage, heating temperature, etc.). The aerosol generated by ENDS is less complex than that produced from combustible cigarettes. Cigarette smoke, which is generated by combustion, contains over 7000 chemicals that include several carcinogens, as well as harmful and potentially harmful constituents (HPHCs) . In contrast, e-liquids generally consist of 4 main ingredients (vegetable glycerin, propylene glycol, nicotine and flavorings). However, the composition of ENDS aerosol varies depending on the flavorings and the device use conditions. For example, several investigators reported that ENDS aerosols contain far less quantities of volatile carbonyls, tobacco specific nitrosamines or polycyclic aromatic hydrocarbons, compared to cigarette smoke . However, some investigators also have reported higher levels of some HPHCs in certain flavors and/or when the devices were operated under sub-ohm conditions .
Consistent with the chemistry of aerosols, ENDS users exhibit significantly lower biomarkers of exposure as shown in product-specific studies [17, 18] or population surveys . Clinical findings of ENDS products have indicated that ENDS users have fewer or distinctly different gene expression profiles compared to cigarette smokers [18, 20, 21]. However, some studies report that former smokers who switched to ENDS vaping express shared, as well as distinct, patterns of gene expression when compared to cigarette smoking . In addition, ENDS use has also shown impact on virus responses, including decreased gene expression .
As alluded to above, ENDS are diverse products, and the toxicogenomic effects of ENDS are still being understood. For rapid evaluation of ENDS effects, suitable in vitro systems would be necessary. The lung is a complex organ that consists of several highly differentiated cell types (e.g. ciliated cells and goblet cells), which facilitates mucociliary clearance and serves as the first-line of defense against pathogens and pollution. Therefore, in vitro models that assess lung function should ideally consist of differentiated cells that capture the functions of airway cell types. Existing monolayer cell culture models, although useful, are not necessarily differentiated; hence they may not recapitulate the functions of lung cell types. The differences in cell culture conditions (monolayer vs organoids), and differentiation status (undifferentiated vs differentiated cells) hence can yield different results.
Given the general complexity of culture conditions, there have been some studies investigating the effects of ENDS exposure in 3D differentiated airway cultures [23,24,25], which are more physiologically relevant than monolayer cultures. From the few 3D organoid studies, ENDS products were reported to alter significantly fewer or distinctly different gene expression profiles compared to cigarette smoke in bronchial and nasal cell cultures, similar to that reported in clinical studies [23,24,25]. To our knowledge, only one study to date used differentiated primary normal human bronchial epithelial (NHBE) cells at air–liquid interface (ALI), and reported differential gene expressions after ENDS exposure compared to cigarette smoke . Other studies investigating the effect of ENDS/vaping on gene expression in differentiated NHBE cultures either investigated specific ingredients/flavorings  or used cell-lines at ALI such as BEAS-2B cell-line . MucilAir™ (Epithelix) primary nasal cells have been directly exposed to ENDS aerosol and 3R4F Kentucky reference cigarette smoke, revealing that ENDS aerosol had reduced effect on gene expression [24, 25]. However, other studies have reported extensive reprogramming of gene expression, suggesting that ENDS use is associated with inflammatory responses . Further investigation into the effects of cigarette and ENDS products on gene expression in fully-differentiated primary NHBE cultures are therefore required.
To address some of the challenges associated with assessing lung function, we have demonstrated the applicability of 3D lung culture models in evaluating the acute and longer-term effects of exposure to cigarette smoke and ENDS aerosol [29, 30]. Fully-differentiated primary NHBE cultures can be used to study the effect of cigarette smoke and ENDS preparations on airway mucociliary function, ion channel function, and physiological changes to the epithelium . In this manuscript, to better understand the acute effects (early changes) of ENDS exposure, we utilized a transcriptomic approach with the fully-differentiated primary NHBE culture model and compared them with the profiles derived from cigarette smoke exposure. These data further support and build on our previous findings that cigarette smoke preparations, not ENDS, alter/disrupt targeted ion channel and mucociliary function, and perturb pathways related to xenobiotic and oxidative stress metabolisms.
Exposure system on fully-differentiated primary NHBE cultures
In an effort to understand how exposures to cigarette or ENDS preparations impact gene expression in fully-differentiated primary NHBE cultures, we first examined whether treatment variables (treatment, concentration, time, donors) affected RNA-seq results (see Methods). Primary NHBE cultures from 4 donors (Additional file 2: Table S1) were fully-differentiated after 4 weeks at ALI, as determined by optimal trans-epithelial electrical resistance (TEER), and presence of ciliated and mucus producing cells . Based on previously published data using lactate dehydrogenase (LDH) release as an indication of cytotoxicity , non-cytotoxic doses of whole-smoke conditioned medium (WS-CM) from cigarette preparations, and aerosol conditioned medium (ACM) from ENDS preparations were chosen. Low, medium, and high doses of equivalent-nicotine units (Eq-Nic.) were designated for WS-CM (3.6, 7.0 and 10.0 µg/mL Eq-Nic.) and for ACM (7.0, 14.0 and 28.0 µg/mL Eq-Nic.). Higher doses of ACM were used because no cytotoxicity was observed even at 28 µg/ml Eq-Nic units. .
Next generation sequencing (NGS) analysis
NGS analysis yielded a dataset of 60,564 gene identifiers that were mapped to the Ingenuity Knowledge Base. Any duplicate identifiers were resolved to their respective genes, resulting in a total of 58,457 genes for the analysis. Samples collected at time 0 h were used for normalization.
Using unsupervised clustering identified treatment period (vehicle 0 h vs. 4 h vs. 24 h) had a strong effect on the gene expression since all the samples were primarily grouped based on exposure time (data not shown). Despite the effect of time, treatments (WS-CM and ACM) exerted an effect on gene expression in this in vitro model system. Therefore, the statistical modeling focused on different treatments/dosages within each time point to identify the differentially expressed genes (DEGs).
WS-CM treatment exerts stronger effect on gene expression, whereas ACM elicited minimal response
WS-CM exerted marked effects on gene expression in a dose- and time-dependent fashion. Compared to vehicle control, treatment with WS-CM at low, medium and high doses of Eq-Nic. units altered expression of 6, 76, and 153 genes, respectively, after 4 h exposure using a criteria cut-off of log2 fold change > 2 or < − 2, and adjusted p-value < 0.05 (Table 1). The number of significantly regulated genes increased after a longer exposure (24 h) of WS-CM, with the low, medium and high doses of WS-CM producing 189, 281 and 388 differentially regulated genes, respectively (Table 1).
In contrast, we did not observe dosage or time effect on gene expression upon treatment with ACM, per the criteria defined for statistical analysis (Table 1). Treatment with ACM resulted in no significant gene expression changes at the 3 doses tested (7.0, 14.0, and 28.0 μg/mL Eq-Nic.), which are similar or higher than the corresponding doses used for WS-CM treatment. However, when the log2 fold change cut off was lowered to > 1.5 or < -1.5 and adjusted p ≤ 0.05, the highest dose of ACM (28 µg/mL Eq-Nic.) at 4 h revealed differential expression of 8 genes (Additional file 2: Table S2). Nevertheless, even under the reduced threshold conditions, no DEGs were detected at 24 h of treatment. Furthermore, only with a consideration of a p-value ≤ 0.05, and no fold-change criteria applied, ACM treatment at 24 h did not result in the detection of any DEGs.
Comparing across the three different doses of WS-CM, there were a few common genes at 4 h (6 genes) exposure, and 24 h (136 genes) (Fig. 1). The six differentially expressed genes (AHRR-196, AHRR-438, ALDH1A3, LPAR6, TRBC1 and TRBC2) after 4 h of low dose WS-CM exposure were also found to be altered at medium and high doses (Additional file 1). Only AHRR-196, ALDH1A3 and LPAR6 were also sustained after 24 h exposure (Additional file 1). If we only compared the medium and high doses at 4 h WS-CM exposure, 69 genes were common. When we compared the two different time-points, 4 h and 24 h WS-CM exposure, 4 and 24 genes were regulated by the medium and high dose at both time-points, respectively (Additional file 1). A set of 21 DEGs were identified in all doses and across both time-points except for low dose at 4 h. Amongst these 21 DEGs, several were increasingly upregulated with higher dose and longer exposure time (e.g. HMOX1, CYP1A1, CYP1B1, CYP1B1-AS1), whilst some others did not vary in expression between the two exposure times (e.g. TIPARP, C5AR2, HILPDA, EGLN3). A couple of genes had decreased regulation at 24 h compared to 4 h, including downregulation of ANGPTL4 and PGF, suggesting that these DEGs were more strongly regulated in a shorter time period (4 h) which has begun to return to base levels by 24 h.
Hierarchical clustering was performed to examine if DEGs would classify samples into vehicle and WS-CM-treated groups. The clustering analyses of DEGs (statistically significant p < 0.05 and ± 2 log2 fold change) allowed a clear separation of vehicle and WS-CM (Fig. 2).
Since most of the doses of ACM tested were higher than used with WS-CM, we compared the effects of WS-CM and ACM at a common dose of 7.0 μg/ml Eq-Nic. At this dose, several genes were significantly altered by WS-CM, but not ACM at 4 h or 24 h (Fig. 3, Additional file 2: Figs. S2B and S2C). WS-CM treatment resulted in 146 DEGs, and in 303 DEGs at 4 h and 24 h, respectively, relative to ACM treatment (Additional file 2: Fig. S2A). A second method of comparison was used to determine the effect of WS-CM compared to ACM using hierarchical clustering (Additional file 2: Fig. S2B). We then attempted to identify the overall effect of treatment (WS-CM vs ACM) with time effect removed. When running full model gene expression differences between WS-CM and ACM, we identified 59 DEGs (adjusted p-value < 1E-5) that could be used as marker genes to differentiate between groups treated with WS-CM or ACM (Additional file 2: Table S3). Hierarchical clustering revealed that while all ACM-treated samples displayed a pattern similar to the control samples, the WS-CM treated samples were distinct and co-clustered (Fig. 3).
Top DEGs identified with WS-CM exposure include those associated with phase I and II enzymes and oxidative stress
Treatment with WS-CM resulted in significant changes in gene expression at 4 h and 24 h, as previously shown in Table 1. The top 20 DEGs resulting from high dose treatment are shown in Fig. 4A. Metallothionein genes were most prominently induced by WS-CM after 24 h (Fig. 4B). Several xenobiotic metabolizing genes, including those responsible for the expression of phase I and phase II enzymes were also prominently induced by WS-CM. The cytochrome P450 genes, CYP1A1, CYP1B1, and AHRR (aryl hydrocarbon receptor repressor) genes were among those significantly upregulated by WS-CM (Fig. 4C). Genes involved in the regulation of oxidative stress such as HMOX1 (hemoxygenase 1), GPX (glutathione peroxidase) and TXN (thioredoxin) were also induced by WS-CM treatment (Fig. 4D). Significant upregulation of HMOX1 was observed at both 4 h and 24 h exposure to WS-CM (Additional file 2: Tables S4 and S5). Furthermore, genes involved in mucociliary function including ion channels responsible for fluid homeostasis and genes regulating differentiation of goblet cells and mucin production were also differentially regulated by WS-CM at 24 h (Fig. 4E). While genes associated with goblet cell hyperplasia and mucus secretion were up-regulated (FOXA3, SPDEF, MUC5AC and MUC2), genes encoding ion channels were down-regulated. Thus, SCNN1G, which encodes for the gamma subunit of the epithelial sodium channel (ENaC), was significantly down-regulated. If the log2 fold change was relaxed to ± 1.5, the beta subunit of ENaC, SCNN1B, was also significantly down-regulated. Another ion channel, cystic fibrosis transmembrane conductance regulator (CFTR) gene was significantly downregulated with cut-off log2 fold change ± 1.5. The changes to the genes associated with mucociliary function were not observed at 4 h exposure to WS-CM preparations.
Top canonical pathways affected by WS-CM
We wanted to identify significantly associated canonical pathways, predicted upstream regulators and the top predicted diseases, and functions associated with each dataset. We used Ingenuity Pathway Analysis (IPA) to identify the differential expression analysis (see Methods for details). For the comparison analyses, the core analyses were collated and common threshold cut-offs of adjusted p-value < 0.05 and a log2 fold change threshold ± 1.5 were used to identify significant differentially expressed genes in each individual analysis. When visualizing the comparison analysis results, a Benjamini–Hochberg adjusted p-value < 0.05 and a Z-score of > ± 2 was used to determine significance.
The top canonical pathways identified with WS-CM (medium and high doses) compared to vehicle control after 4 h in this comparison analysis included two pathways regulating aspects of Melatonin and Nicotine Degradation (Nicotine Degradation II and Nicotine Degradation III) (Fig. 5A). Several Phase I and II enzymes are upregulated in these pathways after 4 h, including CYP1A1, CYP1B1, UTG1A6 and UGT1A7 (Additional file 2: Table S4). A xenobiotic metabolism pathway was identified at both 4 h and 24 h, however there were not enough significantly regulated genes to determine whether this pathway was activated or inhibited.
After 24 h treatment of WS-CM, the top canonical pathways identified in this analysis included the SPINK1 general cancer pathway, nicotine degradation III, melatonin degradation and NRF-2 oxidative stress response (Fig. 5B). SPINK1 cancer pathway was predicted to be inhibited in all 3 comparisons of WS-CM vs vehicle control at various dosages (Fig. 5B). Several metallothionein genes, including MT1E, MT1G and MT1M were upregulated in the SPINK1 general cancer pathway (Additional file 2: Table S5). Other cancer related pathways that were altered included inhibition of gluconeogenesis I, glycolysis I, basal cell carcinoma signaling and small cell lung cancer signaling (activated) (Fig. 5B). The NRF2-mediated oxidative stress pathway was activated in all doses of WS-CM (Fig. 5B). Several oxidative stress related genes were up-regulated, including HMOX1, GPX2, GCLM and TXN, as well as down-regulated genes including MAF and RAP2B (Additional file 2: Table S5). Similar to what was observed at 4 h, nicotine and melatonin degradation pathways were activated after 24 h exposure, with several Phase I and II enzymes upregulated including CYP1A1, CYP1B1, and UGT1A6 (Additional file 2: Table S5).
The objective of this study was to evaluate the early toxicogenomic changes after an acute exposure (4 and 24 h) of cigarette smoke preparations (WS-CM) and ENDS preparations (ACM) on fully differentiated primary human bronchial cultures. There are two major findings to be taken from this study. Firstly, extensive gene expression changes were detected due to WS-CM exposure in a time- and dose-dependent manner. ENDS preparations (ACM), on the other hand, exerted minimal effect, and the changes in gene expression resembled that of the vehicle control. Secondly, several cellular processes were differentially regulated by WS-CM including general cancer pathways and xenobiotic metabolism, oxidative stress and genes associated with mucociliary function. In contrast, no significant changes in cellular pathways were detected with ACM treatment.
As reported in this study, an acute exposure of WS-CM preparations exerted more pronounced transcriptomic changes in fully differentiated primary NHBE cultures compared to ACM. Even if the log2 fold change cut-off was lowered to ± 1.5 in this current study, only eight DEGs were identified at only the highest dose of ACM after 4 h exposure, but were not persistently seen after 24 h exposure. Other studies have demonstrated that ENDS aerosol induces less pronounced gene expression changes than cigarette smoke in primary NHBE cells; genes that were regulated by ENDS aerosol were related to oxidative stress, xenobiotic metabolism, and genes involved in cilia assembly and movement [23, 24, 32]. Two main factors may play a role in the differential effects observed between WS-CM and ACM: the chemical composition present in the preparations, and/or the formation of chemicals from the combustion of tobacco. It has been shown that heating tobacco products (HTPs), which rely on heating tobacco rather than combusting tobacco, have reduced impact on gene expression compared to 3R4F reference cigarette smoke .
Few studies have used fully-differentiated primary NHBE cultures for investigation of transcriptomic changes as a result of cigarette or ENDS exposure. We have expanded the use of in vitro airway models by investigating the global transcriptomic changes after an acute exposure to WS-CM and ACM preparations. Shen et al.  observed expression changes in CYP1A1 and HMOX1 from direct cigarette exposure (1R5F cigarette smoke) in fully-differentiated primary NHBE, similar to what we also observed in this study. Furthermore, other studies have reported regulation of genes and/or proteins by cigarette smoke in in vitro models, including metallothionein (MT1G, MT1M, MT1A) , CYP1B1 , AHRR , GPX2, CFTR [29, 37, 38], FOXA3, SPDEF [39, 40] and MUC5AC [29, 40]. On the other hand, expression from ENDS aerosol was distinctly different from cigarette smoke, with phospholipid and fatty acid triacylglycerol metabolism pathways significantly regulated . Another in vitro study of e-cigarette liquids and aerosol reported induction of GCLM, GCLC, GPX2, NQO1 and HMOX1 . Comparisons to these studies are difficult since we exposed the cells with ENDS preparation (aerosol bubbled through media) rather than aerosol or e-cigarette liquid. However, a couple of in vitro studies with e-cigarette aerosols have reported lower cytotoxicity and minimal changes in the mRNA, microRNA and protein markers [42, 43], supporting our findings with ENDS preparations. Using a less stringent cut-off of log2 FC ± 1.5, we observed a couple of nuclear proteins (Hes Related Family BHLH Transcription Factor with YRPW Motif 1 [HEY1] and Serum Response Factor [SRF]) with ACM (28 µg/mL Eq-Nic.) at 4 h. The other genes appear to code for proteins related to cellular proliferation and structure, including Inertermediate Filament Family Orphan 2 (IFFO2), Cellular Communication Network Factor 2 (CCN2), FosB Proto-oncogene, AP-1 Transcription Factor Subunit (FOSB), and palladin, cytoskeletal associated protein (PALLD). No pathways were identified due to the low number of genes regulated by ACM. Higher concentrations of nicotine in the ACM may have stronger effect on gene expression, and is a question for future studies.
There are several strengths in our present study, including that we used cells from four different donors, whilst other publications use 1–2 donors, reducing the possible variabilities between different donors. Another factor of our study is that we used bubbled aerosol/smoke so that we could control for nicotine concentration for side-by-side comparisons of cigarette and ENDS preparations. Furthermore, the findings from these in vitro fully-differentiated primary NHBE cultures have been reported in smokers. For example, current smokers, compared to never smokers, exhibit significantly differentially regulated genes in airway epithelial cells recovered from bronchoscopy associated with oxidative stress (GPX2 and ALDH3A1), and glutathione and xenobiotic metabolism (CYP1B1 and DBDD) . On the other hand, few clinical studies have investigated the effect of ENDS on gene expression; of those reported, there seems to be fewer gene expression changes as a result of ENDS vaping [20, 45]. One recent investigation observed no significant changes in mRNA or miRNA gene expression in human bronchial epithelial brushings collected from subjects vaping a 50:50 propylene glycol:vegetable glycerin mix for 4 weeks . Another study showed that when cigarette users switched to ENDS, gene expression profiles of the ENDS users were similar to subjects that were former smokers, suggesting that cigarette-induced transcriptional changes revert back to baseline when users switch to ENDS . On the other hand, certain flavoring compounds have been associated with suppressed ciliary beat frequency, mitochondrial respiration, inflammation and immune responses in people that use e-cigarettes [46, 47]. Further investigations into flavoring compounds of ENDS/e-cigarette aerosols is therefore important.
Previously, we have investigated the physiological effects of an acute exposure of WS-CM and ACM on fully-differentiated primary NHBE 3D cultures, including airway mucociliary function, ion channels’ activity, and physiological changes to mucociliary clearance in the airway epithelium [29, 31]. Ion channel function of both ENaC and CFTR was inhibited with cigarette preparations . In this study, WS-CM significantly altered a number of genes associated with regulation of airway mucociliary clearance, including decreased CFTR and ENaC genes (SCNN1B, SCNN1G), increased mucus production (MUC5AC, MUC2B) and regulation of goblet cell differentiation (SPDEF). Altered ion channel function and mucociliary clearance (mucus production) have been strongly linked to lung disease including COPD [48,49,50,51]. FOXA3, another transcription factor, was also upregulated, and has been shown to be co-localized in goblet cells with SPDEF in mouse lung tissue despite that FOXA3 is not expressed abundantly in the normal lung [50, 52]. In this study, we did not observe any genes associated to mucociliary clearance that were altered with even the highest dose (28 µg/mL Eq-Nic.) of ACM preparations, suggesting that the effects observed after cigarette exposure was not similarly observed with ENDS preparations.
Supporting our findings of in vitro fully-differentiated primary NHBE cultures, several clinical studies have also identified several similar pathways that are regulated by cigarette smoke, including xenobiotic and glutathione metabolism, and inflammatory responses [53,54,55,56]. Several canonical pathways were significantly regulated by WS-CM, but not ENDS, in this study. Top pathways affected by WS-CM included upregulation of nicotine degradation and NRF2-mediated oxidative stress response pathways. Oxidative stress is strongly associated with COPD, which in turn can lead to increased risk of lung cancer [53,54,55]. Activation of oxidative stress and xenobiotic metabolism are included as some of the key characteristic end-points caused by agents known to initiate cancer in humans .
Increased expression of metallothionein genes have been associated specifically with acute cigarette smoke exposure, but are suppressed in chronic smokers . Metallothionein, in general, acts as an antioxidant protein. One of the regulators of their expression is Nrf2 signaling pathway (via activation of antioxidant response element), which is a well-known pathway activated by acute cigarette smoke exposure. Nrf2-mediated oxidative stress response pathway was upregulated after 24 h in this study, thus being consistent with metallothionein regulation. A study investigating gene expression in bronchial brushing epithelial cells from individuals who smoked three cigarettes in an acute exposure showed upregulation of several metallothionein genes . In addition, it is known that SPINK1 downregulates metallothionein gene expressions, and therefore metallothionein genes that are involved in the SPINK1 general cancer pathway listed in this study’s IPA canonical pathway analysis. SPINK1 general cancer pathway was significantly inhibited by WS-CM only, with a number of metallothionein genes that were significantly upregulated in this pathway. Metallothionein genes could remain suppressed in former smokers, indicating a persistent risk of lung cancer for these subjects , however we did not observe downregulation of metallothionein in this acute exposure study, suggesting that perhaps longer and/or repeated exposures may be needed to induce downregulation. Taken together, the results observed in this study suggest that upregulation of metallothionein was regulated via Nrf2, and not SPINK1 pathway.
Some limitations of this study needs to be considered. Firstly, no reference ENDS product currently exists. Reported results may vary between studies due to different flavorings, base ingredient composition (e.g. ratio of propylene glycol and vegetable glycerin), nicotine concentrations, and preparations (e.g. whole aerosol, media preparations, e-liquid, etc.). A reference cigarette (3R4F) was considered suitable for cigarette smoke comparisons. Secondly, time of treatment significantly affected DEGs, therefore experimental analysis was required to take this time effect into consideration. Primary NHBE cultures are normally kept at ALI, however WS-CM and ACM preparations were kept on the apical surface of the bronchial epithelium for 4 h and 24 h, to allow longer airway exposure. It is possible that leaving fully-differentiated primary NHBE cultures at liquid/liquid interface caused additional stress which may have masked some of the DEGs changed as a result of WS-CM or ACM exposure. Nevertheless, our study shows more differential regulation of gene expression by cigarette smoke than ENDS. Future studies into long-term gene expression changes in in vitro models should be investigated to determine whether changes are prolonged or observed in a more chronic tobacco exposure model.
We observed a strong differential regulation of gene expression from an acute exposure to cigarette preparations (WS-CM) but not to ENDS preparations (ACM), suggesting a differential response and defense mechanisms towards various tobacco products. Several pathways were regulated by WS-CM including general cancer pathway, oxidative stress, xenobiotic metabolism and mucociliary function. Several gene expression changes related to mucociliary clearance strongly supported physiological changes (decreased CFTR and ENaC ion channels’ activity by WS-CM) we had previously observed in another investigation . The use of fully-differentiated primary NHBE cells for measuring gene expression changes has many advantages since we could use multiple donors (four donors) and achieve a physiological relevant model to study the effects of tobacco products.
Primary human bronchial epithelial cell cultures
Primary NHBE cells were provided by the Nationwide Children’s Hospital Cure Cystic Fibrosis Columbus Epithelial Cell Core, Columbus OH. Cells were isolated from the donor tissues given to the Epithelial Cell Core . Since these cells were provided without identifiers, an exempt status was granted by the Institutional Review Board. The donors were three non-smokers and one with unknown smoking status (see Additional file 2: Table S1 for donor information). Briefly, passage 1 primary NHBE cells were seeded on collagen type IV (0.3 mg/mL; Sigma Aldrich, Saint Louis MO) coated Corning™ Transwells (Fisher Scientific) grown at air–liquid interface (ALI) with PneumaCult™ ALI medium (StemCell Technologies Inc., Tukwila WA) . Medium was changed three times a week and cells were maintained at 37 °C with 5% CO2. Cells were used after 4 weeks at ALI once fully differentiated with presence of ciliated cells, as previously published .
Tobacco cigarette and ENDS preparations
Preparation of whole-smoke conditioned-media (WS-CM) from cigarettes and aerosol conditioned media (ACM) from ENDS vapor has been previously described [29, 58, 59]. Nicotine, tobacco specific nitrosamines, and polycyclic aromatic hydrocarbons from WS-CM and ACM preparations were analyzed as previously described . Briefly, WS-CM was prepared by passing smoke from four 3R4F reference cigarettes through Roswell Park Memorial Institute (RPMI) 1640 medium (without phenol red) using the standard ISO method (35 mL puff volume, 60 s puff interval, 2 s puff duration) . ACM was prepared by passing mainstream aerosol (generated using a refillable 10 watts tank device with a 1.5-Ω coil) in an impinger (25 mL Impinger Midget, Ace Glass) containing RPMI media, using a puff profile of 55 mL puff volume, 30 s puff interval, and 5 s puff duration. The tank was filled with tobacco flavor liquid (1.8% nicotine [weight/volume]). Final nicotine content, measured by gas chromatography-flame ionization detection (GC-FID) of the WS-CM and ACM was used for calculating the exposure of cells, expressed as µg/mL of equivalent-nicotine units (Eq-Nic.) units. The final stock concentrations of nicotine in the WS-CM and ACM was 11.1 µg/mL Eq-Nic. and 116.2 µg/mL Eq-Nic., respectively.
WS-CM and ACM preparations were diluted in Hank’s Balanced Salt Solution (HBSS) containing calcium and magnesium to various concentrations of Eq-Nic. (Table 2). Exposure concentrations were chosen based on our previously reported investigation of acute ENDS and cigarette effects on cytotoxicity in primary NHBE cultures .
In order to standardize the different treatments, preparations contained 3.6–28 µg/ml Eq-Nic. units, as previously reported . Control treatments included diluted vehicle (RPMI 1640) in HBSS. All treatments of cells were performed on separate plates to prevent cross-contamination. Primary NHBE cultures were exposed apically to WS-CM or ACM preparations (100 µL volume) for 0 (no treatment), 4 or 24 h (Fig. 6). All exposures were performed in singlets, with four donors.
Total RNA extraction and purification and sequencing
At treatment endpoints (0, 4 or 24 h), total RNA was isolated using the RNeasy Micro Kit (#74004 Qiagen, Germantown MD). Input RNA quality and quantity was assessed using the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara CA) and Qubit Fluorometer (Thermo Fisher), respectively. Samples with RNA integrity number (RIN) values greater than 7 and RNA concentration greater than 100 ng/µL were sent for sequencing. Sequencing was performed by the Genomics Shared Resources at The Ohio State University, Columbus OH. Messenger RNA (mRNA) sequence libraries were generated with NEBNext® Ultra™ II Directional RNA Library Prep Kit for Illumina (NEB #E7760L) and NEBNext® Poly (A) mRNA Magnetic Isolation Module (NEB #E7490) with an input amount of 200 ng total RNA per sample. Libraries were pooled and sequenced on an Illumina NovaSeq SP flowcell in paired-end 150 bp format (Illumina, San Diego CA) to read yield between 35 and 40 million reads.
Quality control (QC) and alignment
All RNA-seq datasets (72 samples; in fastq.gz format) were sent to Qiagen for data analysis. Raw data QC was performed for each sample using Array Studio (Qiagen, 2019). Statistics performed included total reads, GC%, N%, adapter%, Q20% and Q30%. Human.B38 was chosen as the genome reference and alignment, and the latest Ensembl.R98 as the gene model. All the sample Fastq files were grouped and aligned using STAR  with the default setting. STAR was used in the quantification step to generate the raw count data at the gene level. Normalization was performed using Log Geometric Mean Method .
Principle component analysis (PCA)
PCA was applied to the sample set to identify relationships between different samples and find potential outliers. The first 3 components were used to find potential outliers. PCA was performed based on the normalized count data to detect potential outliers using software Array Studio (Qiagen, 2019). From the PCA, we considered one sample (4 h time-point, vehicle control) as an outlier (Additional file 2: Fig. S1). Therefore, all results described herein were analyzed without the outlier. No particular donor was considered an outlier. None-the-less, any donor variation was accounted for in the downstream statistical model.
Hierarchical clustering of the normalized count data was applied to identify grouping of samples based on their gene expression values, using Array Studio (Qiagen, 2019) and Euclidean distance metric as the variable distance method. Samples that were highly correlated were clustered together. The sample level grouping results were further compared to the true metadata to find hidden patterns between samples’ gene expression values and the corresponding metadata.
Differential expression analysis
DeSeq2 (equivalent to DeSeq V2 of the R package)  analysis was performed to identify Differentially Expressed Genes (DEGs) among different comparisons (e.g. WS-CM vs control, ACM vs control, WS-CM 7 µg/mL Eq-Nic vs. ACM 7 µg/mL Eq-Nic.). DEGs were defined as genes whose adjusted p-value was less than 0.05, and log2 fold change greater than 2 or less than − 2. Statistical significance of DEGs was computed using the Benjamini–Hochberg method  for multiple comparisons, and adjusted from p-values calculated using the Fisher’s Exact Test. A full statistical model to include all relevant factors (time, treatment and dosage) was built to find out the overall DEGs that are affected by treatments. The full statistical model is described as follows:
Ingenuity pathway analysis (IPA; QIAGEN Bioinformatics, Redwood City, CA, USA) was used to analyze the differential expression analysis results calculated from comparisons. IPA was also used to investigate differences between different dosages of the same treatment. The statistical significance of the association between genes in the dataset and canonical pathways was measured using the Benjamini–Hochberg method for multiple hypothesis correction , and adjusted p-value calculated using Fisher’s Exact Test).
Availability of data and materials
The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.
Aerosol conditioned medium
Chronic obstructive pulmonary disease
Cystic fibrosis transmembrane conductance regulator
Sodium epithelial channel
Electronic nicotine delivery system
Equivalent nicotine units
Harmful and potentially harmful constituents
Ingenuity pathway analysis
Next Generation Sequencing
Normal human bronchial epithelial
Trans-epithelial electrical resistance
Whole-smoke conditioned medium
Beane J, Vick J, Schembri F, Anderlind C, Gower A, Campbell J, Luo L, Zhang XH, Xiao J, Alekseyev YO. Characterizing the impact of smoking and lung cancer on the airway transcriptome using RNA-Seq. Cancer Prev Res. 2011;4(6):803–17.
Hübner R-H, Schwartz JD, De BP, Ferris B, Omberg L, Mezey JG, Hackett NR, Crystal RG. Coordinate control of expression of Nrf2-modulated genes in the human small airway epithelium is highly responsive to cigarette smoking. Mol Med. 2009;15(7–8):203–19.
Zhang X, Sebastiani P, Liu G, Schembri F, Zhang X, Dumas YM, Langer EM, Alekseyev Y, O’Connor GT, Brooks DR. Similarities and differences between smoking-related gene expression in nasal and bronchial epithelium. Physiol Gen. 2010;41(1):1–8.
Beineke P, Fitch K, Tao H, Elashoff MR, Rosenberg S, Kraus WE, Wingrove JA, Investigators P. A whole blood gene expression-based signature for smoking status. BMC Med Gen. 2012;5(1):58.
Martin F, Talikka M, Hoeng J, Peitsch M. Identification of gene expression signature for cigarette smoke exposure response—from man to mouse. Human Experiment Toxicol. 2015;34(12):1200–11.
Arimilli S, Schmidt E, Damratski BE, Prasad GL. Role of oxidative stress in the suppression of immune responses in peripheral blood mononuclear cells exposed to combustible tobacco product preparations. Inflammation. 2017;40(5):1622–30.
Arimilli S, Makena P, Liu G, Prasad G. Distinct gene expression changes in human peripheral blood mononuclear cells treated with different tobacco product preparations. Toxicol In Vitro. 2019;57:117–25.
Forey BA, Thornton AJ, Lee PN. Systematic review with meta-analysis of the epidemiological evidence relating smoking to COPD, chronic bronchitis and emphysema. BMC Pulm Med. 2011;11(1):36.
Islami F, Torre LA, Jemal A. Global trends of lung cancer mortality and smoking prevalence. Transl Lung Cancer Res. 2015;4(4):327–38.
Jayes L, Haslam PL, Gratziou CG, Powell P, Britton J, Vardavas C, Jimenez-Ruiz C, Leonardi-Bee J. SmokeHaz: systematic reviews and meta-analysis of the effects of smoking on respiratory health. Chest. 2016;150(1):164–79.
Silvestri M, Franchi S, Pistorio A, Petechia L, Rusconi F. Smoke exposure, wheezing, and asthma development: a systematic review and meta-analysis in unselected birth cohorts. Ped Pulm. 2015;50(4):353–62.
Torre LA, Siegel RL, Jemal A. Lung cancer statistics. Adv Exp Med Biol. 2016;893:1–19.
Rab A, Rowe SM, Raju SV, Bebock Z, Matalon S, Collawn JF. Cigarette smoke and CFTR: implications in the pathogenesis of COPD. AJP Lung Cell Mol Physiol. 2013;305(8):530–41.
Hajek P, Etter JF, Benowitz N, Eissenberg T, McRobbie H. Electronic cigarettes: review of use, content, safety, effects on smokers and potential for harm and benefit. Addiction. 2014;109(11):1801–10.
Harmful and potentially harmful consitituents (HPHCs) https://www.fda.gov/tobacco-products/products-ingredients-components/harmful-and-potentially-harmful-constituents-hphcs
Noël A, Hossain E, Perveen Z, Zaman H, Penn AL. Sub-ohm vaping increases the levels of carbonyls, is cytotoxic, and alters gene expression in human bronchial epithelial cells exposed at the air–liquid interface. Resp Res. 2020;21(1):1–20.
D’Ruiz CD, Graff DW, Robinson E. Reductions in biomarkers of exposure, impacts on smoking urge and assessment of product use and tolerability in adult smokers following partial or complete substitution of cigarettes with electronic cigarettes. BMC Public Health. 2016;16(1):543.
Song M-A, Freudenheim JL, Brasky TM, Mathe EA, McElroy JP, Nickerson QA, Reisinger SA, Smiraglia DJ, Weng DY, Ying KL. Biomarkers of exposure and effect in the lungs of smokers, nonsmokers, and electronic cigarette users. Cancer Epide Biomar. 2020;29(2):443–51.
De Jesús VR, Bhandari D, Zhang L, Reese C, Capella K, Tevis D, Zhu W, Valle-Pinero D, Arseima Y, Lagaud G. Urinary biomarkers of exposure to volatile organic compounds from the population assessment of tobacco and health study wave 1 (2013–2014). Int J Environ Res Pub Health. 2020;17(15):5408.
Corbett SE, Nitzberg M, Moses E, Kleerup E, Wang T, Perdomo C, Perdomo C, Liu G, Xiao X, Liu H. Gene expression alterations in the bronchial epithelium of e-cigarette users. Chest. 2019;156(4):764–73.
Tommasi S, Caliri AW, Caceres A, Moreno DE, Li M, Chen Y, Siegmund KD, Besaratinia A. Deregulation of biologically significant genes and associated molecular pathways in the oral epithelium of electronic cigarette users. Int J Mol Sci. 2019;20(3):738.
Rebuli ME, Glista-Baker E, Hoffman JR, Duffney PF, Robinette C, Speen AM, Pawlak EA, Dhingra R, Noah TL, Jaspers I. E-cigarette use alters nasal mucosal immune response to live-attenuated influenza virus (LAIV). Am J Resp Cell Mol Biol. 2021;64(1):126–37.
Shen Y, Wolkowicz MJ, Kotova T, Fan L, Timko MP. Transcriptome sequencing reveals e-cigarette vapor and mainstream-smoke from tobacco cigarettes activate different gene expression profiles in human bronchial epithelial cells. Sci Rep. 2016;6(4):23984.
Haswell LE, Baxter A, Banerjee A, Verrastro I, Mushonganono J, Adamson J, Thorne D, Gaça M, Minet E. Reduced biological effect of e-cigarette aerosol compared to cigarette smoke evaluated in vitro using normalized nicotine dose and RNA-seq-based toxicogenomics. Sci Rep. 2017;7(1):1–16.
Banerjee A, Haswell LE, Baxter A, Parmar A, Azzopardi D, Corke S, Thorne D, Adamson J, Mushonganono J, Gaca MD, et al. Differential gene expression using RNA sequencing profiling in a reconstituted airway epithelium exposed to conventional cigarette smoke or electronic cigarette aerosols. App In Vitro Toxicol. 2017;3(1):84–98.
Park SJ, Walser TC, Tran LM, Perdomo C, Wang T, Hong L, Pagano P, Li R, Jing Z, Liclican E, et al. The biological impacts of e-cigarettes on airway epithelial cell transformation and gene expression. J Thoracic Oncol. 2016;11(2):35–6.
Antherieu S, Garat A, Beauval N, Soyez M, Allorge D, Garcon G, Lo-Guidice J. Comparison of cellular and transcriptomic effects between electronic cigarette vapor and cigarette smoke in human bronchial epithelial cells. Toxicol In Vitro. 2017;45(3):417–25.
Martin EM, Clapp PW, Rebuli ME, Pawlak EA, Glista-Baker E, Benowitz NL, Fry RC, Jaspers I. E-cigarette use results in suppression of immune and inflammatory-response genes in nasal epithelial cells similar to cigarette smoke. AJP-Lung Cell Mol Physiol. 2016;311(1):L135-144.
Rayner RE, Makena P, Prasad GL, Cormet-Boyaka E. Cigarette and ENDS preparations differentially regulate ion channels and mucociliary clearance in primary normal human bronchial 3D cultures. AJP Lung Mol Cell Phys. 2019;317(2):295–302.
Rayner RE, Makena P, Prasad GL, Cormet-Boyaka E. Cigarette smoke preparations, not electronic nicotine delivery systems (ENDS) preparations, induce features of lung disease in a 3D lung repeat-dose model. AJP-Lung Cell Mol Physiol; 2020.
Rayner RE, Makena P, Prasad GL, Cormet-Boyaka E. Optimization of normal human bronchial epithelial (NHBE) cell 3D culture for in vitro lung model studies. Sci Rep. 2019;9:500.
Moses E, Wang T, Corbett S, Jackson GR, Drizik E, Perdomo C, Perdomo C, Kleerup E, Brooks D, O’Connor G, et al. Molecular impact of electronic cigarette aerosol exposure in human bronchial epithelium. Toxicol Sci. 2017;155(1):248–57.
Haswell LE, Corke S, Verrastro I, Baxter A, Banerjee A, Adamson J, Jaunky T, Proctor C, Gaça M, Minet E. In vitro RNA-seq-based toxicogenomics assessment shows reduced biological effect of tobacco heating products when compared to cigarette smoke. Sci Rep. 2018;8(1):1–18.
Cobb E, Hall J, Palazzolo DL. Induction of metallothionein expression after exposure to conventional cigarette smoke but not electronic cigarette (ECIG)-generated aerosol in caenorhabditis elegans. Front Physiol. 2018;9:426.
Lin J, Peng J, Liu G, Deng L. Overexpression of MECP2 attenuates cigarette smoke extracts induced lung epithelial cell injury by promoting CYP1B1 methylation. J Toxicol Sci. 2020;45(3):177–86.
Iu M, Zago M, de Souza AR, Bouttier M, Pareek S, White JH, Hamid Q, Eidelman DH, Baglole CJ. RelB attenuates cigarette smoke extract-induced apoptosis in association with transcriptional regulation of the aryl hydrocarbon receptor. Free Radical Biol Med. 2017;108:19–31.
Marklew AJ, Patel W, Moore PJ, Tan CD, Smith AJ, Sassano MF, Gray MA, Tarran R. Cigarette smoke exposure induces retrograde trafficking of CFTR to the endoplasmic reticulum. Sci Rep. 2019;9(1):1–17.
Hassan F, Rennolds J, Nuovo G, Ezzie ME, Diaz P, Cormet-Boyaka E. Cigarette smoke regulates the expression of the CFTR chloride channel in bronchial epithelial cells. Am J Resp Crit Care. 2011;183:A4226.
Bodas M, Moore AR, Subramaniyan B, Georgescu C, Wren JD, Freeman WM, Brown BR, Metcalf JP, Walters MS. Cigarette smoke activates NOTCH3 to promote goblet cell differentiation in human airway epithelial cells. Am J Respir Cell Mol Biol. 2021;64(4):426–40.
Hussain SS, George S, Singh S, Jayant R, Hu C-A, Sopori M, Chand HS. A small molecule BH3-mimetic suppresses cigarette smoke-induced mucous expression in airway epithelial cells. Sci Rep. 2018;8(1):1–11.
Solleti SK, Bhattacharya S, Ahmad A, Wang Q, Mereness J, Rangasamy T, Mariani TJ. MicroRNA expression profiling defines the impact of electronic cigarettes on human airway epithelial cells. Sci Rep. 2017;7(1):1–10.
Iskandar AR, Zanetti F, Marescotti D, Titz B, Sewer A, Kondylis A, Leroy P, Belcastro V, Torres LO, Acali S, et al. Application of a multi-layer systems toxicology framework for in vitro assessment of the biological effects of Classic Tobacco e-liquid and its corresponding aerosol using an e-cigarette device with MESH technology. Arch Toxicol. 2019;93:3229–47.
Giralt A, Iskandar AR, Martin F, Moschini E, Serchi T, Kondylis A, Marescotti D, Leroy P, Ortega-Torres L, Majeed S. Comparison of the biological impact of aerosol of e-vapor device with MESH® technology and cigarette smoke on human bronchial and alveolar cultures. Toxicol Lett. 2021;337:98–110.
Spira A, Beane J, Shah V, Liu G, Schembri F, Yang X, Palma J, Brody JS. Effects of cigarette smoke on the human airway epithelial cell transcriptome. PNAS. 2004;101(27):10143–8.
Song M-A, Reisinger SA, Freudenheim JL, Brasky TM, Mathé EA, McElroy JP, Nickerson QA, Weng DY, Wewers MD, Shields PG. Effects of electronic cigarette constituents on the human lung: a pilot clinical trial. Cancer Prev Res. 2020;13(2):145–52.
Clapp PW, Lavrich KS, van Heusden CA, Lazarowski ER, Carson JL, Jaspers I. Cinnamaldehyde in flavored e-cigarette liquids temporarily suppresses bronchial epithelial cell ciliary motility by dysregulation of mitochondrial function. Am J Physiol Lung Cell Mol Physiol. 2019;316(3):L470–86.
Clapp PW, Pawlak EA, Lackey JT, Keating JE, Reeber SL, Glish GL, Jaspers I. Flavored e-cigarette liquids and cinnamaldehyde impair respiratory innate immune cell function. Am J Physiol Lung Cell Mol Physiol 2017.
Fernandez Fernandez E, De Santi C, De Rose V, Greene CM. CFTR dysfunction in cystic fibrosis and chronic obstructive pulmonary disease. Expert Rev Resp Med. 2018;12(6):483–92.
Raju SV, Jackson PL, Courville CA, McNicholas CM, Sloane PA, Sabbatini G, Tidwell S, Tang LP, Liu B, Fortenberry JA. Cigarette smoke induces systemic defects in cystic fibrosis transmembrane conductance regulator function. Am J Resp Crit Care Med. 2013;188(11):1321–30.
Chen G, Korfhagen TR, Xu Y, Kitzmiller J, Wert SE, Maeda Y, Gregorieff A, Clevers H, Whitsett JA. SPDEF is required for mouse pulmonary goblet cell differentiation and regulates a network of genes associated with mucus production. J Clin Invest. 2009;119(10):2914–24.
Chen G, Volmer AS, Wilkinson KJ, Deng Y, Jones LC, Yu D, Bustamante-Marin XM, Burns KA, Grubb BR, O’Neal WK. Role of Spdef in the regulation of Muc5b expression in the airways of naive and mucoobstructed mice. Am J Resp Cell Mol Biol. 2018;59(3):383–96.
Chen G, Korfhagen TR, Karp CL, Impey S, Xu Y, Randell SH, Kitzmiller J, Maeda Y, Haitchi HM, Sridharan A. Foxa3 induces goblet cell metaplasia and inhibits innate antiviral immunity. Am J Resp Crit Care Med. 2014;189(3):301–13.
Hikichi M, Mizumura K, Maruoka S, Gon Y. Pathogenesis of chronic obstructive pulmonary disease (COPD) induced by cigarette smoke. J Thoracic Dis. 2019;11(Suppl 17):S2129.
Parris BA, O’Farrell HE, Fong KM, Yang IA. Chronic obstructive pulmonary disease (COPD) and lung cancer: common pathways for pathogenesis. J Thoracic Dis. 2019;11(Suppl 17):S2155.
Boukhenouna S, Wilson MA, Bahmed K, Kosmider B. Reactive oxygen species in chronic obstructive pulmonary disease. Oxidative Med Cell Longevity 2018, 2018.
Krewski D, Bird M, Al-Zoughool M, Birkett N, Billard M, Milton B, Rice JM, Grosse Y, Cogliano VJ, Hill MA. Key characteristics of 86 agents known to cause cancer in humans. J Toxicol Environ Health B. 2019;22(7–8):244–63.
Billatos E, Faiz A, Gesthalter Y, LeClerc A, Alekseyev Y, Xiao X, Liu G, ten Hacken NH, Heijink I, Timens W. Impact of acute exposure to cigarette smoke on airway gene expression. Physiol Genomics. 2018;50(9):705–13.
Arimilli S, Damratoski BE, Bombick B, Borgerding MF, Prasad GL. Evaluation of cytotoxicity of different tobacco product preparations. Reg Toxicol Pharmacol. 2012;64(3):350–60.
Johnson MD. Evaluation of in vitro assays for assessing the toxicity of cigarette smoke and smokeless tobacco. Cancer Epidemiol Biomarkers Prev. 2009;18(12):3263–304.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M, Gingeras TR. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Gen Biol. 2014;15(12):550.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc. 1995;57(1):289–300.
We thank the Genomics Shared Resources at The Ohio State University Comprehensive Cancer Center, Columbus, OH for the Next-Generation sequencing which was supported by NIH Grant P30 CA016058. Human bronchial epithelial cultures for this work were supplied by the Cure CF Columbus (C3) Epithelial Cell Core at Nationwide Children’s Hospital. C3 is supported by a Research Development Program Grant (MCCOY19R0) from the Cystic Fibrosis Foundation. We thank Neha Jalan from Qiagen Genomic Services (Germantown, MD) for helping with data analysis.
This work was funded by RAI Services Company under a collaborative research agreement with The Ohio State University.
Ethics approval and consent to participate
An exempt status was granted by the OSU Institutional Review Board for the use of primary airway cells since these were provided without identifiers.
Consent for publication
Competing financial interests: P.M. and G.L.P. are employees of RAI Services Company. RAI Services Company is a wholly-owned subsidiary of Reynolds American, Inc., which is a wholly-owned subsidiary of British American Tobacco plc. R.E.R., G.L., and E.C-B. have no competing financial interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1.
showing the genes listed in Figure 1: “Significant DEGs were identified with low, medium and high dose of WS-CM at 4hrs (Tab 1) and 24hrs (Tab 2) exposure. Genes highlighted yellow were consistently observed in all three doses. Genes highlighted in orange were observed in both medium and high doses, but not low dose. Log2 fold change cut-off is set to ±2 and p-value ≤ 0.05. A list of genes identified across both exposure times (4hrs and 24hrs), and within each dose (low, medium and high) is provided in Tab 3.”
Additional file 2.
Supplementary Tables 1–5 and Supplementary Figures 1 and 2.
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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Rayner, R.E., Makena, P., Liu, G. et al. Differential gene expression of 3D primary human airway cultures exposed to cigarette smoke and electronic nicotine delivery system (ENDS) preparations. BMC Med Genomics 15, 76 (2022). https://doi.org/10.1186/s12920-022-01215-x
- Primary Normal Human Bronchial Epithelial (NHBE) cells
- Electronic Nicotine Delivery Systems (ENDS)