An NKX2-1/ERK/WNT feedback loop modulates gastric identity and response to targeted therapy in lung adenocarcinoma

Significance of NKX2-1 as a biomarker for clinical prognosis, immune infiltration, and drug therapy in lung squamous cell carcinoma

Background

This study was performed to determine the biological processes in which NKX2-1 is involved and thus its role in the development of lung squamous cell carcinoma (LUSC) toward improving the prognosis and treatment of LUSC.

Methods

Raw RNA sequencing (RNA-seq) data of LUSC from The Cancer Genome Atlas (TCGA) were used in bioinformatics analysis to characterize NKX2-1 expression levels in tumor and normal tissues. Survival analysis of Kaplan-Meier curve, the time-dependent receiver operating characteristic (ROC) curve, and a nomogram were used to analyze the prognosis value of NKX2-1 for LUSC in terms of overall survival (OS) and progression-free survival (PFS). Then, differentially expressed genes (DEGs) were identified, and Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), and Gene Set Enrichment Analysis (GSEA) were used to clarify the biological mechanisms potentially involved in the development of LUSC. Moreover, the correlation between the NKX2-1 expression level and tumor mutation burden (TMB), tumor microenvironment (TME), and immune cell infiltration revealed that NKX2-1 participates in the development of LUSC. Finally, we studied the effects of NKX2-1 on drug therapy. To validate the protein and gene expression levels of NKX2-1 in LUSC, we employed immunohistochemistry(IHC) datasets, The Gene Expression Omnibus (GEO) database, and qRT-PCR analysis.

Results

NKX2-1 expression levels were significantly lower in LUSC than in normal lung tissue. It significantly differed in gender, stage and N classification. The survival analysis revealed that high expression of NKX2-1 had shorter OS and PFS in LUSC. The multivariate Cox regression hazard model showed the NKX2-1 expression as an independent prognostic factor. Then, the nomogram predicted LUSC prognosis. There are 51 upregulated DEGs and 49 downregulated DEGs in the NKX2-1 high-level groups. GO, KEGG and GSEA analysis revealed that DEGs were enriched in cell cycle and DNA replication.The TME results show that NKX2-1 expression was positively associated with mast cells resting, neutrophils, monocytes, T cells CD4 memory resting, and M2 macrophages but negatively associated with M1 macrophages. The TMB correlated negatively with NKX2-1 expression. The pharmacotherapy had great sensitivity in the NKX2-1 low-level group, the immunotherapy is no significant difference in the NKX2-1 low-level and high-level groups. The analysis of GEO data demonstrated concurrence with TCGA results. IHC revealed NKX2-1 protein expression in tumor tissues of both LUAD and LUSC. Meanwhile qRT-PCR analysis indicated a significantly lower NKX2-1 expression level in LUSC compared to LUAD. These qRT-PCR findings were consistent with co-expression analysis of NKX2-1.

Conclusion

We conclude that NKX2-1 is a potential biomarker for prognosis and treatment LUSC. A new insights of NKX2-1 in LUSC is still needed further research.

Mcam stabilizes luminal progenitor breast cancer phenotypes via Ck2 control and Src/Akt/Stat3 attenuation

Breast cancers are categorized into subtypes with distinctive therapeutic vulnerabilities and prognoses based on their expression of clinically targetable receptors and gene expression patterns mimicking different cell types of the normal gland. Here, we tested the role of Mcam in breast cancer cell state control and tumorigenicity in a luminal progenitor-like murine tumor cell line (Py230) that exhibits lineage and tumor subtype plasticity. Mcam knockdown Py230 cells show augmented Stat3 and Pi3K/Akt activation associated with a lineage state switch away from a hormone-sensing/luminal progenitor state toward alveolar and basal cell related phenotypes that were refractory to growth inhibition by the anti-estrogen therapeutic, tamoxifen. Inhibition of Stat3, or the upstream activator Ck2, reversed these cell state changes. Mcam binds Ck2 and acts as a regulator of Ck2 substrate utilization across multiple mammary tumor cell lines. In Py230 cells this activity manifests as increased mesenchymal morphology, migration, and Src/Fak/Mapk/Paxillin adhesion complex signaling in vitro, in contrast to Mcam's reported roles in promoting mesenchymal phenotypes. In vivo, Mcam knockdown reduced tumor growth and take rate and inhibited cell state transition to Sox10+/neural crest like cells previously been associated with tumor aggressiveness. This contrasts with human luminal breast cancers where MCAM copy number loss is highly coupled to Cyclin D amplification, increased proliferation, and the more aggressive Luminal B subtype. Together these data indicate a critical role for Mcam and its regulation of Ck2 in control of breast cancer cell state plasticity with implications for progression, evasion of targeted therapies and combination therapy design.

A Lung Cancer Mouse Model Database

Lung cancer, the leading cause of cancer mortality, exhibits diverse histological subtypes and genetic complexities. Numerous preclinical models have been developed to study lung cancer, but data from these models are disparate, siloed, and difficult to compare in a centralized fashion. Here we established the Lung Cancer Mouse Model Database (LCMMDB), an extensive repository of 1,354 samples from 77 transcriptomic datasets. Meticulous curation and collaboration with data depositors have produced a robust and comprehensive database, enhancing the fidelity of the genetic landscape it depicts. The LCMMDB aligns 859 tumors from genetically engineered mouse models (GEMMs) with human lung cancer mutations, enabling comparative analysis and revealing a pressing need to broaden the diversity of genetic aberrations modeled in GEMMs. Accompanying this resource, we developed a web application that offers researchers intuitive tools for in-depth gene expression analysis and fostering potential collaborations. With standardized reprocessing of gene expression data, the LCMMDB serves as a powerful platform for cross-study comparison and lays the groundwork for future research, aiming to bridge the gap between mouse models and human lung cancer for improved translational relevance.

Decoding the basis of histological variation in human cancer

Molecular abnormalities that shape human neoplasms dissociate their phenotypic landscape from that of the healthy counterpart. Through the lens of a microscope, tumour pathology optically captures such aberrations projected onto a tissue slide and has categorized human epithelial neoplasms into distinct histological subtypes based on the diverse morphogenetic and molecular programmes that they manifest. Tumour histology often reflects tumour aggressiveness, patient prognosis and therapeutic vulnerability, and thus has been used as a de facto diagnostic tool and for making clinical decisions. However, it remains elusive how the diverse histological subtypes arise and translate into pleiotropic biological phenotypes. Molecular analysis of clinical tumour tissues and their culture, including patient-derived organoids, and add-back genetic reconstruction of tumorigenic pathways using gene engineering in culture models and rodents further elucidated molecular mechanisms that underlie morphological variations. Such mechanisms include genetic mutations and epigenetic alterations in cellular identity codes that erode hard-wired morphological programmes and histologically digress tumours from the native tissues. Interestingly, tumours acquire the ability to grow independently of the niche-driven stem cell ecosystem along with these morphological alterations, providing a biological rationale for histological diversification during tumorigenesis. This Review comprehensively summarizes our current understanding of such plasticity in the histological and lineage commitment fostered cooperatively by molecular alterations and the tumour environment, and describes basic and clinical implications for future cancer therapy.

FoxA1/2-dependent epigenomic reprogramming drives lineage switching in lung adenocarcinoma

The ability of cancer cells to alter their identity is essential for tumor survival and progression. Loss of the pulmonary lineage specifier NKX2-1 within KRAS-driven lung adenocarcinoma (LUAD) enhances tumor progression and results in a pulmonary-to-gastric lineage switch that is dependent upon the activity of pioneer factors FoxA1 and FoxA2; however, the underlying mechanism remains largely unknown. Here, we show that FoxA1/2 reprogram the epigenetic landscape of NKX2-1-negative LUAD to facilitate a gastric identity. After Nkx2-1 deletion, FoxA1/2 mediate demethylation of gastric-defining genes through recruitment of TET3, an enzyme that induces DNA demethylation. H3K27ac ChIP-seq and HiChIP show that FoxA1/2 also control the activity of regulatory elements and their 3D interactions at gastric loci. Furthermore, oncogenic KRAS is required for the FoxA1/2-dependent epigenetic reprogramming. This work demonstrates the role of FoxA1/2 in rewiring the methylation and histone landscape and cis-regulatory dynamics of NKX2-1-negative LUAD to drive cancer cell lineage switching.

HDAC3 is critical in tumor development and therapeutic resistance in Kras-mutant non-small cell lung cancer

HDAC3 is one of the main targets of histone deacetylase (HDAC) inhibitors in clinical development as cancer therapies, yet the in vivo role of HDAC3 in solid tumors is unknown. We identified a critical role for HDAC3 in Kras-mutant lung cancer. Using genetically engineered mouse models (GEMMs), we found that HDAC3 is required for lung tumor growth in vivo. HDAC3 was found to direct and enhance the transcription effects of the lung cancer lineage transcription factor NKX2-1 to mediate expression of a common set of target genes. We identified FGFR1 as a critical previously unidentified target of HDAC3. Leveraging this, we identified that an HDAC3-dependent transcriptional cassette becomes hyperactivated as Kras/LKB1-mutant cells develop resistance to the MEK inhibitor trametinib, and this can be reversed by treatment with the HDAC1/HDAC3 inhibitor entinostat. We found that the combination of entinostat plus trametinib treatment elicits therapeutic benefit in the Kras/LKB1 GEMM.

Organoid modeling reveals the tumorigenic potential of the alveolar progenitor cell state

Alveolar type 2 (AT2) cells, the epithelial progenitor cells of the distal lung, are known to be the prominent cell of origin for lung adenocarcinoma. The regulatory programs that control chromatin and gene expression in AT2 cells during the early stages of tumor initiation are not well understood. Here, we dissected the response of AT2 cells to Kras activation and p53 loss (KP) using combined single cell RNA and ATAC sequencing in an established tumor organoid system. Multi-omic analysis showed that KP tumor organoid cells exhibit two major cellular states: one more closely resembling AT2 cells (SPC-high) and another with loss of AT2 identity (hereafter, Hmga2-high). These cell states are characterized by unique transcription factor (TF) networks, with SPC-high states associated with TFs known to regulate AT2 cell fate during development and homeostasis, and distinct TFs associated with the Hmga2-high state. CD44 was identified as a marker of the Hmga2-high state, and was used to separate organoid cultures for functional comparison of these two cell states. Organoid assays and orthotopic transplantation studies indicated that SPC-high cells have higher tumorigenic capacity in the lung microenvironment compared to Hmga2-high cells. These findings highlight the utility of understanding chromatin regulation in the early oncogenic versions of epithelial cells, which may reveal more effective means to intervene the progression of Kras-driven lung cancer.

Transcriptional Circuitry of NKX2-1 and SOX1 Defines an Unrecognized Lineage Subtype of Small-Cell Lung Cancer

Rationale: The current molecular classification of small-cell lung cancer (SCLC) on the basis of the expression of four lineage transcription factors still leaves its major subtype SCLC-A as a heterogeneous group, necessitating more precise characterization of lineage subclasses. Objectives: To refine the current SCLC classification with epigenomic profiles and to identify features of the redefined SCLC subtypes. Methods: We performed unsupervised clustering of epigenomic profiles on 25 SCLC cell lines. Functional significance of NKX2-1 (NK2 homeobox 1) was evaluated by cell growth, apoptosis, and xenograft using clustered regularly interspaced short palindromic repeats-Cas9 (CRISPR-associated protein 9)-mediated deletion. NKX2-1-specific cistromic profiles were determined using chromatin immunoprecipitation followed by sequencing, and its functional transcriptional partners were determined using coimmunoprecipitation followed by mass spectrometry. Rb1flox/flox; Trp53flox/flox and Rb1flox/flox; Trp53flox/flox; Nkx2-1flox/flox mouse models were engineered to explore the function of Nkx2-1 in SCLC tumorigenesis. Epigenomic landscapes of six human SCLC specimens and 20 tumors from two mouse models were characterized. Measurements and Main Results: We identified two epigenomic subclusters of the major SCLC-A subtype: SCLC-Aα and SCLC-Aσ. SCLC-Aα was characterized by the presence of a super-enhancer at the NKX2-1 locus, which was observed in human SCLC specimens and a murine SCLC model. We found that NKX2-1, a dual lung and neural lineage factor, is uniquely relevant in SCLC-Aα. In addition, we found that maintenance of this neural identity in SCLC-Aα is mediated by collaborative transcriptional activity with another neuronal transcriptional factor, SOX1 (SRY-box transcription factor 1). Conclusions: We comprehensively describe additional epigenomic heterogeneity of the major SCLC-A subtype and define the SCLC-Aα subtype by the core regulatory circuitry of NKX2-1 and SOX1 super-enhancers and their functional collaborations to maintain neuronal linage state.

Transposon Mutagenesis Reveals RBMS3 Silencing as a Promoter of Malignant Progression of BRAFV600E-Driven Lung Tumorigenesis

Mutationally activated BRAF is detected in approximately 7% of human lung adenocarcinomas, with BRAFT1799A serving as a predictive biomarker for treatment of patients with FDA-approved inhibitors of BRAFV600E oncoprotein signaling. In genetically engineered mouse (GEM) models, expression of BRAFV600E in the lung epithelium initiates growth of benign lung tumors that, without additional genetic alterations, rarely progress to malignant lung adenocarcinoma. To identify genes that cooperate with BRAFV600E for malignant progression, we used Sleeping Beauty-mediated transposon mutagenesis, which dramatically accelerated the emergence of lethal lung cancers. Among the genes identified was Rbms3, which encodes an RNA-binding protein previously implicated as a putative tumor suppressor. Silencing of RBMS3 via CRISPR/Cas9 gene editing promoted growth of BRAFV600E lung organoids and promoted development of malignant lung cancers with a distinct micropapillary architecture in BRAFV600E and EGFRL858R GEM models. BRAFV600E/RBMS3Null lung tumors displayed elevated expression of Ctnnb1, Ccnd1, Axin2, Lgr5, and c-Myc mRNAs, suggesting that RBMS3 silencing elevates signaling through the WNT/β-catenin signaling axis. Although RBMS3 silencing rendered BRAFV600E-driven lung tumors resistant to the effects of dabrafenib plus trametinib, the tumors were sensitive to inhibition of porcupine, an acyltransferase of WNT ligands necessary for their secretion. Analysis of The Cancer Genome Atlas patient samples revealed that chromosome 3p24, which encompasses RBMS3, is frequently lost in non-small cell lung cancer and correlates with poor prognosis. Collectively, these data reveal the role of RBMS3 as a lung cancer suppressor and suggest that RBMS3 silencing may contribute to malignant NSCLC progression.

SIGNIFICANCE

Loss of RBMS3 cooperates with BRAFV600E to induce lung tumorigenesis, providing a deeper understanding of the molecular mechanisms underlying mutant BRAF-driven lung cancer and potential strategies to more effectively target this disease.

Modeling diverse genetic subtypes of lung adenocarcinoma with a next-generation alveolar type 2 organoid platform

Lung cancer is the leading cause of cancer-related death worldwide. Lung adenocarcinoma (LUAD), the most common histological subtype, accounts for 40% of all cases. While existing genetically engineered mouse models (GEMMs) recapitulate the histological progression and transcriptional evolution of human LUAD, they are time-consuming and technically demanding. In contrast, cell line transplant models are fast and flexible, but these models fail to capture the full spectrum of disease progression. Organoid technologies provide a means to create next-generation cancer models that integrate the most advantageous features of autochthonous and transplant-based systems. However, robust and faithful LUAD organoid platforms are currently lacking. Here, we describe optimized conditions to continuously expand murine alveolar type 2 (AT2) cells, a prominent cell of origin for LUAD, in organoid culture. These organoids display canonical features of AT2 cells, including marker gene expression, the presence of lamellar bodies, and an ability to differentiate into the AT1 lineage. We used this system to develop flexible and versatile immunocompetent organoid-based models of KRAS, BRAF, and ALK mutant LUAD. Notably, organoid-based tumors display extensive burden and complete penetrance and are histopathologically indistinguishable from their autochthonous counterparts. Altogether, this organoid platform is a powerful, versatile new model system to study LUAD.

NKX2-1 controls lung cancer progression by inducing DUSP6 to dampen ERK activity

The RAS→RAF→MEK→ERK pathway is hyperactivated in the majority of human lung adenocarcinoma (LUAD). However, the initial activating mutations induce homeostatic feedback mechanisms that limit ERK activity. How ERK activation reaches the tumor-promoting levels that overcome the feedback and drive malignant progression is unclear. We show here that the lung lineage transcription factor NKX2-1 suppresses ERK activity. In human tissue samples and cell lines, xenografts, and genetic mouse models, NKX2-1 induces the ERK phosphatase DUSP6, which inactivates ERK. In tumor cells from late-stage LUAD with silenced NKX2-1, re-introduction of NKX2-1 induces DUSP6 and inhibits tumor growth and metastasis. We show that DUSP6 is necessary for NKX2-1-mediated inhibition of tumor progression in vivo and that DUSP6 expression is sufficient to inhibit RAS-driven LUAD. Our results indicate that NKX2-1 silencing, and thereby DUSP6 downregulation, is a mechanism by which early LUAD can unleash ERK hyperactivation for tumor progression.

Epithelial cell plasticity: breaking boundaries and changing landscapes

Epithelial tissues respond to a wide variety of environmental and genotoxic stresses. As an adaptive mechanism, cells can deviate from their natural paths to acquire new identities, both within and across lineages. Under extreme conditions, epithelial tissues can utilize "shape-shifting" mechanisms whereby they alter their form and function at a tissue-wide scale. Mounting evidence suggests that in order to acquire these alternate tissue identities, cells follow a core set of "tissue logic" principles based on developmental paradigms. Here, we review the terminology and the concepts that have been put forward to describe cell plasticity. We also provide insights into various cell intrinsic and extrinsic factors, including genetic mutations, inflammation, microbiota, and therapeutic agents that contribute to cell plasticity. Additionally, we discuss recent studies that have sought to decode the "syntax" of plasticity-i.e., the cellular and molecular principles through which cells acquire new identities in both homeostatic and malignant epithelial tissues-and how these processes can be manipulated for developing novel cancer therapeutics.

An NKX2-1/ERK/WNT feedback loop modulates gastric identity and response to targeted therapy in lung adenocarcinoma

Cancer cells undergo lineage switching during natural progression and in response to therapy. NKX2-1 loss in human and murine lung adenocarcinoma leads to invasive mucinous adenocarcinoma (IMA), a lung cancer subtype that exhibits gastric differentiation and harbors a distinct spectrum of driver oncogenes. In murine BRAFV600E-driven lung adenocarcinoma, NKX2-1 is required for early tumorigenesis, but dispensable for established tumor growth. NKX2-1-deficient, BRAFV600E-driven tumors resemble human IMA and exhibit a distinct response to BRAF/MEK inhibitors. Whereas BRAF/MEK inhibitors drive NKX2-1-positive tumor cells into quiescence, NKX2-1-negative cells fail to exit the cell cycle after the same therapy. BRAF/MEK inhibitors induce cell identity switching in NKX2-1-negative lung tumors within the gastric lineage, which is driven in part by WNT signaling and FoxA1/2. These data elucidate a complex, reciprocal relationship between lineage specifiers and oncogenic signaling pathways in the regulation of lung adenocarcinoma identity that is likely to impact lineage-specific therapeutic strategies.