Effect of active smoking on the human bronchial epithelium transcriptome

Integrative genomics identifies SHPRH as a tumor suppressor gene in lung adenocarcinoma that regulates DNA damage response

BACKGROUND

Identification of driver mutations and development of targeted therapies has considerably improved outcomes for lung cancer patients. However, significant limitations remain with the lack of identified drivers in a large subset of patients. Here, we aimed to assess the genomic landscape of lung adenocarcinomas (LUADs) from individuals without a history of tobacco use to reveal new genetic drivers of lung cancer.

METHODS

Integrative genomic analyses combining whole-exome sequencing, copy number, and mutational information for 83 LUAD tumors was performed and validated using external datasets to identify genetic variants with a predicted functional consequence and assess association with clinical outcomes. LUAD cell lines with alteration of identified candidates were used to functionally characterize tumor suppressive potential using a conditional expression system both in vitro and in vivo.

RESULTS

We identified 21 genes with evidence of positive selection, including 12 novel candidates that have yet to be characterized in LUAD. In particular, SNF2 Histone Linker PHD RING Helicase (SHPRH) was identified due to its frequency of biallelic disruption and location within the familial susceptibility locus on chromosome arm 6q. We found that low SHPRH mRNA expression is associated with poor survival outcomes in LUAD patients. Furthermore, we showed that re-expression of SHPRH in LUAD cell lines with inactivating alterations for SHPRH reduces their in vitro colony formation and tumor burden in vivo. Finally, we explored the biological pathways associated SHPRH inactivation and found an association with the tolerance of LUAD cells to DNA damage.

CONCLUSIONS

These data suggest that SHPRH is a tumor suppressor gene in LUAD, whereby its expression is associated with more favorable patient outcomes, reduced tumor and mutational burden, and may serve as a predictor of response to DNA damage. Thus, further exploration into the role of SHPRH in LUAD development may make it a valuable biomarker for predicting LUAD risk and prognosis.

Smoking modulates different secretory subpopulations expressing SARS-CoV-2 entry genes in the nasal and bronchial airways

SARS-CoV-2 infection and disease severity are influenced by viral entry (VE) gene expression patterns in the airway epithelium. The similarities and differences of VE gene expression (ACE2, TMPRSS2, and CTSL) across nasal and bronchial compartments have not been fully characterized using matched samples from large cohorts. Gene expression data from 793 nasal and 1673 bronchial brushes obtained from individuals participating in lung cancer screening or diagnostic workup revealed that smoking status (current versus former) was the only clinical factor significantly and reproducibly associated with VE gene expression. The expression of ACE2 and TMPRSS2 was higher in smokers in the bronchus but not in the nose. scRNA-seq of nasal brushings indicated that ACE2 co-expressed genes were highly expressed in club and C15orf48+ secretory cells while TMPRSS2 co-expressed genes were highly expressed in keratinizing epithelial cells. In contrast, these ACE2 and TMPRSS2 modules were highly expressed in goblet cells in scRNA-seq from bronchial brushings. Cell-type deconvolution of the gene expression data confirmed that smoking increased the abundance of several secretory cell populations in the bronchus, but only goblet cells in the nose. The association of ACE2 and TMPRSS2 with smoking in the bronchus is due to their high expression in goblet cells which increase in abundance in current smoker airways. In contrast, in the nose, these genes are not predominantly expressed in cell populations modulated by smoking. In individuals with elevated lung cancer risk, smoking-induced VE gene expression changes in the nose likely have minimal impact on SARS-CoV-2 infection, but in the bronchus, smoking may lead to higher viral loads and more severe disease.

Evaluating the Influences of Confounding Variables on Benchmark Dose using a Case Study in the Field of Ionizing Radiation

Purpose A vast amount of data regarding the effects of radiation stressors on transcriptional changes has been produced over the past few decades. These data have shown remarkable consistency across platforms and experimental design, enabling increased understanding of early molecular effects of radiation exposure. However, the value of transcriptomic data in the context of risk assessment is not clear and represents a gap that is worthy of further consideration. Recently, benchmark dose (BMD) modeling has shown promise in correlating a transcriptional point of departure (POD) to that derived using phenotypic outcomes relevant to human health risk assessment. Although frequently applied in chemical toxicity evaluation, our group has recently demonstrated application within the field of radiation research. This approach allows the possibility to quantitatively compare radiation-induced gene and pathway alterations across various datasets using BMD values and derive meaningful biological effects. However, before BMD modeling can confidently be used, an understanding of the impact of confounding variables on BMD outputs is needed. Methods: To this end, BMD modeling was applied to a publicly available microarray dataset (Gene Expression Omnibus #GSE23515) that used peripheral blood ex-vivo gamma-irradiated at 0.82 Gy/min, at doses of 0, 0.1, 0.5 or 2 Gy, and assessed 6 hours post-exposure. The dataset comprised six female smokers (F-S), six female non-smokers (F-NS), six male smokers (M-S), and six male non-smokers (M-NS). Results: A combined total of 412 genes were fit to models and the BMD distribution was noted to be bi-modal across the four groups. A total of 74, 41, 62 and 62 genes were unique to the F-NS, M-NS, F-S and M-S groups. Sixty-two BMD modeled genes and nine pathways were common across all four groups. There were no differential sensitivity of responses in the robust common genes and pathways. Conclusion: For radiation-responsive genes and pathways common across the study groups, the BMD distribution of transcriptional activity was unaltered by sex and smoking status. Although further validation of the data is needed, these initial findings suggest BMD values for radiation relevant genes and pathways are robust and could be explored further in future studies.

Identification of apoptosis-associated protein factors distinctly expressed in cigarette smoke condensate-exposed airway bronchial epithelial cells

Smoking is associated with an increased risk of respiratory diseases, including lung cancer and asthma. However, the mechanisms or diagnostic markers for smoking-related diseases remain largely unknown. Here we investigated the role of cigarette smoke condensate (CSC) in the regulation of human bronchial epithelial cell (BEAS-2B) behavior. We found that exposure to CSC significantly inhibited BEAS-2B cell viability, impaired cell morphology, induced cell apoptosis, triggered oxidative damage, and promoted inflammatory response, which suggests a deleterious effect of CSC on bronchial epithelial cells. In addition, CSC markedly altered the expression of apoptosis-associated protein factors, including p21, soluble tumor necrosis factor receptor 1, and Fas ligand. In sum, our study identified a panel of novel protein factors that may mediate the actions of CSC on bronchial epithelial cells and have a predictive value for the development and progression of smoking-related diseases, thus providing insights into the development of potential diagnostic and therapeutic strategies against these diseases.

Toxicological effects of ambient fine (PM2.5-0.18) and ultrafine (PM0.18) particles in healthy and diseased 3D organo-typic mucocilary-phenotype models

The knowledge of the underlying mechanisms by which particulate matter (PM) exerts its health effects is still incomplete since it may trigger various symptoms as some persons may be more susceptible than others. Detailed studies realized in more relevant in vitro models are highly needed. Healthy normal human bronchial epithelial (NHBE), asthma-diseased human bronchial epithelial (DHBE), and COPD-DHBE cells, differentiated at the air-liquid interface, were acutely or repeatedly exposed to fine (i.e., PM_2.5-0.18, also called FP) and quasi-ultrafine (i.e., PM_0.18, also called UFP) particles. Immunofluorescence labelling of pan-cytokeratin, MUC5AC, and ZO-1 confirmed their specific cell-types. Baselines of the inflammatory mediators secreted by all the cells were quite similar. Slight changes of TNFα, IL-1β, IL-6, IL-8, GM-CSF, MCP-1, and/or TGFα, and of H3K9 histone acetylation supported a higher inflammatory response of asthma- and especially COPD-DHBE cells, after exposure to FP and especially UFP. At baseline, 35 differentially expressed genes (DEG) in asthma-DHBE, and 23 DEG in COPD-DHBE, compared to NHBE cells, were reported. They were involved in biological processes implicated in the development of asthma and COPD diseases, such as cellular process (e.g., PLA2G4C, NLRP1, S100A5, MUC1), biological regulation (e.g., CCNE1), developmental process (e.g., WNT10B), and cell component organization and synthesis (e.g., KRT34, COL6A1, COL6A2). In all the FP or UFP-exposed cell models, DEG were also functionally annotated to the chemical metabolic process (e.g., CYP1A1, CYP1B1, CYP1A2) and inflammatory response (e.g., EREG). Another DEG, FGF-1, was only down-regulated in asthma and specially COPD-DHBE cells repeatedly exposed. While RAB37 could help to counteract the down-regulation of FGF-1 in asthma-DHBE cells, the deregulation of FGR, WNT7B, VIPR1, and PPARGC1A could dramatically contribute to make it worse in COPD-DHBE cells. Taken together, these data contributed to support the highest effects of UFP versus FP and highest sensitivity of asthma- and notably COPD-DHBE versus NHBE cells.

Comparative mechanisms of PAH toxicity by benzo[a]pyrene and dibenzo[def,p]chrysene in primary human bronchial epithelial cells cultured at air-liquid interface

Current assumption for assessing carcinogenic risk of polycyclic aromatic hydrocarbons (PAHs) is that they function through a common mechanism of action; however, recent studies demonstrate that PAHs can act through unique mechanisms potentially contributing to cancer outcomes in a non-additive manner. Using a primary human 3D bronchial epithelial culture (HBEC) model, we assessed potential differences in mechanism of toxicity for two PAHs, benzo[a]pyrene (BAP) and dibenzo[def,p]chrysene (DBC), compared to a complex PAH mixture based on short-term biosignatures identified from transcriptional profiling. Differentiated bronchial epithelial cells were treated with BAP (100-500 μg/ml), DBC (10 μg/ml), and coal tar extract (CTE 500-1500 μg/ml, SRM1597a) for 48 h and gene expression was measured by RNA sequencing or quantitative PCR. Comparison of BAP and DBC gene signatures showed that the majority of genes (~60%) were uniquely regulated by treatment, including signaling pathways for inflammation and DNA damage by DBC and processes for cell cycle, hypoxia and oxidative stress by BAP. Specifically, BAP upregulated targets of AhR, NRF2, and KLF4, while DBC downregulated these same targets, suggesting a chemical-specific pattern in transcriptional regulation involved in antioxidant response, potentially contributing to differences in PAH potency. Other processes were regulated in common by all PAH treatments, BAP, DBC and CTE, including downregulation of genes involved in cell adhesion and reduced functional measurements of barrier integrity. This work supports prior in vivo studies and demonstrates the utility of profiling short-term biosignatures in an organotypic 3D model to identify mechanisms linked to carcinogenic risk of PAHs in humans.

Gene Expression Alterations in the Bronchial Epithelium of e-Cigarette Users

BACKGROUND

Although e-cigarette (ECIG) use has increased in the United States, their potential health effects remain uncertain. Understanding the effects of tobacco cigarette (TCIG) smoke on bronchial airway epithelial gene expression have previously provided insights into tobacco-related disease pathogenesis. Identifying the impact of ECIGs on airway gene expression could provide insights into their potential long-term health effects. We sought to compare the bronchial airway gene-expression profiles of former TCIG smokers now using ECIGs with the profiles of former and current TCIG smokers.

METHODS

We performed gene-expression profiling of bronchial epithelial cells collected from current TCIG smokers (n = 9), current ECIG users who are former TCIG smokers (n = 15), and former TCIG smokers (n = 21). We then compared our findings with previous studies of the effects of TCIG use on bronchial epithelium, as well an in vitro model of ECIG exposure.

RESULTS

Among 3,165 genes whose expression varied between the three study groups (q < 0.05), we identified 468 genes altered in ECIG users relative to former smokers (P < .05). Seventy-nine of these genes were up- or down-regulated concordantly among ECIG and TCIG users. We did not detect ECIG-associated gene-expression changes in known pathways associated with TCIG usage. Genes downregulated in ECIG users are enriched among the genes most downregulated by exposure of airway epithelium to ECIG vapor in vitro.

CONCLUSIONS

ECIGs induce both distinct and shared patterns of gene expression relative to TCIGs in the bronchial airway epithelium. The concordance of the genes altered in ECIG users and in the in vitro study suggests that genes altered in ECIG users are likely to be changed as the direct effect of ECIG exposure.

The effects of cigarette smoking extracts on cell cycle and tumor spread: novel evidence

Cigarette smoking is a major preventable risk factor for lung cancer, contributing to lung cancer progression and metastasis. Moreover, cigarette smoking correlates with increased metastasis frequency of pancreatic, breast and bladder cancer. The aim of this review was to examine the role of cigarette smoke extract in cell cycle and cancer progression. Clinical impact and the effects of cigarette smoke extract on carcinogenesis are discussed. 98 of the over 5000 chemicals in tobacco smoke are known carcinogens that can act on cancer genes such as K-RAS and p53. Through various mechanisms these compounds can activate molecules involved in the cell cycle, such as cyclins, and molecules involved in apoptosis and autophagy, such as Beclin-1 or LC3B. A search of the literature, including in vitro and in vivo studies, was carried out and the results summarized.

There is evidence of cancerogenic effects of cigarette smoke compounds. Cigarette smoke extract is a tobacco condensate obtained by filtration processes. Studies have shown that it can modify the cell cycle, inducing uncontrolled cell proliferation. This effect occurs through activation of genetic and epigenetic pathways and increasing the expression of proteins involved in inflammation. The pathways activated by cigarette smoke extract open up opportunities for researchers to develop new targeted therapies toward the specific molecules involved. Furthermore, the effects exerted by cigarette smoke extract on normal epithelial cells hold potential for use in the development of prevention medicine and early cancer diagnosis.

Persistence of smoking induced non-small cell lung carcinogenesis by decreasing ERBB pathway-related microRNA expression

Background

Tobacco use is responsible for approximately 80–90% of non‐small cell lung cancer cases. A large evidence base has shown that the ERBB pathway is associated with the occurrence of lung cancer. However, the mechanisms of how smoking activates the ERBB pathway have yet to be explained. We hypothesized that microRNAs may induce ERBB pathway activity during the process of lung cancer carcinogenesis.

Methods

We analyzed microRNA array data from the Gene Expression Omnibus and the Kyoto Encyclopedia of Genes and Genomes to determine any associations between genes and smoking in three groups of patients with NSCLC: smokers, former smokers, and non‐smokers.

Results

The interaction network among miRNAs, including hsa‐mir‐185‐3p, hsa‐mir‐4295, hsa‐mir‐4288, and hsa‐mir‐613, promotes lung cancer development by affecting the ERBB pathway.

Conclusion

Our findings provide evidence to explain the mechanism of lung cancer development in smokers.

Ramifications of secreted mucin MUC5AC in malignant journey: a holistic view

Heavily glycosylated secreted mucin MUC5AC, by the virtue of its cysteine-rich repeats, can form inter- and intramolecular disulfide linkages resulting in complex polymers, which in turn craft the framework of the polymeric mucus gel on epithelial cell surfaces. MUC5AC is a molecule with versatile functional implications including barrier functions to epithelial cells, host-pathogen interaction, immune cell attraction to sites of premalignant or malignant lesions and tumor progression in a context-dependent manner. Differential expression, glycosylation and localization of MUC5AC have been associated with a plethora of benign and malignant pathologies. In this era of robust technologies, overexpression strategies and genetically engineered mouse models, MUC5AC is emerging as a potential diagnostic, prognostic and therapeutic target for various malignancies. Considering the clinical relevance of MUC5AC, this review holistically encompasses its genomic organization, domain structure, glycosylation patterns, regulation, functional and molecular connotation from benign to malignant pathologies. Furthermore, we have here explored the incipient and significant experimental tools that are being developed to study this structurally complex and evolutionary conserved gel-forming mucin.

The roles of microRNAs in the pathogenesis of chronic obstructive pulmonary disease

Chronic obstructive pulmonary disease (COPD) is characterized by a progressive and irreversible airflow obstruction, with an abnormal lung function. The etiology of COPD correlates with complex interactions between environmental and genetic determinants. However, the exact pathogenesis of COPD is obscure although it involves multiple aspects including oxidative stress, imbalance between proteolytic and anti-proteolytic activity, immunity and inflammation, apoptosis, and repair and destruction in both airways and lungs. Many genes have been demonstrated to be involved in those pathogenic processes of this disease in patients exposed to harmful environmental factors. Previous reports have investigated promising microRNAs (miRNAs) to disclose the molecular mechanisms for COPD development induced by different environmental exposure and genetic predisposition encounter, and find some potential miRNA biomarkers for early diagnosis and treatment targets of COPD. In this review, we summarized the expression profiles of the reported miRNAs from studies of COPD associated with environmental risk factors including cigarette smoking and air pollution exposures, and provided an overview of roles of those miRNAs in the pathogenesis of the disease. We also highlighted the potential utility and limitations of miRNAs serving as diagnostic biomarkers and therapeutic targets for COPD.

The Airway Transcriptome as a Biomarker for Early Lung Cancer Detection

Lung cancer remains the leading cause of cancer-related death due to its advanced stage at diagnosis. Early detection of lung cancer can be improved by better defining who should be screened radiographically and determining which imaging-detected pulmonary nodules are malignant. Gene expression biomarkers measured in normal-appearing airway epithelium provide an opportunity to use lung cancer-associated molecular changes in this tissue for early detection of lung cancer. Molecular changes in the airway may result from an etiologic field of injury and/or field cancerization. The etiologic field of injury reflects the aberrant physiologic response to carcinogen exposure that creates a susceptible microenvironment for cancer initiation. In contrast, field cancerization reflects effects of "first-hit" mutations in a clone of cells from which the tumor ultimately arises or the effects of the tumor on the surrounding tissue. These fields might have value both for assessing lung cancer risk and diagnosis. Cancer-associated gene expression changes in the bronchial airway have recently been used to develop and validate a 23-gene classifier that improves the diagnostic yield of bronchoscopy for lung cancer among intermediate-risk patients. Recent studies have demonstrated that these lung cancer-related gene expression changes extend to nasal epithelial cells that can be sampled noninvasively. While the bronchial gene expression biomarker is being adopted clinically, further work is necessary to explore the potential clinical utility of gene expression profiling in the nasal epithelium for lung cancer diagnosis, lung cancer risk assessment, and precision medicine for lung cancer treatment and chemoprevention. Clin Cancer Res; 24(13); 2984-92. ©2018 AACR.

E-Cigarette Use Causes a Unique Innate Immune Response in the Lung, Involving Increased Neutrophilic Activation and Altered Mucin Secretion

RATIONALE

E-cigarettes have become increasingly popular and little is known about their potential adverse health effects.

OBJECTIVES

To determine the effects of e-cigarette use on the airways.

METHODS

Induced sputum samples from cigarette smokers, e-cigarette users, and nonsmokers were analyzed by quantitative proteomics, and the total and individual concentrations of mucins MUC5AC and MUC5B were determined by light scattering/refractometry and labeled mass spectrometry, respectively. Neutrophil extracellular trap (NET) formation rates were also determined for the same groups.

MEASUREMENTS AND MAIN RESULTS

E-cigarette users exhibited significant increases in aldehyde-detoxification and oxidative stress-related proteins associated with cigarette smoke compared with nonsmokers. The levels of innate defense proteins associated with chronic obstructive pulmonary disease, such as elastase and matrix metalloproteinase-9, were significantly elevated in e-cigarette users as well. E-cigarette users' sputum also uniquely exhibited significant increases in neutrophil granulocyte-related and NET-related proteins, such as myeloperoxidase, azurocidin, and protein-arginine deiminase 4, despite no significant elevation in neutrophil cell counts. Peripheral neutrophils from e-cigarette users showed increased susceptibility to phorbol 12-myristate 13-acetate-induced NETosis. Finally, a compositional change in the gel-forming building blocks of airway mucus (i.e., an elevated concentration of mucin MUC5AC) was observed in both cigarette smokers and e-cigarette users.

CONCLUSIONS

Together, our results indicate that e-cigarette use alters the profile of innate defense proteins in airway secretions, inducing similar and unique changes relative to cigarette smoking. These data challenge the concept that e-cigarettes are a healthier alternative to cigarettes.

Obstructive sleep apnea and rhonchopathy are associated with downregulation of trefoil factor family peptide 3 (TFF3)-Implications of changes in oral mucus composition

Study objectives

Trefoil factor family (TFF) peptides belong to the family of mucin-associated peptides and are expressed in most mucosal surfaces. TFF peptides carry out functions such as proliferation and migration enhancement, anti-apoptosis, and wound healing. Moreover, TFFs are associated with mucins and interact with them as “linker peptides”, thereby influencing mucus viscosity.

To test the hypothesis that in rhonchopathy and obstructive sleep apnea (OSA) changes occur in the expression of TFF3 and -2 that could contribute to changes in mucus viscosity, leading to an increase in upper airway resistance during breathing.

Methods

RT-PCR, Western-blot, immunohistochemistry and ELISA were performed to detect and quantify TFF3 and -2 in uvula samples. In addition, 99 saliva samples from patients with mild, moderate or severe OSA, as well as samples from rhonchopathy patients and from healthy volunteers, were analyzed by ELISA.

Results

TFF3 was detected in all uvula samples. Immunohistochemistry revealed a subjectively decreasing antibody reactivity of the uvula epithelia with increasing disease severity. ELISA demonstrated significantly higher TFF3 saliva protein concentrations in the healthy control group compared to cases with rhonchopathy and OSA. Predisposing factors of OSA such as BMI or age showed no correlation with TFF3. No significant changes were observed with regard to TFF2.

Conclusions

The results suggest the involvement of TFF3 in the pathogenesis of rhonchopathy and OSA and lead to the hypothesis that reduction of TFF3 production by the epithelium and subepithelial mucous glands of the uvula contribute to an increase in breathing resistance due to a change in mucus organization.

Cigarette smoke differentially affects IL-13-induced gene expression in human airway epithelial cells

Allergic airways inflammation in asthma is characterized by an airway epithelial gene signature composed of POSTN, CLCA1, and SERPINB2 This Th2 gene signature is proposed as a tool to classify patients with asthma into Th2-high and Th2-low phenotypes. However, many asthmatics smoke and the effects of cigarette smoke exposure on the epithelial Th2 gene signature are largely unknown. Therefore, we investigated the combined effect of IL-13 and whole cigarette smoke (CS) on the Th2 gene signature and the mucin-related genes MUC5AC and SPDEF in air-liquid interface differentiated human bronchial (ALI-PBEC) and tracheal epithelial cells (ALI-PTEC). Cultures were exposed to IL-13 for 14 days followed by 5 days of IL-13 with CS exposure. Alternatively, cultures were exposed once daily to CS for 14 days, followed by 5 days CS with IL-13. POSTN, SERPINB2, and CLCA1 expression were measured 24 h after the last exposure to CS and IL-13. In both models POSTN, SERPINB2, and CLCA1 expression were increased by IL-13. CS markedly affected the IL-13-induced Th2 gene signature as indicated by a reduced POSTN, CLCA1, and MUC5AC expression in both models. In contrast, IL-13-induced SERPINB2 expression remained unaffected by CS, whereas SPDEF expression was additively increased. Importantly, cessation of CS exposure failed to restore IL-13-induced POSTN and CLCA1 expression. We show for the first time that CS differentially affects the IL-13-induced gene signature for Th2-high asthma. These findings provide novel insights into the interaction between Th2 inflammation and cigarette smoke that is important for asthma pathogenesis and biomarker-guided therapy in asthma.

Shared Gene Expression Alterations in Nasal and Bronchial Epithelium for Lung Cancer Detection

Background

We previously derived and validated a bronchial epithelial gene expression biomarker to detect lung cancer in current and former smokers. Given that bronchial and nasal epithelial gene expression are similarly altered by cigarette smoke exposure, we sought to determine if cancer-associated gene expression might also be detectable in the more readily accessible nasal epithelium.

Methods

Nasal epithelial brushings were prospectively collected from current and former smokers undergoing diagnostic evaluation for pulmonary lesions suspicious for lung cancer in the AEGIS-1 (n = 375) and AEGIS-2 (n = 130) clinical trials and gene expression profiled using microarrays. All statistical tests were two-sided.

Results

We identified 535 genes that were differentially expressed in the nasal epithelium of AEGIS-1 patients diagnosed with lung cancer vs those with benign disease after one year of follow-up ( P  < .001). Using bronchial gene expression data from the AEGIS-1 patients, we found statistically significant concordant cancer-associated gene expression alterations between the two airway sites ( P  < .001). Differentially expressed genes in the nose were enriched for genes associated with the regulation of apoptosis and immune system signaling. A nasal lung cancer classifier derived in the AEGIS-1 cohort that combined clinical factors (age, smoking status, time since quit, mass size) and nasal gene expression (30 genes) had statistically significantly higher area under the curve (0.81; 95% confidence interval [CI] = 0.74 to 0.89, P  = .01) and sensitivity (0.91; 95% CI = 0.81 to 0.97, P  = .03) than a clinical-factor only model in independent samples from the AEGIS-2 cohort.

Conclusions

These results support that the airway epithelial field of lung cancer-associated injury in ever smokers extends to the nose and demonstrates the potential of using nasal gene expression as a noninvasive biomarker for lung cancer detection.

Detecting the Presence and Progression of Premalignant Lung Lesions via Airway Gene Expression

Purpose: Lung cancer is the leading cause of cancer-related death in the United States. The molecular events preceding the onset of disease are poorly understood, and no effective tools exist to identify smokers with premalignant lesions (PMLs) that will progress to invasive cancer. Prior work identified molecular alterations in the smoke-exposed airway field of injury associated with lung cancer. Here, we focus on an earlier stage in the disease process leveraging the airway field of injury to study PMLs and its utility in lung cancer chemoprevention.Experimental Design: Bronchial epithelial cells from normal appearing bronchial mucosa were profiled by mRNA-Seq from subjects with (n = 50) and without (n = 25) PMLs. Using surrogate variable and gene set enrichment analysis, we identified genes, pathways, and lung cancer-related gene sets differentially expressed between subjects with and without PMLs. A computational pipeline was developed to build and test a chemoprevention-relevant biomarker.Results: We identified 280 genes in the airway field associated with the presence of PMLs. Among the upregulated genes, oxidative phosphorylation was strongly enriched, and IHC and bioenergetics studies confirmed pathway findings in PMLs. The relationship between PMLs and squamous cell carcinomas (SCC) was also confirmed using published lung cancer datasets. The biomarker performed well predicting the presence of PMLs (AUC = 0.92, n = 17), and changes in the biomarker score associated with progression/stability versus regression of PMLs (AUC = 0.75, n = 51).Conclusions: Transcriptomic alterations in the airway field of smokers with PMLs reflect metabolic and early lung SCC alterations and may be leveraged to stratify smokers at high risk for PML progression and monitor outcome in chemoprevention trials. Clin Cancer Res; 23(17); 5091-100. ©2017 AACR.

Reduced biological effect of e-cigarette aerosol compared to cigarette smoke evaluated in vitro using normalized nicotine dose and RNA-seq-based toxicogenomics

Electronic cigarettes (e-cigarettes) use has increased globally and could potentially offer a lower risk alternative to cigarette smoking. Here, we assessed the transcriptional response of a primary 3D airway model acutely exposed to e-cigarette aerosol and cigarette (3R4F) smoke. Aerosols were generated with standard intense smoking regimens with careful consideration for dose by normalizing the exposures to nicotine. Two e-cigarette aerosol dilutions were tested for equivalent and higher nicotine delivery compared to 3R4F. RNA was extracted at 24 hrs and 48 hrs post exposure for RNA-seq. 873 and 205 RNAs were differentially expressed for 3R4F smoke at 24 hrs and 48 hrs using a pFDR < 0.01 and a [fold change] > 2 threshold. 113 RNAs were differentially expressed at the highest dose of e-cigarette aerosol using a looser threshold of pFDR < 0.05, 3 RNAs exceeded a fold change of 2. Geneset enrichment analysis revealed a clear response from lung cancer, inflammation, and fibrosis associated genes after 3R4F smoke exposure. Metabolic/biosynthetic processes, extracellular membrane, apoptosis, and hypoxia were identified for e-cigarette exposures, albeit with a lower confidence score. Based on equivalent or higher nicotine delivery, an acute exposure to e-cigarette aerosol had a reduced impact on gene expression compared to 3R4F smoke exposure in vitro.

Analysis to Estimate Genetic Variations in the Idarubicin-Resistant Derivative MOLT-3

Gene alterations are a well-established mechanism leading to drug resistance in acute leukemia cells. A full understanding of the mechanisms of drug resistance in these cells will facilitate more effective chemotherapy. In this study, we investigated the mechanism(s) of drug resistance in the human acute leukemia cell line MOLT-3 and its idarubicin-resistant derivative MOLT-3/IDR through complete mitochondrial and nuclear DNA analyses. We identified genetic differences between these two cell lines. The ND3 mutation site (p.Thr61Ile) in the mitochondrial DNA sequence was unique to MOLT-3/IDR cells. Moreover, we identified five candidate genes harboring genetic alterations, including GALNT2, via CGH array analysis. Sequencing of the GALNT2 exon revealed a G1716K mutation present within the stop codon in MOLT-3/IDR cells but absent from MOLT-3 cells. This mutation led to an additional 18 amino acids in the protein encoded by GALNT2. Using real-time PCR, we determined an expression value for this gene of 0.35. Protein structure predictions confirmed a structural change in GALNT2 in MOLT-3/IDR cells that corresponded to the site of the mutation. We speculate that this mutation may be related to idarubicin resistance.

Transcriptome sequencing reveals e-cigarette vapor and mainstream-smoke from tobacco cigarettes activate different gene expression profiles in human bronchial epithelial cells

Electronic cigarettes (e-cigarettes) generate an aerosol vapor (e-vapor) thought to represent a less risky alternative to main stream smoke (MSS) of conventional tobacco cigarettes. RNA-seq analysis was used to examine the transcriptomes of differentiated human bronchial epithelial (HBE) cells exposed to air, MSS from 1R5F tobacco reference cigarettes, and e-vapor with and without added nicotine in an in vitro air-liquid interface model for cellular exposure. Our results indicate that while e-vapor does not elicit many of the cell toxicity responses observed in MSS-exposed HBE cells, e-vapor exposure is not benign, but elicits discrete transcriptomic signatures with and without added nicotine. Among the cellular pathways with the most significantly enriched gene expression following e-vapor exposure are the phospholipid and fatty acid triacylglycerol metabolism pathways. Our data suggest that alterations in cellular glycerophopholipid biosynthesis are an important consequences of e-vapor exposure. Moreover, the presence of nicotine in e-vapor elicits a cellular response distinct from e-vapor alone including alterations of cytochrome P450 function, retinoid metabolism, and nicotine catabolism. These studies establish a baseline for future analysis of e-vapor and e-vapor additives that will better inform the FDA and other governmental bodies in discussions of the risks and future regulation of these products.

Persistence of smoking-induced dysregulation of miRNA expression in the small airway epithelium despite smoking cessation

Even after quitting smoking, the risk of the development of chronic obstructive pulmonary disease (COPD) and lung cancer remains significantly higher compared to healthy nonsmokers. Based on the knowledge that COPD and most lung cancers start in the small airway epithelium (SAE), we hypothesized that smoking modulates miRNA expression in the SAE linked to the pathogenesis of smoking-induced airway disease, and that some of these changes persist after smoking cessation. SAE was collected from 10^th to 12^th order bronchi using fiberoptic bronchoscopy. Affymetrix miRNA 2.0 arrays were used to assess miRNA expression in the SAE from 9 healthy nonsmokers and 10 healthy smokers, before and after they quit smoking for 3 months. Smoking status was determined by urine nicotine and cotinine measurement. There were significant differences in the expression of 34 miRNAs between healthy smokers and healthy nonsmokers (p<0.01, fold-change >1.5), with functions associated with lung development, airway epithelium differentiation, inflammation and cancer. After quitting smoking for 3 months, 12 out of the 34 miRNAs did not return to normal levels, with Wnt/β-catenin signaling pathway being the top identified enriched pathway of the target genes of the persistent dysregulated miRNAs. In the context that many of these persistent smoking-dependent miRNAs are associated with differentiation, inflammatory diseases or lung cancer, it is likely that persistent smoking-related changes in SAE miRNAs play a role in the subsequent development of these disorders.

Cytotoxicity of chronic exposure to 4 cigarette smoke condensates in 2 cell lines

Tobacco use is the leading preventable cause of death. The cytotoxicity of cigarette smoke condensate (CSC), the particulate fraction of cigarette smoke without the vapor phase, has mostly been tested in short-term in vitro studies lasting from a few hours to a few days. Here, we assessed the toxicity of CSCs from 2 reference cigarettes, 3R4F and CM6, using a primary human small airway epithelial (PSAE) cell line by quantifying adenosine 5'-triphosphate (ATP), 3-(4,5-dimethylthiazol-2-yl)-5-(3-carboxy-methoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium (MTS), total glutathione (reduced glutathione [GSH] + oxidized glutathione [GSSG]), and lactate dehydrogenase (LDH) release over the course of 28 days. The CSCs, 0.3 to 10 μg/mL, promoted cell proliferation at 120 hours of exposure, but demonstrated cytotoxicity at days 14 and 28. Interestingly, CSCs, 0.3 to 3 μg/mL, showed a cell death effect at day 14 but induced cell proliferation at day 28. Consistently, transformation associated with morphological changes began by day 14 and the transformed cells grew dramatically at day 28. The LDH assay appeared to be sensitive for assessing early cell damage, whereas the ATP, MTS, and GSH assays were more suitable for determining later stage CSCs-induced cytotoxicity. The ATP assay showed greater sensitivity than the MTS and GSH assays. We also assessed the toxicity of CSCs in an human Telomerase Reverse Transcriptase (hTERT)-immortalized Barrett esophagus cell line (CP-C). The CP-C cells demonstrated dose- and time-dependent cytotoxicity over the course of 28 days but displayed higher resistance to CSCs than PSAE cells. This study demonstrates that CSCs cause cytotoxicity and induce transformation related to cell resistance and cell invasion properties.

Sex and Smoking Status Effects on the Early Detection of Early Lung Cancer in High-Risk Smokers Using an Electronic Nose

OBJECTIVE

Volatile organic compounds (VOCs) in exhaled breath as measured by electronic nose (e-nose) have utility as biomarkers to detect subjects at risk of having lung cancer in a screening setting. We hypothesize that breath analysis using an e-nose chemo-resistive sensor array could be used as a screening tool to discriminate patients diagnosed with lung cancer from high-risk smokers.

METHODS

Breath samples from 191 subjects-25 lung cancer patients and 166 high-risk smoker control subjects without cancer-were analyzed. For clinical relevancy, subjects in both groups were matched for age, sex, and smoking histories. Classification and regression trees and discriminant functions classifiers were used to recognize VOC patterns in e-nose data. Cross-validated results were used to assess classification accuracy. Repeatability and reproducibility of e-nose data were assessed by measuring subject-exhaled breath in parallel across two e-nose devices.

RESULTS

e-Nose measurements could distinguish lung cancer patients from high-risk control subjects, with a better than 80% classification accuracy. Subject sex and smoking status impacted classification as area under the curve results (ex-smoker males 0.846, ex-smoker female 0.816, current smoker male 0.745, and current smoker female 0.725) demonstrated. Two e-nose systems could be calibrated to give equivalent readings across subject-exhaled breath measured in parallel.

CONCLUSIONS

e-Nose technology may have significant utility as a noninvasive screening tool for detecting individuals at increased risk for lung cancer.

SIGNIFICANCE

The results presented further the case that VOC patterns could have real clinical utility to screen for lung cancer in the important growing ex-smoker population.

Biomarkers of tobacco smoke exposure

Diseases and death caused by exposure to tobacco smoke have become the single most serious preventable public health concern. Thus, biomarkers that can monitor tobacco exposure and health effects can play a critical role in tobacco product regulation and public health policy. Biomarkers of exposure to tobacco toxicants are well established and have been used in population studies to establish public policy regarding exposure to second-hand smoke, an example being the nicotine metabolite cotinine, which can be measured in urine. Biomarkers of biological response to tobacco smoking range from those indicative of inflammation to mRNA and microRNA patterns related to tobacco use and/or disease state. Biomarkers identifying individuals with an increased risk for a pathological response to tobacco have also been described. The challenge for any novel technology or biomarker is its translation to clinical and/or regulatory application, a process that requires first technical validation of the assay and then careful consideration of the context the biomarker assay may be used in the regulatory setting. Nonetheless, the current efforts to investigate new biomarker of tobacco smoke exposure promise to offer powerful new tools in addressing the health hazards of tobacco product use. This review will examine such biomarkers, albeit with a focus on those related to cigarette smoking.

Smoking status impacts microRNA mediated prognosis and lung adenocarcinoma biology

Background

Cigarette smoke is associated with the majority of lung cancers: however, 25% of lung cancer patients are non-smokers, and half of all newly diagnosed lung cancer patients are former smokers. Lung tumors exhibit distinct epidemiological, clinical, pathological, and molecular features depending on smoking status, suggesting divergent mechanisms underlie tumorigenesis in smokers and non-smokers. MicroRNAs (miRNAs) are integral contributors to tumorigenesis and mediate biological responses to smoking. Based on the hypothesis that smoking-specific miRNA differences in lung adenocarcinomas reflect distinct tumorigenic processes selected by different smoking and non-smoking environments, we investigated the contribution of miRNA disruption to lung tumor biology and patient outcome in the context of smoking status.

Methods

We applied a whole transcriptome sequencing based approach to interrogate miRNA levels in 94 patient-matched lung adenocarcinoma and non-malignant lung parenchymal tissue pairs from current, former and never smokers.

Results

We discovered novel and distinct smoking status-specific patterns of miRNA and miRNA-mediated gene networks, and identified miRNAs that were prognostically significant in a smoking dependent manner.

Conclusions

We conclude that miRNAs disrupted in a smoking status-dependent manner affect distinct cellular pathways and differentially influence lung cancer patient prognosis in current, former and never smokers. Our findings may represent promising biologically relevant markers for lung cancer prognosis or therapeutic intervention.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2407-14-778) contains supplementary material, which is available to authorized users.

Buccal spectral markers for lung cancer risk stratification

Lung cancer remains the leading cause of cancer deaths in the US with >150,000 deaths per year. In order to more effectively reduce lung cancer mortality, more sophisticated screening paradigms are needed. Previously, our group demonstrated the use of low-coherence enhanced backscattering (LEBS) spectroscopy to detect and quantify the micro/nano-architectural correlates of colorectal and pancreatic field carcinogenesis. In the lung, the buccal (cheek) mucosa has been suggested as an excellent surrogate site in the “field of injury”. We, therefore, wanted to assess whether LEBS could similarly sense the presence of lung. To this end, we applied a fiber-optic LEBS probe to a dataset of 27 smokers without diagnosed lung cancer (controls) and 46 with lung cancer (cases), which was divided into a training and a blinded validation set (32 and 41 subjects, respectively). LEBS readings of the buccal mucosa were taken from the oral cavity applying gentle contact. The diagnostic LEBS marker was notably altered in patients harboring lung cancer compared to smoking controls. The prediction rule developed on training set data provided excellent diagnostics with 94% sensitivity, 80% specificity, and 95% accuracy. Applying the same threshold to the blinded validation set yielded 79% sensitivity and 83% specificity. These results were not confounded by patient demographics or impacted by cancer type or location. Moreover, the prediction rule was robust across all stages of cancer including stage I. We envision the use of LEBS as the first part of a two-step paradigm shift in lung cancer screening in which patients with high LEBS risk markers are funnelled into more invasive screening for confirmation.

AP-2α downregulation by cigarette smoke condensate is counteracted by p53 in human lung cancer cells

Cumulative findings have demonstrated that the dysregulation of tumor suppressor genes may be implicated in cigarette smoke-induced carcinogenesis. Activating enhancer-binding protein 2 (AP-2) is a eukaryotic transcriptional factor that plays a significant role in embryonic development and tumorigenesis. The vertebrate AP-2 family consists of AP-2α, AP-2β, AP-2γ, AP-2δ and AP-2ε. Previous studies have suggested that cigarette smoking disrupts AP-2 regulation. In the present study, we investigated the effects of cigarette smoke condensate (CSC) on AP-2α expression in human lung cancer cell lines (NCI-H1299, NCI-H446 and A549), as well as the potential mechanisms involved. Using RT-qPCR, we found that CSC decreased AP-2α expression by suppressing its transcription in human lung cancer cell lines, particularly in p53-deficient NCI-H1299 cells. Western blotting and luciferase assays were implemented and we found that the restoration of p53 expression rescued the NCI-H1299 cells from CSC-induced AP-2α loss, while the silencing of p53 resulted in increased AP-2α loss induced by CSC, suggesting an antagonizing role of p53 in the regulation of AP-2α by CSC. Our results indicate that AP-2α downregulation may be involved in smoke-induced lung carcinogenesis.

Quantification of biological network perturbations for mechanistic insight and diagnostics using two-layer causal models

BACKGROUND

High-throughput measurement technologies such as microarrays provide complex datasets reflecting mechanisms perturbed in an experiment, typically a treatment vs. control design. Analysis of these information rich data can be guided based on a priori knowledge, such as networks or set of related proteins or genes. Among those, cause-and-effect network models are becoming increasingly popular and more than eighty such models, describing processes involved in cell proliferation, cell fate, cell stress, and inflammation have already been published. A meaningful systems toxicology approach to study the response of a cell system, or organism, exposed to bio-active substances requires a quantitative measure of dose-response at network level, to go beyond the differential expression of single genes.

RESULTS

We developed a method that quantifies network response in an interpretable manner. It fully exploits the (signed graph) structure of cause-and-effect networks models to integrate and mine transcriptomics measurements. The presented approach also enables the extraction of network-based signatures for predicting a phenotype of interest. The obtained signatures are coherent with the underlying network perturbation and can lead to more robust predictions across independent studies. The value of the various components of our mathematically coherent approach is substantiated using several in vivo and in vitro transcriptomics datasets. As a proof-of-principle, our methodology was applied to unravel mechanisms related to the efficacy of a specific anti-inflammatory drug in patients suffering from ulcerative colitis. A plausible mechanistic explanation of the unequal efficacy of the drug is provided. Moreover, by utilizing the underlying mechanisms, an accurate and robust network-based diagnosis was built to predict the response to the treatment.

CONCLUSION

The presented framework efficiently integrates transcriptomics data and "cause and effect" network models to enable a mathematically coherent framework from quantitative impact assessment and data interpretation to patient stratification for diagnosis purposes.

MicroRNA 4423 is a primate-specific regulator of airway epithelial cell differentiation and lung carcinogenesis

Smoking is a significant risk factor for lung cancer, the leading cause of cancer-related deaths worldwide. Although microRNAs are regulators of many airway gene-expression changes induced by smoking, their role in modulating changes associated with lung cancer in these cells remains unknown. Here, we use next-generation sequencing of small RNAs in the airway to identify microRNA 4423 (miR-4423) as a primate-specific microRNA associated with lung cancer and expressed primarily in mucociliary epithelium. The endogenous expression of miR-4423 increases as bronchial epithelial cells undergo differentiation into mucociliary epithelium in vitro, and its overexpression during this process causes an increase in the number of ciliated cells. Furthermore, expression of miR-4423 is reduced in most lung tumors and in cytologically normal epithelium of the mainstem bronchus of smokers with lung cancer. In addition, ectopic expression of miR-4423 in a subset of lung cancer cell lines reduces their anchorage-independent growth and significantly decreases the size of the tumors formed in a mouse xenograft model. Consistent with these phenotypes, overexpression of miR-4423 induces a differentiated-like pattern of airway epithelium gene expression and reverses the expression of many genes that are altered in lung cancer. Together, our results indicate that miR-4423 is a regulator of airway epithelium differentiation and that the abrogation of its function contributes to lung carcinogenesis.

Case study: the role of mechanistic network models in systems toxicology

Twenty first century systems toxicology approaches enable the discovery of biological pathways affected in response to active substances. Here, we briefly summarize current network approaches that facilitate the detailed mechanistic understanding of the impact of a given stimulus on a biological system. We also introduce our network-based method with two use cases and show how causal biological network models combined with computational methods provide quantitative mechanistic insights. Our approach provides a robust comparison of the transcriptional responses in different experimental systems and enables the identification of network-based biomarkers modulated in response to exposure. These advances can also be applied to pharmacology, where the understanding of disease mechanisms and adverse drug effects is imperative for the development of efficient and safe treatment options.

Cigarette smoke and CFTR: implications in the pathogenesis of COPD

Chronic obstructive pulmonary disease (COPD) is a progressive respiratory disorder consisting of chronic bronchitis and/or emphysema. COPD patients suffer from chronic infections and display exaggerated inflammatory responses and a progressive decline in respiratory function. The respiratory symptoms of COPD are similar to those seen in cystic fibrosis (CF), although the molecular basis of the two disorders differs. CF is a genetic disease caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene encoding a chloride and bicarbonate channel (CFTR), leading to CFTR dysfunction. The majority of COPD cases result from chronic oxidative insults such as cigarette smoke. Interestingly, environmental stresses including cigarette smoke, hypoxia, and chronic inflammation have also been implicated in reduced CFTR function, and this suggests a common mechanism that may contribute to both the CF and COPD. Therefore, improving CFTR function may offer an excellent opportunity for the development of a common treatment for CF and COPD. In this article, we review what is known about the CF respiratory phenotype and discuss how diminished CFTR expression-associated ion transport defects may contribute to some of the pathological changes seen in COPD.

Cancer testis antigens: a novel target in lung cancer

Lung cancer is the main cause of cancer mortality worldwide. This is mainly due to the fact that it is diagnosed in advanced stage patients, which are no more surgically curable. Consequently, searching for novel treatments and new modalities for early diagnosis offers great promise to improve the clinical outcome. Recently, a new group of antigens, the cancer testis antigens, have been described as possible early diagnostic tools and therapeutic targets in cancer therapy.This review will report emerging evidences of cancer testis antigens deregulation in lung cancer and explore the state of the art of their currently known role and potential as markers for early diagnosis and disease progression and targets of an immunotherapeutic approach aiming to improve the cure rate of this tumor.

Activation of uPAR is required for cigarette smoke extract-induced epithelial-mesenchymal transition in lung epithelial cells

Cigarette smoke is a major risk factor for lung cancer, which may contribute to lung cancer invasion and metastasis. However, the mechanism remains unclear. Epithelial-mesenchymal transition (EMT) is a critical phenotypic alteration of cells that triggers invasion and metastasis. The urokinase-type plasminogen activator receptor (uPAR) is originally thought to assist the directional invasion of migrating cells, and increasing evidences show that overexpression of uPAR in cancer cells promotes EMT. Therefore, we intend to study the role of uPAR in cigarette smoke extract (CSE)-induced EMT in lung epithelial cells. In this study, we showed that lung epithelial cells cultured after CSE treatment demonstrated changes consistent with EMT. E-cadherin was decreased, while vimentin, N-cadherin, and α-SMA expression was increased in both A549 and BEAS-2B cells. Cells acquired a mesenchymal-like morphology and increased cell motility and invasion. In addition, CSE-induced EMT was accompanied by increased expression of uPAR and activation of AKT downstream of uPAR. CSE-induced EMT and activation of AKT were blocked by uPAR gene silencing. Antagonizing PI3K also inhibits development of CSE-induced EMT. We conclude that CSE can induce EMT, and the activity of uPAR-dependent signal pathway in EMT is recapitulated in lung epithelial cells in vitro.

SIRT1 pathway dysregulation in the smoke-exposed airway epithelium and lung tumor tissue

Cigarette smoke produces a molecular field of injury in epithelial cells lining the respiratory tract. However, the specific signaling pathways that are altered in the airway of smokers and the signaling processes responsible for the transition from smoking-induced airway damage to lung cancer remain unknown. In this study, we use a genomic approach to study the signaling processes associated with tobacco smoke exposure and lung cancer. First, we developed and validated pathway-specific gene expression signatures in bronchial airway epithelium that reflect activation of signaling pathways relevant to tobacco exposure, including ATM, BCL2, GPX1, NOS2, IKBKB, and SIRT1. Using these profiles and four independent gene expression datasets, we found that SIRT1 activity is significantly upregulated in cytologically normal bronchial airway epithelial cells from active smokers compared with nonsmokers. In contrast, this activity is strikingly downregulated in non-small cell lung cancer. This pattern of signaling modulation was unique to SIRT1, and downregulation of SIRT1 activity is confined to tumors from smokers. Decreased activity of SIRT1 was validated using genomic analyses of mouse models of lung cancer and biochemical testing of SIRT1 activity in patient lung tumors. Together, our findings indicate a role of SIRT1 in response to smoke and a potential role in repressing lung cancer. Furthermore, our findings suggest that the airway gene expression signatures derived in this study can provide novel insights into signaling pathways altered in the "field of injury" induced by tobacco smoke and thus may impact strategies for prevention of tobacco-related lung cancer.

Genes associated with MUC5AC expression in small airway epithelium of human smokers and non-smokers

Background

Mucus hypersecretion contributes to the morbidity and mortality of smoking-related lung diseases, especially chronic obstructive pulmonary disease (COPD), which starts in the small airways. Despite progress in animal studies, the genes and their expression pattern involved in mucus production and secretion in human airway epithelium are not well understood. We hypothesized that comparison of the transcriptomes of the small airway epithelium of individuals that express high vs low levels of MUC5AC, the major macromolecular component of airway mucus, could be used as a probe to identify the genes related to human small airway mucus production/secretion.

Methods

Flexible bronchoscopy and brushing were used to obtain small airway epithelium (10^th to 12^th order bronchi) from healthy nonsmokers (n=60) and healthy smokers (n=72). Affymetrix HG-U133 plus 2.0 microarrays were used to assess gene expression. Massive parallel sequencing (RNA-Seq) was used to verify gene expression of small airway epithelium from 5 nonsmokers and 6 smokers.

Results

MUC5AC expression varied 31-fold among the healthy nonsmokers. Genome-wide comparison between healthy nonsmokers (n = 60) grouped as “high MUC5AC expressors” vs “low MUC5AC expressors” identified 528 genes significantly up-regulated and 15 genes significantly down-regulated in the high vs low expressors. This strategy identified both mucus production and secretion related genes under control of a network composed of multiple transcription factors. Based on the literature, genes in the up-regulated list were used to identify a 73 “MUC5AC-associated core gene” list with 9 categories: mucus component; mucus-producing cell differentiation-related transcription factor; mucus-producing cell differentiation-related pathway or mediator; post-translational modification of mucin; vesicle transport; endoplasmic reticulum stress-related; secretory granule-associated; mucus secretion-related regulator and mucus hypersecretory-related ion channel. As a validation cohort, we assessed the MUC5AC-associated core gene list in the small airway epithelium of an independent set of healthy smokers (n = 72). There was up-regulation of MUC5AC in the small airway epithelium of smokers (2.3-fold, p < 10^-8) associated with a coordinated up-regulation of MUC5AC-associated core gene expression pattern in the small airway epithelium of smokers (p < 0.01). Deep sequencing confirmed these observations.

Conclusion

The identification of the genes associated with increased airway mucin production in humans should be useful in understanding the pathogenesis of airway mucus hypersecretion and identifying therapeutic targets.

Author summary

Mucus hypersecretion contributes to the morbidity and mortality of smoking-related lung diseases, especially chronic obstructive pulmonary disease (COPD), which starts in the small airways. Little is known about the gene networks associated with the synthesis and secretion of mucins in the human small airway epithelium. Taking advantage of the knowledge that MUC5AC is a major mucin secreted by the small airway epithelium, the expression of MUC5AC in small airway epithelium is highly regulated at the transcriptional level and our observation that healthy nonsmokers have variable numbers of MUC5AC⁺ secretory cells in the human small airway epithelium, we compared genome-wide gene expression of the small airway epithelium of high vs low MUC5AC expressors from 60 nonsmokers to identify the genes associated with MUC5AC expression. This novel strategy enabled identification of a 73 “MUC5AC-associated core gene” list with 9 categories, which control a series of processes from mucin biosynthesis to mucus secretion. The coordinated gene expression pattern of MUC5AC-associated core genes were corroborated in an independent cohort of 72 healthy smokers. Deep sequencing of small airway epithelium RNA confirmed these observations. This finding will be useful in identifying therapeutic targets to treat small airway mucus hypersecretion.

Relation between smoking history and gene expression profiles in lung adenocarcinomas

Background

Lung cancer is the worldwide leading cause of death from cancer. Tobacco usage is the major pathogenic factor, but all lung cancers are not attributable to smoking. Specifically, lung cancer in never-smokers has been suggested to represent a distinct disease entity compared to lung cancer arising in smokers due to differences in etiology, natural history and response to specific treatment regimes. However, the genetic aberrations that differ between smokers and never-smokers’ lung carcinomas remain to a large extent unclear.

Methods

Unsupervised gene expression analysis of 39 primary lung adenocarcinomas was performed using Illumina HT-12 microarrays. Results from unsupervised analysis were validated in six external adenocarcinoma data sets (n=687), and six data sets comprising normal airway epithelial or normal lung tissue specimens (n=467). Supervised gene expression analysis between smokers and never-smokers were performed in seven adenocarcinoma data sets, and results validated in the six normal data sets.

Results

Initial unsupervised analysis of 39 adenocarcinomas identified two subgroups of which one harbored all never-smokers. A generated gene expression signature could subsequently identify never-smokers with 79-100% sensitivity in external adenocarcinoma data sets and with 76-88% sensitivity in the normal materials. A notable fraction of current/former smokers were grouped with never-smokers. Intriguingly, supervised analysis of never-smokers versus smokers in seven adenocarcinoma data sets generated similar results. Overlap in classification between the two approaches was high, indicating that both approaches identify a common set of samples from current/former smokers as potential never-smokers. The gene signature from unsupervised analysis included several genes implicated in lung tumorigenesis, immune-response associated pathways, genes previously associated with smoking, as well as marker genes for alveolar type II pneumocytes, while the best classifier from supervised analysis comprised genes strongly associated with proliferation, but also genes previously associated with smoking.

Conclusions

Based on gene expression profiling, we demonstrate that never-smokers can be identified with high sensitivity in both tumor material and normal airway epithelial specimens. Our results indicate that tumors arising in never-smokers, together with a subset of tumors from smokers, represent a distinct entity of lung adenocarcinomas. Taken together, these analyses provide further insight into the transcriptional patterns occurring in lung adenocarcinoma stratified by smoking history.

Molecular signature of smoking in human lung tissues

Cigarette smoking is the leading risk factor for lung cancer. To identify genes deregulated by smoking and to distinguish gene expression changes that are reversible and persistent following smoking cessation, we carried out genome-wide gene expression profiling on nontumor lung tissue from 853 patients with lung cancer. Gene expression levels were compared between never and current smokers, and time-dependent changes in gene expression were studied in former smokers. A total of 3,223 transcripts were differentially expressed between smoking groups in the discovery set (n = 344, P < 1.29 × 10(-6)). A substantial number of smoking-induced genes also were validated in two replication sets (n = 285 and 224), and a gene expression signature of 599 transcripts consistently segregated never from current smokers across all three sets. The expression of the majority of these genes reverted to never-smoker levels following smoking cessation, although the time course of normalization differed widely among transcripts. Moreover, some genes showed very slow or no reversibility in expression, including SERPIND1, which was found to be the most consistent gene permanently altered by smoking in the three sets. Our findings therefore indicate that smoking deregulates many genes, many of which reverse to normal following smoking cessation. However, a subset of genes remains altered even decades following smoking cessation and may account, at least in part, for the residual risk of lung cancer among former smokers. Cancer Res; 72(15); 3753-63. ©2012 AACR.

Divergent genomic and epigenomic landscapes of lung cancer subtypes underscore the selection of different oncogenic pathways during tumor development

For therapeutic purposes, non-small cell lung cancer (NSCLC) has traditionally been regarded as a single disease. However, recent evidence suggest that the two major subtypes of NSCLC, adenocarcinoma (AC) and squamous cell carcinoma (SqCC) respond differently to both molecular targeted and new generation chemotherapies. Therefore, identifying the molecular differences between these tumor types may impact novel treatment strategy. We performed the first large-scale analysis of 261 primary NSCLC tumors (169 AC and 92 SqCC), integrating genome-wide DNA copy number, methylation and gene expression profiles to identify subtype-specific molecular alterations relevant to new agent design and choice of therapy. Comparison of AC and SqCC genomic and epigenomic landscapes revealed 778 altered genes with corresponding expression changes that are selected during tumor development in a subtype-specific manner. Analysis of >200 additional NSCLCs confirmed that these genes are responsible for driving the differential development and resulting phenotypes of AC and SqCC. Importantly, we identified key oncogenic pathways disrupted in each subtype that likely serve as the basis for their differential tumor biology and clinical outcomes. Downregulation of HNF4α target genes was the most common pathway specific to AC, while SqCC demonstrated disruption of numerous histone modifying enzymes as well as the transcription factor E2F1. In silico screening of candidate therapeutic compounds using subtype-specific pathway components identified HDAC and PI3K inhibitors as potential treatments tailored to lung SqCC. Together, our findings suggest that AC and SqCC develop through distinct pathogenetic pathways that have significant implication in our approach to the clinical management of NSCLC.

RNA-Seq quantification of the human small airway epithelium transcriptome

Background

The small airway epithelium (SAE), the cell population that covers the human airway surface from the 6^th generation of airway branching to the alveoli, is the major site of lung disease caused by smoking. The focus of this study is to provide quantitative assessment of the SAE transcriptome in the resting state and in response to chronic cigarette smoking using massive parallel mRNA sequencing (RNA-Seq).

Results

The data demonstrate that 48% of SAE expressed genes are ubiquitous, shared with many tissues, with 52% enriched in this cell population. The most highly expressed gene, SCGB1A1, is characteristic of Clara cells, the cell type unique to the human SAE. Among other genes expressed by the SAE are those related to Clara cell differentiation, secretory mucosal defense, and mucociliary differentiation. The high sensitivity of RNA-Seq permitted quantification of gene expression related to infrequent cell populations such as neuroendocrine cells and epithelial stem/progenitor cells. Quantification of the absolute smoking-induced changes in SAE gene expression revealed that, compared to ubiquitous genes, more SAE-enriched genes responded to smoking with up-regulation, and those with the highest basal expression levels showed most dramatic changes. Smoking had no effect on SAE gene splicing, but was associated with a shift in molecular pattern from Clara cell-associated towards the mucus-secreting cell differentiation pathway with multiple features of cancer-associated molecular phenotype.

Conclusions

These observations provide insights into the unique biology of human SAE by providing quantit-ative assessment of the global transcriptome under physiological conditions and in response to the stress of chronic cigarette smoking.

Gene expression signatures of radiation exposure in peripheral white blood cells of smokers and non-smokers

PURPOSE

The issue of potential confounding factors is critical to the development of any approach to radiation biodosimetry, and has not been fully addressed for gene expression-based approaches.

MATERIALS AND METHODS

As a step in this direction, we have investigated the effect of smoking on the global radiation gene expression response in ex vivo-irradiated peripheral blood cells using microarray analysis. We also evaluated the ability of gene expression signatures to predict the radiation exposure level of ex vivo-exposed samples from smokers and non-smokers of both genders.

RESULTS

We identified eight genes with a radiation response that was significantly affected by smoking status, and confirmed an effect of smoking on the radiation response of the four and a half LIM domains 2 (FHL2) gene using quantitative real-time polymerase chain reaction. The performance of our previously defined 74-gene signature in predicting the radiation dose to samples in this study was unaffected by differences in gender or smoking status, however, giving 98% correct prediction of dose category. This is the same accuracy as that found in the original study from which the signature was derived, using different donors.

CONCLUSION

The results support the development of peripheral blood gene expression as a viable strategy for radiation biodosimetry.

Multiwalled carbon nanotube-induced gene signatures in the mouse lung: potential predictive value for human lung cancer risk and prognosis

Concerns over the potential for multiwalled carbon nanotubes (MWCNT) to induce lung carcinogenesis have emerged. This study sought to (1) identify gene expression signatures in the mouse lungs following pharyngeal aspiration of well-dispersed MWCNT and (2) determine if these genes were associated with human lung cancer risk and progression. Genome-wide mRNA expression profiles were analyzed in mouse lungs (n = 160) exposed to 0, 10, 20, 40, or 80 μg of MWCNT by pharyngeal aspiration at 1, 7, 28, and 56 d postexposure. By using pairwise statistical analysis of microarray (SAM) and linear modeling, 24 genes were selected, which have significant changes in at least two time points, have a more than 1.5-fold change at all doses, and are significant in the linear model for the dose or the interaction of time and dose. Additionally, a 38-gene set was identified as related to cancer from 330 genes differentially expressed at d 56 postexposure in functional pathway analysis. Using the expression profiles of the cancer-related gene set in 8 mice at d 56 postexposure to 10 μg of MWCNT, a nearest centroid classification accurately predicts human lung cancer survival with a significant hazard ratio in training set (n = 256) and test set (n = 186). Furthermore, both gene signatures were associated with human lung cancer risk (n = 164) with significant odds ratios. These results may lead to development of a surveillance approach for early detection of lung cancer and prognosis associated with MWCNT in the workplace.

Characterizing the impact of smoking and lung cancer on the airway transcriptome using RNA-Seq

Cigarette smoke creates a molecular field of injury in epithelial cells that line the respiratory tract. We hypothesized that transcriptome sequencing (RNA-Seq) will enhance our understanding of the field of molecular injury in response to tobacco smoke exposure and lung cancer pathogenesis by identifying gene expression differences not interrogated or accurately measured by microarrays. We sequenced the high-molecular-weight fraction of total RNA (>200 nt) from pooled bronchial airway epithelial cell brushings (n = 3 patients per pool) obtained during bronchoscopy from healthy never smoker (NS) and current smoker (S) volunteers and smokers with (C) and without (NC) lung cancer undergoing lung nodule resection surgery. RNA-Seq libraries were prepared using 2 distinct approaches, one capable of capturing non-polyadenylated RNA (the prototype NuGEN Ovation RNA-Seq protocol) and the other designed to measure only polyadenylated RNA (the standard Illumina mRNA-Seq protocol) followed by sequencing generating approximately 29 million 36 nt reads per pool and approximately 22 million 75 nt paired-end reads per pool, respectively. The NuGEN protocol captured additional transcripts not detected by the Illumina protocol at the expense of reduced coverage of polyadenylated transcripts, while longer read lengths and a paired-end sequencing strategy significantly improved the number of reads that could be aligned to the genome. The aligned reads derived from the two complementary protocols were used to define the compendium of genes expressed in the airway epithelium (n = 20,573 genes). Pathways related to the metabolism of xenobiotics by cytochrome P450, retinol metabolism, and oxidoreductase activity were enriched among genes differentially expressed in smokers, whereas chemokine signaling pathways, cytokine-cytokine receptor interactions, and cell adhesion molecules were enriched among genes differentially expressed in smokers with lung cancer. There was a significant correlation between the RNA-Seq gene expression data and Affymetrix microarray data generated from the same samples (P < 0.001); however, the RNA-Seq data detected additional smoking- and cancer-related transcripts whose expression was were either not interrogated by or was not found to be significantly altered when using microarrays, including smoking-related changes in the inflammatory genes S100A8 and S100A9 and cancer-related changes in MUC5AC and secretoglobin (SCGB3A1). Quantitative real-time PCR confirmed differential expression of select genes and non-coding RNAs within individual samples. These results demonstrate that transcriptome sequencing has the potential to provide new insights into the biology of the airway field of injury associated with smoking and lung cancer. The measurement of both coding and non-coding transcripts by RNA-Seq has the potential to help elucidate mechanisms of response to tobacco smoke and to identify additional biomarkers of lung cancer risk and novel targets for chemoprevention.

Computational challenges of personal genomics

It is widely predicted that cost and efficiency gains in sequencing will usher in an era of personal genomics and personalized, predictive, preventive, and participatory medicine within a decade. I review the computational challenges ahead and propose general and specific directions for research and development. There is an urgent need to develop semantic ontologies that span genomics, molecular systems biology, and medical data. Although the development of such ontologies would be costly and difficult, the benefits will far outweigh the costs. I argue that availability of such ontologies would allow a revolution in web-services for personal genomics and medicine.

Ultrastructural, immunofluorescence, and RNA evidence support the hypothesis of a "new" virus associated with Kawasaki disease

BACKGROUND

Intracytoplasmic inclusion bodies (ICI) have been identified in ciliated bronchial epithelium of Kawasaki disease (KD) patients using a synthetic antibody derived from acute KD arterial IgA plasma cells; ICI may derive from the KD etiologic agent.

METHODS

Acute KD bronchial epithelium was subjected to immunofluorescence for ICI and cytokeratin, high-throughput sequencing, and transmission electron microscopy (TEM). Interferon pathway gene expression profiling was performed on KD lung.

RESULTS

An intermediate filament cytokeratin "cage" was not observed around KD ICI, making it unlikely that ICI are overproduced or misfolded human protein aggregates. Many interferon-stimulated genes were detected in the bronchial epithelium, and significant modulation of the interferon response pathway was observed in the lung tissue of KD patients. No known virus was identified by sequencing. Aggregates of virus-like particles (VLP) were detected by TEM in all 3 acute KD patients from whom nonembedded formalin-fixed lung tissue was available.

CONCLUSIONS

KD ICI are most likely virus induced; bronchial cells with ICI contain VLP that share morphologic features among several different RNA viral families. Expedited autopsies and tissue fixation from acute KD fatalities are urgently needed to more clearly ascertain the VLP. These findings are compatible with the hypothesis that the infectious etiologic agent of KD may be a "new" RNA virus.