Serial Sampling of the Small Airway Epithelium to Identify Persistent Smoking-dysregulated Genes.

Rationale: The small airway epithelium (beyond the sixth generation), the initiation site of smoking-induced airway disorders, is highly sensitive to the stress of smoking. Because of variations over time in smoking habits, the small airway epithelium transcriptome is dynamic, fluctuating not only among smokers but also within each smoker. Objectives: To perform accurate assessment of the smoking-related dysregulation of the human small airway epithelium despite the variation of smoking within the same individual and of the effects of smoking cessation on the dysregulated transcriptome. Methods: We conducted serial sampling of the same smokers and nonsmoker control subjects over time to identify persistent smoking dysregulation of the biology of the small airway epithelium over 1 year. We conducted serial sampling of smokers who quit smoking, before and after smoking cessation, to assess the effect of smoking cessation on the smoking-dysregulated genes. Measurements and Main Results: Repeated measures ANOVA of the small airway epithelium transcriptome sampled four times in the same individuals over 1 year enabled the identification of 475 persistent smoking-dysregulated genes. Most genes were normalized after 12 months of smoking cessation; however, 53 (11%) genes, including CYP1B1, PIR, ME1, and TRIM16, remained persistently abnormally expressed. Dysregulated pathways enriched with the nonreversible genes included xenobiotic metabolism signaling, bupropion degradation, and nicotine degradation. Conclusions: Analysis of repetitive sampling of the same individuals identified persistent smoking-induced dysregulation of the small airway epithelium transcriptome and the effect of smoking cessation. These results help identify targets for the development of therapies that can be applicable to smoking-related airway diseases.

High apoptotic endothelial microparticle levels measured in asthma with elevated IgE and eosinophils.

While asthma is considered an inflammatory-mediated airway epithelial and smooth muscle disorder, there is increasing evidence of airway capillary endothelial dysfunction associated with vascular remodelling and angiogenesis in some individuals with this condition. The inflammation is typically characterized as type-2 high (eosinophilic) vs type 2-low (neutrophilic and pauci-granulocytic); we hypothesized that the type-2 high group would be more likely to evidence endothelial dysfunction. As a biomarker of these processes, we hypothesized that nonsmokers with allergic asthma may have elevated plasma levels of endothelial microparticles (EMPs), membrane vesicles that are shed when endothelial cells undergo activation or apoptosis. Total and apoptotic circulating EMPs were measured by fluorescence-activated cell analysis in patients with allergic asthma (n = 29) and control subjects (n = 26), all nonsmokers. When the entire group of patients with asthma were compared to the control subjects, there were no differences in total circulating EMPs nor apoptotic EMPs. However, patients with asthma with elevated levels of IgE and eosinophils had higher levels of apoptotic EMPs, compared to patients with asthma with mildly increased IgE and eosinophil levels. This observation is relevant to precision therapies for asthma and highlights the importance of sub-phenotyping in the condition.

Oral epithelium response of electronic cigarette users to electronic cigarette.

BACKGROUND: Electronic cigarettes are increasing in popularity, but there is only little information on their biologic effects on the oral epithelium, the initial site exposed to electronic cigarette smoke. METHODS: We assessed the oral epithelium response to electronic cigarettes by comparing the histology and RNA transcriptome (mRNA and miRNA) of healthy electronic cigarette vapers to nonsmokers. mRNA was assessed based on: (1) genome-wide; (2) genes previously identified as dysregulated in the oral epithelium of electronic cigarette vapers versus nonsmokers; (3) immune and inflammatory-related genes previously identified as dysregulated in the nasal epithelium of electronic cigarette vapers compared to nonsmokers; (4) genes previously identified as dysregulated in the small airway epithelium of nonsmokers following an acute exposure to electronic cigarette; and (5) genes related to the initial steps of COVID-19 infection. In addition, miRNA was assessed genome-wide. Comparisons were performed using analysis of variance, and Benajmini-Hochberg corrected p < 0.05 was considered significant. RESULTS: The histology of the epithelium, lamina propria and basal layer in electronic cigarette vapers appeared normal. Assessment of mRNA and miRNA, based on all gene lists, did not identify any genes significantly modified in the oral epithelium of electronic cigarette vapers in response to electronic cigarette use. CONCLUSION: An average history of 2 years of vaping results in no detectable histologic or transcriptome abnormalities in the buccal mucosa.

Up-regulation of ACE2, the SARS-CoV-2 receptor, in asthmatics on maintenance inhaled corticosteroids.

BACKGROUND: The first step in SARS-CoV-2 infection is binding of the virus to angiotensin converting enzyme 2 (ACE2) on the airway epithelium. Asthma affects over 300 million people world-wide, many of whom may encounter SARS-CoV-2. Epidemiologic data suggests that asthmatics who get infected may be at increased risk of more severe disease. Our objective was to assess whether maintenance inhaled corticosteroids (ICS), a major treatment for asthma, is associated with airway ACE2 expression in asthmatics. METHODS: Large airway epithelium (LAE) of asthmatics treated with maintenance ICS (ICS+), asthmatics not treated with ICS (ICS-), and healthy controls (controls) was analyzed for expression of ACE2 and other coronavirus infection-related genes using microarrays. RESULTS: As a group, there was no difference in LAE ACE2 expression in all asthmatics vs controls. In contrast, subgroup analysis demonstrated that LAE ACE2 expression was higher in asthmatics ICS+ compared to ICS‾ and ACE2 expression was higher in male ICS+ compared to female ICS+ and ICS‾ of either sex. ACE2 expression did not correlate with serum IgE, absolute eosinophil level, or change in FEV1 in response to bronchodilators in either ICS- or ICS+. CONCLUSION: Airway ACE2 expression is increased in asthmatics on long-term treatment with ICS, an observation that should be taken into consideration when assessing the use of inhaled corticosteroids during the pandemic.

Expression of the SARS-CoV-2 ACE2 Receptor in the Human Airway Epithelium.

Rationale: Infection with the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes coronavirus disease (COVID-19), a predominantly respiratory illness. The first step in SARS-CoV-2 infection is binding of the virus to ACE2 (angiotensin-converting enzyme 2) on the airway epithelium.Objectives: The objective was to gain insight into the expression of ACE2 in the human airway epithelium.Methods: Airway epithelia sampled by fiberoptic bronchoscopy of trachea, large airway epithelia (LAE), and small airway epithelia (SAE) of nonsmokers and smokers were analyzed for expression of ACE2 and other coronavirus infection-related genes using microarray, RNA sequencing, and 10x single-cell transcriptome analysis, with associated examination of ACE2-related microRNA.Measurements and Main Results:1) ACE2 is expressed similarly in the trachea and LAE, with lower expression in the SAE; 2) in the SAE, ACE2 is expressed in basal, intermediate, club, mucus, and ciliated cells; 3) ACE2 is upregulated in the SAE by smoking, significantly in men; 4) levels of miR-1246 expression could play a role in ACE2 upregulation in the SAE of smokers; and 5) ACE2 is expressed in airway epithelium differentiated in vitro on air-liquid interface cultures from primary airway basal stem/progenitor cells; this can be replicated using LAE and SAE immortalized basal cell lines derived from healthy nonsmokers.Conclusions:ACE2, the gene encoding the receptor for SARS-CoV-2, is expressed in the human airway epithelium, with variations in expression relevant to the biology of initial steps in SARS-CoV-2 infection.

Cell-specific expression of lung disease risk-related genes in the human small airway epithelium.

BACKGROUND: The human small airway epithelium (SAE) plays a central role in the early events in the pathogenesis of most inherited and acquired lung disorders. Little is known about the molecular phenotypes of the specific cell populations comprising the SAE in humans, and the contribution of SAE specific cell populations to the risk for lung diseases. METHODS: Drop-seq single-cell RNA-sequencing was used to characterize the transcriptome of single cells from human SAE of nonsmokers and smokers by bronchoscopic brushing. RESULTS: Eleven distinct cell populations were identified, including major and rare epithelial cells, and immune/inflammatory cells. There was cell type-specific expression of genes relevant to the risk of the inherited pulmonary disorders, genes associated with risk of chronic obstructive pulmonary disease and idiopathic pulmonary fibrosis and (non-mutated) driver genes for lung cancers. Cigarette smoking significantly altered the cell type-specific transcriptomes and disease risk-related genes. CONCLUSIONS: This data provides new insights into the possible contribution of specific lung cells to the pathogenesis of lung disorders.

Exaggerated BMP4 signalling alters human airway basal progenitor cell differentiation to cigarette smoking-related phenotypes.

Airway remodelling in chronic obstructive pulmonary disease (COPD) originates, in part, from smoking-induced changes in airway basal stem/progenitor cells (BCs). Based on the knowledge that bone morphogenetic protein 4 (BMP4) influences epithelial progenitor function in the developing and adult mouse lung, we hypothesised that BMP4 signalling may regulate the biology of adult human airway BCs relevant to COPD.BMP4 signalling components in human airway epithelium were analysed at the mRNA and protein levels, and the differentiation of BCs was assessed using the BC expansion and air-liquid interface models in the absence/presence of BMP4, BMP receptor inhibitor and/or small interfering RNAs against BMP receptors and downstream signalling.The data demonstrate that in cigarette smokers, BMP4 is upregulated in ciliated and intermediate undifferentiated cells, and expression of the BMP4 receptor BMPR1A is enriched in BCs. BMP4 induced BCs to acquire a smoking-related abnormal phenotype in vitro mediated by BMPR1A/Smad signalling, characterised by decreased capacity to differentiate into normal mucociliary epithelium, while generating squamous metaplasia.Exaggerated BMP4 signalling promotes cigarette smoking-relevant airway epithelial remodelling by inducing abnormal phenotypes in human airway BCs. Targeting of BMP4 signalling in airway BCs may represent a novel target to prevent/treat COPD-associated airway disease.

Ambient Pollution-related Reprogramming of the Human Small Airway Epithelial Transcriptome.

RATIONALE: Epidemiologic studies have demonstrated that exposure to particulate matter ambient pollution has adverse effects on lung health, exacerbated by cigarette smoking. Particulate matter less than or equal to 2.5 µm in aerodynamic diameter (PM2.5) is among the most harmful urban pollutants and is closely linked to respiratory disease. OBJECTIVES: Based on the knowledge that the small airway epithelium (SAE) plays a central role in the pathogenesis of smoking-related lung disease, we hypothesized that elevated PM2.5 levels are associated with dysregulation of SAE gene expression, which may contribute to the development of respiratory disease. METHODS: From 2009 to 2012, healthy nonsmoker (n = 29) and smoker (n = 129) residents of New York City underwent bronchoscopy with SAE brushing (2.6 ± 1.3 samples/subject; total of 405 samples). SAE gene expression was assessed by Affymetrix HG-U133 Plus 2.0 microarray. New York City PM2.5 levels (Environmental Protection Agency data) were averaged for the 30 days before bronchoscopy. A linear mixed model was used to assess PM2.5-related gene dysregulation accounting for multiple clinical and methodologic variables. MEASUREMENTS AND MAIN RESULTS: Thirty-day mean PM2.5 levels varied from 6.2 to 18 µg/m3. In nonsmokers, there was no dysregulation of SAE gene expression associated with ambient PM2.5 levels. In marked contrast, n = 219 genes were significantly dysregulated in association with PM2.5 levels in the SAE of smokers. Many of these genes relate to cell growth and transcription regulation. Interestingly, 11% of genes were mitochondria associated. CONCLUSIONS: PM2.5 exposure contributes to significant dysregulation of the SAE transcriptome of smokers, linking pollution and airway epithelial biology in the risk of development of respiratory disease in susceptible individuals.

Ontogeny and Biology of Human Small Airway Epithelial Club Cells.

RATIONALE: Little is known about human club cells, dome-shaped cells with dense cytoplasmic granules and microvilli that represent the major secretory cells of the human small airways (at least sixth-generation bronchi). OBJECTIVES: To define the ontogeny and biology of the human small airway epithelium club cell. METHODS: The small airway epithelium was sampled from the normal human lung by bronchoscopy and brushing. Single-cell transcriptome analysis and air-liquid interface culture were used to assess club cell ontogeny and biology. MEASUREMENTS AND MAIN RESULTS: We identified the club cell population by unbiased clustering using single-cell transcriptome sequencing. Principal component gradient analysis uncovered an ontologic link between KRT5 (keratin 5)+ basal cells and SCGB1A1 (secretoglobin family 1A member 1)+ club cells, a hypothesis verified by demonstrating in vitro that a pure population of human KRT5+ SCGB1A1- small airway epithelial basal cells differentiate into SCGB1A1+KRT5- club cells on air-liquid interface culture. Using SCGB1A1 as the marker of club cells, the single-cell analysis identified novel roles for these cells in host defense, xenobiotic metabolism, antiprotease, physical barrier function, monogenic lung disorders, and receptors for human viruses. CONCLUSIONS: These observations provide novel insights into the molecular phenotype and biologic functions of the human club cell population and identify basal cells as the human progenitor cells for club cells.

POU2AF1 Functions in the Human Airway Epithelium To Regulate Expression of Host Defense Genes.

In the process of seeking novel lung host defense regulators by analyzing genome-wide RNA sequence data from normal human airway epithelium, we detected expression of POU domain class 2-associating factor 1 (POU2AF1), a known transcription cofactor previously thought to be expressed only in lymphocytes. Lymphocyte contamination of human airway epithelial samples obtained by bronchoscopy and brushing was excluded by immunohistochemistry staining, the observation of upregulation of POU2AF1 in purified airway basal stem/progenitor cells undergoing differentiation, and analysis of differentiating single basal cell clones. Lentivirus-mediated upregulation of POU2AF1 in airway basal cells induced upregulation of host defense genes, including MX1, IFIT3, IFITM, and known POU2AF1 downstream genes HLA-DRA, ID2, ID3, IL6, and BCL6. Interestingly, expression of these genes paralleled changes of POU2AF1 expression during airway epithelium differentiation in vitro, suggesting POU2AF1 helps to maintain a host defense tone even in pathogen-free condition. Cigarette smoke, a known risk factor for airway infection, suppressed POU2AF1 expression both in vivo in humans and in vitro in human airway epithelial cultures, accompanied by deregulation of POU2AF1 downstream genes. Finally, enhancing POU2AF1 expression in human airway epithelium attenuated the suppression of host defense genes by smoking. Together, these findings suggest a novel function of POU2AF1 as a potential regulator of host defense genes in the human airway epithelium.

Smoking-Dependent Distal-to-Proximal Repatterning of the Adult Human Small Airway Epithelium.

RATIONALE: Small airways are the primary site of pathologic changes in chronic obstructive pulmonary disease (COPD), the major smoking-induced lung disorder. OBJECTIVES: On the basis of the concept of proximal-distal patterning that determines regional specialization of the airway epithelium during lung development, we hypothesized that a similar program operates in the adult human lung being altered by smoking, leading to decreased regional identity of the small airway epithelium (SAE). METHODS: The proximal and distal airway signatures were identified by comparing the transcriptomes of large and small airway epithelium samples obtained by bronchoscopy from healthy nonsmokers. The expression of these signatures was evaluated in the SAE of healthy smokers and smokers with COPD compared with that of healthy nonsmokers. The capacity of airway basal stem cells (BCs) to maintain region-associated phenotypes was evaluated using the air-liquid interface model. MEASUREMENTS AND MAIN RESULTS: The distal and proximal airway signatures, containing 134 and 233 genes, respectively, were identified. These signatures included known developmental regulators of airway patterning, as well as novel regulators such as epidermal growth factor receptor, which was associated with the proximal airway phenotype. In the SAE of smokers with COPD, there was a dramatic smoking-dependent loss of the regional transcriptome identity with concomitant proximalization. This repatterning phenotype was reproduced by stimulating SAE BCs with epidermal growth factor, which was up-regulated in the SAE of smokers, during differentiation of SAE BCs in vitro. CONCLUSIONS: Smoking-induced global distal-to-proximal reprogramming of the SAE represents a novel pathologic feature of COPD and is mediated by exaggerated epidermal growth factor/epidermal growth factor receptor signaling in SAE BCs.

Pulmonary Abnormalities in Young, Light-Use Waterpipe (Hookah) Smokers.

RATIONALE: Waterpipes, also called hookahs, are currently used by millions of people worldwide. Despite the increasing use of waterpipe smoking, there is limited data on the health effects of waterpipe smoking and there are no federal regulations regarding its use. OBJECTIVES: To assess the effects of waterpipe smoking on the human lung using clinical and biological parameters in young, light-use waterpipe smokers. METHODS: We assessed young, light-use, waterpipe-only smokers in comparison with lifelong nonsmokers using clinical parameters of cough and sputum scores, lung function, and chest high-resolution computed tomography as well as biological parameters of lung epithelial lining fluid metabolome, small airway epithelial (SAE) cell differential and transcriptome, alveolar macrophage transcriptome, and plasma apoptotic endothelial cell microparticles. MEASUREMENTS AND MAIN RESULTS: Compared with nonsmokers, waterpipe smokers had more cough and sputum as well as a lower lung diffusing capacity, abnormal epithelial lining fluid metabolome profile, increased proportions of SAE secretory and intermediate cells, reduced proportions of SAE ciliated and basal cells, markedly abnormal SAE and alveolar macrophage transcriptomes, and elevated levels of apoptotic endothelial cell microparticles. CONCLUSIONS: Young, light-use, waterpipe-only smokers have a variety of abnormalities in multiple lung-related biological and clinical parameters, suggesting that even limited waterpipe use has broad consequences on human lung biology and health. We suggest that large epidemiological studies should be initiated to investigate the harmful effects of waterpipe smoking.

Smoking-Associated Disordering of the Airway Basal Stem/Progenitor Cell Metabotype.

The airway epithelium is a complex pseudostratified multicellular layer lining the tracheobronchial tree, functioning as the primary defense against inhaled environmental contaminants. The major cell types of the airway epithelium include basal, intermediate columnar, ciliated, and secretory. Basal cells (BCs) are the proliferating stem/progenitor population that differentiate into the other specialized cell types of the airway epithelium during normal turnover and repair. Given that cigarette smoke delivers thousands of xenobiotics and high levels of reactive molecules to the lung epithelial surface, we hypothesized that cigarette smoke broadly perturbs BC metabolism. To test this hypothesis, primary airway BCs were isolated from healthy nonsmokers (n = 11) and healthy smokers (n = 7) and assessed by global metabolic profiling by liquid chromatography-mass spectrometry. The analysis identified 52 significant metabolites in BCs differentially expressed between smokers and nonsmokers (P < 0.05). These changes included metabolites associated with redox pathways, energy production, and inflammatory processes. Notably, BCs from smokers exhibited altered levels of the key enzyme cofactors/substrates nicotinamide adenine dinucleotide, flavin adenine dinucleotide, acetyl coenzyme A, and membrane phospholipid levels. Consistent with the high burden of oxidants in cigarette smoke, glutathione levels were diminished, whereas 3-nitrotyrosine levels were increased, suggesting that protection of airway epithelial cells against oxidative and nitrosative stress is significantly compromised in smoker BCs. It is likely that this altered metabotype is a reflection of, and likely contributes to, the disordered biology of airway BCs consequent to the stress cigarette smoking puts on the airway epithelium.

Persistence of circulating endothelial microparticles in COPD despite smoking cessation.

INTRODUCTION: Increasing evidence links COPD pathogenesis with pulmonary capillary apoptosis. We previously demonstrated that plasma levels of circulating microparticles released from endothelial cells (EMPs) due to apoptosis are elevated in smokers with normal spirometry but low diffusion capacity, that is, with early evidence of lung destruction. We hypothesised that pulmonary capillary apoptosis persists with the development of COPD and assessed its reversibility in healthy smokers and COPD smokers following smoking cessation. METHODS: Pulmonary function and high-resolution CT (HRCT) were assessed in 28 non-smokers, 61 healthy smokers and 49 COPD smokers; 17 healthy smokers and 18 COPD smokers quit smoking for 12 months following the baseline visit. Total EMP (CD42b-CD31+), pulmonary capillary EMP (CD42b-CD31+ACE+) and apoptotic EMP (CD42b-CD62E+/CD42b-CD31+) levels were quantified by flow cytometry. RESULTS: Compared with non-smokers, healthy smokers and COPD smokers had elevated levels of circulating EMPs due to active pulmonary capillary endothelial apoptosis. Levels remained elevated over 12 months in healthy smokers and COPD smokers who continued smoking, but returned to non-smoker levels in healthy smokers who quit. In contrast, levels remained significantly abnormal in COPD smokers who quit. CONCLUSIONS: Pulmonary capillary apoptosis is reversible in healthy smokers who quit, but continues to play a role in COPD pathogenesis in smokers who progressed to airflow obstruction despite smoking cessation. TRIAL REGISTRATION NUMBER: NCT00974064; NCT01776398.

Cigarette smoking induces small airway epithelial epigenetic changes with corresponding modulation of gene expression.

The small airway epithelium (SAE), the first site of smoking-induced lung pathology, exhibits genome-wide changes in gene expression in response to cigarette smoking. Based on the increasing evidence that the epigenome can respond to external stimuli in a rapid manner, we assessed the SAE of smokers for genome-wide DNA methylation changes compared with nonsmokers, and whether changes in SAE DNA methylation were linked to the transcriptional output of these cells. Using genome-wide methylation analysis of SAE DNA of nonsmokers and smokers, the data identified 204 unique genes differentially methylated in SAE DNA of smokers compared with nonsmokers, with 67% of the regions with differential methylation occurring within 2 kb of the transcriptional start site. Among the genes with differential methylation were those related to metabolism, transcription, signal transduction and transport. For the differentially methylated genes, 35 exhibited a correlation with gene expression, 54% with an inverse correlation of DNA methylation with gene expression and 46% a direct correlation. These observations provide evidence that cigarette smoking alters the DNA methylation patterning of the SAE and that, for some genes, these changes are associated with the smoking-related changes in gene expression.

Risk of COPD with obstruction in active smokers with normal spirometry and reduced diffusion capacity.

Smokers are assessed for chronic obstructive pulmonary disease (COPD) using spirometry, with COPD defined by the Global Initiative for Chronic Obstructive Lung Disease (GOLD) as airflow limitation that is not fully reversible with bronchodilators. There is a subset of smokers with normal spirometry (by GOLD criteria), who have a low diffusing capacity of the lung for carbon monoxide (DLCO), a parameter linked to emphysema and small airway disease. The natural history of these "normal spirometry/low DLCO" smokers is unknown.From a cohort of 1570 smokers in the New York City metropolitian area, all of whom had normal spirometry, two groups were randomly selected for lung function follow-up: smokers with normal spirometry/normal DLCO (n=59) and smokers with normal spirometry/low DLCO (n=46). All had normal history, physical examination, complete blood count, urinalysis, HIV status, α1-antitrypsin level, chest radiography, forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC), FEV1/FVC ratio and total lung capacity. Throughout the study, all continued to be active smokers.In the normal spirometry/normal DLCO group assessed over 45±20 months, 3% developed GOLD-defined COPD. In contrast, in the normal spirometry/low DLCO group, followed over 41±31 months, 22% developed GOLD-defined COPD.Despite appearing "normal" according to GOLD, smokers with normal spirometry but low DLCO are at significant risk of developing COPD with obstruction to airflow.

HEFT: eQTL analysis of many thousands of expressed genes while simultaneously controlling for hidden factors.

MOTIVATION: Identification of expression Quantitative Trait Loci (eQTL), the genetic loci that contribute to heritable variation in gene expression, can be obstructed by factors that produce variation in expression profiles if these factors are unmeasured or hidden from direct analysis. METHODS: We have developed a method for Hidden Expression Factor analysis (HEFT) that identifies individual and pleiotropic effects of eQTL in the presence of hidden factors. The HEFT model is a combined multivariate regression and factor analysis, where the complete likelihood of the model is used to derive a ridge estimator for simultaneous factor learning and detection of eQTL. HEFT requires no pre-estimation of hidden factor effects; it provides P-values and is extremely fast, requiring just a few hours to complete an eQTL analysis of thousands of expression variables when analyzing hundreds of thousands of single nucleotide polymorphisms on a standard 8 core 2.6 G desktop. RESULTS: By analyzing simulated data, we demonstrate that HEFT can correct for an unknown number of hidden factors and significantly outperforms all related hidden factor methods for eQTL analysis when there are eQTL with univariate and multivariate (pleiotropic) effects. To demonstrate a real-world application, we applied HEFT to identify eQTL affecting gene expression in the human lung for a study that included presumptive hidden factors. HEFT identified all of the cis-eQTL found by other hidden factor methods and 91 additional cis-eQTL. HEFT also identified a number of eQTLs with direct relevance to lung disease that could not be found without a hidden factor analysis, including cis-eQTL for GTF2H1 and MTRR, genes that have been independently associated with lung cancer. AVAILABILITY: Software is available at http://mezeylab.cb.bscb.cornell.edu/Software.aspx. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Airway basal cells of healthy smokers express an embryonic stem cell signature relevant to lung cancer.

Activation of the human embryonic stem cell (hESC) signature genes has been observed in various epithelial cancers. In this study, we found that the hESC signature is selectively induced in the airway basal stem/progenitor cell population of healthy smokers (BC-S), with a pattern similar to that activated in all major types of human lung cancer. We further identified a subset of 6 BC-S hESC genes, whose coherent overexpression in lung adenocarcinoma (AdCa) was associated with reduced lung function, poorer differentiation grade, more advanced tumor stage, remarkably shorter survival, and higher frequency of TP53 mutations. BC-S shared with hESC and a considerable subset of lung carcinomas a common TP53 inactivation molecular pattern which strongly correlated with the BC-S hESC gene expression. These data provide transcriptome-based evidence that smoking-induced reprogramming of airway BC toward the hESC-like phenotype might represent a common early molecular event in the development of aggressive lung carcinomas in humans.