Bronchial epithelial gene expression and interstitial lung abnormalities.

INTRODUCTION: Interstitial lung abnormalities (ILA) often represent early fibrotic changes that can portend a progressive fibrotic phenotype. In particular, the fibrotic subtype of ILA is associated with increased mortality and rapid decline in lung function. Understanding the differential gene expression that occurs in the lungs of participants with fibrotic ILA may provide insight into development of a useful biomarker for early detection and therapeutic targets for progressive pulmonary fibrosis. METHODS: Measures of ILA and gene expression data were available in 213 participants in the Detection of Early Lung Cancer Among Military Personnel (DECAMP1 and DECAMP2) cohorts. ILA was defined using Fleischner Society guidelines and determined by sequential reading of computed tomography (CT) scans. Primary analysis focused on comparing gene expression in ILA with usual interstitial pneumonia (UIP) pattern with those with no ILA. RESULTS: ILA was present in 51 (24%) participants, of which 16 (7%) were subtyped as ILA with a UIP pattern. One gene, pro platelet basic protein (PPBP) and seventeen pathways (e.g. TNF-α signalling) were significantly differentially expressed between those with a probable or definite UIP pattern of ILA compared to those without ILA. 16 of these 17 pathways, but no individual gene, met significance when comparing those with ILA to those without ILA. CONCLUSION: Our study demonstrates that abnormal inflammatory processes are apparent in the bronchial airway gene expression profiles of smokers with and without lung cancer with ILA. Future studies with larger and more diverse populations will be needed to confirm these findings.

Single-Cell Characterization of Pulmonary Nodules Implicates Suppression of Immunosurveillance across Early Stages of Lung Adenocarcinoma.

A greater understanding of molecular, cellular, and immunological changes during the early stages of lung adenocarcinoma development could improve diagnostic and therapeutic approaches in patients with pulmonary nodules at risk for lung cancer. To elucidate the immunopathogenesis of early lung tumorigenesis, we evaluated surgically resected pulmonary nodules representing the spectrum of early lung adenocarcinoma as well as associated normal lung tissues using single-cell RNA sequencing and validated the results by flow cytometry and multiplex immunofluorescence (MIF). Single-cell transcriptomics revealed a significant decrease in gene expression associated with cytolytic activities of tumor-infiltrating natural killer and natural killer T cells. This was accompanied by a reduction in effector T cells and an increase of CD4+ regulatory T cells (Treg) in subsolid nodules. An independent set of resected pulmonary nodules consisting of both adenocarcinomas and associated premalignant lesions corroborated the early increment of Tregs in premalignant lesions compared with the associated normal lung tissues by MIF. Gene expression analysis indicated that cancer-associated alveolar type 2 cells and fibroblasts may contribute to the deregulation of the extracellular matrix, potentially affecting immune infiltration in subsolid nodules through ligand-receptor interactions. These findings suggest that there is a suppression of immune surveillance across the spectrum of early-stage lung adenocarcinoma. SIGNIFICANCE: Analysis of a spectrum of subsolid pulmonary nodules by single-cell RNA sequencing provides insights into the immune regulation and cell-cell interactions in the tumor microenvironment during early lung tumor development.

Update on Biomarkers for the Stratification of Indeterminate Pulmonary Nodules.

Lung cancer is the leading cause of cancer-related deaths. Early detection and diagnosis are critical, as survival decreases with advanced stages. Approximately 1.6 million nodules are incidentally detected every year on chest CT scan images in the United States. This number of nodules identified is likely much larger after accounting for screening-detected nodules. Most of these nodules, whether incidentally or screening detected, are benign. Despite this, many patients undergo unnecessary invasive procedures to rule out cancer because our current stratification approaches are suboptimal, particularly for intermediate probability nodules. Thus, noninvasive strategies are urgently needed. Biomarkers have been developed to assist through the continuum of lung cancer care and include blood protein-based biomarkers, liquid biopsies, quantitative imaging analysis (radiomics), exhaled volatile organic compounds, and bronchial or nasal epithelium genomic classifiers, among others. Although many biomarkers have been developed, few have been integrated into clinical practice as they lack clinical utility studies showing improved patient-centered outcomes. Rapid technologic advances and large network collaborative efforts will continue to drive the discovery and validation of many novel biomarkers. Ultimately, however, randomized clinical utility studies showing improved patient outcomes will be required to bring biomarkers into clinical practice.

Convergence of YAP/TAZ, TEAD and TP63 activity is associated with bronchial premalignant severity and progression.

BACKGROUND: Bronchial premalignant lesions (PMLs) are composed primarily of cells resembling basal epithelial cells of the airways, which through poorly understood mechanisms have the potential to progress to lung squamous cell carcinoma (LUSC). Despite ongoing efforts that have mapped gene expression and cell diversity across bronchial PML pathologies, signaling and transcriptional events driving malignancy are poorly understood. Evidence has suggested key roles for the Hippo pathway effectors YAP and TAZ and associated TEAD and TP63 transcription factor families in bronchial basal cell biology and LUSC. In this study we examine the functional association of YAP/TAZ, TEADs and TP63 in bronchial epithelial cells and PMLs. METHODS: We performed RNA-seq in primary human bronchial epithelial cells following small interfering RNA (siRNA)-mediated depletion of YAP/TAZ, TEADs or TP63, and combined these data with ChIP-seq analysis of these factors. Directly activated or repressed genes were identified and overlapping genes were profiled across gene expression data obtained from progressive or regressive human PMLs and across lung single cell RNA-seq data sets. RESULTS: Analysis of genes regulated by YAP/TAZ, TEADs, and TP63 in human bronchial epithelial cells revealed a converged transcriptional network that is strongly associated with the pathological progression of bronchial PMLs. Our observations suggest that YAP/TAZ-TEAD-TP63 associate to cooperatively promote basal epithelial cell proliferation and repress signals associated with interferon responses and immune cell communication. Directly repressed targets we identified include the MHC Class II transactivator CIITA, which is repressed in progressive PMLs and associates with adaptive immune responses in the lung. Our findings provide molecular insight into the control of gene expression events driving PML progression, including those contributing to immune evasion, offering potential new avenues for lung cancer interception. CONCLUSIONS: Our study identifies important gene regulatory functions for YAP/TAZ-TEAD-TP63 in the early stages of lung cancer development, which notably includes immune-suppressive roles, and suggest that an assessment of the activity of this transcriptional complex may offer a means to identify immune evasive bronchial PMLs and serve as a potential therapeutic target.

Systematic dissection, preservation, and multiomics in whole human and bovine hearts.

OBJECTIVES: We sought to develop a rigorous, systematic protocol for the dissection and preservation of human hearts for biobanking that expands previous success in postmortem transcriptomics to multiomics from paired tissue. BACKGROUND: Existing cardiac biobanks consist largely of biopsy tissue or explanted hearts in select diseases and are insufficient for correlating whole organ phenotype with clinical data. METHODS: We demonstrate optimal conditions for multiomics interrogation (ribonucleic acid (RNA) sequencing, untargeted metabolomics) in hearts by evaluating the effect of technical variables (storage solution, temperature) and simulated postmortem interval (PMI) on RNA and metabolite stability. We used bovine (n=3) and human (n=2) hearts fixed in PAXgene or snap-frozen with liquid nitrogen. RESULTS: Using a paired Wald test, only two of the genes assessed were differentially expressed between left ventricular samples from bovine hearts stored in PAXgene at 0 and 12 hours PMI (FDR q<0.05). We obtained similar findings in human left ventricular samples, suggesting stability of RNA transcripts at PMIs up to 12 hours. Different library preparation methods (mRNA poly-A capture vs. rRNA depletion) resulted in similar quality metrics with both library preparations achieving >95% of reads properly aligning to the reference genomes across all PMIs for bovine and human hearts. PMI had no effect on RNA Integrity Number or quantity of RNA recovered at the time points evaluated. Of the metabolites identified (855 total) using untargeted metabolomics of human left ventricular tissue, 503 metabolites remained stable across PMIs (0, 4, 8, 12 hours). Most metabolic pathways retained several stable metabolites. CONCLUSIONS: Our data demonstrate a technically rigorous, reproducible protocol that will enhance cardiac biobanking practices and facilitate novel insights into human CVD. CONDENSED ABSTRACT: Cardiovascular disease (CVD) is the leading cause of mortality worldwide. Current biobanking practices insufficiently capture both the diverse array of phenotypes present in CVDs and the spatial heterogeneity across cardiac tissue sites. We have developed a rigorous and systematic protocol for the dissection and preservation of human cardiac biospecimens to enhance the availability of whole organ tissue for multiple applications. When combined with longitudinal clinical phenotyping, our protocol will enable multiomics in hearts to deepen our understanding of CVDs.

Bronchial gene expression alterations associated with radiological bronchiectasis.

OBJECTIVES: Discovering airway gene expression alterations associated with radiological bronchiectasis may improve the understanding of the pathobiology of early-stage bronchiectasis. METHODS: Presence of radiological bronchiectasis in 173 individuals without a clinical diagnosis of bronchiectasis was evaluated. Bronchial brushings from these individuals were transcriptomically profiled and analysed. Single-cell deconvolution was performed to estimate changes in cellular landscape that may be associated with early disease progression. RESULTS: 20 participants have widespread radiological bronchiectasis (three or more lobes). Transcriptomic analysis reflects biological processes associated with bronchiectasis including decreased expression of genes involved in cell adhesion and increased expression of genes involved in inflammatory pathways (655 genes, false discovery rate <0.1, log2 fold-change >0.25). Deconvolution analysis suggests that radiological bronchiectasis is associated with an increased proportion of ciliated and deuterosomal cells, and a decreased proportion of basal cells. Gene expression patterns separated participants into three clusters: normal, intermediate and bronchiectatic. The bronchiectatic cluster was enriched by participants with more lobes of radiological bronchiectasis (p<0.0001), more symptoms (p=0.002), higher SERPINA1 mutation rates (p=0.03) and higher computed tomography derived bronchiectasis scores (p<0.0001). CONCLUSIONS: Genes involved in cell adhesion, Wnt signalling, ciliogenesis and interferon-γ pathways had altered expression in the bronchus of participants with widespread radiological bronchiectasis, possibly associated with decreased basal and increased ciliated cells. This gene expression pattern is not only highly enriched among individuals with radiological bronchiectasis, but also associated with airway-related symptoms in those without discernible radiological bronchiectasis, suggesting that it reflects a bronchiectasis-associated, but non-bronchiectasis-specific lung pathophysiological process.

Inactivation of LATS1/2 drives luminal-basal plasticity to initiate basal-like mammary carcinomas.

Basal-like breast cancers, an aggressive breast cancer subtype that has poor treatment options, are thought to arise from luminal mammary epithelial cells that undergo basal plasticity through poorly understood mechanisms. Using genetic mouse models and ex vivo primary organoid cultures, we show that conditional co-deletion of the LATS1 and LATS2 kinases, key effectors of Hippo pathway signaling, in mature mammary luminal epithelial cells promotes the development of Krt14 and Sox9-expressing basal-like carcinomas that metastasize over time. Genetic co-deletion experiments revealed that phenotypes resulting from the loss of LATS1/2 activity are dependent on the transcriptional regulators YAP/TAZ. Gene expression analyses of LATS1/2-deleted mammary epithelial cells notably revealed a transcriptional program that associates with human basal-like breast cancers. Our study demonstrates in vivo roles for the LATS1/2 kinases in mammary epithelial homeostasis and luminal-basal fate control and implicates signaling networks induced upon the loss of LATS1/2 activity in the development of basal-like breast cancer.

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.

Bronchial gene expression signature associated with rate of subsequent FEV1 decline in individuals with and at risk of COPD.

BACKGROUND: COPD is characterised by progressive lung function decline. Leveraging prior work demonstrating bronchial airway COPD-associated gene expression alterations, we sought to determine if there are alterations associated with differences in the rate of FEV1 decline. METHODS: We examined gene expression among ever smokers with and without COPD who at baseline had bronchial brushings profiled by Affymetrix microarrays and had longitudinal lung function measurements (n=134; mean follow-up=6.38±2.48 years). Gene expression profiles associated with the rate of FEV1 decline were identified by linear modelling. RESULTS: Expression differences in 171 genes were associated with rate of FEV1 decline (false discovery rate <0.05). The FEV1 decline signature was replicated in an independent dataset of bronchial biopsies from patients with COPD (n=46; p=0.018; mean follow-up=6.76±1.32 years). Genes elevated in individuals with more rapid FEV1 decline are significantly enriched among the genes altered by modulation of XBP1 in two independent datasets (Gene Set Enrichment Analysis (GSEA) p<0.05) and are enriched in mucin-related genes (GSEA p<0.05). CONCLUSION: We have identified and replicated an airway gene expression signature associated with the rate of FEV1 decline. Aspects of this signature are related to increased expression of XBP1-regulated genes, a transcription factor involved in the unfolded protein response, and genes related to mucin production. Collectively, these data suggest that molecular processes related to the rate of FEV1 decline can be detected in airway epithelium, identify a possible indicator of FEV1 decline and make it possible to detect, in an early phase, ever smokers with and without COPD most at risk of rapid FEV1 decline.

Pathologic and gene expression comparison of CT- screen detected and routinely detected stage I/0 lung adenocarcinoma in NCCN risk-matched cohorts.

INTRODUCTION: Although three randomized control trials have proven mortality benefit of CT lung cancer screening (CTLS), <5% of eligible US smokers are screened. Some attribute this to fear of harm conveyed at shared decision visits, including the harm of overdiagnosis/overtreatment of indolent BAC-like adenocarcinoma. METHODS: Since the frequency of indolent cancers has not been compared between CTLS and routinely detected cohorts, we compare pathology and RNA expression of 86 NCCN high-risk CTLS subjects to 83 high-risk (HR-R) and 51 low-risk (LR-R) routinely detected patients. Indolent adenocarcinoma was defined as previously described for low malignant potential (LMP) adenocarcinoma along with AIS/MIA. Exome RNA sequencing was performed on a subset of high-risk (CTLS and HR-R) FFPE tumor samples. RESULTS: Indolent adenocarcinoma (AIS, MIA, and LMP) showed 100% disease-specific survival (DSS) with similar frequency in CTLS (18%) and HR-R (20%) which were comparatively lower than LR-R (33%). Despite this observation, CTLS exhibited intermediate DSS between HR-R and LR-R (5-year DSS: 88% CTLS, 82% HR-R, & 95% LR-R, p = 0.047), possibly reflecting a 0.4 cm smaller median tumor size and lower frequency of tumor necrosis compared to HR-R. WGCNA gene modules derived from TCGA lung adenocarcinoma correlated with aggressive histologic patterns, mitotic activity, and tumor invasive features, but no significant differential expression between CTLS and HR-R was observed. CONCLUSION: CTLS subjects are at no greater risk of overdiagnosis from indolent adenocarcinoma (AIS, MIA, and LMP) than risk-matched patients whose cancers are discovered in routine clinical practice. Improved outcomes likely reflect detection and treatment at smaller size.

Elevated T cell repertoire diversity is associated with progression of lung squamous cell premalignant lesions.

OBJECTIVE: The immune response to invasive carcinoma has been the focus of published work, but little is known about the adaptive immune response to bronchial premalignant lesions (PMLs), precursors of lung squamous cell carcinoma. This study was designed to characterize the T cell receptor (TCR) repertoire in PMLs and its association with clinical, pathological, and molecular features. METHODS: Endobronchial biopsies (n=295) and brushings (n=137) from high-risk subjects (n=50), undergoing lung cancer screening at approximately 1-year intervals via autofluorescence bronchoscopy and CT, were profiled by RNA-seq. We applied the TCR Repertoire Utilities for Solid Tissue/Tumor tool to the RNA-seq data to identify TCR CDR3 sequences across all samples. In the biopsies, we measured the correlation of TCR diversity with previously derived immune-associated PML transcriptional signatures and PML outcome. We also quantified the spatial and temporal distribution of shared and clonally expanded TCRs. Using the biopsies and brushes, the ratio of private (ie, found in one patient only) and public (ie, found in two or more patients) TCRs was quantified, and the CDR3 sequences were compared with those found in curated databases with known antigen specificities. RESULTS: We detected 39,303 unique TCR sequences across all samples. In PML biopsies, TCR diversity was negatively associated with a transcriptional signature of T cell mediated immune activation (p=4e-4) associated with PML outcome. Additionally, in lesions of the proliferative molecular subtype, TCR diversity was decreased in regressive versus progressive/persistent PMLs (p=0.045). Within each patient, TCRs were more likely to be shared between biopsies sampled at the same timepoint than biopsies sampled at the same anatomic location at different times. Clonally expanded TCRs, within a biopsied lesion, were more likely to be expanded at future time points than non-expanded clones. The majority of TCR sequences were found in a single sample, with only 3396 (8.6%) found in more than one sample and 1057 (2.7%) found in two or more patients (ie, public); however, when compared with a public database of CDR3 sequences, 4543 (11.6%) of TCRs were identified as public. TCRs with known antigen specificities were enriched among public TCRs (p<0.001). CONCLUSIONS: Decreased TCR diversity may reflect nascent immune responses that contribute to PML elimination. Further studies are needed to explore the potential for immunoprevention of PMLs.

Yap/Taz inhibit goblet cell fate to maintain lung epithelial homeostasis.

Proper lung function relies on the precise balance of specialized epithelial cells that coordinate to maintain homeostasis. Herein, we describe essential roles for the transcriptional regulators YAP/TAZ in maintaining lung epithelial homeostasis, reporting that conditional deletion of Yap and Wwtr1/Taz in the lung epithelium of adult mice results in severe defects, including alveolar disorganization and the development of airway mucin hypersecretion. Through in vivo lineage tracing and in vitro molecular experiments, we reveal that reduced YAP/TAZ activity promotes intrinsic goblet transdifferentiation of secretory airway epithelial cells. Global gene expression and chromatin immunoprecipitation sequencing (ChIP-seq) analyses suggest that YAP/TAZ act cooperatively with TEA domain (TEAD) transcription factors and the NuRD complex to suppress the goblet cell fate program, directly repressing the SPDEF gene. Collectively, our study identifies YAP/TAZ as critical factors in lung epithelial homeostasis and offers molecular insight into the mechanisms promoting goblet cell differentiation, which is a hallmark of many lung diseases.

The transition from normal lung anatomy to minimal and established fibrosis in idiopathic pulmonary fibrosis (IPF).

BACKGROUND: The transition from normal lung anatomy to minimal and established fibrosis is an important feature of the pathology of idiopathic pulmonary fibrosis (IPF). The purpose of this report is to examine the molecular and cellular mechanisms associated with this transition. METHODS: Pre-operative thoracic Multidetector Computed Tomography (MDCT) scans of patients with severe IPF (n = 9) were used to identify regions of minimal(n = 27) and established fibrosis(n = 27). MDCT, Micro-CT, quantitative histology, and next-generation sequencing were used to compare 24 samples from donor controls (n = 4) to minimal and established fibrosis samples. FINDINGS: The present results extended earlier reports about the transition from normal lung anatomy to minimal and established fibrosis by showing that there are activations of TGFBI, T cell co-stimulatory genes, and the down-regulation of inhibitory immune-checkpoint genes compared to controls. The expression patterns of these genes indicated activation of a field immune response, which is further supported by the increased infiltration of inflammatory immune cells dominated by lymphocytes that are capable of forming lymphoid follicles. Moreover, fibrosis pathways, mucin secretion, surfactant, TLRs, and cytokine storm-related genes also participate in the transitions from normal lung anatomy to minimal and established fibrosis. INTERPRETATION: The transition from normal lung anatomy to minimal and established fibrosis is associated with genes that are involved in the tissue repair processes, the activation of immune responses as well as the increased infiltration of CD4, CD8, B cell lymphocytes, and macrophages. These molecular and cellular events correlate with the development of structural abnormality of IPF and probably contribute to its pathogenesis.

Distinguishing Smoking-Related Lung Disease Phenotypes Via Imaging and Molecular Features.

BACKGROUND: Chronic tobacco smoke exposure results in a broad range of lung pathologies including emphysema, airway disease and parenchymal fibrosis as well as a multitude of extra-pulmonary comorbidities. Prior work using CT imaging has identified several clinically relevant subgroups of smoking related lung disease, but these investigations have generally lacked organ specific molecular correlates. RESEARCH QUESTION: Can CT imaging be used to identify clinical phenotypes of smoking related lung disease that have specific bronchial epithelial gene expression patterns to better understand disease pathogenesis? STUDY DESIGN AND METHODS: Using K-means clustering, we clustered participants from the COPDGene study (n = 5,273) based on CT imaging characteristics and then evaluated their clinical phenotypes. These clusters were replicated in the Detection of Early Lung Cancer Among Military Personnel (DECAMP) cohort (n = 360), and were further characterized using bronchial epithelial gene expression. RESULTS: Three clusters (preserved, interstitial predominant and emphysema predominant) were identified. Compared to the preserved cluster, the interstitial and emphysema clusters had worse lung function, exercise capacity and quality of life. In longitudinal follow-up, individuals from the emphysema group had greater declines in exercise capacity and lung function, more emphysema, more exacerbations, and higher mortality. Similarly, genes involved in inflammatory pathways (tumor necrosis factor-α, interferon-ß) are more highly expressed in bronchial epithelial cells from individuals in the emphysema cluster, while genes associated with T-cell related biology are decreased in these samples. Samples from individuals in the interstitial cluster generally had intermediate levels of expression of these genes. INTERPRETATION: Using quantitative CT imaging, we identified three groups of individuals in older ever-smokers that replicate in two cohorts. Airway gene expression differences between the three groups suggests increased levels of inflammation in the most severe clinical phenotype, possibly mediated by the tumor necrosis factor-α and interferon-ß pathways. CLINICAL TRIAL REGISTRATION: COPDGene (NCT00608764), DECAMP-1 (NCT01785342), DECAMP-2 (NCT02504697).

Improving lung cancer risk stratification leveraging whole transcriptome RNA sequencing and machine learning across multiple cohorts.

BACKGROUND: Bronchoscopy for suspected lung cancer has low diagnostic sensitivity, rendering many inconclusive results. The Bronchial Genomic Classifier (BGC) was developed to help with patient management by identifying those with low risk of lung cancer when bronchoscopy is inconclusive. The BGC was trained and validated on patients in the Airway Epithelial Gene Expression in the Diagnosis of Lung Cancer (AEGIS) trials. A modern patient cohort, the BGC Registry, showed differences in key clinical factors from the AEGIS cohorts, with less smoking history, smaller nodules and older age. Additionally, we discovered interfering factors (inhaled medication and sample collection timing) that impacted gene expressions and potentially disguised genomic cancer signals. METHODS: In this study, we leveraged multiple cohorts and next generation sequencing technology to develop a robust Genomic Sequencing Classifier (GSC). To address demographic composition shift and interfering factors, we synergized three algorithmic strategies: 1) ensemble of clinical dominant and genomic dominant models; 2) development of hierarchical regression models where the main effects from clinical variables were regressed out prior to the genomic impact being fitted in the model; and 3) targeted placement of genomic and clinical interaction terms to stabilize the effect of interfering factors. The final GSC model uses 1232 genes and four clinical covariates - age, pack-years, inhaled medication use, and specimen collection timing. RESULTS: In the validation set (N = 412), the GSC down-classified low and intermediate pre-test risk subjects to very low and low post-test risk with a specificity of 45% (95% CI 37-53%) and a sensitivity of 91% (95%CI 81-97%), resulting in a negative predictive value of 95% (95% CI 89-98%). Twelve percent of intermediate pre-test risk subjects were up-classified to high post-test risk with a positive predictive value of 65% (95%CI 44-82%), and 27% of high pre-test risk subjects were up-classified to very high post-test risk with a positive predictive value of 91% (95% CI 78-97%). CONCLUSIONS: The GSC overcame the impact of interfering factors and achieved consistent performance across multiple cohorts. It demonstrated diagnostic accuracy in both down- and up-classification of cancer risk, providing physicians actionable information for many patients with inconclusive bronchoscopy.

Characterizing smoking-induced transcriptional heterogeneity in the human bronchial epithelium at single-cell resolution.

The human bronchial epithelium is composed of multiple distinct cell types that cooperate to defend against environmental insults. While studies have shown that smoking alters bronchial epithelial function and morphology, its precise effects on specific cell types and overall tissue composition are unclear. We used single-cell RNA sequencing to profile bronchial epithelial cells from six never and six current smokers. Unsupervised analyses led to the characterization of a set of toxin metabolism genes that localized to smoker ciliated cells, tissue remodeling associated with a loss of club cells and extensive goblet cell hyperplasia, and a previously unidentified peri-goblet epithelial subpopulation in smokers who expressed a marker of bronchial premalignant lesions. Our data demonstrate that smoke exposure drives a complex landscape of cellular alterations that may prime the human bronchial epithelium for disease.

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.

Tobacco-Related Alterations in Airway Gene Expression are Rapidly Reversed Within Weeks Following Smoking-Cessation.

The physiologic response to tobacco smoke can be measured by gene-expression profiling of the airway epithelium. Temporal resolution of kinetics of gene-expression alterations upon smoking-cessation might delineate distinct biological processes that are activated during recovery from tobacco smoke exposure. Using whole genome gene-expression profiling of individuals initiating a smoking-cessation attempt, we sought to characterize the kinetics of gene-expression alterations in response to short-term smoking-cessation in the nasal epithelium. RNA was extracted from the nasal epithelial of active smokers at baseline and at 4, 8, 16, and 24-weeks after smoking-cessation and put onto Gene ST arrays. Gene-expression levels of 119 genes were associated with smoking-cessation (FDR < 0.05, FC ≥1.7) with a majority of the changes occurring by 8-weeks and a subset changing by 4-weeks. Genes down-regulated by 4- and 8-weeks post-smoking-cessation were involved in xenobiotic metabolism and anti-apoptotic functions respectively. These genes were enriched among genes previously found to be induced in smokers and following short-term in vitro exposure of airway epithelial cells to cigarette smoke (FDR < 0.05). Our findings suggest that the nasal epithelium can serve as a minimally-invasive tool to measure the reversible impact of smoking and broadly, may serve to assess the physiological impact of changes in smoking behavior.

Molecular subtyping reveals immune alterations associated with progression of bronchial premalignant lesions.

Bronchial premalignant lesions (PMLs) are precursors of lung squamous cell carcinoma, but have variable outcome, and we lack tools to identify and treat PMLs at risk for progression to cancer. Here we report the identification of four molecular subtypes of PMLs with distinct differences in epithelial and immune processes based on RNA-Seq profiling of endobronchial biopsies from high-risk smokers. The Proliferative subtype is enriched with bronchial dysplasia and exhibits up-regulation of metabolic and cell cycle pathways. A Proliferative subtype-associated gene signature identifies subjects with Proliferative PMLs from normal-appearing uninvolved large airway brushings with high specificity. In progressive/persistent Proliferative lesions expression of interferon signaling and antigen processing/presentation pathways decrease and immunofluorescence indicates a depletion of innate and adaptive immune cells compared with regressive lesions. Molecular biomarkers measured in PMLs or the uninvolved airway can enhance histopathological grading and suggest immunoprevention strategies for intercepting the progression of PMLs to lung cancer.

Effect of long-term corticosteroid treatment on microRNA and gene-expression profiles in COPD.

The aim was to investigate whether microRNA (miRNA) expression is modulated by inhaled corticosteroid (ICS) treatmentWe performed genome-wide miRNA analysis on bronchial biopsies of 69 moderate/severe chronic obstructive pulmonary disease (COPD) patients at baseline and after 6- and 30-month treatment with the ICS fluticasone propionate or placebo. The effect of ICS on miRNA expression was validated in differentiated primary bronchial epithelial cultures, and functional studies were conducted in BEAS-2B cells. MiRNAs affected by ICS and their predicted targets were compared to an independent miRNA dataset of bronchial brushings from COPD patients and healthy controls.Treatment with ICS for both 6 and 30 months significantly altered the expression of four miRNAs, including miR-320d, which was increased during ICS treatment compared with placebo. The ICS-induced increase of miR-320d was confirmed in primary airway epithelial cells. MiR-320d negatively correlated targets were enriched for pro-inflammatory genes and were increased in the bronchial brushes of patients with lower lung function in the independent dataset. Overexpression of miR-320d in BEAS-2B cells dampened cigarette smoke extract-induced pro-inflammatory activity via inhibition of nuclear factor-κB.Collectively, we identified miR-320d as a novel mediator of ICS, regulating the pro-inflammatory response of the airway epithelium.