An integrative approach unveils FOSL1 as an oncogene vulnerability in KRAS-driven lung and pancreatic cancer

FRA1 controls acinar cell plasticity during murine KrasG12D-induced pancreatic acinar to ductal metaplasia

Acinar cells have been proposed as a cell-of-origin for pancreatic ductal adenocarcinoma (PDAC) after undergoing acinar-to-ductal metaplasia (ADM). ADM can be triggered by pancreatitis, causing acinar cells to de-differentiate to a ductal-like state. We identify FRA1 (gene name Fosl1) as the most active transcription factor during KrasG12D acute pancreatitis-mediated injury, and we have elucidated a functional role of FRA1 by generating an acinar-specific Fosl1 knockout mouse expressing KrasG12D. Using a gene regulatory network and pseudotime trajectory inferred from single-nuclei ATAC-seq and bulk RNA sequencing (RNA-seq), we hypothesized a regulatory model of the acinar-ADM-pancreatic intraepithelial neoplasia (PanIN) continuum and experimentally validated that Fosl1 knockout mice are delayed in the onset of ADM and neoplastic transformation. Our study also identifies that pro-inflammatory cytokines, such as granulocyte colony stimulating factor (G-CSF), can regulate FRA1 activity to modulate ADM. Our findings identify that FRA1 is a mediator of acinar cell plasticity and is critical for acinar cell de-differentiation and transformation.

CRISPR-Cas9 Screening Identifies KRAS-Induced COX2 as a Driver of Immunotherapy Resistance in Lung Cancer

Oncogenic KRAS impairs antitumor immune responses. As effective strategies to combine KRAS inhibitors and immunotherapies have so far proven elusive, a better understanding of the mechanisms by which oncogenic KRAS drives immune evasion is needed to identify approaches that could sensitize KRAS-mutant lung cancer to immunotherapy. In vivo CRISPR-Cas9 screening in an immunogenic murine lung cancer model identified mechanisms by which oncogenic KRAS promotes immune evasion, most notably via upregulation of immunosuppressive COX2 in cancer cells. Oncogenic KRAS potently induced COX2 in both mouse and human lung cancer, which was suppressed using KRAS inhibitors. COX2 acted via prostaglandin E2 (PGE2) to promote resistance to immune checkpoint blockade (ICB) in lung adenocarcinoma. Targeting COX2/PGE2 remodeled the tumor microenvironment by inducing proinflammatory polarization of myeloid cells and influx of activated cytotoxic CD8+ T cells, which increased the efficacy of ICB. Restoration of COX2 expression contributed to tumor relapse after prolonged KRAS inhibition. These results provide the rationale for testing COX2/PGE2 pathway inhibitors in combination with KRASG12C inhibition or ICB in patients with KRAS-mutant lung cancer. Significance: COX2 signaling via prostaglandin E2 is a major mediator of immune evasion driven by oncogenic KRAS that promotes immunotherapy and KRAS-targeted therapy resistance, suggesting effective combination treatments for KRAS-mutant lung cancer.

Interaction of NF-κB and FOSL1 drives glioma stemness

Glioblastoma multiforme (GBM) is the most common and malignant primary brain tumor; GBM's inevitable recurrence suggests that glioblastoma stem cells (GSC) allow these tumors to persist. Our previous work showed that FOSL1, transactivated by the STAT3 gene, functions as a tumorigenic gene in glioma pathogenesis and acts as a diagnostic marker and potential drug target in glioma patients. Accumulating evidence shows that STAT3 and NF-κB cooperate to promote the development and progression of various cancers. The link between STAT3 and NF-κB suggests that NF-κB can also transcriptionally regulate FOSL1 and contribute to gliomagenesis. To investigate downstream molecules of FOSL1, we analyzed the transcriptome after overexpressing FOSL1 in a PDX-L14 line characterized by deficient FOSL1 expression. We then conducted immunohistochemical staining for FOSL1 and NF-κB p65 using rabbit polyclonal anti-FOSL1 and NF-κB p65 in glioma tissue microarrays (TMA) derived from 141 glioma patients and 15 healthy individuals. Next, mutants of the human FOSL1 promoter, featuring mutations in essential binding sites for NF-κB were generated using a Q5 site-directed mutagenesis kit. Subsequently, we examined luciferase activity in glioma cells and compared it to the wild-type FOSL1 promoter. Then, we explored the mutual regulation between NF-κB signaling and FOSL1 by modulating the expression of NF-κB or FOSL1. Subsequently, we assessed the activity of FOSL1 and NF-κB. To understand the role of FOSL1 in cell growth and stemness, we conducted a CCK-8 assay and cell cycle analysis, assessing apoptosis and GSC markers, ALDH1, and CD133 under varying FOSL1 expression conditions. Transcriptome analyses of downstream molecules of FOSL1 show that NF-κB signaling pathway is regulated by FOSL1. NF-κB p65 protein expression correlates to the expression of FOSL1 in glioma patients, and both are associated with glioma grades. NF-κB is a crucial transcription factor activating the FOSL1 promoter in glioma cells. Mutual regulation between NF-κB and FOSL1 contributes to glioma tumorigenesis and stemness through promoting G1/S transition and inhibiting apoptosis. Therefore, the FOSL1 molecular pathway is functionally connected to NF-κB activation, enhances stemness, and is indicative that FOSL1 may potentially be a novel GBM drug target.

Defining the KRAS- and ERK-dependent transcriptome in KRAS-mutant cancers

How the KRAS oncogene drives cancer growth remains poorly understood. Therefore, we established a systemwide portrait of KRAS- and extracellular signal-regulated kinase (ERK)-dependent gene transcription in KRAS-mutant cancer to delineate the molecular mechanisms of growth and of inhibitor resistance. Unexpectedly, our KRAS-dependent gene signature diverges substantially from the frequently cited Hallmark KRAS signaling gene signature, is driven predominantly through the ERK mitogen-activated protein kinase (MAPK) cascade, and accurately reflects KRAS- and ERK-regulated gene transcription in KRAS-mutant cancer patients. Integration with our ERK-regulated phospho- and total proteome highlights ERK deregulation of the anaphase promoting complex/cyclosome (APC/C) and other components of the cell cycle machinery as key processes that drive pancreatic ductal adenocarcinoma (PDAC) growth. Our findings elucidate mechanistically the critical role of ERK in driving KRAS-mutant tumor growth and in resistance to KRAS-ERK MAPK targeted therapies.

Determining the ERK-regulated phosphoproteome driving KRAS-mutant cancer

To delineate the mechanisms by which the ERK1 and ERK2 mitogen-activated protein kinases support mutant KRAS-driven cancer growth, we determined the ERK-dependent phosphoproteome in KRAS-mutant pancreatic cancer. We determined that ERK1 and ERK2 share near-identical signaling and transforming outputs and that the KRAS-regulated phosphoproteome is driven nearly completely by ERK. We identified 4666 ERK-dependent phosphosites on 2123 proteins, of which 79 and 66%, respectively, were not previously associated with ERK, substantially expanding the depth and breadth of ERK-dependent phosphorylation events and revealing a considerably more complex function for ERK in cancer. We established that ERK controls a highly dynamic and complex phosphoproteome that converges on cyclin-dependent kinase regulation and RAS homolog guanosine triphosphatase function (RHO GTPase). Our findings establish the most comprehensive molecular portrait and mechanisms by which ERK drives KRAS-dependent pancreatic cancer growth.

Trametinib sensitizes KRAS-mutant lung adenocarcinoma tumors to PD-1/PD-L1 axis blockade via Id1 downregulation

BACKGROUND

The identification of novel therapeutic strategies to overcome resistance to the MEK inhibitor trametinib in mutant KRAS lung adenocarcinoma (LUAD) is a challenge. This study analyzes the effects of trametinib on Id1 protein, a key factor involved in the KRAS oncogenic pathway, and investigates the role of Id1 in the acquired resistance to trametinib as well as the synergistic anticancer effect of trametinib combined with immunotherapy in KRAS-mutant LUAD.

METHODS

We evaluated the effects of trametinib on KRAS-mutant LUAD by Western blot, RNA-seq and different syngeneic mouse models. Genetic modulation of Id1 expression was performed in KRAS-mutant LUAD cells by lentiviral or retroviral transductions of specific vectors. Cell viability was assessed by cell proliferation and colony formation assays. PD-L1 expression and apoptosis were measured by flow cytometry. The anti-tumor efficacy of the combined treatment with trametinib and PD-1 blockade was investigated in KRAS-mutant LUAD mouse models, and the effects on the tumor immune infiltrate were analyzed by flow cytometry and immunohistochemistry.

RESULTS

We found that trametinib activates the proteasome-ubiquitin system to downregulate Id1 in KRAS-mutant LUAD tumors. Moreover, we found that Id1 plays a major role in the acquired resistance to trametinib treatment in KRAS-mutant LUAD cells. Using two preclinical syngeneic KRAS-mutant LUAD mouse models, we found that trametinib synergizes with PD-1/PD-L1 blockade to hamper lung cancer progression and increase survival. This anti-tumor activity depended on trametinib-mediated Id1 reduction and was associated with a less immunosuppressive tumor microenvironment and increased PD-L1 expression on tumor cells.

CONCLUSIONS

Our data demonstrate that Id1 expression is involved in the resistance to trametinib and in the synergistic effect of trametinib with anti-PD-1 therapy in KRAS-mutant LUAD tumors. These findings suggest a potential therapeutic approach for immunotherapy-refractory KRAS-mutant lung cancers.

High somatic mutations in circulating tumor DNA predict response of metastatic pancreatic ductal adenocarcinoma to first-line nab-paclitaxel plus S-1: prospective study

AIMS

We previously showed that the nab-paclitaxel plus S-1 (NPS) regimen had promising effects against metastatic pancreatic ducal adenocarcinoma (mPDAC), whose efficacy however could not be precisely predicted by routine biomarkers. This prospective study aimed to investigate the values of mutations in circulating tumor DNA (ctDNA) and their dynamic changes in predicting response of mPDAC to NPS chemotherapy.

METHODS

Paired tumor tissue and blood samples were prospectively collected from patients with mPDAC receiving first-line NPS chemotherapy, and underwent next-generation sequencing with genomic profiling of 425 genes for ctDNA. High mutation allelic frequency (MAF) was defined as ≥ 30% and ≥ 5% in tumor tissue and blood, respectively. Kappa statistics were used to assess agreement between mutant genes in tumor and ctDNA. Associations of mutations in ctDNA and their dynamic changes with tumor response, overall survival (OS), and progression-free survival (PFS) were assessed using the Kaplan-Meier method, multivariable-adjusted Cox proportional hazards regression, and longitudinal data analysis.

RESULTS

147 blood samples and 43 paired tumor specimens from 43 patients with mPDAC were sequenced. The most common driver genes with high MAF were KRAS (tumor, 35%; ctDNA, 37%) and TP53 (tumor, 37%; ctDNA, 33%). Mutation rates of KRAS and TP53 in ctDNA were significantly higher in patients with liver metastasis, with baseline CA19-9 ≥ 2000 U/mL, and/or without an early CA19-9 response. κ values for the 5 most commonly mutated genes between tumor and ctDNA ranged from 0.48 to 0.76. MAFs of the genes mostly decreased sequentially during subsequent measurements, which significantly correlated with objective response, with an increase indicating cancer progression. High mutations of KRAS and ARID1A in both tumor and ctDNA, and of TP53, CDKN2A, and SMAD4 in ctDNA but not in tumor were significantly associated with shorter survival. When predicting 6-month OS, AUCs for the 5 most commonly mutated genes in ctDNA ranged from 0.59 to 0.84, larger than for genes in tumor (0.56 to 0.71) and for clinicopathologic characteristics (0.51 to 0.68). Repeated measurements of mutations in ctDNA significantly differentiated survival and tumor response. Among the 31 patients with ≥ 2 ctDNA tests, longitudinal analysis of changes in gene MAF showed that ctDNA progression was 60 and 58 days ahead of radiologic and CA19-9 progression for 48% and 42% of the patients, respectively.

CONCLUSIONS

High mutations of multiple driving genes in ctDNA and their dynamic changes could effectively predict response of mPDAC to NPS chemotherapy, with promising reliable predictive performance superior to routine clinicopathologic parameters. Inspiringly, longitudinal ctDNA tracking could predict disease progression about 2 months ahead of radiologic or CA19-9 evaluations, with the potential to precisely devise individualized therapeutic strategies for mPDAC.

Cost and time-efficient construction of a 3'-end mRNA library from unpurified bulk RNA in a single tube

The major drawbacks of RNA sequencing (RNA-seq), a remarkably accurate transcriptome profiling method, is its high cost and poor scalability. Here, we report a highly scalable and cost-effective method for transcriptomics profiling called Bulk transcriptOme profiling of cell Lysate in a single poT (BOLT-seq), which is performed using unpurified bulk 3'-end mRNA in crude cell lysates. During BOLT-seq, RNA/DNA hybrids are directly subjected to tagmentation, and second-strand cDNA synthesis and RNA purification are omitted, allowing libraries to be constructed in 2 h of hands-on time. BOLT-seq was successfully used to cluster small molecule drugs based on their mechanisms of action and intended targets. BOLT-seq competes effectively with alternative library construction and transcriptome profiling methods.

TEAD Inhibition Overcomes YAP1/TAZ-Driven Primary and Acquired Resistance to KRASG12C Inhibitors

Primary/intrinsic and treatment-induced acquired resistance limit the initial response rate to and long-term efficacy of direct inhibitors of the KRASG12C mutant in cancer. To identify potential mechanisms of resistance, we applied a CRISPR/Cas9 loss-of-function screen and observed loss of multiple components of the Hippo tumor suppressor pathway, which acts to suppress YAP1/TAZ-regulated gene transcription. YAP1/TAZ activation impaired the antiproliferative and proapoptotic effects of KRASG12C inhibitor (G12Ci) treatment in KRASG12C-mutant cancer cell lines. Conversely, genetic suppression of YAP1/WWTR1 (TAZ) enhanced G12Ci sensitivity. YAP1/TAZ activity overcame KRAS dependency through two distinct TEAD transcription factor-dependent mechanisms, which phenocopy KRAS effector signaling. First, TEAD stimulated ERK-independent transcription of genes normally regulated by ERK (BIRC5, CDC20, ECT2, FOSL1, and MYC) to promote progression through the cell cycle. Second, TEAD caused activation of PI3K-AKT-mTOR signaling to overcome apoptosis. G12Ci treatment-induced acquired resistance was also caused by YAP1/TAZ-TEAD activation. Accordingly, concurrent treatment with pharmacologic inhibitors of TEAD synergistically enhanced KRASG12C inhibitor antitumor activity in vitro and prolonged tumor suppression in vivo. In summary, these observations reveal YAP1/TAZ-TEAD signaling as a crucial driver of primary and acquired resistance to KRAS inhibition and support the use of TEAD inhibitors to enhance the antitumor efficacy of KRAS-targeted therapies.

SIGNIFICANCE

YAP1/TAZ-TEAD activation compensates for loss of KRAS effector signaling, establishing a mechanistic basis for concurrent inhibition of TEAD to enhance the efficacy of KRASG12C-selective inhibitor treatment of KRASG12C-mutant cancers. See related commentary by Johnson and Haigis, p. 4005.

PARP7-mediated ADP-ribosylation of FRA1 promotes cancer cell growth by repressing IRF1- and IRF3-dependent apoptosis

PARP7 was reported to promote tumor growth in a cell-autonomous manner and by repressing the antitumor immune response. Nevertheless, the molecular mechanism of how PARP7-mediated ADP-ribosylation exerts these effects in cancer cells remains elusive. Here, we identified PARP7 as a nuclear and cysteine-specific mono-ADP-ribosyltransferase that modifies targets critical for regulating transcription, including the AP-1 transcription factor FRA1. Loss of FRA1 ADP-ribosylation via PARP7 inhibition by RBN-2397 or mutation of the ADP-ribosylation site C97 increased FRA1 degradation by the proteasome via PSMC3. The reduction in FRA1 protein levels promoted IRF1- and IRF3-dependent cytokine as well as proapoptotic gene expression, culminating in CASP8-mediated apoptosis. Furthermore, high PARP7 expression was indicative of the PARP7 inhibitor response in FRA1-positive lung and breast cancer cells. Collectively, our findings highlight the connected roles of PARP7 and FRA1 and emphasize the clinical potential of PARP7 inhibitors for FRA1-driven cancers.

Constructing an APOBEC-related gene signature with predictive value in the overall survival and therapeutic sensitivity in lung adenocarcinoma

Background

APOBEC family play an important role in cancer mutagenesis and tumor development. The role of APOBEC family in lung adenocarcinoma (LUAD) has not been studied comprehensively.

Materials and methods

The expression data of pan-cancer as well as LUAD was obtained from public databases. The expression level of APOBEC family genes was analyzed in different normal and cancer tissues. APOBEC mutagenesis enrichment score (AMES) was utilized to evaluate the APOBEC-induced mutations and the relation of APOBEC with genomic instability. Gene set enrichment analysis was used to identify differentially enriched pathways. Univariate Cox regression and Lasso regression were applied to screen key prognostic genes. The immune cell infiltration was estimated by CIBERSORT. RT-qPCR assay, CCK-8 and Transwell assay were conducted to explore gene expression and lung cancer cell invasion.

Results

Cancer tissues had significantly altered expression of APOBEC family genes and the expression patterns of APOBEC family were different in different cancer types. APOBEC3B was the most aberrantly expressed in most cancer types. In LUAD, we observed a significantly positive correlation of AMES with intratumor heterogeneity (ITH), tumor neoantigen burden (TNB), and tumor mutation burden (TMB). High AMES group had high mutation counts of DNA damage repair pathways, and high enrichment of cell cycle and DNA repair pathways. We identified four prognostic genes (LYPD3, ANLN, MUC5B, and FOSL1) based on AMES, and constructed an AMES-related gene signature. The expressions of four genes were enhanced and accelerated the invasion ability and viability of lung cancer cells. Furthermore, we found that high group increased oxidative stress level.

Conclusions

APOBEC family was associated with genomic instability, DNA damage-related pathways, and cell cycle in LUAD. The AMES-related gene signature had a great potential to indicate the prognosis and guide immunotherapy/chemotherapy for patients suffering from LUAD.

Signature-driven repurposing of Midostaurin for combination with MEK1/2 and KRASG12C inhibitors in lung cancer

Drug combinations are key to circumvent resistance mechanisms compromising response to single anti-cancer targeted therapies. The implementation of combinatorial approaches involving MEK1/2 or KRASG12C inhibitors in the context of KRAS-mutated lung cancers focuses fundamentally on targeting KRAS proximal activators or effectors. However, the antitumor effect is highly determined by compensatory mechanisms arising in defined cell types or tumor subgroups. A potential strategy to find drug combinations targeting a larger fraction of KRAS-mutated lung cancers may capitalize on the common, distal gene expression output elicited by oncogenic KRAS. By integrating a signature-driven drug repurposing approach with a pairwise pharmacological screen, here we show synergistic drug combinations consisting of multi-tyrosine kinase PKC inhibitors together with MEK1/2 or KRASG12C inhibitors. Such combinations elicit a cytotoxic response in both in vitro and in vivo models, which in part involves inhibition of the PKC inhibitor target AURKB. Proteome profiling links dysregulation of MYC expression to the effect of both PKC inhibitor-based drug combinations. Furthermore, MYC overexpression appears as a resistance mechanism to MEK1/2 and KRASG12C inhibitors. Our study provides a rational framework for selecting drugs entering combinatorial strategies and unveils MEK1/2- and KRASG12C-based therapies for lung cancer.

Actionable Genetic Screens Unveil Targeting of AURKA, MEK, and Fatty Acid Metabolism as an Alternative Therapeutic Approach for Advanced Melanoma

Despite the remarkable improvements achieved in the management of metastatic melanoma, there are still unmet clinical needs. A considerable fraction of patients does not respond to immune and/or targeted therapies owing to primary and acquired resistance, high-grade immune-related adverse events, and a lack of alternative treatment options. To design effective combination therapies, we set up a functional ex vivo preclinical assay on the basis of a drop-out genetic screen in metastatic melanoma patient-derived xenografts. We showed that this approach can be used to isolate actionable vulnerabilities predictive of drug efficacy. In particular, we highlighted that the dual targeting of AURKA and MAPK/extracellular signal-regulated kinase kinase employing the combination of alisertib and trametinib is highly effective in a cohort of metastatic melanoma patient-derived xenografts, both ex vivo and in vivo. Alisertib and trametinib combination therapy outperforms standard-of-care therapy in both BRAF-mutant patient-derived xenografts and targeted therapy-resistant models. Furthermore, alisertib and trametinib treatment modulates several critical cancer pathways, including an early metabolic reprogramming that leads to the transcriptional upregulation of the fatty acid oxidation pathway. This acquired trait unveiled an additional point of intervention for pharmacological targeting, and indeed, the triple combination of alisertib and trametinib with the fatty acid oxidation inhibitor etomoxir proved to be further beneficial, inducing tumor regression and remarkably prolonging the overall survival of the mice.

Shear Stress Drives the Cleavage Activation of Protease-Activated Receptor 2 by PRSS3/Mesotrypsin to Promote Invasion and Metastasis of Circulating Lung Cancer Cells

When circulating tumor cells (CTCs) travel in circulation, they can be killed by detachment-induced anoikis and fluidic shear stress (SS)-mediated apoptosis. Circulatory treatment, which can make CTCs detached but also generate SS, can increase metastasis of cancer cells. To identify SS-specific mechanosensors without detachment impacts, a microfluidic circulatory system is used to generate arteriosus SS and compare transcriptome profiles of circulating lung cancer cells with suspended cells. Half of the cancer cells can survive SS damage and show higher invasion ability. Mesotrypsin (PRSS3), protease-activated receptor 2 (PAR2), and the subunit of activating protein 1, Fos-related antigen 1 (FOSL1), are upregulated by SS, and their high expression is responsible for promoting invasion and metastasis. SS triggers PRSS3 to cleave the N-terminal inhibitory domain of PAR2 within 2 h. As a G protein-coupled receptor, PAR2 further activates the Gαi protein to turn on the Src-ERK/p38/JNK-FRA1/cJUN axis to promote the expression of epithelial-mesenchymal transition markers, and also PRSS3, which facilitates metastasis. Enriched PRSS3, PAR2, and FOSL1 in human tumor samples and their correlations with worse outcomes reveal their clinical significance. PAR2 may serve as an SS-specific mechanosensor cleavable by PRSS3 in circulation, which provides new insights for targeting metastasis-initiating CTCs.

scKINETICS: inference of regulatory velocity with single-cell transcriptomics data

MOTIVATION

Transcriptional dynamics are governed by the action of regulatory proteins and are fundamental to systems ranging from normal development to disease. RNA velocity methods for tracking phenotypic dynamics ignore information on the regulatory drivers of gene expression variability through time.

RESULTS

We introduce scKINETICS (Key regulatory Interaction NETwork for Inferring Cell Speed), a dynamical model of gene expression change which is fit with the simultaneous learning of per-cell transcriptional velocities and a governing gene regulatory network. Fitting is accomplished through an expectation-maximization approach designed to learn the impact of each regulator on its target genes, leveraging biologically motivated priors from epigenetic data, gene-gene coexpression, and constraints on cells' future states imposed by the phenotypic manifold. Applying this approach to an acute pancreatitis dataset recapitulates a well-studied axis of acinar-to-ductal transdifferentiation whilst proposing novel regulators of this process, including factors with previously appreciated roles in driving pancreatic tumorigenesis. In benchmarking experiments, we show that scKINETICS successfully extends and improves existing velocity approaches to generate interpretable, mechanistic models of gene regulatory dynamics.

AVAILABILITY AND IMPLEMENTATION

All python code and an accompanying Jupyter notebook with demonstrations are available at http://github.com/dpeerlab/scKINETICS.

Stepwise activation of p63 and the MEK/ERK pathway induces the expression of ARL4C to promote oral squamous cell carcinoma cell proliferation

Carcinogenesis is a multistep process wherein cells accumulate multiple genetic alterations and progress to a more malignant phenotype. It has been proposed that sequential accumulation of gene abnormalities in specific genes drives the transition from non-tumorous epithelia through a preneoplastic lesion/benign tumor to cancer. Histologically, oral squamous cell carcinoma (OSCC) progresses in multiple ordered steps that begin with mucosal epithelial cell hyperplasia, which is followed by dysplasia, carcinoma in situ and invasive carcinoma. It is therefore hypothesized that genetic alteration-mediated multistep carcinogenesis would be involved in the development of OSCC; however, the detailed molecular mechanisms are unknown. We clarified the comprehensive gene expression patterns and carried out an enrichment analysis using DNA microarray data from a pathological specimen of OSCC (including a non-tumor region, carcinoma in situ lesion and invasive carcinoma lesion). The expression of numerous genes and signal activation were altered in the development of OSCC. Among these, the p63 expression was increased and the MEK/ERK-MAPK pathway was activated in carcinoma in situ lesion and in invasive carcinoma lesion. Immunohistochemical analyses revealed that p63 was initially upregulated in carcinoma in situ and ERK was sequentially activated in invasive carcinoma lesions in OSCC specimens. ADP-ribosylation factor (ARF)-like 4c (ARL4C), the expression of which is reportedly induced by p63 and/or the MEK/ERK-MAPK pathway in OSCC cells, has been shown to promote tumorigenesis. Immunohistochemically, in OSCC specimens, ARL4C was more frequently detected in tumor lesions, especially in invasive carcinoma lesions, than in carcinoma in situ lesions. Additionally, ARL4C and phosphorylated ERK were frequently merged in invasive carcinoma lesions. Loss-of-function experiments using inhibitors and siRNAs revealed that p63 and MEK/ERK-MAPK cooperatively induce the expression of ARL4C and cell growth in OSCC cells. These results suggest that the stepwise activation of p63 and MEK/ERK-MAPK contributes to OSCC tumor cell growth through regulation of ARL4C expression.

Erianin inhibits the growth and metastasis through autophagy-dependent ferroptosis in KRASG13D colorectal cancer

Colorectal cancer (CRC) is the third most common cause of cancer mortality worldwide. Approximately 40% of CRC patients are KRAS sequence variation, including KRAS G13D mutation (KRAS^G13D) CRC patients, accounting for approximately 8% of all KRAS mutations in CRC patients and showing little benefit from anti-EGFR therapy. Therefore, there is an urgent need for new and effective anticancer agents in patients with KRAS^G13D CRC. Here, we identified a natural product, erianin, that directly interacted with purified recombinant human KRAS^G13D with a Kd of 1.1163 μM, which also significantly improve the thermal stability of KRAS^G13D. The cell viability assay showed that KRAS^G13D cells were more sensitive to erianin than KRAS^WT or KRAS^G12V cells. In vitro, results showed that erianin suppressed the migration, invasion and epithelial-mesenchymal transition (EMT) of KRAS^G13D CRC cells. Furthermore, erianin induced ferroptosis, as evidenced by the accumulation of Fe²⁺ and reactive oxygen species (ROS), lipid peroxidation, and changes in the mitochondrial morphology of KRAS^G13D CRC cells. Interestingly, we also found that erianin-induced ferroptosis was accompanied by autophagy. Moreover, the occurrence of erianin-induced ferroptosis is reversed by autophagy inhibitors (NH4Cl and Bafilomycin A1) and ATG5 knockdown, suggesting that erianin-induced ferroptosis is autophagy-dependent. In addition, we evaluated the inhibition of tumor growth and metastasis by erianin in vivo using a subcutaneous tumor model and a spleen-liver metastasis model, respectively. Collectively, these data provide novel insights into the anticancer activity of erianin, which is valuable for the further discussion and investigation of the use of erianin in clinical anticancer chemotherapy for KRAS^G13D CRC.

The phospholipid transporter PITPNC1 links KRAS to MYC to prevent autophagy in lung and pancreatic cancer

BACKGROUND

The discovery of functionally relevant KRAS effectors in lung and pancreatic ductal adenocarcinoma (LUAD and PDAC) may yield novel molecular targets or mechanisms amenable to inhibition strategies. Phospholipids availability has been appreciated as a mechanism to modulate KRAS oncogenic potential. Thus, phospholipid transporters may play a functional role in KRAS-driven oncogenesis. Here, we identified and systematically studied the phospholipid transporter PITPNC1 and its controlled network in LUAD and PDAC.

METHODS

Genetic modulation of KRAS expression as well as pharmacological inhibition of canonical effectors was completed. PITPNC1 genetic depletion was performed in in vitro and in vivo LUAD and PDAC models. PITPNC1-deficient cells were RNA sequenced, and Gene Ontology and enrichment analyses were applied to the output data. Protein-based biochemical and subcellular localization assays were run to investigate PITPNC1-regulated pathways. A drug repurposing approach was used to predict surrogate PITPNC1 inhibitors that were tested in combination with KRASG12C inhibitors in 2D, 3D, and in vivo models.

RESULTS

PITPNC1 was increased in human LUAD and PDAC, and associated with poor patients' survival. PITPNC1 was regulated by KRAS through MEK1/2 and JNK1/2. Functional experiments showed PITPNC1 requirement for cell proliferation, cell cycle progression and tumour growth. Furthermore, PITPNC1 overexpression enhanced lung colonization and liver metastasis. PITPNC1 regulated a transcriptional signature which highly overlapped with that of KRAS, and controlled mTOR localization via enhanced MYC protein stability to prevent autophagy. JAK2 inhibitors were predicted as putative PITPNC1 inhibitors with antiproliferative effect and their combination with KRASG12C inhibitors elicited a substantial anti-tumour effect in LUAD and PDAC.

CONCLUSIONS

Our data highlight the functional and clinical relevance of PITPNC1 in LUAD and PDAC. Moreover, PITPNC1 constitutes a new mechanism linking KRAS to MYC, and controls a druggable transcriptional network for combinatorial treatments.

Comprehensive Analyses Reveal Effects on Tumor Immune Infiltration and Immunotherapy Response of APOBEC Mutagenesis and Its Molecular Mechanisms in Esophageal Squamous Cell Carcinoma

The apolipoprotein B mRNA editing enzyme catalytic polypeptide (APOBEC) mutagenesis is prevalent in esophageal squamous cell carcinoma (ESCC). However, the functional role of APOBEC mutagenesis has yet to be fully delineated. To address this, we collect matched multi-omics data of 169 ESCC patients and evaluate characteristics of immune infiltration using multiple bioinformatic approaches based on bulk and single-cell RNA sequencing (scRNA-seq) data and verified by functional assays. We find that APOBEC mutagenesis prolongs overall survival (OS) of ESCC patients. The reason for this outcome is probably due to high anti-tumor immune infiltration, immune checkpoints expression and immune related pathway enrichment, such as interferon (IFN) signaling, innate and adaptive immune system. The elevated AOBEC3A (A3A) activity paramountly contributes to the footprints of APOBEC mutagenesis and is first discovered to be transactivated by FOSL1. Mechanistically, upregulated A3A exacerbates cytosolic double-stranded DNA (dsDNA) accumulation, thus stimulating cGAS-STING pathway. Simultaneously, A3A is associated with immunotherapy response which is predicted by TIDE algorithm, validated in a clinical cohort and further confirmed in mouse models. These findings systematically elucidate the clinical relevance, immunological characteristics, prognostic value for immunotherapy and underlying mechanisms of APOBEC mutagenesis in ESCC, which demonstrate great potential in clinical utility to facilitate clinical decisions.

FRA-1 as a Regulator of EMT and Metastasis in Breast Cancer

Among FOS-related components of the dimeric AP-1 transcription factor, the oncoprotein FRA-1 (encoded by FOSL1) is a key regulator of invasion and metastasis. The well-established FRA-1 pro-invasive activity in breast cancer, in which FOSL1 is overexpressed in the TNBC (Triple Negative Breast Cancer)/basal subtypes, correlates with the FRA-1-dependent transcriptional regulation of EMT (Epithelial-to-Mesenchymal Transition). After summarizing the major findings on FRA-1 in breast cancer invasiveness, we discuss the FRA-1 mechanistic links with EMT and cancer cell stemness, mediated by transcriptional and posttranscriptional interactions between FOSL1/FRA-1 and EMT-regulating transcription factors, miRNAs, RNA binding proteins and cytokines, along with other target genes involved in EMT. In addition to the FRA-1/AP-1 effects on the architecture of target promoters, we discuss the diagnostic and prognostic significance of the EMT-related FRA-1 transcriptome, along with therapeutic implications. Finally, we consider several novel perspectives regarding the less explored roles of FRA-1 in the tumor microenvironment and in control of the recently characterized hybrid EMT correlated with cancer cell plasticity, stemness, and metastatic potential. We will also examine the application of emerging technologies, such as single-cell analyses, along with animal models of TNBC and tumor-derived CTCs and PDXs (Circulating Tumor Cells and Patient-Derived Xenografts) for studying the FRA-1-mediated mechanisms in in vivo systems of EMT and metastasis.

TPX2 expression as a negative predictor of gemcitabine efficacy in pancreatic cancer

BACKGROUND

Targeting protein for Xenopus kinesin-like protein 2 (TPX2) overexpression in human tumours is associated with increased malignancy. Its effect on gemcitabine resistance in pancreatic ductal adenocarcinoma (PDAC) has not been studied yet.

METHODS

The prognostic impact of TPX2 expression was examined in the tumour tissue of 139 patients with advanced PDAC (aPDAC) treated within the AIO-PK0104 trial or translational trials and of 400 resected PDAC (rPDAC) patients. The findings were validated using RNAseq data of 149 resected PDAC patients.

RESULTS

In the aPDAC cohorts, 13.7% of all samples showed high TPX2 expression, conferring significantly shorter progression-free survival (PFS, HR 5.25, P < 0.001) and overall survival times (OS, HR 4.36, P < 0.001) restricted to gemcitabine-based treated patients (n = 99). In the rPDAC cohort, 14.5% of all samples showed high TPX2 expression, conferring significantly shorter disease-free survival times (DFS, HR 2.56, P < 0.001) and OS times (HR 1.56, P = 0.04) restricted to patients treated with adjuvant gemcitabine. RNAseq data from the validation cohort confirmed the findings.

CONCLUSIONS

High TPX2 expression may serve as a negative predictor of gemcitabine-based palliative and adjuvant chemotherapy in PDAC and could be used to inform clinical therapy decisions.

CLINICAL TRIAL REGISTRY

The clinical trial registry identifier is NCT00440167.

FOSL1 transcriptionally regulates PHLDA2 to promote 5-FU resistance in colon cancer cells

BACKGROUND

Tumor drug resistance is a leading cause of tumor treatment failure. To date, the association between FOS-Like antigen-1 (FOSL1) and chemotherapy sensitivity in colon cancer is unclear. The present study investigated the molecular mechanism of FOSL1 regulating 5-Fluorouracil (5-FU) resistance in colon cancer.

METHODS

FOSL1 expression in colon cancer was analyzed by bioinformatics methods, and its downstream regulatory factors were predicted. Pearson correlation analyzed the expression of FOSL1 and downstream regulatory gene. Meanwhile, the expression of FOSL1 and its downstream factor Pleckstrin Homology-Like Domain Family A Member 2 (PHLDA2) in colon cancer cell lines was measured by qRT-PCR and western blot. The regulatory relationship between FOSL1 and PHLDA2 was verified by chromatin immunoprecipitation (ChIP) assay and dual-luciferase reporter assay. The effects of the FOSL1/PHLDA2 axis on the resistance in colon cancer cells to 5-FU were analyzed by cell experiments.

RESULTS

FOSL1 expression was evidently up-regulated in colon cancer and 5-FU resistant cells. FOSL1 was positively correlated with PHLDA2 in colon cancer. In vitro cell assays showed that low expression of FOSL1 significantly enhanced 5-FU sensitivity in colon cancer cells, significantly suppressed the proliferation of cancer cells, and induced apoptosis. Overexpression of FOSL1 presented the opposite regulatory trend. Mechanistically, FOSL1 activated PHLDA2 and up-regulated its expression. Moreover, by activating glycolysis, PHLDA2 promoted 5-Fu resistance and cell proliferation, and reduced cell apoptosis in colon cancer.

CONCLUSION

Down-regulated FOSL1 expression could enhance the 5-FU sensitivity of colon cancer cells, and FOSL1/PHLDA2 axis may be an effective target for overcoming chemotherapy resistance in colon cancer.

Upregulated GPRC5A disrupting the Hippo pathway promotes the proliferation and migration of pancreatic cancer cells via the cAMP-CREB axis

BACKGROUND

Pancreatic cancer has a high mortality rate worldwide, and is predicted to be third leading cause of death in the near future. However, the regulatory mechanisms that inhibit the progression of pancreatic cancer remain elusive. Currently, exploring the function and mechanisms of GPCRs (G-protein coupled receptors) is an important way to discover promising therapeutic targets for cancer.

METHODS

GPRC5A expression was measured using real-time quantitative PCR, immunohistochemistry and western blot assays. Cell proliferation and migration were assessed using CCK-8, clone formation, wound-healing and transwell assays. A cytosolic/nuclear distribution experiment was used to detect the protein location transfer. A xenograft model of pancreatic cancer was established to explore the role of GPRC5A in vivo.

RESULTS

GPRC5A expression was increased in pancreatic cancer, and disruption of GPRC5A expression inhibited tumor growth in vivo. Mechanistically, GPRC5A positively regulated the transcription of YAP1 through cAMP-CREB signaling. Moreover, we show that the proliferation and migration induced by GPRC5A in pancreatic cancer could be rescued by inhibiting YAP1 expression.

CONCLUSIONS

GPRC5A interacts with the Hippo pathway to promote the progression of pancreatic cancer. These findings reveal an important crosstalk model and provide potential targets for pancreatic cancer therapy.

Identification of two molecular subtypes and a novel prognostic model of lung adenocarcinoma based on a cuproptosis-associated gene signature

Lung adenocarcinoma is the most common subtype of lung cancer clinically, with high mortality and poor prognosis. Cuproptosis present a newly discovered mode of cell death characterized by aggregation of fatty acylated proteins, depletion of iron-sulfur clusterin, triggering of HSP70, and induction of intracellular toxic oxidative stress. However, the impact of cuproptosis on lung adenocarcinoma development, prognosis, and treatment has not been elucidated. By systematically analyzing the genetic alterations of 10 cuproptosis-related genes in lung adenocarcinoma, we found that CDKN2A, DLAT, LIAS, PDHA1, FDX1, GLS, and MTF1 were differentially expressed between lung cancer tissues and adjacent tissues. Based on the expression levels of 10 cuproptosis-related genes, we classified lung adenocarcinoma patients into two molecular subtypes using the Consensus clustering method, of which subtype 2 had a worse prognosis. Differential expression genes associated with prognosis between the two subtypes were obtained by differential analysis and survival analysis, and cox lasso regression was applied to construct a cuproptosis-related prognostic model. Its survival predicting ability was validated in three extrinsic validation cohorts. The results of multivariate cox analysis indicated that cuproptosis risk score was an independent prognostic predictor, and the mixed model formed by cupproptosis prognostic model combined with stage had more robust prognostic prediction accuracy. We found the differences in cell cycle, mitosis, and p53 signaling pathways between high- and low-risk groups according to GO and KEGG enrichment analysis. The results of immune microenvironment analysis showed that the enrichment score of activated dendritic cells, mast cells, and type 2 interferon response were down-regulated in the high-risk group, while the fraction of neutrophils and M0 macrophages were upregulated in the high-risk group. Compared with the high-risk group, subjects in the low-risk group had higher Immunophenoscore and may be more sensitive to immunotherapy. We identified seven chemotherapy agents may improve the curative effect in LUAD samples with higher risk score. Overall, we discovered that cuproptosis is closely related to the occurrence, prognosis, and treatment of lung adenocarcinoma. The cuproptosis prognostic model is a potential prognostic predictor and may provide new strategies for precision therapy in lung adenocarcinoma.

AP1/Fra1 confers resistance to MAPK cascade inhibition in pancreatic cancer

Targeting KRAS downstream signaling remains an important therapeutic approach in pancreatic cancer. We used primary pancreatic ductal epithelial cells and mouse models allowing the conditional expression of oncogenic Kras^G12D, to investigate KRAS signaling integrators. We observed that the AP1 family member FRA1 is tightly linked to the KRAS signal and expressed in pre-malignant lesions and the basal-like subtype of pancreatic cancer. However, genetic-loss-of-function experiments revealed that FRA1 is dispensable for Kras^G12D-induced pancreatic cancer development in mice. Using FRA1 gain- and loss-of-function models in an unbiased drug screen, we observed that FRA1 is a modulator of the responsiveness of pancreatic cancer to inhibitors of the RAF-MEK-ERK cascade. Mechanistically, context-dependent FRA1-associated adaptive rewiring of oncogenic ERK signaling was observed and correlated with sensitivity to inhibitors of canonical KRAS signaling. Furthermore, pharmacological-induced degradation of FRA1 synergizes with MEK inhibitors. Our studies establish FRA1 as a part of the molecular machinery controlling sensitivity to MAPK cascade inhibition allowing the development of mechanism-based therapies.

A Comprehensive Analysis of the Effects of Key Mitophagy Genes on the Progression and Prognosis of Lung Adenocarcinoma

The aim of our study was to perform a comprehensive analysis of the gene expression, copy number variation (CNV) and mutation of key mitophagy genes in the progression and prognosis of lung adenocarcinoma (LUAD). We obtained the data from The Cancer Genome Atlas (TCGA). Clustering analysis was performed to stratify the mitophagy related groups. The least absolute shrinkage and selection operator (LASSO) based cox model was used to select hub survival genes. An independent validation cohort was retrieved from Gene Expression Omnibus database. We found 24 out of 27 mitophagy genes were aberrantly expressed between tumor and normal samples. CNV gains were associated with higher expression of mitophagy genes in 23 of 27 mitophagy genes. The clustering analysis identified high and low risk mitophagy groups with distinct survival differences. The high risk mitophagy groups had higher tumor mutation burden, stemness phenotype, total CNVs and lower CD4+ T cells infiltration. Drugs targeted to high risk mitophagy groups were identified including the PI3K/AKT/mTOR inhibitor, HDAC inhibitor and chemotherapy agents such as cisplatin and gemcitabine. In addition, the differentially expressed genes (DEGs) were identified between mitophagy groups. Further univariate Cox analysis of each DEG and subsequent LASSO-based Cox model revealed a mitophagy-related prognostic signature. The risk score model of this signature showed a strong ability to predict the overall survival of LUAD patients in training and validation datasets. In conclusion, the mitophagy genes played an important role in the progression and prognosis of LUAD, which might provide useful information for the treatment of LUAD.

Brain metastatic outgrowth and osimertinib resistance are potentiated by RhoA in EGFR-mutant lung cancer

The brain is a major sanctuary site for metastatic cancer cells that evade systemic therapies. Through pre-clinical pharmacological, biological, and molecular studies, we characterize the functional link between drug resistance and central nervous system (CNS) relapse in Epidermal Growth Factor Receptor- (EGFR-) mutant non-small cell lung cancer, which can progress in the brain when treated with the CNS-penetrant EGFR inhibitor osimertinib. Despite widespread osimertinib distribution in vivo, the brain microvascular tumor microenvironment (TME) is associated with the persistence of malignant cell sub-populations, which are poised to proliferate in the brain as osimertinib-resistant lesions over time. Cellular and molecular features of this poised state are regulated through a Ras homolog family member A (RhoA) and Serum Responsive Factor (SRF) gene expression program. RhoA potentiates the outgrowth of disseminated tumor cells on osimertinib treatment, preferentially in response to extracellular laminin and in the brain. Thus, we identify pre-existing and adaptive features of metastatic and drug-resistant cancer cells, which are enhanced by RhoA/SRF signaling and the brain TME during the evolution of osimertinib-resistant disease.

DSTYK inhibition increases the sensitivity of lung cancer cells to T cell-mediated cytotoxicity

Lung cancer remains the leading cause of cancer-related death worldwide. We identify DSTYK, a dual serine/threonine and tyrosine non-receptor protein kinase, as a novel actionable target altered in non-small cell lung cancer (NSCLC). We also show DSTYK's association with a lower overall survival (OS) and poorer progression-free survival (PFS) in multiple patient cohorts. Abrogation of DSTYK in lung cancer experimental systems prevents mTOR-dependent cytoprotective autophagy, impairs lysosomal biogenesis and maturation, and induces accumulation of autophagosomes. Moreover, DSTYK inhibition severely affects mitochondrial fitness. We demonstrate in vivo that inhibition of DSTYK sensitizes lung cancer cells to TNF-α-mediated CD8+-killing and immune-resistant lung tumors to anti-PD-1 treatment. Finally, in a series of lung cancer patients, DSTYK copy number gain predicts lack of response to the immunotherapy. In summary, we have uncovered DSTYK as new therapeutic target in lung cancer. Prioritization of this novel target for drug development and clinical testing may expand the percentage of NSCLC patients benefiting from immune-based treatments.

Pharmacologic Targeting of TFIIH Suppresses KRAS-Mutant Pancreatic Ductal Adenocarcinoma and Synergizes with TRAIL

Pancreatic ductal adenocarcinoma (PDAC) typically presents as metastatic disease at diagnosis and remains refractory to treatment. Next-generation sequencing efforts have described the genomic landscape, classified molecular subtypes, and confirmed frequent alterations in major driver genes, with coexistent alterations in KRAS and TP53 correlating with the highest metastatic burden and poorest outcomes. However, translating this information to guide therapy remains a challenge. By integrating genomic analysis with an arrayed RNAi druggable genome screen and drug profiling of a KRAS/TP53 mutant PDAC cell line derived from a patient-derived xenograft (PDCL), we identified numerous targetable vulnerabilities that reveal both known and novel functional aspects of pancreatic cancer biology. A dependence on the general transcription and DNA repair factor TFIIH complex, particularly the XPB subunit and the CAK complex (CDK7/CyclinH/MAT1), was identified and further validated utilizing a panel of genomically subtyped KRAS mutant PDCLs. TFIIH function was inhibited with a covalent inhibitor of CDK7/12/13 (THZ1), a CDK7/CDK9 kinase inhibitor (SNS-032), and a covalent inhibitor of XPB (triptolide), which led to disruption of the protein stability of the RNA polymerase II subunit RPB1. Loss of RPB1 following TFIIH inhibition led to downregulation of key transcriptional effectors of KRAS-mutant signaling and negative regulators of apoptosis, including MCL1, XIAP, and CFLAR, initiating caspase-8 dependent apoptosis. All three drugs exhibited synergy in combination with a multivalent TRAIL, effectively reinforcing mitochondrial-mediated apoptosis. These findings present a novel combination therapy, with direct translational implications for current clinical trials on metastatic pancreatic cancer patients. Significance: This study utilizes functional genetic and pharmacological profiling of KRAS-mutant pancreatic adenocarcinoma to identify therapeutic strategies and finds that TFIIH inhibition synergizes with TRAIL to induce apoptosis in KRAS-driven pancreatic cancer.

FOSL1 promotes tumor growth and invasion in ameloblastoma

Background

FOSL1, a key component of the Activating protein-1 (AP-1) transcriptional complex, plays an important role in cancer cell migration, invasion, and proliferation. However, the impact of FOSL1 in ameloblastoma (AM) has not been clarified. Herein, we aimed to assess the expression of FOSL1 and investigate its functional role in AM.

Methods

The expression of FOSL1 was examined based on an immunohistochemistry analysis of 96 AM samples. Cell proliferation, migration, invasion, and tumorigenesis were assessed using Cell Counting Kit-8 (CCK-8), colony formation, Transwell, and sphere formation assays. RNA sequencing (RNA-seq) was employed to investigate the molecular alterations of AM cells upon FOSL depletion. Microarrays of AMs were downloaded from the Gene Expression Omnibus (GEO) database for bioinformatics analysis. In addition, patient-derived AM organoids were used to evaluate the therapeutic value of the AP-1 inhibitor.

Results

FOSL1 was detected in the nuclei of AMs and upregulated in conventional AMs compared to unicystic AMs and normal oral epithelium. Compared with primary AM, FOSL1 expression was significantly increased in recurrent AM. Genetic knockdown of FOSL1 suppressed the proliferation, migration, invasion, and sphere formation of AMs. Similar results were also observed by pharmacological inhibition of AP-1 activity. Moreover, the AP-1 inhibitor T5224 impeded the growth of organoids derived from AM patients. Mechanistically, our Ingenuity Pathway Analysis (IPA) and gene set enrichment analysis (GSEA) results revealed that depletion of FOSL1 inactivated kinetochore metaphase signaling and the epithelial–mesenchymal transition pathway and then impaired the aggressiveness of AM cells accordingly.

Conclusion

FOSL1 promotes tumor recurrence and invasive growth in AM by modulating kinetochore metaphase signaling and the epithelial–mesenchymal transition pathway; thus, it represents a promising therapeutic target for AM treatment.

MiR-33a targets FOSL1 and EN2 as a clinical prognostic marker for sarcopenia by glioma

To determine the relationship between glioma and muscle aging and to predict prognosis by screening for co-expressed genes, this study examined the relationship between glioma and sarcopenia. The study identified eight co-downregulated miRNAs, three co-upregulated miRNAs, and seven genes associated with overall glioma survival, namely, KRAS, IFNB1, ALCAM, ERBB2, STAT3, FOSL1, and EN2. With a multi-factor Cox regression model incorporating FOSL1 and EN2, we obtained ROC curves of 0.702 and 0.709, respectively, suggesting that glioma prognosis can be predicted by FOSL1 and EN2, which are differentially expressed in both cancer and aged muscle. FOSL1 and EN2 were analyzed using Gene Set Enrichment Analysis to identify possible functional pathways. RT-qPCR and a dual-luciferase reporter gene system verified that hsa-miR-33a targets FOSL1 and EN2. We found that hsa-mir-33a co-targeting FOSL1 and EN2 has a good predictive value for glioblastoma and skeletal muscle reduction.

The Fra-1: Novel role in regulating extensive immune cell states and affecting inflammatory diseases

Fra-1(Fos-related antigen1), a member of transcription factor activator protein (AP-1), plays an important role in cell proliferation, apoptosis, differentiation, inflammation, oncogenesis and tumor metastasis. Accumulating evidence suggest that the malignancy and invasive ability of tumors can be significantly changed by directly targeting Fra-1. Besides, the effects of Fra-1 are gradually revealed in immune and inflammatory settings, such as arthritis, pneumonia, psoriasis and cardiovascular disease. These regulatory mechanisms that orchestrate immune and non-immune cells underlie Fra-1 as a potential therapeutic target for a variety of human diseases. In this review, we focus on the current knowledge of Fra-1 in immune system, highlighting its unique importance in regulating tissue homeostasis. In addition, we also discuss the possible critical intervention strategy in diseases, which also outline future research and development avenues.

Therapeutic targeting of ATR in alveolar rhabdomyosarcoma

Despite advances in multi-modal treatment approaches, clinical outcomes of patients suffering from PAX3-FOXO1 fusion oncogene-expressing alveolar rhabdomyosarcoma (ARMS) remain dismal. Here we show that PAX3-FOXO1-expressing ARMS cells are sensitive to pharmacological ataxia telangiectasia and Rad3 related protein (ATR) inhibition. Expression of PAX3-FOXO1 in muscle progenitor cells is not only sufficient to increase sensitivity to ATR inhibition, but PAX3-FOXO1-expressing rhabdomyosarcoma cells also exhibit increased sensitivity to structurally diverse inhibitors of ATR. Mechanistically, ATR inhibition leads to replication stress exacerbation, decreased BRCA1 phosphorylation and reduced homologous recombination-mediated DNA repair pathway activity. Consequently, ATR inhibitor treatment increases sensitivity of ARMS cells to PARP1 inhibition in vitro, and combined treatment with ATR and PARP1 inhibitors induces complete regression of primary patient-derived ARMS xenografts in vivo. Lastly, a genome-wide CRISPR activation screen (CRISPRa) in combination with transcriptional analyses of ATR inhibitor resistant ARMS cells identifies the RAS-MAPK pathway and its targets, the FOS gene family, as inducers of resistance to ATR inhibition. Our findings provide a rationale for upcoming biomarker-driven clinical trials of ATR inhibitors in patients suffering from ARMS.

Transcriptome profile of spleen tissues from locally-adapted Kenyan pigs (Sus scrofa) experimentally infected with three varying doses of a highly virulent African swine fever virus genotype IX isolate: Ken12/busia.1 (ken-1033)

Background

African swine fever (ASF) is a lethal hemorrhagic disease affecting domestic pigs resulting in up to 100% mortality rates caused by the ASF virus (ASFV). The locally-adapted pigs in South-western Kenya have been reported to be resilient to disease and harsh climatic conditions and tolerate ASF; however, the mechanisms by which this tolerance is sustained remain largely unknown. We evaluated the gene expression patterns in spleen tissues of these locally-adapted pigs in response to varying infective doses of ASFV to elucidate the virus-host interaction dynamics.

Methods

Locally adapted pigs (n = 14) were experimentally infected with a high dose (1x10⁶HAD₅₀), medium dose (1x10⁴HAD₅₀), and low dose (1x10²HAD₅₀) of the highly virulent genotype IX ASFV Ken12/busia.1 (Ken-1033) isolate diluted in PBS and followed through the course of infection for 29 days. The in vivo pig host and ASFV pathogen gene expression in spleen tissues from 10 pigs (including three from each infective group and one uninfected control) were analyzed in a dual-RNASeq fashion. We compared gene expression between three varying doses in the host and pathogen by contrasting experiment groups against the naïve control.

Results

A total of 4954 differentially expressed genes (DEGs) were detected after ASFV Ken12/1 infection, including 3055, 1771, and 128 DEGs in the high, medium, and low doses, respectively. Gene ontology and KEGG pathway analysis showed that the DEGs were enriched for genes involved in the innate immune response, inflammatory response, autophagy, and apoptosis in lethal dose groups. The surviving low dose group suppressed genes in pathways of physiopathological importance. We found a strong association between severe ASF pathogenesis in the high and medium dose groups with upregulation of proinflammatory cytokines and immunomodulation of cytokine expression possibly induced by overproduction of prostaglandin E synthase (4-fold; p < 0.05) or through downregulation of expression of M1-activating receptors, signal transductors, and transcription factors. The host-pathogen interaction resulted in induction of expression of immune-suppressive cytokines (IL-27), inactivation of autophagy and apoptosis through up-regulation of NUPR1 [5.7-fold (high dose) and 5.1-fold (medium dose) [p < 0.05] and IL7R expression. We detected repression of genes involved in MHC class II antigen processing and presentation, such as cathepsins, SLA-DQB1, SLA-DOB, SLA-DMB, SLA-DRA, and SLA-DQA in the medium and high dose groups. Additionally, the host-pathogen interaction activated the CD8⁺ cytotoxicity and neutrophil machinery by increasing the expression of neutrophils/CD8⁺ T effector cell-recruiting chemokines (CCL2, CXCL2, CXCL10, CCL23, CCL4, CXCL8, and CXCL13) in the lethal high and medium dose groups. The recovered pigs infected with ASFV at a low dose significantly repressed the expression of CXCL10, averting induction of T lymphocyte apoptosis and FUNDC1 that suppressed neutrophilia.

Conclusions

We provide the first in vivo gene expression profile data from locally-adapted pigs from south-western Kenya following experimental infection with a highly virulent ASFV genotype IX isolate at varying doses that mimic acute and mild disease. Our study showed that the locally-adapted pigs induced the expression of genes associated with tolerance to infection and repression of genes involved in inflammation at varying levels depending upon the ASFV dose administered.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08754-8.

Deeper insights into long-term survival heterogeneity of pancreatic ductal adenocarcinoma (PDAC) patients using integrative individual- and group-level transcriptome network analyses

Pancreatic ductal adenocarcinoma (PDAC) is categorized as the leading cause of cancer mortality worldwide. However, its predictive markers for long-term survival are not well known. It is interesting to delineate individual-specific perturbed genes when comparing long-term (LT) and short-term (ST) PDAC survivors and integrate individual- and group-based transcriptome profiling. Using a discovery cohort of 19 PDAC patients from CHU-Liège (Belgium), we first performed differential gene expression analysis comparing LT to ST survivor. Second, we adopted systems biology approaches to obtain clinically relevant gene modules. Third, we created individual-specific perturbation profiles. Furthermore, we used Degree-Aware disease gene prioritizing (DADA) method to develop PDAC disease modules; Network-based Integration of Multi-omics Data (NetICS) to integrate group-based and individual-specific perturbed genes in relation to PDAC LT survival. We identified 173 differentially expressed genes (DEGs) in ST and LT survivors and five modules (including 38 DEGs) showing associations to clinical traits. Validation of DEGs in the molecular lab suggested a role of REG4 and TSPAN8 in PDAC survival. Via NetICS and DADA, we identified various known oncogenes such as CUL1 and TGFB1. Our proposed analytic workflow shows the advantages of combining clinical and omics data as well as individual- and group-level transcriptome profiling.

Combined Inhibition of FOSL-1 and YAP Using siRNA-Lipoplexes Reduces the Growth of Pancreatic Tumor

Pancreatic cancer evades most of the current therapies and there is an urgent need for new treatments that could efficiently eliminate this aggressive tumor, such as the blocking of routes driving cell proliferation. In this work, we propose the use of small interfering RNA (siRNA) to inhibit the combined expression of FOSL-1 and YAP, two signaling proteins related with tumor cell proliferation and survival. To improve the efficacy of cell transfection, DODAB:MO (1:2) liposomes were used as siRNA nanocarriers, forming a complex denominated siRNA-lipoplexes. Liposomes and lipoplexes (carrying two siRNA for each targeted protein, or the combination of four siRNAs) were physico-chemically and biologically characterized. They showed very good biocompatibility and stability. The efficient targeting of FOSL-1 and YAP expression at both mRNA and protein levels was first proved in vitro using mouse pancreatic tumoral cell lines (KRASG12V and p53 knockout), followed by in vivo studies using subcutaneous allografts on mice. The peri-tumoral injection of lipoplexes lead to a significant decrease in the tumor growth in both Athymic Nude-Foxn1nu and C57BL/6 mice, mainly in those receiving the combination of four siRNAs, targeting both YAP and FOSL-1. These results open a new perspective to overcome the fast tumor progression in pancreatic cancer.

PDGFA-associated protein 1 is a novel target of c-Myc and contributes to colorectal cancer initiation and progression

BACKGROUND

The mechanism underlying colorectal cancer (CRC) initiation and progression remains elusive, and overall survival is far from satisfactory. Previous studies have shown that PDGFA-associated protein 1 (PDAP1) is upregulated in several cancers including CRC. Here, we aimed to identify the cause and consequence of PDAP1 dysregulation in CRC and evaluate its role as a potential therapeutic target.

METHODS

Multi-omics data analysis was performed to identify potential key players in CRC initiation and progression. Immunohistochemistry (IHC) staining was applied to determine the expression pattern of PDAP1 in CRC tissues. Pdap1 conditional knockout mice were used to establish colitis and CRC mouse models. RNA sequencing, a phosphoprotein antibody array, western blotting, histological analysis, 5-bromo-2'-deoxyuridine (BrdU) incorporation assay, and interactome analysis were applied to identify the underlying mechanisms of PDAP1. A human patient-derived xenograft (PDX) model was used to assess the potential of PDAP1 as a therapeutic target.

RESULTS

PDAP1 was identified as a potential key player in CRC development using multi-omics data analysis. PDAP1 was overexpressed in CRC cells and correlated with reduced overall survival. Further investigation showed that PDAP1 was critical for the regulation of cell proliferation, migration, invasion, and metastasis. Significantly, depletion of Pdap1 in intestinal epithelial cells impaired mucosal restitution in dextran sulfate sodium salt-induced colitis and inhibited tumor initiation and growth in colitis-associated cancers. Mechanistic studies showed that c-Myc directly transactivated PDAP1, which contributed to the high PDAP1 expression in CRC cells. PDAP1 interacted with the juxtamembrane domain of epidermal growth factor receptor (EGFR) and facilitated EGFR-mitogen-activated protein kinase (MAPK) signaling activation, which resulted in FOS-related antigen 1 (FRA-1) expression, thereby facilitating CRC progression. Notably, silencing of PDAP1 could hinder the growth of patient-derived xenografts that sustain high PDAP1 levels.

CONCLUSIONS

PDAP1 facilitates mucosal restitution and carcinogenesis in colitis-associated cancer. c-Myc-driven upregulation of PDAP1 promotes proliferation, migration, invasion, and metastasis of CRC cells via the EGFR-MAPK-FRA-1 signaling axis. These findings indicated that PDAP1 inhibition is warranted for CRC patients with PDAP1 overexpression.

LncRNA HOTAIR sponges miR-301a-3p to promote glioblastoma proliferation and invasion through upregulating FOSL1

Glioblastoma, one of the most fatal brain tumors, is associated with a dismal prognosis and an extremely short overall survival. We previously reported that the overexpressed transient receptor potential channel TRPM7 is an essential glioblastoma regulator. Accumulating evidence suggests that long noncoding RNAs (lncRNAs) play an important role in glioma's initiation and progression. However, the function of lncRNA, HOX transcript antisense intergenic RNA (HOTAIR) mediated by TRPM7 in glioma remains unclear. In this study, HOTAIR expression was found to be positively regulated by TRPM7, significantly upregulated in glioma tissues, and is a poor prognosis factor for glioma patients. Moreover, reduced HOTAIR expression impeded the proliferation and invasion of glioma cells. Mechanistically, HOTAIR directly interacted with miR-301a-3p, and downregulation of miR-301a-3p efficiently reversed FOSL1 suppression induced by siRNA HOTAIR, which implied that HOTAIR positively regulated FOSL1 level through sponging miR-301a-3p and played an oncogenic role in glioma progression. In contrast to HOTAIR's role, miR-301a-3p alone served as a tumor suppressor to decrease glioma cell viability and migration/invasion. In agreement with HOTAIR's role, FOSL1 functioned as a tumorigenic gene in glioma pathogenesis, which was highly expressed in glioma tissues, and was shown to be an unfavorable prognostic factor for glioma patients. Mechanically, FOSL1 inhibition by siRNA FOSL1 efficiently rescued the oncogenic-like phenotypes caused by the miR-301a-3p inhibitor in glioma pathogenesis. SIGNIFICANCE: Our study elucidated the role of TRPM7-mediated HOTAIR as a miRNA sponge to target downstream FOSL1 oncogene and therefore consequently contribute to gliomagenesis, which shed new light on TRPM7/lncRNA-directed diagnostic and therapeutic approach in glioma.

The Fra-1/AP-1 Oncoprotein: From the "Undruggable" Transcription Factor to Therapeutic Targeting

The genetic and epigenetic changes affecting transcription factors, coactivators, and chromatin modifiers are key determinants of the hallmarks of cancer. The acquired dependence on oncogenic transcriptional regulators, representing a major determinant of cancer cell vulnerability, points to transcription factors as ideal therapeutic targets. However, given the unavailability of catalytic activities or binding pockets for small-molecule inhibitors, transcription factors are generally regarded as undruggable proteins. Among components of the AP-1 complex, the FOS-family transcription factor Fra-1, encoded by FOSL1, has emerged as a prominent therapeutic target. Fra-1 is overexpressed in most solid tumors, in response to the BRAF-MAPK, Wnt-beta-catenin, Hippo-YAP, IL-6-Stat3, and other major oncogenic pathways. In vitro functional analyses, validated in onco-mouse models and corroborated by prognostic correlations, show that Fra-1-containing dimers control tumor growth and disease progression. Fra-1 participates in key mechanisms of cancer cell invasion, Epithelial-to-Mesenchymal Transition, and metastatic spreading, by driving the expression of EMT-inducing transcription factors, cytokines, and microRNAs. Here we survey various strategies aimed at inhibiting tumor growth, metastatic dissemination, and drug resistance by interfering with Fra-1 expression, stability, and transcriptional activity. We summarize several tools aimed at the design and tumor-specific delivery of Fra-1/AP-1-specific drugs. Along with RNA-based therapeutics targeting the FOSL1 gene, its mRNA, or cognate regulatory circRNAs, we will examine the exploitation of blocking peptides, small molecule inhibitors, and innovative Fra-1 protein degraders. We also consider the possible caveats concerning Fra-1 inhibition in specific therapeutic contexts. Finally, we discuss a recent suicide gene therapy-based approach, aimed at selectively killing the Fra-1-overexpressing neoplastic cells.

Inhibiting mutant KRAS G12D gene expression using novel peptide nucleic acid‑based antisense: A potential new drug candidate for pancreatic cancer

KRAS mutations, which are the main cause of the pathogenesis of lethal pancreatic adenocarcinomas, impair the functioning of the GTPase subunit, thus rendering it constitutively active and signaling intracellular pathways that end with cell transformation. In the present study, the AsPC-1 cell line, which has a G12D-mutated KRAS gene sequence, was utilized as a cellular model to test peptide nucleic acid-based antisense technology. The use of peptide nucleic acids (PNAs) that are built to exhibit improved hybridization specificity and have an affinity for complementary RNA and DNA sequences, as well as a simple chemical structure and high biological stability that affords resistance to nucleases and proteases, enabled targeting of the KRAS-mutated gene to inhibit its expression at the translation level. Because PNA-based antisense molecules should be capable of binding to KRAS mRNA sequences, PNAs were utilized to target the mRNA of the mutated KRAS gene, a strategy that could lead to the development of a novel drug for pancreatic cancer. Moreover, it was demonstrated that introducing new PNA to cells inhibited the growth of cancer cells and induced apoptotic death and, notably, that it can inhibit G12D-mutated KRAS gene expression, as demonstrated by RT-PCR and western blotting. Altogether, these data strongly suggest that the use of PNA-based antisense agents is an attractive therapeutic approach to treating KRAS-driven cancers and may lead to the development of novel drugs that target the expression of other mutated genes.

ALKBH5 promotes the progression of infantile hemangioma through regulating the NEAT1/miR-378b/FOSL1 axis

Our work aims to investigate long non-coding RNA (lncRNA) N6-methyladenosine (m⁶A) modification and its role in infantile hemangioma (IH). The mRNA and protein expression levels were assessed using quantitative real-time polymerase chain reaction, western blot and immunohistochemistry. Me-RIP assay was performed to evaluate lncRNA NEAT1 m⁶A levels. Cell proliferation, migration and invasion were evaluated using cell counting kit-8 assay, transwell migration and invasion assay, respectively. Photo-activatable ribonucleoside-enhanced crosslinking and immunoprecipitation assay was conducted to verify the binding relationship between lncRNA nuclear paraspeckle assembly transcript 1 (NEAT1) and ALKBH5 (an RNA demethylase). The binding relationship between lncRNA NEAT1, microRNA (miR)-378b and FOS-like antigen 1 (FOSL1) was verified using dual-luciferase reporter gene assay and/or RNA immunoprecipitation assay. ALKBH5, lncRNA NEAT1 and FOLS1 expression was elevated in IH tissues, while miR-378b was downregulated. ALKBH5 knockdown suppressed cell proliferation, migration and invasion of IH cells, while promoting cell apoptosis. ALKBH5 promoted lncRNA NEAT1 expression by reducing the m⁶A modification of lncRNA NEAT1. In addition, miR-378b was the target of lncRNA NEAT1, and its overexpression reversed the promotion effect of lncRNA NEAT1 overexpression on IH cell tumor-like behaviors. Moreover, FOLS1 was the target of miR-378b, and its overexpression reversed the inhibitory effect of miR-378b overexpression on IH cell tumor-like behaviors in vitro. ALKBH5 might have great potential as therapeutic target for IH, since ALKBH5 silencing suppressed IH progression by regulation of the NEAT1/miR-378b/FOSL1 axis.

Integrated single-cell multiomics analysis reveals novel candidate markers for prognosis in human pancreatic ductal adenocarcinoma

The epigenomic abnormality of pancreatic ductal adenocarcinoma (PDAC) has rarely been investigated due to its strong heterogeneity. Here, we used single-cell multiomics sequencing to simultaneously analyze the DNA methylome, chromatin accessibility and transcriptome in individual tumor cells of PDAC patients. We identified normal epithelial cells in the tumor lesion, which have euploid genomes, normal patterns of DNA methylation, and chromatin accessibility. Using all these normal epithelial cells as controls, we determined that DNA demethylation in the cancer genome was strongly enriched in heterochromatin regions but depleted in euchromatin regions. There were stronger negative correlations between RNA expression and promoter DNA methylation in cancer cells compared to those in normal epithelial cells. Through in-depth integrated analyses, a set of novel candidate biomarkers were identified, including ZNF667 and ZNF667-AS1, whose expressions were linked to a better prognosis for PDAC patients by affecting the proliferation of cancer cells. Our work systematically revealed the critical epigenomic features of cancer cells in PDAC patients at the single-cell level.

Epigenetic basis of oncogenic-Kras-mediated epithelial-cellular proliferation and plasticity

Oncogenic Kras induces a hyper-proliferative state that permits cells to progress to neoplasms in diverse epithelial tissues. Depending on the cell of origin, this also involves lineage transformation. Although a multitude of downstream factors have been implicated in these processes, the precise chronology of molecular events controlling them remains elusive. Using mouse models, primary human tissues, and cell lines, we show that, in Kras-mutant alveolar type II cells (AEC2), FOSL1-based AP-1 factor guides the mSWI/SNF complex to increase chromatin accessibility at genomic loci controlling the expression of genes necessary for neoplastic transformation. We identified two orthogonal processes in Kras-mutant distal airway club cells. The first promoted their transdifferentiation into an AEC2-like state through NKX2.1, and the second controlled oncogenic transformation through the AP-1 complex. Our results suggest that neoplasms retain an epigenetic memory of their cell of origin through cell-type-specific transcription factors. Our analysis showed that a cross-tissue-conserved AP-1-dependent chromatin remodeling program regulates carcinogenesis.

ERK inhibitor ASN007 effectively overcomes acquired resistance to EGFR inhibitor in non-small cell lung cancer

The emergence of acquired resistance limits the long-term efficacy of EGFR tyrosine kinase inhibitors (EGFR TKIs). Thus, development of effective strategies to overcome resistance to EGFR TKI is urgently needed. Multiple mechanisms to reactivate ERK signaling have been successfully demonstrated in acquired resistance models. We found that in EGFR mutant non-small cell lung cancer (NSCLC) patients, acquired resistance to EGFR TKIs was accompanied by increased activation of ERK. Increased ERK activation was also found in in vitro models of acquired EGFR TKI resistance. ASN007 is a potent selective ERK1/2 inhibitor with promising antitumor activity in cancers with BRAF and RAS mutations. ASN007 treatment impeded tumor cell growth and the cell cycle in EGFR TKI-resistant cells. In addition, combination treatment with ASN007 and EGFR TKIs significantly decreased the survival of resistant cells, enhanced induction of apoptosis, and effectively inhibited the growth of erlotinib-resistant xenografts, providing the preclinical rationale for testing combinations of ASN007 and EGFR TKIs in EGFR-mutated NSCLC patients. This study emphasizes the importance of targeting ERK signaling in maintaining the long-term benefits of EGFR TKIs by overcoming acquired resistance.

Histone deacetylases: A novel class of therapeutic targets for pancreatic cancer

Pancreatic cancer is the seventh leading cause of cancer death worldwide, with a low 5-year survival rate. Novel agents are urgently necessary to treat the main pathological type, known as pancreatic ductal carcinoma (PDAC). The dysregulation of histone deacetylases (HDACs) has been identified in association with PDAC, which can be more easily targeted by small molecular inhibitors than gene mutations and may represent a therapeutic breakthrough for PDAC. However, the contributions of HDACs to PDAC remain controversial, and pharmacokinetic challenges have limited the application of HDAC inhibitors (HDACis) in PDAC. This review summarizes the mechanisms associated with success and failure of HDACis in PDAC and discusses the recent progress made in HDACi development and application, such as combination therapies designed to enhance efficacy. More precise strategies involving HDACis might eventually improve the outcomes of PDAC treatment.

Master lineage transcription factors anchor trans mega transcriptional complexes at highly accessible enhancer sites to promote long-range chromatin clustering and transcription of distal target genes

The term 'super enhancers' (SE) has been widely used to describe stretches of closely localized enhancers that are occupied collectively by large numbers of transcription factors (TFs) and co-factors, and control the transcription of highly-expressed genes. Through integrated analysis of >600 DNase-seq, ChIP-seq, GRO-seq, STARR-seq, RNA-seq, Hi-C and ChIA-PET data in five human cancer cell lines, we identified a new class of autonomous SEs (aSEs) that are excluded from classic SE calls by the widely used Rank Ordering of Super-Enhancers (ROSE) method. TF footprint analysis revealed that compared to classic SEs and regular enhancers, aSEs are tightly bound by a dense array of master lineage TFs, which serve as anchors to recruit additional TFs and co-factors in trans. In addition, aSEs are preferentially enriched for Cohesins, which likely involve in stabilizing long-distance interactions between aSEs and their distal target genes. Finally, we showed that aSEs can be reliably predicted using a single DNase-seq data or combined with Mediator and/or P300 ChIP-seq. Overall, our study demonstrates that aSEs represent a unique class of functionally important enhancer elements that distally regulate the transcription of highly expressed genes.

Proteome Landscape of Epithelial-to-Mesenchymal Transition (EMT) of Retinal Pigment Epithelium Shares Commonalities With Malignancy-Associated EMT

Stress and injury to the retinal pigment epithelium (RPE) often lead to dedifferentiation and epithelial-to-mesenchymal transition (EMT). These processes have been implicated in several retinal diseases, including proliferative vitreoretinopathy, diabetic retinopathy, and age-related macular degeneration. Despite the importance of RPE-EMT and the large body of data characterizing malignancy-related EMT, comprehensive proteomic studies to define the protein changes and pathways underlying RPE-EMT have not been reported. This study sought to investigate the temporal protein expression changes that occur in a human-induced pluripotent stem cell–based RPE-EMT model. We utilized multiplexed isobaric tandem mass tag labeling followed by high-resolution tandem MS for precise and in-depth quantification of the RPE-EMT proteome. We have identified and quantified 7937 protein groups in our tandem mass tag–based MS analysis. We observed a total of 532 proteins that are differentially regulated during RPE-EMT. Furthermore, we integrated our proteomic data with prior transcriptomic (RNA-Seq) data to provide additional insights into RPE-EMT mechanisms. To validate these results, we have performed a label-free single-shot data-independent acquisition MS study. Our integrated analysis indicates both the commonality and uniqueness of RPE-EMT compared with malignancy-associated EMT. Our comparative analysis also revealed that multiple age-related macular degeneration–associated risk factors are differentially regulated during RPE-EMT. Together, our integrated dataset provides a comprehensive RPE-EMT atlas and resource for understanding the molecular signaling events and associated biological pathways that underlie RPE-EMT onset. This resource has already facilitated the identification of chemical modulators that could inhibit RPE-EMT, and it will hopefully aid in ongoing efforts to develop EMT inhibition as an approach for the treatment of retinal disease.

•

Proteomics data were integrated with prior transcriptomic (RNA-Seq) data on RPE-EMT.

•

Dysregulated RPE-EMT proteome shares commonality with malignancy-associated EMT.

•

Altered RPE-EMT proteome signatures correlated with known AMD-associated risk factors.

•

Protein kinases and phosphatases crosstalk modulate RPE-EMT.

Immune-Omics Networks of CD27, PD1, and PDL1 in Non-Small Cell Lung Cancer

To date, there are no prognostic/predictive biomarkers to select chemotherapy, immunotherapy, and radiotherapy in individual non-small cell lung cancer (NSCLC) patients. Major immune-checkpoint inhibitors (ICIs) have more DNA copy number variations (CNV) than mutations in The Cancer Genome Atlas (TCGA) NSCLC tumors. Nevertheless, CNV-mediated dysregulated gene expression in NSCLC is not well understood. Integrated CNV and transcriptional profiles in NSCLC tumors (n = 371) were analyzed using Boolean implication networks for the identification of a multi-omics CD27, PD1, and PDL1 network, containing novel prognostic genes and proliferation genes. A 5-gene (EIF2AK3, F2RL3, FOSL1, SLC25A26, and SPP1) prognostic model was developed and validated for patient stratification (p < 0.02, Kaplan-Meier analyses) in NSCLC tumors (n = 1163). A total of 13 genes (COPA, CSE1L, EIF2B3, LSM3, MCM5, PMPCB, POLR1B, POLR2F, PSMC3, PSMD11, RPL32, RPS18, and SNRPE) had a significant impact on proliferation in 100% of the NSCLC cell lines in both CRISPR-Cas9 (n = 78) and RNA interference (RNAi) assays (n = 92). Multiple identified genes were associated with chemoresponse and radiotherapy response in NSCLC cell lines (n = 117) and patient tumors (n = 966). Repurposing drugs were discovered based on this immune-omics network to improve NSCLC treatment.

Signature identification of relapse-related overall survival of early lung adenocarcinoma after radical surgery

Background

The widespread use of low-dose chest CT screening has improved the detection of early lung adenocarcinoma. Radical surgery is the best treatment strategy for patients with early lung adenocarcinoma; however, some patients present with postoperative recurrence and poor prognosis. Through this study, we hope to establish a model that can identify patients that are prone to recurrence and have poor prognosis after surgery for early lung adenocarcinoma.

Materials and Methods

We screened prognostic and relapse-related genes using The Cancer Genome Atlas (TCGA) database and the GSE50081 dataset from the Gene Expression Omnibus (GEO) database. The GSE30219 dataset was used to further screen target genes and construct a risk prognosis signature. Time-dependent ROC analysis, calibration degree analysis, and DCA were used to evaluate the reliability of the model. We validated the TCGA dataset, GSE50081, and GSE30219 internally. External validation was conducted in the GSE31210 dataset.

Results

A novel four-gene signature (INPP5B, FOSL2, CDCA3, RASAL2) was established to predict relapse-related survival outcomes in patients with early lung adenocarcinoma after surgery. The discovery of these genes may reveal the molecular mechanism of recurrence and poor prognosis of early lung adenocarcinoma. In addition, ROC analysis, calibration analysis and DCA were used to verify the genetic signature internally and externally. Our results showed that our gene signature had a good predictive ability for recurrence and prognosis.

Conclusions

We established a four-gene signature and predictive model to predict the recurrence and corresponding survival rates in patients with early lung adenocarcinoma after surgery. These may be helpful for reforumulating post-operative consolidation treatment strategies.

NF1 regulates mesenchymal glioblastoma plasticity and aggressiveness through the AP-1 transcription factor FOSL1

The molecular basis underlying glioblastoma (GBM) heterogeneity and plasticity is not fully understood. Using transcriptomic data of human patient-derived brain tumor stem cell lines (BTSCs), classified based on GBM-intrinsic signatures, we identify the AP-1 transcription factor FOSL1 as a key regulator of the mesenchymal (MES) subtype. We provide a mechanistic basis to the role of the neurofibromatosis type 1 gene (NF1), a negative regulator of the RAS/MAPK pathway, in GBM mesenchymal transformation through the modulation of FOSL1 expression. Depletion of FOSL1 in NF1-mutant human BTSCs and Kras-mutant mouse neural stem cells results in loss of the mesenchymal gene signature and reduction in stem cell properties and in vivo tumorigenic potential. Our data demonstrate that FOSL1 controls GBM plasticity and aggressiveness in response to NF1 alterations.