The evolution of lung adenocarcinoma precursors is associated with chromosomal instability and transition from innate to adaptive immune response/evasion.

Studying lung adenocarcinoma (LUAD) early carcinogenesis is challenging, primarily due to the lack of LUAD precursors specimens. We amassed multi-omics data from 213 LUAD and LUAD precursors to identify molecular features underlying LUAD precancer evolution. We observed progressively increasing mutations, chromosomal aberrations, whole genome doubling and genomic instability from precancer to invasive LUAD, indicating aggravating chromosomal instability (CIN). Telomere shortening, a crucial genomic alteration linked to CIN, emerged at precancer stage. Moreover, later-stage lesions demonstrated increasing cancer stemness and decreasing alveolar identity, suggesting epithelial de-differentiation during early LUAD carcinogenesis. The innate immune cells progressively diminished from precancer to invasive LUAD, concomitant with a gradual recruitment of adaptive immune cells (except CD8+ and gamma-delta T cells that decreased in later stages) and upregulation of numerous immune checkpoints, suggesting LUAD precancer evolution is associated with a shift from innate to adaptive immune response and immune evasion mediated by various mechanisms.

An atlas of epithelial cell states and plasticity in lung adenocarcinoma.

Understanding the cellular processes that underlie early lung adenocarcinoma (LUAD) development is needed to devise intervention strategies1. Here we studied 246,102 single epithelial cells from 16 early-stage LUADs and 47 matched normal lung samples. Epithelial cells comprised diverse normal and cancer cell states, and diversity among cancer cells was strongly linked to LUAD-specific oncogenic drivers. KRAS mutant cancer cells showed distinct transcriptional features, reduced differentiation and low levels of aneuploidy. Non-malignant areas surrounding human LUAD samples were enriched with alveolar intermediate cells that displayed elevated KRT8 expression (termed KRT8+ alveolar intermediate cells (KACs) here), reduced differentiation, increased plasticity and driver KRAS mutations. Expression profiles of KACs were enriched in lung precancer cells and in LUAD cells and signified poor survival. In mice exposed to tobacco carcinogen, KACs emerged before lung tumours and persisted for months after cessation of carcinogen exposure. Moreover, they acquired Kras mutations and conveyed sensitivity to targeted KRAS inhibition in KAC-enriched organoids derived from alveolar type 2 (AT2) cells. Last, lineage-labelling of AT2 cells or KRT8+ cells following carcinogen exposure showed that KACs are possible intermediates in AT2-to-tumour cell transformation. This study provides new insights into epithelial cell states at the root of LUAD development, and such states could harbour potential targets for prevention or intervention.

A Nasal Swab Classifier to Evaluate the Probability of Lung Cancer in Patients With Pulmonary Nodules.

BACKGROUND: Accurate assessment of the probability of lung cancer (pCA) is critical in patients with pulmonary nodules (PNs) to help guide decision-making. We sought to validate a clinical-genomic classifier developed using whole-transcriptome sequencing of nasal epithelial cells from patients with a PN ≤ 30 mm who smoke or have previously smoked. RESEARCH QUESTION: Can the pCA in individuals with a PN and a history of smoking be predicted by a classifier that uses clinical factors and genomic data from nasal epithelial cells obtained by cytologic brushing? STUDY DESIGN AND METHODS: Machine learning was used to train a classifier using genomic and clinical features on 1,120 patients with PNs labeled as benign or malignant established by a final diagnosis or a minimum of 12 months of radiographic surveillance. The classifier was designed to yield low-, intermediate-, and high-risk categories. The classifier was validated in an independent set of 312 patients, including 63 patients with a prior history of cancer (other than lung cancer), comparing the classifier prediction with the known clinical outcome. RESULTS: In the primary validation set, sensitivity and specificity for low-risk classification were 96% and 42%, whereas sensitivity and specificity for high-risk classification was 58% and 90%, respectively. Sensitivity was similar across stages of non-small cell lung cancer, independent of subtype. Performance compared favorably with clinical-only risk models. Analysis of 63 patients with prior cancer showed similar performance as did subanalyses of patients with light vs heavy smoking burden and those eligible for lung cancer screening vs those who were not. INTERPRETATION: The nasal classifier provides an accurate assessment of pCA in individuals with a PN ≤ 30 mm who smoke or have previously smoked. Classifier-guided decision-making could lead to fewer diagnostic procedures in patients without cancer and more timely treatment in patients with lung cancer.

Integrative genetic and genomic networks identify microRNA associated with COPD and ILD.

Chronic obstructive pulmonary disease (COPD) and interstitial lung disease (ILD) are clinically and molecularly heterogeneous diseases. We utilized clustering and integrative network analyses to elucidate roles for microRNAs (miRNAs) and miRNA isoforms (isomiRs) in COPD and ILD pathogenesis. Short RNA sequencing was performed on 351 lung tissue samples of COPD (n = 145), ILD (n = 144) and controls (n = 64). Five distinct subclusters of samples were identified including 1 COPD-predominant cluster and 2 ILD-predominant clusters which associated with different clinical measurements of disease severity. Utilizing 262 samples with gene expression and SNP microarrays, we built disease-specific genetic and expression networks to predict key miRNA regulators of gene expression. Members of miR-449/34 family, known to promote airway differentiation by repressing the Notch pathway, were among the top connected miRNAs in both COPD and ILD networks. Genes associated with miR-449/34 members in the disease networks were enriched among genes that increase in expression with airway differentiation at an air-liquid interface. A highly expressed isomiR containing a novel seed sequence was identified at the miR-34c-5p locus. 47% of the anticorrelated predicted targets for this isomiR were distinct from the canonical seed sequence for miR-34c-5p. Overexpression of the canonical miR-34c-5p and the miR-34c-5p isomiR with an alternative seed sequence down-regulated NOTCH1 and NOTCH4. However, only overexpression of the isomiR down-regulated genes involved in Ras signaling such as CRKL and GRB2. Overall, these findings elucidate molecular heterogeneity inherent across COPD and ILD patients and further suggest roles for miR-34c in regulating disease-associated gene-expression.

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.

Single-molecule genome-wide mutation profiles of cell-free DNA for non-invasive detection of cancer.

Somatic mutations are a hallmark of tumorigenesis and may be useful for non-invasive diagnosis of cancer. We analyzed whole-genome sequencing data from 2,511 individuals in the Pan-Cancer Analysis of Whole Genomes (PCAWG) study as well as 489 individuals from four prospective cohorts and found distinct regional mutation type-specific frequencies in tissue and cell-free DNA from patients with cancer that were associated with replication timing and other chromatin features. A machine-learning model using genome-wide mutational profiles combined with other features and followed by CT imaging detected >90% of patients with lung cancer, including those with stage I and II disease. The fixed model was validated in an independent cohort, detected patients with cancer earlier than standard approaches and could be used to monitor response to therapy. This approach lays the groundwork for non-invasive cancer detection using genome-wide mutation features that may facilitate cancer screening and monitoring.

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.

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.

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.

The Single-Cell Immunogenomic Landscape of B and Plasma Cells in Early-Stage Lung Adenocarcinoma.

Tumor-infiltrating B and plasma cells (TIB) are prevalent in lung adenocarcinoma (LUAD); however, they are poorly characterized. We performed paired single-cell RNA and B-cell receptor (BCR) sequencing of 16 early-stage LUADs and 47 matching multiregion normal tissues. By integrative analysis of ∼50,000 TIBs, we define 12 TIB subsets in the LUAD and adjacent normal ecosystems and demonstrate extensive remodeling of TIBs in LUADs. Memory B cells and plasma cells (PC) were highly enriched in tumor tissues with more differentiated states and increased frequencies of somatic hypermutation. Smokers exhibited markedly elevated PCs and PCs with distinct differentiation trajectories. BCR clonotype diversity increased but clonality decreased in LUADs, smokers, and with increasing pathologic stage. TIBs were mostly localized within CXCL13+ lymphoid aggregates, and immune cell sources of CXCL13 production evolved with LUAD progression and included elevated fractions of CD4 regulatory T cells. This study provides a spatial landscape of TIBs in early-stage LUAD. SIGNIFICANCE: While TIBs are highly enriched in LUADs, they are poorly characterized. This study provides a much-needed understanding of the transcriptional, clonotypic states and phenotypes of TIBs, unraveling their potential roles in the immunopathology of early-stage LUADs and constituting a road map for the development of TIB-targeted immunotherapies for the treatment of this morbid malignancy. This article is highlighted in the In This Issue feature, p. 2483.

Clinical validation and utility of Percepta GSC for the evaluation of lung cancer.

The Percepta Genomic Sequencing Classifier (GSC) was developed to up-classify as well as down-classify the risk of malignancy for lung lesions when bronchoscopy is non-diagnostic. We evaluated the performance of Percepta GSC in risk re-classification of indeterminate lung lesions. This multicenter study included individuals who currently or formerly smoked undergoing bronchoscopy for suspected lung cancer from the AEGIS I/ II cohorts and the Percepta Registry. The classifier was measured in normal-appearing bronchial epithelium from bronchial brushings. The sensitivity, specificity, and predictive values were calculated using predefined thresholds. The ability of the classifier to decrease unnecessary invasive procedures was estimated. A set of 412 patients were included in the validation (prevalence of malignancy was 39.6%). Overall, 29% of intermediate-risk lung lesions were down-classified to low-risk with a 91.0% negative predictive value (NPV) and 12.2% of intermediate-risk lesions were up-classified to high-risk with a 65.4% positive predictive value (PPV). In addition, 54.5% of low-risk lesions were down-classified to very low risk with >99% NPV and 27.3% of high-risk lesions were up-classified to very high risk with a 91.5% PPV. If the classifier results were used in nodule management, 50% of patients with benign lesions and 29% of patients with malignant lesions undergoing additional invasive procedures could have avoided these procedures. The Percepta GSC is highly accurate as both a rule-out and rule-in test. This high accuracy of risk re-classification may lead to improved management of lung lesions.

Clinical Study of Aspirin and Zileuton on Biomarkers of Tobacco-Related Carcinogenesis in Current Smokers.

The chemopreventive effect of aspirin and other non-steroidal anti-inflammatory drugs (NSAIDs) on lung cancer risk is supported by epidemiologic and preclinical studies. Zileuton, a 5-lipoxygenase inhibitor, has additive activity with NSAIDs against tobacco carcinogenesis in preclinical models. We hypothesized that cyclooxygenase plus 5-lipoxygenase inhibition would be more effective than a placebo in modulating the nasal epithelium gene signatures of tobacco exposure and lung cancer. We conducted a randomized, double-blinded study of low-dose aspirin plus zileuton vs. double placebo in current smokers to compare the modulating effects on nasal gene expression and arachidonic acid metabolism. In total, 63 participants took aspirin 81 mg daily plus zileuton (Zyflo CR) 600 mg BID or the placebo for 12 weeks. Nasal brushes from the baseline, end-of-intervention, and one-week post intervention were profiled via microarray. Aspirin plus zilueton had minimal effects on the modulation of the nasal or bronchial gene expression signatures of smoking, lung cancer, and COPD but favorably modulated a bronchial gene expression signature of squamous dysplasia. Aspirin plus zileuton suppressed urinary leukotriene but not prostaglandin E2, suggesting shunting through the cyclooxygenase pathway when combined with 5-lipoxygenase inhibition. Continued investigation of leukotriene inhibitors is needed to confirm these findings, understand the long-term effects on the airway epithelium, and identify the safest, optimally dosed agents.

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.

The molecular and cellular mechanisms associated with the destruction of terminal bronchioles in COPD.

RATIONALE: Peripheral airway obstruction is a key feature of chronic obstructive pulmonary disease (COPD), but the mechanisms of airway loss are unknown. This study aims to identify the molecular and cellular mechanisms associated with peripheral airway obstruction in COPD. METHODS: Ten explanted lung specimens donated by patients with very severe COPD treated by lung transplantation and five unused donor control lungs were sampled using systematic uniform random sampling (SURS), resulting in 240 samples. These samples were further examined by micro-computed tomography (CT), quantitative histology and gene expression profiling. RESULTS: Micro-CT analysis showed that the loss of terminal bronchioles in COPD occurs in regions of microscopic emphysematous destruction with an average airspace size of ≥500 and <1000 µm, which we have termed a "hot spot". Based on microarray gene expression profiling, the hot spot was associated with an 11-gene signature, with upregulation of pro-inflammatory genes and downregulation of inhibitory immune checkpoint genes, indicating immune response activation. Results from both quantitative histology and the bioinformatics computational tool CIBERSORT, which predicts the percentage of immune cells in tissues from transcriptomic data, showed that the hot spot regions were associated with increased infiltration of CD4 and CD8 T-cell and B-cell lymphocytes. INTERPRETATION: The reduction in terminal bronchioles observed in lungs from patients with COPD occurs in a hot spot of microscopic emphysema, where there is upregulation of IFNG signalling, co-stimulatory immune checkpoint genes and genes related to the inflammasome pathway, and increased infiltration of immune cells. These could be potential targets for therapeutic interventions in COPD.

Smoking Modulates Different Secretory Subpopulations Expressing SARS-CoV-2 Entry Genes in the Nasal and Bronchial Airways.

Background : SARS-CoV-2 infection and disease severity are influenced by viral entry (VE) gene expression patterns in airway epithelium. The similarities and differences of VE gene expression (ACE2, TMPRSS2, and CTSL) across nasal and bronchial compartments has not been fully characterized using matched samples from large cohorts. Results : Gene expression data from 793 nasal and 1,673 bronchial brushes obtained from individuals participating in lung cancer screening or diagnostic workup revealed that smoking was the only clinical factor significantly and reproducibly associated with VE gene expression. ACE2 and TMPRSS2 expression were 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 RNA-seq confirmed that smoking increased the abundance of several secretory cell populations in the bronchus, but only goblet cells in the nose. Conclusions : 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. Smoking-induced VE gene expression changes in the nose likely has minimal impact on SARS-CoV-2 infection, but in the bronchus, smoking may lead to higher viral loads and more severe disease.

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.

Integrated Biomarkers for the Management of Indeterminate Pulmonary Nodules.

Rationale: Patients with indeterminate pulmonary nodules (IPNs) at risk of cancer undergo high rates of invasive, costly, and morbid procedures. Objectives: To train and externally validate a risk prediction model that combined clinical, blood, and imaging biomarkers to improve the noninvasive management of IPNs. Methods: In this prospectively collected, retrospective blinded evaluation study, probability of cancer was calculated for 456 patient nodules using the Mayo Clinic model, and patients were categorized into low-, intermediate-, and high-risk groups. A combined biomarker model (CBM) including clinical variables, serum high sensitivity CYFRA 21-1 level, and a radiomic signature was trained in cohort 1 (n = 170) and validated in cohorts 2-4 (total n = 286). All patients were pooled to recalibrate the model for clinical implementation. The clinical utility of the CBM compared with current clinical care was evaluated in 2 cohorts. Measurements and Main Results: The CBM provided improved diagnostic accuracy over the Mayo Clinic model with an improvement in area under the curve of 0.124 (95% bootstrap confidence interval, 0.091-0.156; P < 2 × 10-16). Applying 10% and 70% risk thresholds resulted in a bias-corrected clinical reclassification index for cases and control subjects of 0.15 and 0.12, respectively. A clinical utility analysis of patient medical records estimated that a CBM-guided strategy would have reduced invasive procedures from 62.9% to 50.6% in the intermediate-risk benign population and shortened the median time to diagnosis of cancer from 60 to 21 days in intermediate-risk cancers. Conclusions: Integration of clinical, blood, and image biomarkers improves noninvasive diagnosis of patients with IPNs, potentially reducing the rate of unnecessary invasive procedures while shortening the time to diagnosis.

Resolving the Spatial and Cellular Architecture of Lung Adenocarcinoma by Multiregion Single-Cell Sequencing.

Little is known of the geospatial architecture of individual cell populations in lung adenocarcinoma (LUAD) evolution. Here, we perform single-cell RNA sequencing of 186,916 cells from five early-stage LUADs and 14 multiregion normal lung tissues of defined spatial proximities from the tumors. We show that cellular lineages, states, and transcriptomic features geospatially evolve across normal regions to LUADs. LUADs also exhibit pronounced intratumor cell heterogeneity within single sites and transcriptional lineage-plasticity programs. T regulatory cell phenotypes are increased in normal tissues with proximity to LUAD, in contrast to diminished signatures and fractions of cytotoxic CD8+ T cells, antigen-presenting macrophages, and inflammatory dendritic cells. We further find that the LUAD ligand-receptor interactome harbors increased expression of epithelial CD24, which mediates protumor phenotypes. These data provide a spatial atlas of LUAD evolution, and a resource for identification of targets for its treatment. SIGNIFICANCE: The geospatial ecosystem of the peripheral lung and early-stage LUAD is not known. Our multiregion single-cell sequencing analyses unravel cell populations, states, and phenotypes in the spatial and ecologic evolution of LUAD from the lung that comprise high-potential targets for early interception.This article is highlighted in the In This Issue feature, p. 2355.