RBM10 Deficiency Is Associated With Increased Immune Activity in Lung Adenocarcinoma

Harnessing DNA replication stress to target RBM10 deficiency in lung adenocarcinoma

The splicing factor RNA-binding motif protein 10 (RBM10) is frequently mutated in lung adenocarcinoma (LUAD) (9-25%). Most RBM10 cancer mutations are loss-of-function, correlating with increased tumorigenesis and limiting the efficacy of current LUAD targeted therapies. Remarkably, therapeutic strategies leveraging RBM10 deficiency remain unexplored. Here, we conduct a CRISPR-Cas9 synthetic lethality (SL) screen and identify ~60 RBM10 SL genes, including WEE1 kinase. WEE1 inhibition sensitizes RBM10-deficient LUAD cells in-vitro and in-vivo. Mechanistically, we identify a splicing-independent role of RBM10 in regulating DNA replication fork progression and replication stress response, which underpins RBM10-WEE1 SL. Additionally, RBM10 interacts with active DNA replication forks, relying on DNA Primase Subunit 1 (PRIM1) that synthesizes Okazaki RNA primers. Functionally, we demonstrate that RBM10 serves as an anchor for recruiting Histone Deacetylase 1 (HDAC1) to facilitate H4K16 deacetylation and R-loop homeostasis to maintain replication fork stability. Collectively, our data reveal a role of RBM10 in fine-tuning DNA replication and provide therapeutic arsenal for targeting RBM10-deficient tumors.

Investigation of RBM10 mutation and its associations with clinical and molecular characteristics in EGFR-mutant and EGFR-wildtype lung adenocarcinoma

Background

RBM10 is commonly mutated in lung adenocarcinoma (LUAD). However, its role in the pathogenesis of LUAD remains undefined. EGFR-mutant LUAD represents a distinct subset of non-small cell lung cancer (NSCLC). The function of RBM10 in tumor pathogenesis is supposed to differ between EGFR-mutant and EGFR-wt LUAD. This study aimed to interrogate the prevalence of RBM10 mutation in a large cohort of Chinese patients with LUAD and investigate the association of RBM10 mutation with clinical and molecular characteristics of EGFR-mutant and EGFR-wt LUAD.

Methods

Tumor sequencing data from 2848 Chinese patients with LUAD were retrospectively reviewed and analyzed. The prevalence of RBM10 was also compared with other three cohorts: OrigMed (n = 1222), MSKCC (n = 1267), and TCGA (n = 566). The associations of RBM10 mutation with clinical and molecular characteristics were assessed. An external cohort of 182 patients with LUAD who received PD-1 inhibitor were used to investigate the association of RBM10 mutation with clinical outcomes upon immunotherapy.

Results

Our cohort showed a higher prevalence of RBM10 in EGFR-mutant LUAD than in EGFR-wt LUAD (14.8 % vs. 6.5 %, p < 0.001). The enrichment of RBM10 mutations in EGFR-mutant LUAD was also seen in another Chinese cohort (OrigMed: 14.9 % vs. 7.8 %, p < 0.001), but not in the two western cohorts (MSKCC: 7.4 % vs. 9.5 %, p = 0.272; TCGA: 8.1 % vs. 6.7 %, p = 0.624). RBM10 mutations co-occurred more frequently with EGFR L858R mutations (23.7 %) than with other types of EGFR mutations (19 del: 7.7 %; other: 7.1 % in others, p < 0.001). In EGFR-mutant LUAD, RBM10 mutations were more commonly found in stage I (18.2 %) and II (21.8 %) vs. stage III (9.4 %) and IV (11.3 %) tumors (p < 0.001). The proportion of PD-L1 positive expression in EGFR-mutant LUAD with concomitant RBM10 mutation was not different from that those without RBM10 mutations (41.8 % vs. 47.9 %, p = 0.566). In contrast, RBM10 mutation occurred more frequently in EGFR-wt LUAD at stage II-IV (stage II: 12.0 %, stage III: 8.7 %, stage IV: 6.6 %) than at stage I (2.8 %). EGFR-wt LUAD with concomitant RBM10 mutations had higher proportions of PD-L1 expression positivity (78.9 % vs. 61.9 %, p = 0.014) and higher tumor mutational load (8.97 vs. 2.99 muts/Mb, p < 0.001) than those without. Patients with EGFR-wt LUAD who also harbored RBM10 loss of function (LOF) mutations had a longer median PFS upon immunotherapy than those with RBM10 non-LOF mutations (7.15 m vs. 2.60 m, HR = 4.83 [1.30-17.94], p = 0.010).

Conclusion

We comprehensively investigated RBM10 mutations in a Chinese cohort with LUAD. Compared to western cohorts, a significant enrichment of RBM10 mutations in EGFR-mutant LUAD compared to EGFR-wildtype LUAD in the Chinese population. RBM10 mutation shows different associations with clinical and molecular characteristics between EGFR-mutant and EGFR-wt LUAD, suggesting a divergent mechanism between these two subsets via which RBM10 deficiency contribute to tumor pathogenesis. The findings contribute to our understanding of the molecular landscape of LUAD and highlight the importance of considering population-specific factors in cancer genomics research.

To develop a prognostic model for neoadjuvant immunochemotherapy efficacy in esophageal squamous cell carcinoma by analyzing the immune microenvironment

Objective

The choice of neoadjuvant therapy for esophageal squamous cell carcinoma (ESCC) is controversial. This study aims to provide a basis for clinical treatment selection by establishing a predictive model for the efficacy of neoadjuvant immunochemotherapy (NICT).

Methods

A retrospective analysis of 30 patients was conducted, divided into Response and Non-response groups based on whether they achieved major pathological remission (MPR). Differences in genes and immune microenvironment between the two groups were analyzed through next-generation sequencing (NGS) and multiplex immunofluorescence (mIF). Variables most closely related to therapeutic efficacy were selected through LASSO regression and ROC curves to establish a predictive model. An additional 48 patients were prospectively collected as a validation set to verify the model's effectiveness.

Results

NGS suggested seven differential genes (ATM, ATR, BIVM-ERCC5, MAP3K1, PRG, RBM10, and TSHR) between the two groups (P < 0.05). mIF indicated significant differences in the quantity and location of CD3+, PD-L1+, CD3+PD-L1+, CD4+PD-1+, CD4+LAG-3+, CD8+LAG-3+, LAG-3+ between the two groups before treatment (P < 0.05). Dynamic mIF analysis also indicated that CD3+, CD8+, and CD20+ all increased after treatment in both groups, with a more significant increase in CD8+ and CD20+ in the Response group (P < 0.05), and a more significant decrease in PD-L1+ (P < 0.05). The three variables most closely related to therapeutic efficacy were selected through LASSO regression and ROC curves: Tumor area PD-L1+ (AUC= 0.881), CD3+PD-L1+ (AUC= 0.833), and CD3+ (AUC= 0.826), and a predictive model was established. The model showed high performance in both the training set (AUC= 0.938) and the validation set (AUC= 0.832). Compared to the traditional CPS scoring criteria, the model showed significant improvements in accuracy (83.3% vs 70.8%), sensitivity (0.625 vs 0.312), and specificity (0.937 vs 0.906).

Conclusion

NICT treatment may exert anti-tumor effects by enriching immune cells and activating exhausted T cells. Tumor area CD3+, PD-L1+, and CD3+PD-L1+ are closely related to therapeutic efficacy. The model containing these three variables can accurately predict treatment outcomes, providing a reliable basis for the selection of neoadjuvant treatment plans.

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

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

RNA-binding motif protein 10 inactivates c-Myc by partnering with ribosomal proteins uL18 and uL5

RNA-binding motif protein 10 (RBM10) is a frequently mutated tumor suppressor in lung adenocarcinoma (LUAD). Yet, it remains unknown whether cancer-derived mutant RBM10 compromises its tumor suppression function and, if so, the molecular insight of the underlying mechanisms. Here, we show that wild-type RBM10 suppresses lung cancer cell growth and proliferation by inactivating c-Myc that is essential for cancer cell survival. RBM10 directly binds to c-Myc and promotes c-Myc's ubiquitin-dependent degradation, while RBM10 knockdown leads to the induction of c-Myc level and activity. This negative action on c-Myc is further boosted by ribosomal proteins (RPs) uL18 (RPL5) and uL5 (RPL11) via their direct binding to RBM10. Cancer-derived mutant RBM10-I316F fails to bind to uL18 and uL5 and to inactivate c-Myc, thus incapable of suppressing tumorigenesis. Our findings uncover RBM10 as a pivotal c-Myc repressor by cooperating with uL18 and uL5 in lung cancer cells, as its failure to do so upon mutation favors tumorigenesis.

RNA splicing alterations in lung cancer pathogenesis and therapy

RNA splicing alterations are widespread and play critical roles in cancer pathogenesis and therapy. Lung cancer is highly heterogeneous and causes the most cancer-related deaths worldwide. Large-scale multi-omics studies have not only characterized the mutational landscapes but also discovered a plethora of transcriptional and post-transcriptional changes in lung cancer. Such resources have greatly facilitated the development of new diagnostic markers and therapeutic options over the past two decades. Intriguingly, altered RNA splicing has emerged as an important molecular feature and therapeutic target of lung cancer. In this review, we provide a brief overview of splicing dysregulation in lung cancer and summarize the recent progress on key splicing events and splicing factors that contribute to lung cancer pathogenesis. Moreover, we describe the general strategies targeting splicing alterations in lung cancer and highlight the potential of combining splicing modulation with currently approved therapies to combat this deadly disease. This review provides new mechanistic and therapeutic insights into splicing dysregulation in cancer.

Multiple Omics Analysis of the Role of RBM10 Gene Instability in Immune Regulation and Drug Sensitivity in Patients with Lung Adenocarcinoma (LUAD)

OBJECTIVE

The RNA-binding protein RBM10 can regulate apoptosis during the proliferation and migration of pancreatic cancer, endometrial cancer, and osteosarcoma cells; however, the molecular mechanism underlying lung adenocarcinoma is rarely reported. Recent studies have detected multiple truncated and missense mutations in RBM10 in lung adenocarcinoma, but the role of RBM10 in lung adenocarcinoma is unclear. This study mainly explored the immune regulation mechanism of RBM10 in the development of lung adenocarcinoma and its influence on sensitivity to targeted therapy drugs.

METHODS

The transcriptome data of CGAP were used to analyze the RNA-seq data of lung adenocarcinoma patients from different subgroups by using the CIBERSORT algorithm to infer the relative proportion of various immune infiltrating cells, and Spearman correlation analysis was performed to determine the gene expression and immune cell content. In addition, this study utilized drug trial data from the GDSC database. The IC50 estimates for each specific targeted therapy were obtained by using a regression method, and the regression and prediction accuracy were tested via ten cross-validations with the GDSC training set. An immunohistochemical test was performed on the samples of 20 patients with lung adenocarcinoma in the subcomponent analysis of immune cells, and the protein expression of RBM10 in lung adenocarcinoma tissues was verified by cellular immunofluorescence assays. Nucleic acids were extracted at low temperatures, and qRT-PCR was used to verify the expression levels of the mRNA of RBM10 in lung adenocarcinoma tissues and normal tissues (p < 0.05).

RESULTS

After screening and inclusion using a machine language, the results showed that RBM10 was significantly highly expressed in the lung adenocarcinoma tissues. The related signaling pathways were mainly concentrated in ncRNA processing, rRNA metabolic processes, ribosome biogenesis, and the regulation of translation. The qRT-PCR for 20 lung adenocarcinoma tissues showed that the expression of RBM10 in these tissues was significantly different from that in normal tissues (p = 0.0255). Immunohistochemistry analysis and cell immunofluorescence staining also confirmed that RBM10 was involved in the immune regulation of lung adenocarcinoma tissues, and the number of immune cell aggregations was significantly higher than that of the control group. RBM10 regulates B cell memory-CIBERSORT (p = 0.042) and B cell memory-CIBERSOTRT-abs (p = 0.027), cancer-associated fibroblast-EPIC (p = 0.001), cancer-associated fibroblast- MCPCounter (p = 0.0037), etc. The risk score was significantly associated with the sensitivity of patients to lapatinib (p = 0.049), nilotinib (p = 0.015), pazopanib (p = 0.001), and sorafenib (p = 0.048).

CONCLUSIONS

RBM10 can inhibit the proliferation and invasion of lung adenocarcinoma cells through negative regulation and promote the apoptosis of lung adenocarcinoma cells through immunomodulatory mechanisms. The expression level of RBM10 affects the efficacy of targeted drug therapy and the survival prognosis of lung adenocarcinoma patients, which has a certain guiding significance for the clinical treatment of these patients.

RBM10 Is a Biomarker Associated with Pan-Cancer Prognosis and Immune Infiltration: System Analysis Combined with In Vitro and Vivo Experiments

RNA binding motif protein 10 (RBM10) is a splicing factor that has been reported to be involved in the occurrence and progression of multiple malignancies. However, the RBM10 involvement in pan-cancer immunotherapy is not clear. In here, we provide the first comprehensive assessment of the prognostic value and immunological function of RBM10 in human pan-cancer utilizing multiple public databases. Data reveal the aberrant RBM10 expression in most tumors, and its expression is positively or negatively linked with the clinical prognosis of various cancers, depending on the different types and subtypes of cancers. In most tumors, RBM10 mutations are frequently occurred, which is closely related to tumor progression. Moreover, our results also show that RBM10 is considerably linked with most of the immune checkpoint genes, tumor immune cell infiltration, tumor mutation burden, and microsatellite instability. Additionally, RBM10 is significantly positively correlated with the sensitivity of trametinib, 17-AAG, PD-0325901, RDEA119, cetuximab, and afatinib, indicating potential antagonism between RBM10 inhibitors and these antitumor drugs, and more likely to develop drug resistance. We also verify that downregulation of RBM10 enhances the malignant phenotype of lung adenocarcinoma cells using in vitro cell experiments, and in vivo animal experiments show that the overexpression of RBM10 reduces the growth of tumors. Furthermore, upregulating RBM10 greatly reduces the PD-L1 protein levels, while silencing RBM10 considerably enhances PD-L1 protein levels. Moreover, the overexpression of RBM10 decreases the protein stability of PD-L1. To sum up, our pan-cancer analysis indicates that RBM10 is a promising biomarker for prognosis and immunotherapy, which provides a new insight for cancer immunotherapy.

Immune Subtypes in LUAD Identify Novel Tumor Microenvironment Profiles With Prognostic and Therapeutic Implications

The six transcriptomic immune subtypes (ISs) (C1 - C6) were reported to have complex and different interplay between TME and cancer cells in TCGA (The Cancer Genome Atlas) pan-cancer cohort. Our study specifically explored how the consequence of interplay determines the prognosis and the response to therapy in LUAD cohorts. Clinical and molecular information of LUAD patients were from TCGA and Gene Expression Omnibus (GEO). The immune cell populations and gene/pathway enrichment analysis were performed to explore the molecular differences among the C3 IS and other ISs in the LUAD population. The proportion of C3 inflammatory IS was identified as the most common IS in both TCGA (N = 457) and GEO (N = 901) cohorts. The C3 IS was also found to be the most accurate prognostic subtype, which was associated with significantly longer OS (p <0.001) and DFS (p <0.001). The C3 IS presented higher levels of CD8 T, M1 macrophage, and myeloid dendritic cells, while lower levels of M2 macrophages and cancer-associated fibroblast cells. Moreover, the C3 subtype was enriched in the antigen process and presenting, interferon-gamma response, T cell receptor signaling, and natural killer cell-mediated cytotoxicity pathways than C1/C2. In contrast, the C1/C2 presented greater activation of pathways related to the cell cycles, DNA repair, and p53 signaling pathways. The immune-related C3 IS had a great ability to stratify the prognosis of LUAD, providing clues for further pathogenic research. This classification might help direct precision medicine screenings of LUAD patients, thus possibly improving their prognoses.

The genomic landscape of metastatic clear cell renal cell carcinoma after systemic therapy

Primary clear cell renal cell carcinoma (ccRCC) has been previously characterized, but the genomic landscape of metastatic ccRCC is largely unexplored. Here, we performed whole exome sequencing (WES) in 68 samples from 44 patients with ccRCC, including 52 samples from a metastatic site. SETD2, PBRM1, APC and VHL were the most frequently mutated genes in the metastatic ccRCC cohort. RBM10 and FBXW7 were also among the 10 most frequently mutated genes in metastatic tissues. Recurrent somatic copy number variations (CNV) were observed at the previously identified regions 3p25, 9p21 and 14q25, but also at 6p21 (CDKN1A) and 13q14 (RB1). No statistically significant differences were found between samples from therapy‐naïve and pretreated patients. Clonal evolution analyses with multiple samples from 13 patients suggested that early appearance of CNVs at 3p25, 9p21 and 14q25 may be associated with rapid clinical progression. Overall, the genomic landscapes of primary and metastatic ccRCC seem to share frequent CNVs at 3p25, 9p21 and 14q25. Future work will clarify the implication of RBM10 and FBXW7 mutations and 6p21 and 13q14 CNVs in metastatic ccRCC.