Multi-region exome sequencing reveals genomic evolution from preneoplasia to lung adenocarcinoma

Tumor cell heterogeneity drives spatial organization of the intratumoral immune response

Intratumoral heterogeneity (ITH)-defined as genetic and cellular diversity within a tumor-is linked to failure of immunotherapy and an inferior anti-tumor immune response. We modeled heterogeneous tumors comprised of "hot" and "cold" tumor populations (giving rise to T cell-rich and T cell-poor tumors, respectively) and introduced fluorescent labels to enable precise spatial tracking. We found the cold tumor cell population exerted a "dominant cold" effect in mixed tumors. Strikingly, spatial analysis revealed that the tumor cells themselves created distinct local microenvironments within heterogeneous tumors: regions occupied by cold tumor cells showed pronounced immunosuppression, harboring increased CD206Hi macrophages and diminished local T cell function. This inferior T cell activity in cold regions persisted even after immunotherapy and mechanistically was mediated by CX3CL1 produced by the cold tumor cells. An immune cold tumor population within a heterogeneous tumor thus impairs tumor immunity on both a tumor-wide and a highly localized spatial scale.

Spatial and multiomics analysis of human and mouse lung adenocarcinoma precursors reveals TIM-3 as a putative target for precancer interception

How tumor microenvironment shapes lung adenocarcinoma (LUAD) precancer evolution remains poorly understood. Spatial immune profiling of 114 human LUAD and LUAD precursors reveals a progressive increase of adaptive response and a relative decrease of innate immune response as LUAD precursors progress. The immune evasion features align the immune response patterns at various stages. TIM-3-high features are enriched in LUAD precancers, which decrease in later stages. Furthermore, single-cell RNA sequencing (scRNA-seq) and spatial immune and transcriptomics profiling of LUAD and LUAD precursor specimens from 5 mouse models validate high TIM-3 features in LUAD precancers. In vivo TIM-3 blockade at precancer stage, but not at advanced cancer stage, decreases tumor burden. Anti-TIM-3 treatment is associated with enhanced antigen presentation, T cell activation, and increased M1/M2 macrophage ratio. These results highlight the coordination of innate and adaptive immune response/evasion during LUAD precancer evolution and suggest TIM-3 as a potential target for LUAD precancer interception.

Translating premalignant biology to accelerate non-small-cell lung cancer interception

Over the past decade, substantial progress has been made in the development of targeted and immune-based therapies for patients with advanced non-small-cell lung cancer. To further improve outcomes for patients with lung cancer, identifying and intercepting disease at the earliest and most curable stages are crucial next steps. With the recent implementation of low-dose computed tomography scan screening in populations at high risk, there is an emerging unmet need for new diagnostic, prognostic and therapeutic tools to help treat patients suspected of harbouring premalignant lesions and minimally invasive non-small-cell lung cancer. Continued advances in the identification of the earliest drivers of lung carcinogenesis are poised to address these unmet needs. Employing multimodal approaches to chart the temporal and spatial maps of the molecular events driving lung premalignant lesion progression will refine our understanding of early carcinogenesis. Elucidating the molecular drivers of premalignancy is critical to the development of biomarkers to detect those incubating a premalignant lesion, to stratify risk for progression to invasive cancer and to identify novel therapeutic targets to intercept that process. In this Review, we summarize emerging insights into the earliest cellular and molecular events associated with lung squamous and adenocarcinoma carcinogenesis and highlight the growing opportunity for translating these insights into clinical tools for early detection and disease interception to transform the outcomes for those at risk for lung cancer.

The Current Roadmap of Lung Cancer Biology, Genomics and Racial Disparity

Globally, lung cancer is the most prevalent cause of cancer-related death. There are two large histological groups of lung cancer: small-cell lung cancer (SCLC) and non-small-cell lung cancer (NSCLC). Based on histopathological and molecular features, adenocarcinoma (ADC) and squamous cell carcinoma (SCC) are the two major histologic subtypes of NSCLC. Various epidemiological and environmental factors are linked with an increased risk of lung cancer. However, these risk factors show disparities in patients with divergent racial and ethnic backgrounds. Interestingly, different populations were found to harbor distinct molecular features as evidenced by variations in genetic mutation profiles. Moreover, diverse histological and molecular progression patterns are identified in lung cancer, which could be crucial in improving diagnosis, prognosis, and therapeutic planning. In concert with a plethora of nuclear genetic alterations, mitochondrial alteration, epigenetic reprogramming, microbial dysbiosis, and immune alteration signatures have been identified in various lung cancer types. This review article provides a comprehensive overview of screening tests and the treatment strategies for NSCLC and SCLC, including surgery, radiation therapy, chemotherapy, targeted therapies, and immunotherapies. Through the unification of these diverse aspects, this review article aspires to a complete understanding of lung cancer's genomics, biology, microbial landscapes, and racial disparity and seeks to understand the essential role of racial and ethnic factors in lung cancer occurrence and treatment.

Study on the regulatory mechanism of luteolin inhibiting WDR72 on the proliferation and metastasis of non small cell lung cancer

Non-small cell lung cancer (NSCLC) is a major cause of cancer-related mortality worldwide. Understanding molecular mechanisms and identifying potential therapeutic targets are crucial for improving treatment outcomes. This study aims to explore the effect of luteolin on NSCLC progression by regulating WDR72 and to investigate the related molecular mechanisms using cellular and animal models. The study employed a comprehensive set of experiments to evaluate the impact of luteolin and WDR72 on NSCLC cell proliferation and metastasis. Techniques included the CCK- 8 assay, colony formation assay, scratch test, and Transwell assay. Molecular docking experiments were performed to validate the binding interaction between luteolin and WDR72. Experimental groups included OE-WDR72, OE-WDR72 + Luteolin, Control, Control + Luteolin, and sh-WDR72. Western blot analysis was used to examine protein expression related to apoptosis, epithelial-mesenchymal transition (EMT), AKT signaling, and other markers. Additionally, a nude mouse subcutaneous tumor model was established to assess the in vivo tumor-forming ability of NSCLC cells under different treatments. Luteolin significantly inhibited the proliferation, invasion, and migration of NSCLC cell lines (H1299 and A549) and reduced tumor formation in nude mice. Molecular docking demonstrated strong binding affinity between luteolin and WDR72. Overexpression of WDR72 promoted NSCLC cell proliferation and migration, while WDR72 silencing showed the opposite effects. Western blot analysis revealed that WDR72 overexpression increased phosphorylated AKT and Bcl- 2 levels while decreasing caspase- 3. In contrast, silencing WDR72 reduced these protein levels. Luteolin treatment in WDR72-overexpressing cells resulted in decreased phosphorylated AKT, increased apoptosis, and suppressed EMT. Tumor transplantation experiments indicated that tumors in the OE-WDR72 group exhibited the fastest growth, while the sh-WDR72 group showed the slowest growth. Luteolin treatment significantly reduced WDR72 expression, suggesting a regulatory role in NSCLC progression. Luteolin effectively inhibits EMT, invasion, and migration of NSCLC cells by modulating WDR72. WDR72 plays a pivotal role in stimulating the proliferation and metastasis of NSCLC cells. By downregulating WDR72, luteolin suppresses NSCLC progression, potentially through modulation of the PI3 K/AKT/EMT signaling pathway. These findings highlight luteolin as a promising therapeutic agent for NSCLC treatment.

Current insights into the clinico-pathologic characteristics of lung cancer in women

INTRODUCTION

Lung cancer is responsible for premature cancer deaths in women and is the first cause of cancer deaths in women in many countries. The problem of lung cancer in women seems to be underestimated in many aspects, including low participation in clinical trials and screening tests.

AREAS COVERED

Current research progress has contributed to a better understanding of the issue and makes it possible to describe the problem in a new light. In our paper, the problem of lung cancer in women was discussed in a broad aspect, taking into account women's health, the harmful effects of smoking and the current diagnostic and treatment process. The results of treatment also differ in relation to sex. All these aspects of the diversity of women's lung cancer were presented on the basis of newest and most comprehensive literature.

EXPERT OPINION

Lung cancer in women is and will remain an important health problem worldwide, which is justified by epidemiological data, basic research and treatment results.

Prognostic implications, genomic and immune characteristics of lung adenocarcinoma with lepidic growth pattern

AIMS

Conflicting data were provided regarding the prognostic impact and genomic features of lung adenocarcinoma (LUAD) with lepidic growth pattern (LP+A). Delineation of the genomic and immune characteristics of LP+A could provide deeper insights into its prognostic implications and treatment determination.

METHODS

We conducted a search of articles in PubMed, EMBASE and the Cochrane Library from inception to January 2024. A domestic cohort consisting of 52 LUAD samples was subjected to whole-exome sequencing as internal validation. Data from The Cancer Genomic Atlas and the Gene Expression Omnibus datasets were obtained to characterise the genomic and immune profiles of LP+A. Pooled HRs and rates were calculated.

RESULTS

The pooled results indicated that lepidic growth pattern was either predominant (0.35, 95% CI 0.22 to 0.56, p<0.01) or minor (HR 0.50, 95% CI 0.36 to 0.70, p<0.01) histological subtype was associated with favourable disease-free survival. Pooled gene mutation rates suggested higher EGFR mutation (0.55, 95% CI 0.46 to 0.64, p<0.01) and lower KRAS mutation (0.14, 95% CI 0.02 to 0.25, p=0.02) in lepidic-predominant LUAD. Lepidic-predominant LUAD had lower tumour mutation burden and pooled positive rate of PD-L1 expression compared with other subtypes. LP+A was characterised by abundance in resting CD4+memory T cells, monocytes and γδ T cells, as well as scarcity of cancer-associated fibroblasts.

CONCLUSIONS

LP+A was a unique histological subtype with a higher EGFR mutation rate, lower tumour mutation burden and immune checkpoint expression levels. Our findings suggested potential benefits from targeted therapy over immunotherapy in LP+A.

Proteogenomic characterization reveals tumorigenesis and progression of lung cancer manifested as subsolid nodules

Lung adenocarcinoma (LUAD) radiologically displayed as subsolid nodules (SSNs) is prevalent. Nevertheless, the precise clinical management of SSNs necessitates a profound understanding of their tumorigenesis and progression. Here, we analyze 66 LUAD displayed as SSNs covering 3 histological stages including adenocarcinoma in situ (AIS), minimally invasive adenocarcinoma (MIA) and invasive adenocarcinoma (IAC) by incorporating genomics, proteomics, phosphoproteomics and glycoproteomics. Intriguingly, cholesterol metabolism is aberrantly regulated in the preneoplastic AIS stage. Importantly, target ablation of proprotein convertase subtilisin/kexin type 9 (PCSK9) promotes the initiation of LUAD. Furthermore, sustained endoplasmic reticulum stress is demonstrated to be a hallmark and a reliable biomarker of AIS progression to IAC. Consistently, target promotion of ER stress profoundly retards LUAD progression. Our study provides comprehensive proteogenomic landscape of SSNs, sheds lights on the tumorigenesis and progression of SSNs and suggests preventive and therapeutic strategies for LUAD.

CENPF (+) cancer cells promote malignant progression of early-stage TP53 mutant lung adenocarcinoma

The prevention and precise treatment of early-stage lung adenocarcinoma (LUAD) characterized by small nodules (stage IA) remains a significant challenge for clinicians, which is due largely to the limited understanding of the oncogenic mechanisms spanning from preneoplasia to invasive adenocarcinoma. Our study highlights the pivotal role of cancer cells exhibiting high expression of centromere protein F (CENPF), driven by TP53 mutations, which become increasingly prevalent during the transition from preneoplasia to invasive LUAD. Biologically, cancer cells (CENPF+) exhibited robust proliferative and stem-like capabilities, thereby propelling the malignant progression of early-stage LUAD. Clinically, autoantibodies against CENPF in the serum and elevated cancer cells (CENPF+) in tissue correlated positively with the progression of early-stage LUAD, especially those in stage IA. Our findings suggest that cancer cells (CENPF+) play a central role in orchestrating the malignant evolution of LUAD and hold potential as a novel biomarker for early-stage detection and management of the disease.

Multiomics Analysis Reveals Molecular Changes during Early Progression of Precancerous Lesions to Lung Adenocarcinoma in Never-Smokers

Lung cancer is the most common cause of cancer mortality globally, and the prevalence of lung adenocarcinoma, the most common lung cancer subtype, has increased sharply in East Asia. Early diagnosis leads to better survival rates, but this requires an improved understanding of the molecular changes during early tumorigenesis, particularly in nonsmokers. In this study, we performed whole-exome sequencing and RNA sequencing of samples from 94 East Asian patients with precancerous lesions [25 with atypical adenomatous hyperplasia (AAH); 69 with adenocarcinoma in situ (AIS)] and 73 patients with early invasive lesions [minimally invasive adenocarcinoma (MIA)]. Cellular analysis revealed that the activities of endothelial and stromal cells could be used to categorize tumors into molecular subtypes within pathologically defined types of lesions. The subtypes were linked with the radiologically defined type of lesions and corresponded to immune cell infiltration throughout the early progression of lung adenocarcinoma. Spatial transcriptomic analysis revealed the distribution of epithelial cells, endothelial cells, fibroblasts, and plasma cells within MIA samples. Characterization of the molecular lesion subtypes identified positively selected mutational patterns and suggested that angiogenesis in the late-stage AIS type potentially contributes to tissue invasion of the MIA type. This study offers a resource that may help improve early diagnosis and patient prognosis, and the findings suggest possible approaches for early disease interception. Significance: Integrative analysis of multiomics data revealed coordination between immune and nonimmune cells during early progression of precancerous lesions to lung adenocarcinomas and shed light on the molecular characteristics of clinically defined subtypes.

Highly proliferating cancer cells function as novel prognostic biomarkers for lung adenocarcinoma with particular usefulness for stage IA risk stratification

BACKGROUND

The refinement of risk stratification in lung adenocarcinoma (LUAD) plays a pivotal role in advancing precision medicine; however, the current staging classification falls short of comprehensiveness, particularly in the case of stage IA patients. We aimed to molecularly stratify LUAD patients especially for stage IA.

METHODS

We analysed tumour heterogeneity and identified highly proliferating cancer cells (HPCs) in LUAD by performing single-cell RNA sequencing (scRNA-seq) analysis, immunohistochemical (IHC) staining using a tissue microarray, flow cytometry and biological experiments. Then, we quantified the content of HPCs in nine LUAD datasets by single-sample gene set enrichment analysis and evaluated the relationship between the percentage of HPCs and overall survival (OS). Next, we analysed the OS predictive effect of HPCs at different LUAD stages, especially for stage I risk stratification. Furthermore, we established a prognostic prediction model based on HPC-associated genes for clinical application. The above findings were validated in another five LUAD datasets. Finally, we explored the relationship between HPCs and the progressive pathological evolution of early-stage LUAD and the driving mutations by scRNA-seq, bulk RNA-seq and IHC staining.

RESULTS

LUAD tissues carry a small proportion of HPCs, which show potential for malignant proliferation and intense interactions with the microenvironment. A high HPC content is an independent risk factor for OS in LUAD patients, even in stage IA patients. HPCs can be used to establish a cut-off point for the prognosis of stage IA disease, with patients with a higher risk showing a prognosis similar to that of patients with stage IB disease. We built an R package (HSurADs) based on HPC-associated genes, which exhibited good efficacy for the prognostic prediction of LUAD. HPCs gradually increase with the pathological evolution of early-stage LUAD, which may be affected by TP53 mutations.

CONCLUSION

The HPC content can be used as a novel prognostic factor for LUAD, especially for stage IA risk stratification.

Acquisition of discrete immune suppressive barriers contributes to the initiation and progression of preinvasive to invasive human lung cancer

Computerized chest tomography (CT)-guided screening in populations at risk for lung cancer has increased the detection of preinvasive subsolid nodules, which progress to solid invasive adenocarcinoma. Despite the clinical significance, there is a lack of effective therapies for intercepting the progression of preinvasive to invasive adenocarcinoma. To uncover determinants of early disease emergence and progression, we used integrated single-cell approaches, including scRNA-seq, multiplexed imaging mass cytometry and spatial transcriptomics, to construct the first high-resolution map of the composition, lineage/functional states, developmental trajectories and multicellular crosstalk networks from microdissected non-solid (preinvasive) and solid compartments (invasive) of individual part-solid nodules. We found that early disease initiation and subsequent progression are associated with the evolution of immune-suppressive cellular phenotypes characterized by decreased cytotoxic CD8 T and NK cells, increased T cell exhaustion and accumulation of immunosuppressive regulatory T cells (Tregs) and M2-like macrophages expressing TREM2. Within Tregs, we identified a unique population of 4-1BB+ Treg subset enriched for the IL2-STAT5 suppressive pathway with transcription profiles supporting discrete metabolic alterations. Spatial analysis showed increased density of suppressive immune cells around tumor cells, increased exhaustion phenotype of both CD4 and CD8 T cells expressing chemokine CXCL13, and spatial microcomplex of endothelial and lymphocyte interactions within tertiary lymphoid structures. The single-cell architecture identifies determinants of early disease emergence and progression, which may be developed not only as diagnostic/prognostic biomarkers but also as targets for disease interception. Additionally, our dataset constitutes a valuable resource for the preinvasive lung cancer research community.

Genomic Landscape Features of Minimally Invasive Adenocarcinoma and Invasive Lung Adenocarcinoma

Background  The widespread implementation of computed tomography has significantly increased the detection of small pulmonary nodules, including atypical adenomatous hyperplasia, minimally invasive adenocarcinoma (MIA), and invasive adenocarcinoma (IAC). Few studies have focused on the genomic differences between MIA and IAC. Methods  We retrospectively analyzed patients with lung adenocarcinoma (LUAD) who underwent surgery from January 2020 to December 2023. Patients were categorized into MIA and IAC groups. The mutation status of common driver genes was assessed using next-generation sequencing. Results  A total of 422 LUAD patients were included in the study, comprising 119 MIA cases and 303 IAC cases. MIA patients were younger and predominantly female compared with IAC patients. EGFR mutations were detected in 251 patients (59.5%), with the frequency of EGFR mutations increasing from 37.0% in MIA to 68.3% in IAC ( p  < 0.001). TP53 mutations were found in 108 patients (25.6%), with 7 patients (5.9%) in MIA and 101 patients (33.3%) in IAC ( p  < 0.001). ERBB2 mutations were identified in 23 MIA patients (19.3%) and 20 IAC patients (6.6%) ( p  < 0.001). Additionally, CDKN2A mutations were detected in 23 IAC patients (7.6%), while no mutations in this gene were found in the MIA group. Moreover, ALK and RET gene fusions were identified in 11 patients, respectively. Conclusion  ERBB2 mutations and RET fusions are early genomic events in LUAD, while TP53 and CDKN2A mutations and ALK fusions occur later. Genomic intratumor heterogeneity likely arises early, before invasive characteristics develop.

Applications of single-cell analysis in immunotherapy for lung cancer: Current progress, new challenges and expectations

BACKGROUND

Lung cancer is a prevalent form of cancer worldwide, presenting a substantial risk to human well-being. Lung cancer is classified into two main types: non-small cell lung cancer (NSCLC) and small cell lung cancer (SCLC). The advancement of tumor immunotherapy, specifically immune checkpoint inhibitors and adaptive T-cell therapy, has encountered substantial obstacles due to the rapid progression of SCLC and the metastasis, recurrence, and drug resistance of NSCLC. These challenges are believed to stem from the tumor heterogeneity of lung cancer within the tumor microenvironment.

AIM OF REVIEW

This review aims to comprehensively explore recent strides in single-cell analysis, a robust sequencing technology, concerning its application in the realm of tumor immunotherapy for lung cancer. It has been effectively integrated with transcriptomics, epigenomics, genomics, and proteomics for various applications. Specifically, these techniques have proven valuable in mapping the transcriptional activity of tumor-infiltrating lymphocytes in patients with NSCLC, identifying circulating tumor cells, and elucidating the heterogeneity of the tumor microenvironment.

KEY SCIENTIFIC CONCEPTS OF REVIEW

The review emphasizes the paramount significance of single-cell analysis in mapping the immune cells within NSCLC patients, unveiling circulating tumor cells, and elucidating the tumor microenvironment heterogeneity. Notably, these advancements highlight the potential of single-cell analysis to revolutionize lung cancer immunotherapy by characterizing immune cell fates, improving therapeutic strategies, and identifying promising targets or prognostic biomarkers. It is potential to unravel the complexities within the tumor microenvironment and enhance treatment strategies marks a significant step towards more effective therapies and improved patient outcomes.

Super multiple primary lung cancers harbor high-frequency BRAF and low-frequency EGFR mutations in the MAPK pathway

The incidence of multiple primary lung cancer (MPLC) is increasing, with some of our surgical patients exhibiting numerous lesions. We defined lung cancer with five or more primary lesions as super MPLCs. Elucidating the genomic characteristics of this special MPLC subtype can help reduce disease burden and understand tumor evolution. In our cohort of synchronous super early-stage MPLCs (PUMCH-ssesMPLC), whole-exome sequencing on 130 resected malignant specimens from 18 patients provided comprehensive super-MPLC genomic landscapes. Mutations are enriched in PI3k-Akt and MAPK pathways. Their BRAF mutation frequency (31.5%) is significantly higher than MPLC with fewer lesions and early-stage single-lesion cancer, while EGFR mutations are significantly fewer (13.8%). As lesion counts increase, BRAF mutations gradually become dominant. Also, invasive lesions more tend to have classic super-MPLC mutation patterns. High-frequency BRAF mutations, especially Class II, and low-frequency EGFR mutations could be a reason for the limited effectiveness of targeted therapy in super-MPLC patients.

Multiregional transcriptomic profiling provides improved prognostic insight in localized non-small cell lung cancer

Lung Cancer remains the leading cause of cancer deaths in the USA and worldwide. Non-small cell lung cancer (NSCLC) harbors high transcriptomic intratumor heterogeneity (RNA-ITH) that limits the reproducibility of expression-based prognostic models. In this study, we used multiregional RNA-seq data (880 tumor samples from 350 individuals) from both public (TRACERx) and internal (MDAMPLC) cohorts to investigate the effect of RNA-ITH on prognosis in localized NSCLC at the gene, signature, and tumor microenvironment levels. At the gene level, the maximal expression of hazardous genes (expression negatively associated with survival) but the minimal expression of protective genes (expression positively associated with survival) across different regions within a tumor were more prognostic than the average expression. Following that, we examined whether multiregional expression profiling can improve the performance of prognostic signatures. We investigated 11 gene signatures collected from previous publications and one signature developed in this study. For all of them, the prognostic prediction accuracy can be significantly improved by converting the regional expression of signature genes into sample-specific expression with a simple function-taking the maximal expression of hazardous genes and the minimal expression of protective genes. In the tumor microenvironment, we found a similar rule also seems applicable to immune ITH. We calculated the infiltration levels of major immune cell types in each region of a sample based on expression deconvolution. Prognostic analysis indicated that the region with the lowest infiltration level of protective or highest infiltration level of hazardous immune cells determined the prognosis of NSCLC patients. Our study highlighted the impact of RNA-ITH on the prognostication of NSCLC, which should be taken into consideration to optimize the design and application of expression-based prognostic biomarkers and models. Multiregional assays have the great potential to significantly improve their applications to prognostic stratification.

Non-small-cell lung cancer

Non-small-cell lung cancer (NSCLC) is one of the most frequent cancer types and is responsible for the majority of cancer-related deaths worldwide. The management of NSCLC has improved considerably, especially in the past 10 years. The systematic screening of populations at risk with low-dose CT, the implementation of novel surgical and radiotherapeutic techniques and a deeper biological understanding of NSCLC that has led to innovative systemic treatment options have improved the prognosis of patients with NSCLC. In non-metastatic NSCLC, the combination of various perioperative strategies and adjuvant immunotherapy in locally advanced disease seem to enhance cure rates. In metastatic NSCLC, the implementation of novel drugs might prolong disease control together with preserving quality of life. The further development of predictive clinical and genetic markers will be essential for the next steps in individualized treatment concepts.

Genomic and immune heterogeneity of multiple synchronous lung adenocarcinoma at different developmental stages

Multiple synchronous lung cancers (MSLCs) constitute a unique subtype of lung cancer. To explore the genomic and immune heterogeneity across different pathological stages of MSLCs, we analyse 16 MSLCs from 8 patients using single-cell RNA-seq, single-cell TCR sequencing, and bulk whole-exome sequencing. Our investigation indicates clonally independent tumours with convergent evolution driven by shared driver mutations. However, tumours from the same individual exhibit few shared mutations, indicating independent origins. During the transition from pre-invasive to invasive adenocarcinoma, we observe a shift in T cell phenotypes characterized by increased Treg cells and exhausted CD8+ T cells, accompanied by diminished cytotoxicity. Additionally, invasive adenocarcinomas exhibit greater neoantigen abundance and a more diverse TCR repertoire, indicating heightened heterogeneity. In summary, despite having a common genetic background and environmental exposure, our study emphasizes the individuality of MSLCs at different stages, highlighting their unique genomic and immune characteristics.

Premalignant Progression in the Lung: Knowledge Gaps and Novel Opportunities for Interception of Non-Small Cell Lung Cancer. An Official American Thoracic Society Research Statement

Rationale: Despite significant advances in precision treatments and immunotherapy, lung cancer is the most common cause of cancer death worldwide. To reduce incidence and improve survival rates, a deeper understanding of lung premalignancy and the multistep process of tumorigenesis is essential, allowing timely and effective intervention before cancer development. Objectives: To summarize existing information, identify knowledge gaps, formulate research questions, prioritize potential research topics, and propose strategies for future investigations into the premalignant progression in the lung. Methods: An international multidisciplinary team of basic, translational, and clinical scientists reviewed available data to develop and refine research questions pertaining to the transformation of premalignant lung lesions to advanced lung cancer. Results: This research statement identifies significant gaps in knowledge and proposes potential research questions aimed at expanding our understanding of the mechanisms underlying the progression of premalignant lung lesions to lung cancer in an effort to explore potential innovative modalities to intercept lung cancer at its nascent stages. Conclusions: The identified gaps in knowledge about the biological mechanisms of premalignant progression in the lung, together with ongoing challenges in screening, detection, and early intervention, highlight the critical need to prioritize research in this domain. Such focused investigations are essential to devise effective preventive strategies that may ultimately decrease lung cancer incidence and improve patient outcomes.

Multi-omics analysis unveils immunosuppressive microenvironment in the occurrence and development of multiple pulmonary lung cancers

Multiple pulmonary lung cancers (MPLCs) are frequently encountered on computed tomography (CT) scanning of chest, yet their intrinsic characteristics associated with genomic features and radiological or pathological textures that may lead to distinct clinical outcomes remain largely unexplored. A total of 27 pulmonary nodules covering different radiological or pathological textures as well as matched adjacent normal tissues and blood samples were collected from patients diagnosed with MPLCs. Whole-exome sequencing (WES) and whole-transcriptome sequencing were performed. The molecular and immune features of MPLCs associated with distinct radiological or pathological textures were comprehensively investigated. Genomics analysis unveiled the distinct branches of pulmonary nodules originating independently within the same individual. EGFR and KRAS mutations were found to be prevalent in MPLCs, exhibiting mutual exclusivity. The group with KRAS mutations exhibited stronger immune signatures compared to the group with EGFR mutations. Additionally, MPLCs exhibited a pronounced immunosuppressive microenvironment, which was particularly distinct when compared with normal tissues. The expression of the FDSCP gene was specifically observed in MPLCs. When categorizing MPLCs based on radiological or pathological characteristics, a progressive increase in mutation accumulation was observed, accompanied by heightened chromatin-level instability as ground-glass opacity component declined or invasive progression occurred. A close association with the immunosuppressive microenvironment was also observed during the progression of pulmonary nodules. Notably, the upregulation of B cell and regulatory T cell marker genes occurred progressively. Immune cell abundance analysis further demonstrated a marked increase in exhausted cells and regulatory T cells during the progression of pulmonary nodules. These results were further validated by independent datasets including nCounter RNA profiling, single-cell RNA sequencing, and spatial transcriptomic datasets. Our study provided a comprehensive representation of the diverse landscape of MPLCs originating within the same individual and emphasized the significant influence of the immunosuppressive microenvironment in the occurrence and development of pulmonary nodules. These findings hold great potential for enhancing the clinical diagnosis and treatment strategies for MPLCs.

Single-cell spatial multiomics reveals tumor microenvironment vulnerabilities in cancer resistance to immunotherapy

Heterogeneous resistance to immunotherapy remains a major challenge in cancer treatment, often leading to disease progression and death. Using CITE-seq and matched 40-plex PhenoCycler tissue imaging, we performed longitudinal multimodal single-cell analysis of tumors from metastatic melanoma patients with innate resistance, acquired resistance, or response to immunotherapy. We established the multimodal integration toolkit to align transcriptomic features, cellular epitopes, and spatial information to provide deeper insights into the tumors. With longitudinal analysis, we identified an "immune-striving" tumor microenvironment marked by peri-tumor lymphoid aggregates and low infiltration of T cells in the tumor and the emergence of MITF+SPARCL1+ and CENPF+ melanoma subclones after therapy. The enrichment of B cell-associated signatures in the molecular composition of lymphoid aggregates was associated with better survival. These findings provide further insights into the establishment of microenvironmental cell interactions and molecular composition of spatial structures that could inform therapeutic intervention.

Integrating machine learning and single-cell analysis to uncover lung adenocarcinoma progression and prognostic biomarkers

The progression of lung adenocarcinoma (LUAD) from atypical adenomatous hyperplasia (AAH) to invasive adenocarcinoma (IAC) involves a complex evolution of tumour cell clusters, the mechanisms of which remain largely unknown. By integrating single-cell datasets and using inferCNV, we identified and analysed tumour cell clusters to explore their heterogeneity and changes in abundance throughout LUAD progression. We applied gene set variation analysis (GSVA), pseudotime analysis, scMetabolism, and Cytotrace scores to study biological functions, metabolic profiles and stemness traits. A predictive model for prognosis, based on key cluster marker genes, was developed using CoxBoost and plsRcox (CPM), and validated across multiple cohorts for its prognostic prediction capabilities, tumour microenvironment characterization, mutation landscape and immunotherapy response. We identified nine distinct tumour cell clusters, with Cluster 6 indicating an early developmental stage, high stemness and proliferative potential. The abundance of Clusters 0 and 6 increased from AAH to IAC, correlating with prognosis. The CPM model effectively distinguished prognosis in immunotherapy cohorts and predicted genomic alterations, chemotherapy drug sensitivity, and immunotherapy responsiveness. Key gene S100A16 in the CPM model was validated as an oncogene, enhancing LUAD cell proliferation, invasion and migration. The CPM model emerges as a novel biomarker for predicting prognosis and immunotherapy response in LUAD patients, with S100A16 identified as a potential therapeutic target.

Tumor initiation and early tumorigenesis: molecular mechanisms and interventional targets

Tumorigenesis is a multistep process, with oncogenic mutations in a normal cell conferring clonal advantage as the initial event. However, despite pervasive somatic mutations and clonal expansion in normal tissues, their transformation into cancer remains a rare event, indicating the presence of additional driver events for progression to an irreversible, highly heterogeneous, and invasive lesion. Recently, researchers are emphasizing the mechanisms of environmental tumor risk factors and epigenetic alterations that are profoundly influencing early clonal expansion and malignant evolution, independently of inducing mutations. Additionally, clonal evolution in tumorigenesis reflects a multifaceted interplay between cell-intrinsic identities and various cell-extrinsic factors that exert selective pressures to either restrain uncontrolled proliferation or allow specific clones to progress into tumors. However, the mechanisms by which driver events induce both intrinsic cellular competency and remodel environmental stress to facilitate malignant transformation are not fully understood. In this review, we summarize the genetic, epigenetic, and external driver events, and their effects on the co-evolution of the transformed cells and their ecosystem during tumor initiation and early malignant evolution. A deeper understanding of the earliest molecular events holds promise for translational applications, predicting individuals at high-risk of tumor and developing strategies to intercept malignant transformation.

Integrative genomic and transcriptomic profiling of pulmonary sarcomatoid carcinoma identifies molecular subtypes associated with distinct immune features and clinical outcomes

Background

Pulmonary sarcomatoid carcinoma (PSC) is a rare and aggressive subtype of non-small cell lung cancer (NSCLC), characterized by the presence of epithelial and sarcoma-like components. The molecular and immune landscape of PSC has not been well defined.

Methods

Multiomics profiling of 21 pairs of PSCs with matched normal lung tissues was performed through targeted high-depth DNA panel, whole-exome, and RNA sequencing. We describe molecular and immune features that define subgroups of PSC with disparate genomic and immunogenic features as well as distinct clinical outcomes.

Results

In total, 27 canonical cancer gene mutations were identified, with TP53 the most frequently mutated gene, followed by KRAS. Interestingly, most TP53 and KRAS mutations were earlier genomic events mapped to the trunks of the tumors, suggesting branching evolution in most PSC tumors. We identified two distinct molecular subtypes of PSC, driven primarily by immune infiltration and signaling. The Immune High (IM-H) subtype was associated with superior survival, highlighting the impact of immune infiltration on the biological and clinical features of localized PSCs.

Conclusions

We provided detailed insight into the mutational landscape of PSC and identified two molecular subtypes associated with prognosis. IM-H tumors were associated with favorable recurrence-free survival and overall survival, highlighting the importance of tumor immune infiltration in the biological and clinical features of PSCs.

Multidimensional biological characteristics of ground glass nodules

The detection rate of ground glass nodules (GGNs) has increased in recent years because of their malignant potential but relatively indolent biological behavior; thus, correct GGN recognition and management has become a research focus. Many scholars have explored the underlying mechanism of the indolent progression of GGNs from several perspectives, such as pathological type, genomic mutational characteristics, and immune microenvironment. GGNs have different major mutated genes at different stages of development; EGFR mutation is the most common mutation in GGNs, and p53 mutation is the most abundant mutation in the invasive stage of GGNs. Pure GGNs have fewer genomic alterations and a simpler genomic profile and exhibit a gradually evolving genomic mutation profile as the pathology progresses. Compared to advanced lung adenocarcinoma, GGN lung adenocarcinoma has a higher immune cell percentage, is under immune surveillance, and has less immune escape. However, as the pathological progression and solid component increase, negative immune regulation and immune escape increase gradually, and a suppressive immune environment is established gradually. Currently, regular computer tomography monitoring and surgery are the main treatment strategies for persistent GGNs. Stereotactic body radiotherapy and radiofrequency ablation are two local therapeutic alternatives, and systemic therapy has been progressively studied for lung cancer with GGNs. In the present review, we discuss the characterization of the multidimensional molecular evolution of GGNs that could facilitate more precise differentiation of such highly heterogeneous lesions, laying a foundation for the development of more effective individualized treatment plans.

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.

Immunologic tumor microenvironment modulators for turning cold tumors hot

Tumors can be classified into distinct immunophenotypes based on the presence and arrangement of cytotoxic immune cells within the tumor microenvironment (TME). Hot tumors, characterized by heightened immune activity and responsiveness to immune checkpoint inhibitors (ICIs), stand in stark contrast to cold tumors, which lack immune infiltration and remain resistant to therapy. To overcome immune evasion mechanisms employed by tumor cells, novel immunologic modulators have emerged, particularly ICIs targeting cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) and programmed cell death protein 1/programmed death-ligand 1(PD-1/PD-L1). These agents disrupt inhibitory signals and reactivate the immune system, transforming cold tumors into hot ones and promoting effective antitumor responses. However, challenges persist, including primary resistance to immunotherapy, autoimmune side effects, and tumor response heterogeneity. Addressing these challenges requires innovative strategies, deeper mechanistic insights, and a combination of immune interventions to enhance the effectiveness of immunotherapies. In the landscape of cancer medicine, where immune cold tumors represent a formidable hurdle, understanding the TME and harnessing its potential to reprogram the immune response is paramount. This review sheds light on current advancements and future directions in the quest for more effective and safer cancer treatment strategies, offering hope for patients with immune-resistant tumors.

Clonal expansion of shared T cell receptors reveals the existence of immune commonality among different lesions of synchronous multiple primary lung cancer

The increase in the detection rate of synchronous multiple primary lung cancer (MPLC) has posed remarkable clinical challenges due to the limited understanding of its pathogenesis and molecular features. Here, comprehensive comparisons of genomic and immunologic features between MPLC and solitary lung cancer nodule (SN), as well as different lesions of the same patient, were performed. Compared with SN, MPLC displayed a lower rate of EGFR mutation but higher rates of BRAF, MAP2K1, and MTOR mutation, which function exactly in the upstream and downstream of the same signaling pathway. Considerable heterogeneity in T cell receptor (TCR) repertoire exists among not only different patients but also among different lesions of the same patient. Invasive lesions of MPLC exhibited significantly higher TCR diversity and lower TCR expansion than those of SN. Intriguingly, different lesions of the same patient always shared a certain proportion of TCR clonotypes. Significant clonal expansion could be observed in shared TCR clonotypes, particularly in those existing in all lesions of the same patient. In conclusion, this study provided evidences of the distinctive mutational landscape, activation of oncogenic signaling pathways, and TCR repertoire in MPLC as compared with SN. The significant clonal expansion of shared TCR clonotypes demonstrated the existence of immune commonality among different lesions of the same patient and shed new light on the individually tailored precision therapy for MPLC.

Multi-omics Analysis Reveals Immune Features Associated with Immunotherapy Benefit in Patients with Squamous Cell Lung Cancer from Phase III Lung-MAP S1400I Trial

PURPOSE

Identifying molecular and immune features to guide immune checkpoint inhibitor (ICI)-based regimens remains an unmet clinical need.

EXPERIMENTAL DESIGN

Tissue and longitudinal blood specimens from phase III trial S1400I in patients with metastatic squamous non-small cell carcinoma (SqNSCLC) treated with nivolumab monotherapy (nivo) or nivolumab plus ipilimumab (nivo+ipi) were subjected to multi-omics analyses including multiplex immunofluorescence (mIF), nCounter PanCancer Immune Profiling Panel, whole-exome sequencing, and Olink.

RESULTS

Higher immune scores from immune gene expression profiling or immune cell infiltration by mIF were associated with response to ICIs and improved survival, except regulatory T cells, which were associated with worse overall survival (OS) for patients receiving nivo+ipi. Immune cell density and closer proximity of CD8+GZB+ T cells to malignant cells were associated with superior progression-free survival and OS. The cold immune landscape of NSCLC was associated with a higher level of chromosomal copy-number variation (CNV) burden. Patients with LRP1B-mutant tumors had a shorter survival than patients with LRP1B-wild-type tumors. Olink assays revealed soluble proteins such as LAMP3 increased in responders while IL6 and CXCL13 increased in nonresponders. Upregulation of serum CXCL13, MMP12, CSF-1, and IL8 were associated with worse survival before radiologic progression.

CONCLUSIONS

The frequency, distribution, and clustering of immune cells relative to malignant ones can impact ICI efficacy in patients with SqNSCLC. High CNV burden may contribute to the cold immune microenvironment. Soluble inflammation/immune-related proteins in the blood have the potential to monitor therapeutic benefit from ICI treatment in patients with SqNSCLC.

Advancements in the diagnosis and treatment of sub‑centimeter lung cancer in the era of precision medicine (Review)

Lung cancer is the malignancy with the highest global mortality rate and imposes a substantial burden on society. The increasing popularity of lung cancer screening has led to increasing number of patients being diagnosed with pulmonary nodules due to their potential for malignancy, causing considerable distress in the affected population. However, the diagnosis and treatment of sub-centimeter grade pulmonary nodules remain controversial. The evolution of genetic detection technology and the development of targeted drugs have positioned the diagnosis and treatment of lung cancer in the precision medicine era, leading to a marked improvement in the survival rate of patients with lung cancer. It has been established that lung cancer driver genes serve a key role in the development and progression of sub-centimeter lung cancer. The present review aimed to consolidate the findings on genes associated with sub-centimeter lung cancer, with the intent of serving as a reference for future studies and the personalized management of sub-centimeter lung cancer through genetic testing.

Targeted next-generation sequencing of 491 lung cancers in clinical practice: Implications for future detection strategy and targeted therapy

Although lung cancer remains the most common cause of global cancer-related mortality, the identification of oncogenic driver alterations and the development of targeted drugs has dramatically altered the therapeutic landscape. In this retrospective study, we found that 97.7% samples carried at least one mutation in the 25 genes tested in our cohort. 53.6% samples were positive for EGFR mutations, followed by TP53 (41.1%), KRAS (11.8%), ERBB2 (4.3%). EGFR mutations were mainly found in female adenocarcinomas, while TP53 was mainly found in male non-adenocarcinomas. Significant differences can be found in the mutation rate of EGFR (60.9% vs 11.9%), KRAS (12.2% vs 25.0%), STK11 (1.5% vs 11.9%), FGFR3 (2.4% vs 0.0%) and ERBB4 (1.2% vs 6.1%) between adenocarcinoma in our cohort and TCGA-LUAD data (all p < 0.001). What's more, we found that the mutation of EGFR increased significantly from adenocarcinomas in situ (AIS, 21.4%) to microinvasive adenocarcinomas (MIA, 52.4%) and invasive adenocarcinomas (IA, 61.1%), while the mutation of ERBB2 dropped markedly from AIS (21.4%) to MIA (9.5%) and IA (4.1%). At last, comparations between targeted NGS and ARMS-based single gene test in the detection of EGFR showed a 94.6% consistence. In conclusion, targeted NGS can provide a comprehensive mutational profile of lung cancer. Considering the high mutation rate of EGFR in NSCLC of Asian populations, a specialized detection strategy should be conducted.

Mutations, Bottlenecks, and Clonal Sweeps: How Environmental Carcinogens and Genomic Changes Shape Clonal Evolution during Tumor Progression

The transition from a single, initiated cell to a full-blown malignant tumor involves significant genomic evolution. Exposure to carcinogens-whether directly mutagenic or not-can drive progression toward malignancy, as can stochastic acquisition of cancer-promoting genetic events. Mouse models using both carcinogens and germline genetic manipulations have enabled precise inquiry into the evolutionary dynamics that take place as a tumor progresses from benign to malignant to metastatic stages. Tumor progression is characterized by changes in somatic point mutations and copy-number alterations, even though any single tumor can itself have a high or low burden of genomic alterations. Further, lineage-tracing, single-cell analyses and CRISPR barcoding have revealed the distinct clonal dynamics within benign and malignant tumors. Application of these tools in a range of mouse models can shed unique light on the patterns of clonal evolution that take place in both mouse and human tumors.

A gene-based score for the risk stratification of stage IA lung adenocarcinoma

OBJECTIVE

We aim to molecularly stratify stage IA lung adenocarcinoma (LUAD) for precision medicine.

METHODS

Twelve multi-institution datasets (837 cases of IA) were used to classify the high- and low-risk types (based on survival status within 5 years), and the biological differences were compared. Then, a gene-based classifying score (IA score) was trained, tested and validated by several machine learning methods. Furthermore, we estimated the significance of the IA score in the prognostic assessment, chemotherapy prediction and risk stratification of stage IA LUAD. We also developed an R package for the clinical application. The SEER database (15708 IA samples) and TCGA Pan-Cancer (1881 stage I samples) database were used to verify clinical significance.

RESULTS

Compared with the low-risk group, the high-risk group of stage IA LUAD has obvious enrichment of the malignant pathway and more driver mutations and copy number variations. The effect of the IA score on the classification of high- and low-risk stage IA LUAD was much better than that of classical clinicopathological factors (training set: AUC = 0.9, validation set: AUC = 0.7). The IA score can significantly predict the prognosis of stage IA LUAD and has a prognostic effect for stage I pancancer. The IA score can effectively predict chemotherapy sensitivity and occult metastasis or invasion in stage IA LUAD. The R package IAExpSuv has a good risk probability prediction effect for both groups and single stages of IA LUAD.

CONCLUSIONS

The IA score can effectively stratify the risk of stage IA LUAD, offering good assistance in precision medicine.

Whole-genome sequencing reveals the molecular implications of the stepwise progression of lung adenocarcinoma

The mechanism underlying the development of tumors, particularly at early stages, still remains mostly elusive. Here, we report whole-genome long and short read sequencing analysis of 76 lung cancers, focusing on very early-stage lung adenocarcinomas such as adenocarcinoma in situ (AIS) and minimally invasive adenocarcinoma. The obtained data is further integrated with bulk and spatial transcriptomic data and epigenomic data. These analyses reveal key events in lung carcinogenesis. Minimal somatic mutations in pivotal driver mutations and essential proliferative factors are the only detectable somatic mutations in the very early-stage of AIS. These initial events are followed by copy number changes and global DNA hypomethylation. Particularly, drastic changes are initiated at the later AIS stage, i.e., in Noguchi type B tumors, wherein cancer cells are exposed to the surrounding microenvironment. This study sheds light on the pathogenesis of lung adenocarcinoma from integrated pathological and molecular viewpoints.

Susceptibility Genes Associated with Multiple Primary Cancers

With advancements in treatment and screening techniques, we have been witnessing an era where more cancer survivors harbor multiple primary cancers (MPCs), affecting approximately one in six patients. Identifying MPCs is crucial for tumor staging and subsequent treatment choices. However, the current clinicopathological criteria for clinical application are limited and insufficient, making it challenging to differentiate them from recurrences or metastases. The emergence of next-generation sequencing (NGS) technology has provided a genetic perspective for defining multiple primary cancers. Researchers have found that, when considering multiple tumor pairs, it is crucial not only to examine well-known essential mutations like MLH1/MSH2, EGFR, PTEN, BRCA1/2, CHEK2, and TP53 mutations but also to explore certain pleiotropic loci. Moreover, specific deleterious mutations may serve as regulatory factors in second cancer development following treatment. This review aims to discuss these susceptibility genes and provide an explanation of their functions based on the signaling pathway background. Additionally, the association network between genetic signatures and different tumor pairs will be summarized.

Pathomic Features Reveal Immune and Molecular Evolution From Lung Preneoplasia to Invasive Adenocarcinoma

Recent statistics on lung cancer, including the steady decline of advanced diseases and the dramatically increasing detection of early-stage diseases and indeterminate pulmonary nodules, mark the significance of a comprehensive understanding of early lung carcinogenesis. Lung adenocarcinoma (ADC) is the most common histologic subtype of lung cancer, and atypical adenomatous hyperplasia is the only recognized preneoplasia to ADC, which may progress to adenocarcinoma in situ (AIS) and minimally invasive adenocarcinoma (MIA) and eventually to invasive ADC. Although molecular evolution during early lung carcinogenesis has been explored in recent years, the progress has been significantly hindered, largely due to insufficient materials from ADC precursors. Here, we employed state-of-the-art deep learning and artificial intelligence techniques to robustly segment and recognize cells on routinely used hematoxylin and eosin histopathology images and extracted 9 biology-relevant pathomic features to decode lung preneoplasia evolution. We analyzed 3 distinct cohorts (Japan, China, and United States) covering 98 patients, 162 slides, and 669 regions of interest, including 143 normal, 129 atypical adenomatous hyperplasia, 94 AIS, 98 MIA, and 205 ADC. Extracted pathomic features revealed progressive increase of atypical epithelial cells and progressive decrease of lymphocytic cells from normal to AAH, AIS, MIA, and ADC, consistent with the results from tissue-consuming and expensive molecular/immune profiling. Furthermore, pathomics analysis manifested progressively increasing cellular intratumor heterogeneity along with the evolution from normal lung to invasive ADC. These findings demonstrated the feasibility and substantial potential of pathomics in studying lung cancer carcinogenesis directly from the low-cost routine hematoxylin and eosin staining.

Comparative immunological landscape between pre- and early-stage LUAD manifested as ground-glass nodules revealed by scRNA and scTCR integrated analysis

BACKGROUND

Mechanism underlying the malignant progression of precancer to early-stage lung adenocarcinoma (LUAD) as well as their indolence nature remains elusive.

METHODS

Single-cell RNA sequencing (scRNA) with simultaneous T cell receptor (TCR) sequencing on 5 normal lung tissues, 3 precancerous and 4 early-stage LUAD manifested as pulmonary ground-glass nodules (GGNs) were performed.

RESULTS

Through this integrated analysis, we have delineated five key modules that drive the malignant progression of early-stage LUAD in a disease stage-dependent manner. These modules are related to cell proliferation and metabolism, immune response, mitochondria, cilia, and cell adhesion. We also find that the tumor micro-environment (TME) of early-stage LUAD manifested as GGN are featured with regulatory T (Tregs) cells accumulation with three possible origins, and loss-functional state (decreased clonal expansion and cytotoxicity) of CD8 + T cells. Instead of exhaustion, the CD8 + T cells are featured with a shift to memory phenotype, which is significantly different from the late stage LUAD. Furthermore, we have identified monocyte-derived macrophages that undergo a lipid-phenotype transition and may contribute to the suppressive TME. Intense interaction between stromal cells, myeloid cells including lipid associated macrophages and LAMP3 + DCs, and lymphocytes were also characterized.

CONCLUSIONS

Our work provides new insight into the molecular and cellular mechanism underlying malignant progression of LUAD manifested as GGN, and pave way for novel immunotherapies for GGN. Video Abstract.

Spatial downregulation of CD74 signatures may drive invasive component development in part-solid lung adenocarcinoma

Pulmonary nodules with part-solid imaging features manifest during the progression from preinvasive to invasive lung adenocarcinoma. To define the spatial composition and evolutionary trajectories of early-stage lung adenocarcinoma, we combined spatial transcriptomics (ST) and pathological annotations from 20 part-solid nodules (PSNs), four of which were matched with single-cell RNA sequencing. Two malignant cell populations (MC1 and MC2) were identified, and a linear evolutionary relationship was observed. Compared to MC2, the pre-existing malignant MC1 exhibited a lower metastatic signature, corresponding to the preinvasive component (lepidic) on pathology and the ground glass component on PSN imaging. Higher immune infiltration was observed among MC1 regions in ST profiles, and further analysis revealed that macrophages may be involved in this process through the CD74 axis. This work provides deeper insights into the evolutionary process and spatial immune cell composition behind PSNs and highlights the mechanisms of immune escape behind this adenocarcinoma trajectory.

Spatial heterogeneity of T cell repertoire across NSCLC tumors, tumor edges, adjacent and distant lung tissues

Background

A better understanding of T cells in lung cancer and their distribution across tumor-adjacent lungs and peripheral blood is needed to improve efficacy and minimize toxicity from immunotherapy to lung cancer patients.

Methods

Here, we performed CDR3β TCR sequencing of 136 samples from 20 patients with early-stage NSCLC including peripheral blood mononuclear cells, tumors, tumor edges (<1 cm from tumor), as well as adjacent lungs 1 cm, 2 cm, 5 cm, and 10 cm away from the tumor to gain insight into the spatial heterogeneity of T cells across the lungs in patients with NSCLC. PD-L1, CD4, and CD8 expression was assessed using immunohistochemical staining, and genomic features were derived by targeted sequencing of 1,021 cancer-related genes. Multiplex immunohistochemistry against PD-1, CTLA4, LAG3, and TIM3 was performed on four patients to assess T cell exhaustion.

Results

Our study reveals a decreasing gradient in TIL Tumor Infiltrating Lymphocytes homology with tumor edge, adjacent lungs, and peripheral blood but no discernible distance-associated patterns of T cell trafficking within the adjacent lung itself. Furthermore, we show a decrease in pathogen-specific TCRs in regions with high T cell clonality and PD-L1 expression.

Conclusions

Exclusion in T exhaustion cells at play across the lungs of patients with NSCLC may potentially be the mechanism for lung cancer occurrence.

Immune-active tumor-adjacent tissues are associated with favorable prognosis in stage I lung squamous cell carcinoma

The immunogenomic features of tumor-adjacent lungs (TALs) in stage I lung squamous cell carcinoma (LUSC) are not clear. Multiomics analyses of tumor tissues and paired TALs from 59 stage I LUSC patients were performed. Compared to tumors, TALs exhibited a better-preserved immune contexture indicated by upregulation of immune pathways, increased immune infiltration, and higher expression of immune effector molecules. Notably, TALs had no mutations in PTEN and KEAP1, a lower incidence of human leukocyte antigen (HLA) loss and higher expression of HLA class I genes, major histocompatibility complex (MHC) I chaperones, and interferon (IFN)-γ-associated genes. Digital spatial profiling validated the generally higher immune infiltration in TALs and revealed a higher level of immune heterogeneity in LUSC tumors. Importantly, patients with higher immune infiltration in TALs had significantly longer survival, while high immune heterogeneity was associated with inferior patient survival. Our work can be considered in the selection of patients for adjuvant therapy, especially immunotherapy.

Diagnostic performance of RASSF1A and SHOX2 methylation combined with EGFR mutations for differentiation between small pulmonary nodules

BACKGROUND AND AIM

Aberrant methylation of Ras association domain family 1, isoform A (RASSF1A), and short-stature homeobox gene 2 (SHOX2) promoters has been validated as a pair of valuable biomarkers for diagnosing early lung adenocarcinomas (LUADs). Epidermal growth factor receptor (EGFR) is the key driver mutation in lung carcinogenesis. This study aimed to investigate the aberrant promoter methylation of RASSF1A and SHOX2, and the genetic mutation of EGFR in 258 specimens of early LUADs.

METHODS

We retrospectively selected 258 paraffin-embedded samples of pulmonary nodules measuring 2 cm or less in diameter and evaluated the diagnostic performance of individual biomarker assays and multiple panels between noninvasive (group 1) and invasive lesions (groups 2A and 2B). Then, we investigated the interaction between genetic and epigenetic alterations.

RESULTS

The degree of RASSF1A and SHOX2 promoter methylation and EGFR mutation was significantly higher in invasive lesions than in noninvasive lesions. The three biomarkers distinguished between noninvasive and invasive lesions with reliable sensitivity and specificity: 60.9% sensitivity [95% confidence interval (CI) 52.41-68.78] and 80.0% specificity (95% CI 72.14-86.07). The novel panel biomarkers could further discriminate among three invasive pathological subtypes (area under the curve value > 0.6). The distribution of RASSF1A methylation and EGFR mutation was considerably exclusive in early LUAD (P = 0.002).

CONCLUSION

DNA methylation of RASSF1A and SHOX2 is a pair of promising biomarkers, which may be used in combination with other driver alterations, such as EGFR mutation, to support the differential diagnosis of LUADs, especially for stage I.

Single-cell profiling reveals the trajectory of FOLR2-expressing tumor-associated macrophages to regulatory T cells in the progression of lung adenocarcinoma

An immunosuppressive microenvironment enriched with regulatory CD4+ T lymphocytes (Tregs) facilitates the progression of lung adenocarcinoma (LUAD). This study aims to investigate the cellular mechanism underlying the formation of the immunosuppressive microenvironment in LUAD. LUAD samples (n = 12) and normal lung samples (n = 3) were obtained from patients with different pathological stages of LUAD. Single-cell RNA sequencing was performed to classify cellular components and analyze the transcriptomes, including transcription factors/targets and chemokine ligands/receptors, followed by bioinformatics study such as pseudotime analysis. Myeloid cells and T cells were the most abundant cell types in tumors and normal lung tissues, while tumor-associated macrophage-folate receptor 2 (TAM-FOLR2) and CD4+ nuclear receptor subfamily 4 group A member 3 (NR4A3) exhibited sharp increases in invasive adenocarcinoma (IA). The enrichment of TAM-FOLR2 in IA might result from alveolar resident macrophage-resistin (ARM-RETN) transformation and recruitment of dendritic cells (DCs) and other TAMs, as evidenced by temporal trajectories and differential expression profiles of chemokine ligands/receptors versus those in the early stages of tumors. High expression of CCL17/19/22 was observed in IA as well as in DCs, along with the strong interaction of TAM-FOLR2 with DCs. The results of pseudotime analysis suggested that CD4+NR4A3 might potentially convert to CD4+FOXP3, further supported by the high expression of NR4A3 target genes in CD4+FOXP3 cells. This study provides a single-cell transcriptome atlas from preinvasive to invasive LUAD and reveals a potential ARM-RETN/TAM-FOLR2/DCs/CD4+NR4A3/CD4+FOXP3 trajectory in shaping the immune suppressive microenvironment along the pathogenesis of LUAD.

Advances in Mapping Tumor Progression from Precancer Atlases

Tissue profiling technologies present opportunities for understanding transition from precancerous lesions to malignancy, which may impact risk stratification, prevention, and even cancer treatment. A human precancer atlas building effort is ongoing to tackle the significant challenge of decoding the heterogeneity among cells, specimens, and patients. Here, we discuss the findings resulting from atlases built across precancer types, including those found in colon, breast, lung, stomach, cervix, and skin, using bulk, single-cell, and spatial profiling strategies. We highlight two main themes that emerge across precancer types: the ordering of molecular events that occur during tumor progression and the fluctuation of microenvironmental response during precancer progression. We further highlight the key challenges of data integration across large cohorts of patients, and the need for computational tools to reliably annotate and quality control high-volume, high-dimensional data.

Differences of molecular events driving pathological and radiological progression of lung adenocarcinoma

BACKGROUND

Ground-glass opacity (GGO)-like lung adenocarcinoma (LUAD) has been detected increasingly in the clinic and its inert property and superior survival indicate unique biological characteristics. However, we do not know much about them, which hampers identification of key reasons for the inert property of GGO-like LUAD.

METHODS

Using whole-exome sequencing and RNA sequencing, taking into account both radiological and pathological classifications of the same 197 patients concomitantly, we systematically interrogate genes driving the progression from GGO to solid nodule and potential reasons for the inertia of GGO. Using flow cytometry and IHC, we validated the abundance of immune cells and activity of cell proliferation.

FINDINGS

Identifying the differences between GGO and solid nodule, we found adenocarcinoma in situ/minimally invasive adenocarcinoma (AIS/MIA) and GGO-like LUAD exhibited lower TP53 mutation frequency and less active cell proliferation-related pathways than solid nodule in LUAD. Identifying the differences in GGO between AIS/MIA and LUAD, we noticed that EGFR mutation frequency and CNV load were significantly higher in LUAD than in AIS/MIA. Regulatory T cell was also higher in LUAD, while CD8+ T cell decreased from AIS/MIA to LUAD. Finally, we constructed a transcriptomic signature to quantify the development from GGO to solid nodule, which was an independent predictor of patients' prognosis in 11 external LUAD datasets.

INTERPRETATION

Our results provide deeper insights into the indolent nature of GGO and provide a molecular basis for the treatment of GGO-like LUAD.

FUNDING

This study was supported in part by the National Natural Science Foundation of China (32170657), the National Natural Science Foundation of China (82203037), and Shanghai Sailing Program (22YF1408900).

Tumor cell heterogeneity drives spatial organization of the intratumoral immune response in squamous cell skin carcinoma

Intratumoral heterogeneity (ITH)-defined as genetic and cellular diversity within a tumor-is linked to failure of immunotherapy and an inferior anti-tumor immune response. The underlying mechanism of this association is unknown. To address this question, we modeled heterogeneous tumors comprised of a pro-inflammatory ("hot") and an immunosuppressive ("cold") tumor population, labeled with YFP and RFP tags respectively to enable precise spatial tracking. The resulting mixed-population tumors exhibited distinct regions comprised of YFP+ (hot) cells, RFP+ (cold) cells, or a mixture. We found that tumor regions occupied by hot tumor cells (YFP+) harbored more total T cells and a higher frequency of Th1 cells and IFNγ+ CD8 T cells compared to regions occupied by cold tumor cells (RFP+), whereas immunosuppressive macrophages showed the opposite spatial pattern. We identified the chemokine CX3CL1, produced at higher levels by our cold tumors, as a mediator of intratumoral macrophage accumulation, particularly immunosuppressive CD206Hi macrophages. Furthermore, we examined the response of heterogeneous tumors to a therapeutic combination of PD-1 blockade and CD40 agonist on a region-by-region basis. While the combination successfully increases Th1 abundance in "cold" tumor regions, it fails to bring overall T cell activity to the same level as seen in "hot" regions. The presence of the "cold" cells thus ultimately leads to a failure of the therapy to induce tumor rejection. Collectively, our results demonstrate that the organization of heterogeneous tumor cells has a profound impact on directing the spatial organization and function of tumor-infiltrating immune cells as well as on responses to immunotherapy.

B4GALT1 promotes immune escape by regulating the expression of PD-L1 at multiple levels in lung adenocarcinoma

BACKGROUND

Invasive adenocarcinoma (IAC), which is typically preceded by minimally invasive adenocarcinoma (MIA), is the dominant pathological subtype of early-stage lung adenocarcinoma (LUAD). Identifying the molecular events underlying the progression from MIA to IAC may provide a crucial perspective and boost the exploration of novel strategies for early-stage LUAD diagnosis and treatment.

METHODS

Transcriptome sequencing of four pairs of MIA and IAC tumours obtained from four multiple primary lung cancer patients was performed to screen out beta-1,4-galactosyltransferase1 (B4GALT1). Function and mechanism experiments in vitro and in vivo were performed to explore the regulatory mechanism of B4GALT1-mediated immune evasion by regulating programmed cell death ligand 1 (PD-L1).

RESULTS

B4GALT1, a key gene involved in N-glycan biosynthesis, was highly expressed in IAC samples. Further experiments revealed that B4GALT1 regulated LUAD cell proliferation and invasion both in vitro and in vivo and was related to the impaired antitumour capacity of CD8 + T cells. Mechanistically, B4GALT1 directly mediates the N-linked glycosylation of PD-L1 protein, thus preventing PD-L1 degradation at the posttranscriptional level. In addition, B4GALT1 stabilized the TAZ protein via glycosylation, which activated CD274 at the transcriptional level. These factors lead to lung cancer immune escape. Importantly, inhibition of B4GALT1 increased CD8 + T-cell abundance and activity and enhanced the antitumour immunity of anti-PD-1 therapy in vivo.

CONCLUSION

B4GALT1 is a critical molecule in the development of early-stage LUAD and may be a novel target for LUAD intervention and immunotherapy.

Bootstrap confidence for molecular evolutionary estimates from tumor bulk sequencing data

Bulk sequencing is commonly used to characterize the genetic diversity of cancer cell populations in tumors and the evolutionary relationships of cancer clones. However, bulk sequencing produces aggregate information on nucleotide variants and their sample frequencies, necessitating computational methods to predict distinct clone sequences and their frequencies within a sample. Interestingly, no methods are available to measure the statistical confidence in the variants assigned to inferred clones. We introduce a bootstrap resampling approach that combines clone prediction and statistical confidence calculation for every variant assignment. Analysis of computer-simulated datasets showed the bootstrap approach to work well in assessing the reliability of predicted clones as well downstream inferences using the predicted clones (e.g., mapping metastatic migration paths). We found that only a fraction of inferences have good bootstrap support, which means that many inferences are tentative for real data. Using the bootstrap approach, we analyzed empirical datasets from metastatic cancers and placed bootstrap confidence on the estimated number of mutations involved in cell migration events. We found that the numbers of driver mutations involved in metastatic cell migration events sourced from primary tumors are similar to those where metastatic tumors are the source of new metastases. So, mutations with driver potential seem to keep arising during metastasis. The bootstrap approach developed in this study is implemented in software available at https://github.com/SayakaMiura/CloneFinderPlus.

Evolution of lung adenocarcinoma from preneoplasia to invasive adenocarcinoma

OBJECTIVE

Mutations in driver genes contribute to the development and progression of lung adenocarcinoma (LUAD). However, in the dynamic evolutionary process from adenocarcinoma in situ (AIS) to minimally invasive adenocarcinoma (MIA) and eventually to invasive adenocarcinoma (IAC), the role of driver genes is currently unclear. This study aimed to analyse the role of driver gene status in the progression of LUAD from preneoplasia to IAC.

METHODS

Patients with LUAD who underwent surgery in our centre from March 2015 to December 2019 were retrospectively analysed, and LUAD patients with tumour sizes ≤3.0 cm and pN0 were included in the final analysis. The mutation status of common driver genes, including EGFR, ALK and ROS1, was detected. According to the pathological characteristics, the patients were divided into three stages: AIS, MIA and IAC. We analysed the distribution of driver gene mutation frequencies across three stages of LUAD. In addition, we performed univariate and multivariate analyses of IAC patients to screen for relevant variables (driver genes and clinicopathological features) affecting their prognosis.

RESULTS

Ultimately, 759 patients with LUAD were enrolled, including 135, 130, and 494 cases of AIS, MIA, and IAC, respectively. EGFR mutations were identified in 359 (61.8%) patients, and with the transition from AIS to MIA, the frequency of EGFR mutations increased from 33.3% to 50.8%, p = 0.004, whereas the frequency of EGFR mutations was comparable for MIA and IAC (50.8% vs. 50.2%, p = 0.922). Moreover, ALK and ROS1 gene fusions were identified in 17 cases (2.2%) and 2 cases (3.0‰) respectively. For AIS, neither ALK gene nor ROS1 gene fusions were observed. When the tumour progressed to MIA, the ALK fusion frequency was 2.3% (3/130), which was basically consistent with the ALK fusion frequency of 2.8% in IAC, p = 0.143. For IAC, fusions of ROS1 fell into this category. In addition, we found that 40 patients (5.3%) developed metastasis/recurrence, and 14 patients (1.8%) died of cancer-specific related diseases. Notably, for AIS, there were no recurrences and no deaths, and for MIA, only 1 patient died with LUAD. Finally, survival analysis was performed in patients with stage IA invasive adenocarcinoma, and EGFR-mutant patients showed better DFS than EGFR-wild-type patients (p = 0.036). Conversely, patients with ALK fusions showed worse DFS than those with ALK wild-type (p = 0.004), and the same results were found in OS analysis.

CONCLUSIONS

The accumulation of EGFR driver gene mutation frequencies mediates the progression of LUAD from AIS to MIA. When the tumour progresses to stage IA invasive adenocarcinoma, multivariate analysis based on driver gene status can be used as a pivotal prognostic factor.

Diagnosis, Monitoring, and Prognosis of Liquid Biopsy in Cancer Immunotherapy

Liquid biopsy (LB), as a minimally invasive method of gleaning insight into the dynamics of diseases through a patient fluid sample, represents an interesting tool that can advise in disease monitoring, treatment selection, early diagnosis, evaluation of the response, and prognosis. Cancer immunotherapy is a breakthrough in cancer treatment, which is now recognized as the "fourth pillar" of cancer treatment, after surgery, chemotherapy, and radiotherapy. Liquid biopsy offers a different befalling for beneath invasive diagnosis, real-time accommodating monitoring, and analysis options, involving the isolation of circulating biomarkers, such as cell-free DNA (cfDNA), circulating tumor cells (CTCs), exosomes, and microRNAs (miRNAs). The biomarkers herein have great potential to allow the realization of liquid biopsy for predicting the immunotherapy response and precision medicine. Liquid biopsy offers an alternative, less invasive approach to select cancer patients who would benefit from immunotherapy and to monitor patients during their disease course. This review focuses on the use of liquid biopsy in the immunotherapy treatment of patients with cancer. In this review, we addressed the different promising liquid biopsy-based biomarkers in cancer patients that enable the selection of patients who benefit from immunotherapy and the monitoring of patients during this therapy.

Comprehensive characterization of genomic and radiologic features reveals distinct driver patterns of RTK/RAS pathway in ground-glass opacity pulmonary nodules

Ground-glass opacity (GGO)-associated pulmonary nodules have been known as a radiologic feature of early-stage lung cancers and exhibit an indolent biological behavior. However, the correlation between driver genes and radiologic features as well as the immune microenvironment remains poorly understood. We performed a custom 1021-gene panel sequencing of 334 resected pulmonary nodules presenting as GGO from 262 Chinese patients. A total of 130 multiple pulmonary nodules were sampled from 58 patients. Clinical-pathologic and radiologic parameters of these pulmonary nodules were collected. Immunohistochemistry (IHC) and multiplex immunofluorescent staining (mIF) were applied to analyze proliferation and immune cell markers of GGO-associated pulmonary nodules. Compared with pure GGO nodules, mixed GGO nodules were enriched for invasive adenocarcinoma (IAC) (182/216 vs 73/118, P &lt; .001). Eighty-eight percent (294/334) of GGO-associated nodules carried at least one mutation in EGFR/ERBB2/BRAF/KRAS/MAP2K1 of the RTK/RAS signaling pathway, and the alterations in these driver genes were mutually exclusive. The analysis of multifocal pulmonary nodules from the same patient revealed evidence of functional convergence on RTK/RAS pathways. Nodules with ERBB2/BRAF/MAP2K1 mutations tended to be more indolent than those with EGFR and KRAS mutations. IHC and mIF staining showed that KRAS-mutant GGO nodules displayed higher infiltration of CD4+ T cell and CD8+ T cell as well as stronger proliferation and immune inhibitory signals. Our study demonstrates a driver landscape of radiologically detectable GGO-associated pulmonary nodules in Chinese patients and supports that different driver patterns in RTK/RAS pathway are corresponding to different radiologic features.