- Primary research
- Open Access
Pan-cancer analysis identifies LMNB1 as a target to redress Th1/Th2 imbalance and enhance PARP inhibitor response in human cancers
Cancer Cell International volume 22, Article number: 101 (2022)
Emerging evidence suggests that LMNB1 is involved in the development of multiple cancer types. However, there is no study reporting the potential role of LMNB1 in a systematic pan-cancer manner.
The gene expression level and potential oncogenic roles of LMNB1 in The Cancer Genome Atlas (TCGA) database were analyzed with Tumor Immune Estimation Resource version 2 (TIMER2.0), Gene Expression Profiling Interactive Analysis version 2 (GEPIA2), UALCAN and Sangerbox tools. Pathway enrichment analysis was carried out to explore the possible mechanism of LMNB1 on tumorigenesis and tumor progression. The therapeutic effects of LMNB1 knockdown combined with PARP inhibition on human cancers were further investigated in vitro.
LMNB1 upregulation is generally observed in the tumor tissues of most TCGA cancer types, and is verified in kidney renal clear cell carcinoma using clinical specimens of our institute. High level of LMNB1 expression usually predicts poor overall survival and disease free survival for patients with tumors. Mechanically, LMNB1 level is positively correlated with CD4+ Th2 cell infiltration and DNA homologous recombination repair gene expression. In vitro experiments reveal that targeting LMNB1 has a synergistic effect on prostate cancer with PARP inhibitor treatment.
LMNB1 is a biomarker of CD4+ Th2 cell infiltration and DNA homologous recombination repair in human cancers. Blockage of LMNB1 combined with PARP inhibitor treatment could be a promising therapeutic strategy for patients with cancers.
Cancer is a complicated disease caused and influenced by genetic risks and environment factors. It is generally recognized that adverse genetic changes such as p53 gene mutations are significant risk factors associated with nearly all kinds of human tumors. With the bioinformatic advancements based on public cancer datasets and repositories, such as The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO) and Clinical Proteomic Tumor Analysis Consortium (CPTAC), pan-cancer analysis of interested oncogenes or tumor suppressive genes can be performed to explore the universal mechanisms of tumorigenesis and provide potential targets of general strategies for cancer therapy.
Lamin B1 is one of the two B-type lamin proteins [1, 2], which is encoded by LMNB1 gene (NCBI Entrez Gene: 4001). Lamin proteins are evolutionarily conserved among vertebrates , and located between the inner nuclear membrane and the peripheral heterochromatin, playing important roles in making up nuclear envelope and maintaining nuclear stability . Lamins are found implicated in the regulation of in a wide variety of nuclear processes, including DNA replication and repair, mitotic spindle assembly, chromosome distribution, oxidative stress response, and gene expression [5,6,7,8,9,10]. Recently, increased studies have revealed that LMNB1 is abnormally expressed in a lot of human cancer types, such as upregulated in lung adenocarcinoma , prostate cancer [12, 13], cervical cancer , pancreatic cancer , liver cancer , and downregulated in breast cancer , gastric cancer  and two subtypes of lung cancer . Though previous evidences showed that LMNB1 expression affected the clinical behavior and patient outcome of the above cancer types, the precise mechanism underlying carcinogenesis of LMNB1 has not been extensively explored.
Hence, this pan-cancer analysis was designed to investigate the precise role and possible mechanism of LMNB1 in all human cancer types. In brief, we comprehensively analyzed the mRNA and protein levels of LMNB1, and its association with prognostic landscape across all TCGA cancer types. Considering that malignant tumors often show a unique immunosuppressive microenvironment mediated by varieties of inflammation-related stromal cells and tumor-infiltrating immune cells, which results in poor response to immunotherapy and worse prognosis [20,21,22,23], and previous research indicated a significant correlation between LMNB1 expression and several kinds of tumor infiltrating lymphocytes in renal cell carcinoma , so we further investigated whether LMNB1 is involved in immune signaling pathways in a systematic pan-cancer manner, aiming to identify novel therapeutic strategy. Moreover, pathway enrichment analysis led us to focus attention on the relevance between LMNB1 expression and DNA homologous recombination repair (HRR), which plays a critical role in DNA repair and genome integrity. Loss of function alterations in HRR-associated genes such BRCA1 and BRCA2, resulted in better response to poly ADP-ribose polymerase (PARP) inhibitor (PARPi) treatment in cancer patients [25, 26]. We also investigated the effects of LMNB1 knockdown on PARPi treatment response in human prostate cancer.
To the best of our knowledge, this is the first pan-cancer analysis of LMNB1 gene, which showed that LMNB1 is not only a biomarker of CD4+ Th2 cell infiltration and DNA homologous recombination repair, but also a promising therapeutic target to redress Th1/Th2 imbalance and enhance PARPi response.
Materials and methods
Cell lines and cell culture
Human immortalized normal prostate epithelial cell line RWPE-1 and human prostate cancer cell line C4-2 were purchased from the American Type Culture Collect (Manassas, USA). Human benign prostate hyperplasia cell line BPH-1 were obtained from the Leibniz Institute DSMZ (Braunschweig, Germany). Other human prostate cancer cell lines namely LNCaP, VCaP, 22Rv1, DU145 and PC3 were purchased from the Cell Bank of Type Culture Collection, Chinese Academy of Science (Shanghai, China). RWPE-1, LNCaP, 22Rv1, DU145 and PC3 cell lines have been authenticated using STR profiling within the last 3 years. RWPE-1 cells were maintained in Keratinocyte Serum Free Medium (K-SFM, Gibco, 17005042, Grand Island, NY, USA) supplemented with 0.05 mg/ml bovine pituitary extract (Gibco, 17005042) and 5 ng/ml human recombinant EGF (Gibco, 17005042). BPH-1, LNCaP, C4-2 and 22Rv1 cells were maintained in RPMI-1640 (Gibco, 11875093) supplemented with 10% fetal bovine serum (FBS, Gibco, 10100147). VCaP and DU145 cells were maintained in DMEM (Gibco, 11995065) supplemented with 10% FBS. PC3 cells were maintained in F12K (Gibco, 21127022) supplemented with 10% FBS. All complete medium were supplemented with 100 units/ml penicillin and 100 μg/ml streptomycin. All cell experiments were performed with mycoplasma-free cells.
Overexpression/shRNA lentivirus of LMNB1 and transduction
Overexpression and shRNA lentivirus of human LMNB1 were designed and provided by GenePharma (Shanghai, China). PC3 cells were transduced with the overexpression lentivirus (multiplicity of infection (MOI): 50) to establish LMNB1 stably overexpressing cells, while 22Rv1 cells were transduced with the shRNA lentivirus (MOI: 10) to acquire LMNB1 stably knockdown cells. The transformed cells were selected by 1 μg/ml puromycin (Selleck, S7417, Shanghai, China) for 3 days. Overexpression or knockdown efficiency was determined by qRT-PCR.
Cell viability assay
6 × 103 cells suspended in 100 μl complete medium were seeded in 96-well plates per well. 12 h later, olaparib (Selleck, S1060) with a series of concentrations was added to the medium. After incubating at 37 ℃ for 72 h, the viability of prostate cancer cells was evaluated with Cell Counting Kit 8 (CCK-8) (Vazyme Biotech, A311, Nanjing, China) and detected at OD450 nm with a microplate reader (TECAN Group Ltd., infinite M200pro, Männedorf, CH/CHE). The half maximal inhibitory concentration (IC50) was acquired with the cell viability assay.
Patients and clinical specimens
This study was approved by the Ethics Committee of Drum Tower Hospital, Medical School of Nanjing University (Nanjing, China) and conducted in accordance with the Declaration of Helsinki principles. Clinical renal and prostate specimens of radical surgery along with corresponding clinicopathological characteristics were collected from June 2011 to June 2016. Neither chemotherapy nor radiotherapy was performed to the patients before surgery. In total, 25 pairs of frozen KIRC tissues and corresponding adjacent normal tissues were prepared for qRT-PCR analysis, and nine pairs for western blotting experiment. 130 formalin-fixated KIPC tissues of which 46 cases were accompanied by corresponding noncancerous kidney tissues, and 81 formalin-fixated PRAD tissues were collected for immunohistochemistry staining.
RNA isolation and quantitative real-time PCR (qRT-PCR)
Total RNA of cancer cells or surgical specimens was extracted using TRIzol® reagent (Invitrogen, 15596018, Carlsbad, CA, USA) and 1 μg total RNA was reversed to complementary DNA in 20 μl reaction system with PrimeScript™ RT Master Mix kit (TaKaRa Biotechnology (Dalian), RR036A, Dalian, China). qRT-PCR was performed by using ChamQ Universal SYBR qPCR Master Mix (Vazyme Biotech, Q711) with the QuantStudio™ 6 Flex Real-Time PCR System (Applied Biosystems, 4485692, Foster City, CA, USA). The primer concentration used for the qRT-PCR is 10 μM. Relative mRNA levels of target genes were normalized to ACTB by the 2−△△CT algorithm. Specially, when qRT-PCR was performed with clinical samples, both ACTB and GAPDH were set as endogenous control to avoid inaccurate interruption of results. Every qRT-PCR experiment was performed in triplicate and repeated three times independently. The primer sequences used in the study were listed in Additional file 2: Table S1.
Protein extraction and western blotting were performed as described previously . 20 μg total proteins in the lysates were loaded in and separated by SDS-PAGE. Primary antibodies for lamin B1 (1:5000, abcam, ab133741, Cambridge, MA, USA), cleaved PARP (Asp214) (1:1000, Cell Signaling Technology (CST), #5625, Danvers, MA, USA), cleaved caspase-3 (Asp175) (1:1000, CST, #9661), β-actin (1:1000, CST,#4970) and GAPDH (1:1000, CST, #5174) were used to bind the corresponding targets. After incubating with horseradish peroxidase-conjugated secondary antibodies, the protein bands were visualized with standard chemical luminescence methodology.
IHC staining was carried out on 3 µm thick paraffin-embedded sections with primary antibodies against lamin B1 (1:3000, abcam, ab133741) and BRCA1 (1:1000, abcam, ab16780). All the sections of IHC staining were scored by two well-trained pathologists independently. The staining intensity was scored as 0 (negative), 1 (weak), 2 (moderate), and 3 (strong). The staining range in the tumor area was scored as 0 (0%), 1 (1–25%), 2 (26–50%), 3 (51–75%) and 4 (75–100%). The final score was obtained by multiplying the intensity score and the staining range, which was defined as negative or low expression when ≤ 6 or as high expression when > 6. Five fields from each section were randomly chosen and evaluated.
Gene and protein expression analysis
To determine the expression pattern of LMNB1 in healthy human beings, we explored The Human Protein Atlas (HPA) (https://www.proteinatlas.org/). Consensus normalized expression (NX) was obtained by combining the data of several sources, including HPA, The Genotype-Tissue Expression (GTEx), The Functional Annotation of Mammalian Genomes 5 (FANTOM5), Single Cell Expression Atlas, the Human Cell Atlas, GEO, and the European Genome-phenome Archive. The protein concentrations in blood were determined by mass spectrometry-based plasma proteomics and estimated from spectral counts in the Peptide Atlas. “Enhanced” was defined as NX levels of a group (of 1–5 tissues or 1–10 cell types) at least 4 times the mean of other tissue or cell types. “Low specificity” was defined as NX levels not less than 1 in at least one tissue or cell type but not elevated in any tissue or cell type.
Tumor Immune Estimation Resource version 2 (TIMER2.0) (http://timer.cistrome.org/) [28,29,30] and Gene Expression Profiling Interactive Analysis version 2 (GEPIA2) (http://gepia2.cancer-pku.cn/#index)  provided powerful assistance in studying the differential expression of LMNB1 in TCGA database. We employed “Gene_DE” module of TIMER2.0 to obtain the box-whisker plots of LMNB1 expression between tumor and adjacent normal tissues across all TCGA tumors, and “Stage Plot” module of GEPIA2 to acquire the violin plots of LMNB1 expression grouped in different pathological stages of the TCGA tumor types whose data of pathological stage characteristic were available. The log2TPM transformed expression data were used for plotting.
Differential protein expression of LMNB1 was analyzed with UALCAN (http://ualcan.path.uab.edu/index.html) , which is a comprehensive webserver providing protein expression analysis based on the data from CPTAC .
Gene correlation analysis
We computed the correlation of LMNB1 and other genes or signatures of interest in the “Correlation Analysis” pane of GEPIA2. According to the user tutorial of GEPIA2, the non-log scale was used for calculation and the log-scale axis for visualization, and the mean value of the log2TPM of the genes included in a signature was used as the signature score . The heatmap of correlation between LMNB1 gene and target signature across all TCGA tumors was drawn by R language with “heatmap” package. Furthermore, we used “Gene_Corr” module of TIMER2.0 to get a detailed heatmap of the correlation between LMNB1 and every single gene in the specific signature in various TCGA cancer types.
We performed survival analysis with GEPIA2 tools. Heatmap of overall survival (OS) and disease free survival (DFS) based on LMNB1 expression across all TCGA tumor types were obtained in the “Survival Map” module, then Kaplan–Meier curves of the tumor types with logrank p value < 0.05 were drawn in the “Survival analysis” module. The group cutoff was set as 50% of LMNB1 expression.
Immune infiltration analysis
Immune pathways analysis was performed using the Sangerbox tools, a free online platform for TCGA data analysis (http://www.sangerbox.com/tool). Then we employed the “Immune-Gene” module of TIMER2.0 to investigate the details of immune cell infiltration levels of the selected immune pathways (e.g. CD4+ T cells) in a variety of TCGA tumors. The results were displayed in a heatmap with colors indicating the purity-adjusted Spearman’s rho coefficient across all TCGA cancers, and in scatter plots showing the relationship between LMNB1 expression and immune infiltration estimation with detailed Spearman’s rho coefficient and p value of interested cancer types.
Function enrichment analysis
In the “Similar Genes Detection” pane of GEPIA2, we obtained the top 100 genes that presented a similar expression pattern with LMNB1 in the combined data of all the TCGA tumor tissues. On the other hand, we acquired another top 50 genes whose corresponding proteins were capable of binding to lamin B1 in STRING website (https://string-db.org/). The parameter of active interaction sources was set as “experiments”, minimum required interaction score as “low confidence (0.150)”, max number of interactors to show as “no more than 50 interactors (in 1st shell)”, and the other checkboxes were kept default. We combined the two screened gene lists to performed Kyoto encyclopedia of genes and genomes (KEGG) and Gene ontology (GO) enrichment analysis. Briefly, the above gene aggregate was uploaded to Database for Annotation, Visualization, and Integrated Discovery (DAVID) website (https://david.ncifcrf.gov/) to acquire the enriched pathways associated with LMNB1 gene in homo sapiens, which were then visualized by R language with “ggplot2” package. GO term enrichment analysis was also conducted in R language with “clusterProfiler” package.
Gene Set Enrichment Analysis (GSEA) of TCGA datasets was performed with the “GSEA4.0” software (http://www.gsea-msigdb.org/gsea/index.jsp) according to the manufacture’s instruction [34, 35]. The mRNA expression data of TCGA cohort which was downloaded using Sangerbox tools was uploaded to the software, then analyzed with the “h.all.v7.4.symbols.gmt [Hallmarks]”, “c2.cp.kegg.v7.4.symbols.gmt [Curated]” and “c5.go.bp.v7.4.symbols.gmt [Gene ontology]” gene sets of the Molecular Signatures Database (MSigDB), by using LMNB1 gene expression as the phenotype.
R programming language (version 4.0.4), IBM SPSS Statistics (Version 22.0) and GraphPad Prism (version 8.0.2) were used for statistical analysis. Continuous normally distributed variables were presented as mean ± standard deviation, and compared using unpaired Student’s t test. Categorical variables were analyzed by Chi-square test or Fisher exact test. Pearson and Spearman correlation analyze was employed to evaluate the correlation between two continuous normally distributed factors or categorical variables, respectively. p < 0.05 was considered statistically significant.
LMNB1 upregulation was observed in various cancer types and validated in kidney renal clear cell carcinoma
The expression and distribution of LMNB1 in different human tissues and cell types under physiological conditions was explored in the HPA website. LMNB1 was widely expressed in almost all the tissue types while especially enhanced in the lymphoid tissue, including thymus, appendix, lymph node, tonsil and bone marrow (Additional file 1: Fig. S1A). Based on single cell RNAseq of human issues and blood cells, LMNB1 showed highest expression in granulocytes, monocytes, spermatids and extravillous trophoblasts (Additional file 1: Fig. S1B), whereas low expression specificity in blood cell types (Additional file 1: Fig. S1C). Based on the mass spectrometry of human plasma in the publicly available Peptide Atlas, the protein concentration of lamin B1 was estimated as 3.3 μg/l (Additional file 1: Fig. S1D).
The TIMER2.0 webserver was applied to study the differential expression of LMNB1 between tumor and adjacent normal tissues across 33 human cancer types in TCGA database. As shown in Fig. 1A, in the 21 distinct tumor types of which normal tissue data are available, LMNB1 expression in the tumor tissues is upregulated compared to the corresponding adjacent tissues in 19 kinds of tumors, including BLCA (Bladder Urothelial Carcinoma), BRCA (Breast invasive carcinoma), CESC (Cervical squamous cell carcinoma and endocervical adenocarcinoma), CHOL (Cholangio carcinoma), COAD (Colon adenocarcinoma), ESCA (Esophageal carcinoma), GBM (Glioblastoma multiforme), HNSC (Head and Neck squamous cell carcinoma), KIRC (Kidney renal clear cell carcinoma), KIRP (Kidney renal papillary cell carcinoma), LIHC (Liver hepatocellular carcinoma), LUAD (Lung adenocarcinoma), LUSC (Lung squamous cell carcinoma), PCPG (Pheochromocytoma and Paraganglioma), PRAD (Prostate adenocarcinoma), READ (Rectum adenocarcinoma), STAD (Stomach adenocarcinoma), THCA (Thyroid carcinoma), UCEC (Uterine Corpus Endometrial Carcinoma) (Fig. 1A). In contrast, only in KICH (Kidney Chromophobe), LMNB expression in the tumor tissues is downregulated compared to the corresponding control tissues. Besides, LMNB1 expression in the metastasis tissues of SKCM (Skin Cutaneous Melanoma) is much higher than the primary tumor tissues. In addition to transcriptome analysis, we further compared the protein level of lamin B1 in different tumor kinds through UALCAN resource which provides protein expression analysis using data from CPTAC. The total protein expression of lamin B1 was found significantly higher in the tumor tissues of all the six tumor kinds (BRCA, COAD, KIRC, LUAD, OV (Ovarian serous cystadenocarcinoma) and UCEC) whose normal tissue data were available in CPTAC, than the corresponding normal tissues (Fig. 1B). Furthermore, we applied GEPIA2 approach to explore the relationship between LMNB1 mRNA expression and the pathological stages across all TCGA tumors. The results indicated that when the pathological stage increased, LMNB1 expression showed a trend of gradual increase in ACC (Adrenocortical carcinoma), KIRC, LUAD, TGCT (Testicular Germ Cell Tumors), and a trend of gradual decrease in OV, while no obvious gradual trend but with stage-specific expression difference in KICH, KIRP, LIHC and SKCM (Fig. 1C).
Since analysis of TCGA and CPTAC datasets showed that in KIRC and LUAD tumor tissues, LMNB1 expression was upregulated both in mRNA and protein levels, as well as increased gradually with elevated pathological stages, we decided to validate the LMNB1 expression status in specimens of radical nephrectomy in our center. Firstly, we detected LMNB1 mRNA expression in 25 pairs of KIRC tumor tissues and adjacent normal tissues by qRT-PCR. Compared to the corresponding non-cancerous tissues, 21 of 25 (84.0%) cancerous tissues showed higher LMNB1 mRNA expression (Fig. 2A, B). Furthermore, western blotting and IHC staining confirmed that lamin B1 protein level in tumor specimens was much higher than the corresponding normal tissues (Fig. 2C–F). Consistent with TCGA analysis, IHC staining showed that KIRC with higher T stage had higher expression of lamin B1 (Fig. 2F, G), and KIRC with advanced Fuhrman nuclear grade (G3 and G4) also showed higher lamin B1 expression than low grade (G1 and G2) (Fig. 2H). Moreover, we explored the correlation between the clinicopathological features and lamin B1 expression in surgical resected specimens from a cohort of 130 consecutive KIRC patients by IHC staining. The results showed that high lamin B1 IHC staining in KIRC patients was positively associated with male gender (p = 0.001), pathological T stage (p < 0.001), distant metastasis (p = 0.008), Fuhrman grade (p = 0.002) and microvascular invasion (p = 0.001) (Table 1). These results indicated that LMNB1 upregulation was generally correlated with tumor initiation and progression.
LMNB1 upregulation was universally correlated with poor survival in human cancers
Since high expression of LMNB1 was found in the overwhelming majority of cancer types and correlated with cancer aggressiveness, we wonder whether it affects the prognosis of human cancers. The “Survival analysis” module of GEPIA2 tool was applied to compare the OS and DFS based on the expression status of LMNB1 in TCGA datasets. As shown in Fig. 3A, LMNB1 expression higher than median was significantly associated with poor prognosis of OS for cancers of ACC (p < 0.0001), LGG (Brain Lower Grade Glioma) (p = 0.0005), PAAD (p = 0.0076), KIRP (p = 0.0021), LIHC (p = 0.0034), MESO (Mesothelioma) (p = 0.047), SARC (Sarcoma) (p = 0.01), while with favorable prognosis of OS for cancers of LUSC (p = 0.038), THYM (Thymoma) (p = 0.009) (Fig. 3A). On the other hand, DFS analysis revealed that highly expressed LMNB1 was notably correlated with poor prognosis for ACC (p < 0.0001), KIRP (p = 0.0001), LIHC (p < 0.0001), PRAD (p = 0.013), UVM (Uveal Melanoma) (p = 0.012), ESCA (p = 0.04), LGG (p = 0.025), PAAD (p = 0.035), SARC (p = 0.0004), while not correlated with favorable prognosis for any human cancer type (Fig. 3B).
LMNB1 was identified as a biomarker of CD4+ Type 2 T helper cell infiltration in human cancers
For the past few years, immune cells infiltrated in the tumor microenvironment (TME) has been found playing key roles in the tumorigenesis and tumor progression [36, 37]. Previous study has demonstrated that a gene hub including LMNB1 in KIRC was positively correlated with multiple kinds of tumor infiltrating lymphocytes, such as activated CD4+ T cells, CD8+ T cells, regulatory T cells and follicular helper T cells, but negatively correlated with resting mast cells, resting NK cells and activated NK cells . Herein, we aimed to provide a comprehensive research on the potential relationship between LMNB1 expression and immune infiltration levels across all TCGA cancers. Firstly, Sangerbox tools helped us identified that activated CD4+ T cells and type 2 T helper cells were the main immune pathways which were significantly correlated with high expression of LMNB1 in most of human cancer types, especially ACC, KIRP, KIRC, LIHC, LGG, PRAD, THCA and UVM (Fig. 4A). Then we focus our attention on the details of CD4+ T cell infiltration levels. With the aid of TIMER2.0 resource, we employed the EPIC, TIMER, QUANTISEQ, XCELL, CIBERSORT, CIBERSORT-ABS methods to analyze the correlation between the infiltration levels of CD4+ T cell subtypes and the mRNA expression of LMNB1 across different TCGA tumor kinds. We found that all the TCGA cancers excluding UCS (Uterine Carcinosarcoma) showed statistical positive correlation between high LMNB1 expression and the immune infiltration of CD4+ Th2 cells based on XCELL algorithm (Fig. 4B). In contrast, the infiltration level of CD4+ central memory T cells and CD4+ effector memory T cells was found negatively correlated high LMNB1 expression in more than half of the TCGA cancer kinds (Fig. 4B). The scatter plots of the cancer types which occupied the top 10 purity-adjusted Spearman’s rho coefficient were showed (Fig. 4C, Additional file 1: Fig. S2). Since CD4+ Th2 cells and Th2 cytokines such as interleukin-4 (IL-4), IL-5, IL-6, IL10 and IL13 were thought to be associated with immunosuppressive contexture and tumor-promoting effects [38,39,40,41,42,43], we believed that LMNB1 could be considered as a biomarker of immunosuppressive microenvironment.
Pathways of cell cycle and nuclear division were involved in the effects of LMNB1 on tumor pathogenesis
In order to further illustrate the molecular mechanism of LMNB1 in tumor initiation and progression, we screened out LMNB1-related genes for subsequent pathway enrichment analysis. By using GEPIA2 tool, we acquired the top 100 genes which had similar expression patterns with LMNB1 in the combined data of all TCGA cancer types, and the most six correlated genes are TMPO (thymopoietin) (R = 0.81), KIF11 (kinesin family member 11) (R = 0.79), NUSAP1 (nucleolar and spindle associated protein 1) (R = 0.79), KIF15 (kinesin family member 11) (R = 0.78), MCM6 (minichromosome maintenance complex component 6) (R = 0.77), PLK4 (R = 0.77) (Fig. 5A). The corresponding heatmap also verified that the above top six related genes showed statistical positive correlation with LMNB1 in almost every cancer type (Additional file 1: Fig. S3). On the other hand, with STRING website we obtained another top 50 genes which encode the proteins that had been experimentally proved to physically bind to lamin B1 (Fig. 5B), such as RPA1 (replication protein A1) (interaction score = 0.873) and RPA3 (replication protein A3) (interaction score = 0.835). The venn diagram of the above two gene sets showed an intersection consisting of two genes, namely, TMPO and ZWINT (ZW10 interacting kinetochore protein) (Fig. 5C). Then, we carried out KEGG and GO enrichment analysis with the union of the above two gene sets. As shown in Fig. 5D, KEGG analysis revealed that “cell cycle” and “DNA replication” pathways appeared to play important roles in the influence of LMNB1 on the tumorigenesis and development (Fig. 5D). GO term analysis showed that “nuclear division”, “organelle fission”, “mitotic nuclear division” and “chromosome segregation” in biological process (Fig. 5E, F), “chromosomal region” and “spindle” in cellular component (Additional file 1: Fig. S4A, C), “tubulin binding” and “microtubule binding” in molecular function (Additional file 1: Fig. S4B, D) were most involved in the effects of LMNB1 on human tumors. Overall, we hold the opinion that LMNB1 played an important role in tumor growth and cell mitosis.
LMNB1 was required for DNA homologous recombination repair and contributed to PARPi resistance
Interestingly, we found that several genes in the gene aggregate which encodes lamin B1-binding proteins or shows similar expression pattern with LMNB1 were also associated with DNA HRR, such as BRCA1, CHEK1, RAD51, RAD54L and BRIP1 (Fig. 5F, Additional file 1: Fig. S4C, D). To further explore the detailed role of LMNB1 on the DNA HRR, we used GEPIA2 to analyze the correlation between the expression of LMNB1 and a self-defined gene set named as PROfound signature, which consists of 15 HRR-associated genes prespecified by Johann de Bono, et al. in their phase 3 PROfound clinical trials . The heatmap showed that LMNB1 expression had a strong positive correlation with the PROfound signature in the vast majority of TCGA cancer types, and the corresponding scatter plots of the cancer types whose Pearson correlation coefficient was not less than 0.80 were displayed (Fig. 6A). In addition, TIMER2.0 assisted us to show the detail of the correlation between LMNB1 and every single gene in the PROfound signature across all TCGA tumor kinds (Fig. 6B).
PROfound trials revealed that metastatic castration-resistant prostate cancer (CRPC) patients with at least one alteration in the PROfound signature genes would benefit for progression-free survival from olaparib (a selective PARP inhibitor) treatment . Since LMNB1 was revealed to be closely related to the PROfound signature, we speculated that LMNB1 might influence the treatment effects of olaparib. Firstly, in the PRAD cohort of TCGA database, LMNB1 appeared to be strongly correlated with the PROfound signature (R = 0.75, p = 5.6E−90) (Fig. 6C). Then we stably overexpressed LMNB1 gene in a CRPC cell line PC3, and qRT-PCR analysis showed that forced expression of LMNB1 increased the mRNA level of BRCA1, BRCA2, CHEK1, CHEK2 and ATM (Fig. 6D). To further validate the relationship between the expression of LMNB1 and HRR genes in PRAD, lamin B1 and BRCA1, which was selected as a representative of HRR proteins were examined in 81 paraffin-embedded prostate cancer tissues by IHC staining (Fig. 6E). As expected, we found that the tumor epithelial areas in the PRAD specimens displaying strong staining of lamin B1 also had heavy signals of BRCA1 (Fig. 6E). Statistically, a positive relationship between lamin B1 and BRCA1 was observed (Spearman r = 0.5407, p < 0.0001) (Fig. 6F). In addition, Pearson correlation analysis also indicated that LMNB1 mRNA level was significantly associated with BRCA1 mRNA expression in a panel of prostate cell lines (Pearson r = 0.8977, p = 0.0025) (Fig. 6G).
In order to investigate the effects of LMNB1 on the PARPi therapy, we stably knocked know LMNB1 in another PRAD cell line 22Rv1, by infecting with two lentiviruses containing specific short hairpin RNAs (shRNAs) targeting different regions of LMNB1 gene (Fig. 7A, Additional file 2: Table S2). The cell viability assay showed that silencing of LMNB1 in 22Rv1 cells resulted in a prominent decreased IC50 of olaparib (97.1 μM for shLMNB1#1 and 108.8 μM for shLMNB1#2 versus 177.6 μM for shNC) (Fig. 7B). Western blotting also demonstrated that LMNB1 knockdown caused more serious cell apoptosis characterized by cleaved PARP and cleaved caspase3 in the condition of olaparib incubation (Fig. 7C). Finally, GSEA of PRAD cohort in TCGA database showed that the pathways of “DNA repair” in hallmark gene set (Fig. 7D, E, Additional file 3: Table S3), “homologous recombination”, “mismatch repair”, “nucleotide excision repair”, “base excision repair” in KEGG gene set (Additional file 1: Fig. S5A, B, Additional file 4: Table S4), and “double strand break repair” in GO biological process gene set (Additional file 1: Fig. S5C, D, Additional file 5: Table S5) were positively associated with high LMNB1 expression, indicating that LMNB1 did play a crucial role on DNA repair and PARPi therapy.
LMNB1 is a major constituent of the nuclear lamina and plays important roles in diverse chromatin-associated processes. It has been recognized that duplication of LMNB1 gene or mutation in its promoter leading to increased expression of the encoded protein is a genetic cause of adult-onset autosomal dominant leukodystrophy (ADLD), which is a fatal progressive neurological disease characterized by widespread central nervous system demyelination [2, 44, 45]. The role of LMNB1 in tumor initiation and progression still remain elusive. Based on the data of TCGA, we found that LMNB1 was highly expressed in not only lung adenocarcinoma , prostate cancer [12, 13], cervical cancer  and liver cancer , but also most other human cancer types. Though previous reports indicated that LMNB1 was downregulated in breast cancer by qRT-PCR , gastric cancer , lung adenocarcinoma and squamous cell carcinoma by IHC staining , TCGA analysis indicated LMNB1 expression in the above 4 tumor kinds was upregulated at least in the mRNA level. Besides, proteomic data from CPTAC also demonstrated that lamin B1 expression in the protein level was upregulated in breast cancer and lung adenocarcinoma. Consistent with bioinformatics analysis, a moderate cohort of consecutive KIRC patients in our center also confirmed that LMNB1 was highly expressed in cancerous tissues and positively correlated with cancer aggressiveness. Further research based on multicenter design and relatively large sample size will help to draw a solid conclusion of LMNB1 expression in tumor tissues.
Increasing evidence indicates that LMNB1 is associated with tumor proliferation and metastasis in several types of human cancer, but as far as we know, there is little investigation focusing on the relationship between LMNB1 expression and immune cell infiltration. In this study, we first demonstrated a statistical positive correlation between the mRNA expression of LMNB1 and the infiltration level of CD4+ Th2 cells in all the TCGA tumor types excluding UCS. Th1 and Th2 cells are two distinct subsets of CD4+ T lymphocytes, and Th1/Th2 ratio is in balance under normal conditions . During extracellular pathogen infection, naïve Th cells proliferate and differentiate toward the Th2 subtype, which functions through secreting various effector cytokines, including IL-4, IL-5, IL-6, IL-10, IL-13. Previous studies indicate that Th1/Th2 imbalance has been observed with an elevation of Th2 cells and Th2-released cytokines in several kinds of human cancers [38,39,40,41,42,43]. Furthermore, patients with Th2 dominant response in TME have poorer prognosis than those with a balance Th1/Th2 infiltration [43, 47]. Therefore, development of new strategies targeting LMNB1 or Th2 cells may help to redress Th1/Th2 imbalance and improve the prognosis of tumor patients with high LMNB1 expression.
In our study, KEGG enrichment analysis using LMNB1-related partners based on TCGA datasets indicated that “cell cycle” and “DNA replication” are concerned in the effects of LMNB1 on tumor pathogenesis. These findings have been evidenced in previous studies that lamin B1 is a key component of the nuclear lamina providing a framework for the nuclear envelope and is required for the initiation phase or the elongation phase of DNA replication . Butin-Israeli et al. found that LMNB1 knockdown by shRNA in cancer cells slowed the cell cycle with the S phase delay, which was accompanied by the stalling and collapse of replication forks . Furthermore, the replication fork collapse resulted in the increase of double-strand DNA breaks (DSB), which were inefficiently repaired in LMNB1 deficient cells . Their results suggested that the maintain of LMNB1 expression was also required for chromatin stability and DNA repair. Consistently, in the LMNB1-related cellular process in cancer cells including “nuclear division”, “organelle fission”, “mitotic nuclear division” and “chromosome segregation”, which were enriched by GO pathway analysis using TCGA data, we found several HRR genes such as BRCA1, CHEK1, RAD51, RAD54L and BRIP1 were implicated. Further investigation confirmed a strong positive relevance between LMNB1 expression and the 15 HRR genes in the PROfound trial in almost all TCGA cancer types. Practically, IHC staining of prostate cancer tissues in our center and qRT-PCR analysis of prostate cell lines also showed a close correlation between LMNB1 and BRCA1 expression. Forced expression of LMNB1 in prostate cancer cells increased the mRNA level of several HRR genes, especially BRCA1. These findings were inconsistent with Butin-Israeli V’s research that BRCA1 expression was increased in LMNB1 knockdown cells by shRNA in osteosarcoma cell line U-2-OS and colorectal carcinoma cell line HCT116 . One possible explanation for the inconsistent findings may be that the mechanism of LMNB1 on regulating BRCA1 is distinct in different tumor types. On the other hand, the off-target effects in previous study could not be completely ruled out since only one shRNA was used to silence LMNB1 expression .
Tumor cells with deficiency of HRR genes, especially BRCA1 and BRCA2, are sensitive to PARPi through the mechanism of synthetic lethality [49, 50]. Since LMNB1 is tightly correlated with HRR genes and plays an important role in DNA repair, we speculate that targeting LMNB1 could synergistically promote the effects of PARPi on cancer cells. Cell viability assay and western blotting experiment validate our thoughts, and these findings suggest that combination therapy of LMNB1 knockdown and PARPi could improve the prognosis of cancer patients without loss-of-function alteration of HRR genes.
In conclusion, our pan-cancer analysis provides a comprehensive overview of the oncogenic roles of LMNB1 in human cancers. Overexpression of LMNB1 generally predicted poor prognosis for cancer patients and suggested high level of CD4+ Th2 cell infiltration. Moreover, LMNB1 was universally correlated with HRR genes, which may serve as a potential therapeutic target to enhance the antitumor effect effects of PARPi in human cancers.
Availability of data and materials
The public datasets involved in this study are available in online repositories. Original contributions presented in this study are included in the manuscript or additional files. Further inquiries can be directed to the corresponding author.
The Cancer Genome Atlas
Tumor Immune Estimation Resource version 2
Gene Expression Profiling Interactive Analysis version 2
Gene Expression Omnibus
Clinical Proteomic Tumor Analysis Consortium
Poly ADP-ribose polymerase
Cell Counting Kit 8
Quantitative real-time PCR
The Human Protein Atlas
The Genotype-Tissue Expression
The Functional Annotation of Mammalian Genomes 5
Disease free survival
Kyoto encyclopedia of genes and genomes
Database for Annotation, Visualization, and Integrated Discovery
Gene Set Enrichment Analysis
The Molecular Signatures Database
Bladder Urothelial Carcinoma
Breast invasive carcinoma
Cervical squamous cell carcinoma and endocervical adenocarcinoma
Head and neck squamous cell carcinoma
Kidney renal clear cell carcinoma
Kidney renal papillary cell carcinoma
Liver hepatocellular carcinoma
Lung squamous cell carcinoma
Pheochromocytoma and Paraganglioma
Uterine Corpus Endometrial Carcinoma
Skin Cutaneous Melanoma
Ovarian serous cystadenocarcinoma
Testicular Germ Cell Tumors
Brain Lower Grade Glioma
The tumor microenvironment
Homologous recombination repair
Castration-resistant prostate cancer
Short hairpin RNAs
Adult-onset autosomal dominant leukodystrophy
Double-strand DNA breaks
Dittmer TA, Misteli T. The lamin protein family. Genome Biol. 2011;12(5):222.
Padiath QS. Autosomal dominant leukodystrophy: a disease of the nuclear lamina. Front Cell Dev Biol. 2019;7:41.
Worman HJ. Nuclear lamins and laminopathies. J Pathol. 2012;226(2):316–25.
Bukata L, Parker SL, D’Angelo MA. Nuclear pore complexes in the maintenance of genome integrity. Curr Opin Cell Biol. 2013;25(3):378–86.
Dechat T, Adam SA, Taimen P, Shimi T, Goldman RD. Nuclear lamins. Cold Spring Harbor Perspect Biol. 2010;2(11):a000547.
Malhas AN, Lee CF, Vaux DJ. Lamin B1 controls oxidative stress responses via Oct-1. J Cell Biol. 2009;184(1):45–55.
Butin-Israeli V, Adam SA, Jain N, Otte GL, Neems D, Wiesmüller L, et al. Role of lamin b1 in chromatin instability. Mol Cell Biol. 2015;35(5):884–98.
Liu NA, Sun J, Kono K, Horikoshi Y, Ikura T, Tong X, et al. Regulation of homologous recombinational repair by lamin B1 in radiation-induced DNA damage. FASEB J. 2015;29(6):2514–25.
Tang CW, Maya-Mendoza A, Martin C, Zeng K, Chen S, Feret D, et al. The integrity of a lamin-B1-dependent nucleoskeleton is a fundamental determinant of RNA synthesis in human cells. J Cell Sci. 2008;121(Pt 7):1014–24.
Dobrzynska A, Gonzalo S, Shanahan C, Askjaer P. The nuclear lamina in health and disease. Nucleus (Austin, Tex). 2016;7(3):233–48.
Li W, Li X, Li X, Li M, Yang P, Wang X, et al. Lamin B1 overexpresses in lung adenocarcinoma and promotes proliferation in lung cancer cells via AKT pathway. Onco Targets Ther. 2020;13:3129–39.
Luo F, Han J, Chen Y, Yang K, Zhang Z, Li J. Lamin B1 promotes tumor progression and metastasis in primary prostate cancer patients. Future Oncol (Lond, Engl). 2021;17(6):663–73.
Song ZY, Chao F, Zhuo Z, Ma Z, Li W, Chen G. Identification of hub genes in prostate cancer using robust rank aggregation and weighted gene co-expression network analysis. Aging. 2019;11(13):4736–56.
Yang Z, Sun Q, Guo J, Wang S, Song G, Liu W, et al. GRSF1-mediated MIR-G-1 promotes malignant behavior and nuclear autophagy by directly upregulating TMED5 and LMNB1 in cervical cancer cells. Autophagy. 2019;15(4):668–85.
Li L, Du Y, Kong X, Li Z, Jia Z, Cui J, et al. Lamin B1 is a novel therapeutic target of betulinic acid in pancreatic cancer. Clin Cancer Res. 2013;19(17):4651–61.
Sun S, Xu MZ, Poon RT, Day PJ, Luk JM. Circulating Lamin B1 (LMNB1) biomarker detects early stages of liver cancer in patients. J Proteome Res. 2010;9(1):70–8.
Wazir U, Ahmed MH, Bridger JM, Harvey A, Jiang WG, Sharma AK, et al. The clinicopathological significance of lamin A/C, lamin B1 and lamin B receptor mRNA expression in human breast cancer. Cell Mol Biol Lett. 2013;18(4):595–611.
Yu ZY, Jiang XY, Zhao RR, Luo CJ, Ren YX, Ma ZJ, et al. Lamin B1 deficiency promotes malignancy and predicts poor prognosis in gastric cancer. Neoplasma. 2020;67(6):1303–13.
Jia Y, Vong JS, Asafova A, Garvalov BK, Caputo L, Cordero J, et al. Lamin B1 loss promotes lung cancer development and metastasis by epigenetic derepression of RET. J Exp Med. 2019;216(6):1377–95.
Prieto J, Melero I, Sangro B. Immunological landscape and immunotherapy of hepatocellular carcinoma. Nat Rev Gastroenterol Hepatol. 2015;12(12):681–700.
Krishnamoorthy M, Gerhardt L, Maleki VS. Immunosuppressive effects of myeloid-derived suppressor cells in cancer and immunotherapy. Cells. 2021;10(5):1170.
Costa A, Kieffer Y, Scholer-Dahirel A, Pelon F, Bourachot B, Cardon M, et al. Fibroblast heterogeneity and immunosuppressive environment in human breast cancer. Cancer Cell. 2018;33(3):463-79.e10.
Groth C, Hu X, Weber R, Fleming V, Altevogt P, Utikal J, et al. Immunosuppression mediated by myeloid-derived suppressor cells (MDSCs) during tumour progression. Br J Cancer. 2019;120(1):16–25.
Lv D, Wu X, Wang M, Chen W, Yang S, Liu Y, et al. Functional assessment of four novel immune-related biomarkers in the pathogenesis of clear cell renal cell carcinoma. Front Cell Dev Biol. 2021;9:621618.
de Bono J, Mateo J, Fizazi K, Saad F, Shore N, Sandhu S, et al. Olaparib for metastatic castration-resistant prostate cancer. N Engl J Med. 2020;382(22):2091–102.
Golan T, Hammel P, Reni M, Van Cutsem E, Macarulla T, Hall MJ, et al. Maintenance olaparib for germline BRCA-mutated metastatic pancreatic cancer. N Engl J Med. 2019;381(4):317–27.
Qin H, Yang Y, Jiang B, Pan C, Chen W, Diao W, et al. SOX9 in prostate cancer is upregulated by cancer-associated fibroblasts to promote tumor progression through HGF/c-Met-FRA1 signaling. The FEBS journal. 2021;288(18):5406–29.
Li T, Fu J, Zeng Z, Cohen D, Li J, Chen Q, et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020;48(W1):W509–14.
Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77(21):e108–10.
Li B, Severson E, Pignon JC, Zhao H, Li T, Novak J, et al. Comprehensive analyses of tumor immunity: implications for cancer immunotherapy. Genome Biol. 2016;17(1):174.
Tang Z, Kang B, Li C, Chen T, Zhang Z. GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res. 2019;47(W1):W556–60.
Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi B, et al. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia (New York, NY). 2017;19(8):649–58.
Chen F, Chandrashekar DS, Varambally S, Creighton CJ. Pan-cancer molecular subtypes revealed by mass-spectrometry-based proteomic characterization of more than 500 human cancers. Nat Commun. 2019;10(1):5679.
Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, et al. PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003;34(3):267–73.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102(43):15545–50.
Kirtane K, Elmariah H, Chung CH, Abate-Daga D. Adoptive cellular therapy in solid tumor malignancies: review of the literature and challenges ahead. J Immunother Cancer. 2021;9(7):e002723.
Hiltbrunner S, Mannarino L, Kirschner MB, Opitz I, Rigutto A, Laure A, et al. Tumor immune microenvironment and genetic alterations in mesothelioma. Front Oncol. 2021;11:660039.
DeNardo DG, Barreto JB, Andreu P, Vasquez L, Tawfik D, Kolhatkar N, et al. CD4(+) T cells regulate pulmonary metastasis of mammary carcinomas by enhancing protumor properties of macrophages. Cancer Cell. 2009;16(2):91–102.
Mahata B, Zhang X, Kolodziejczyk AA, Proserpio V, Haim-Vilmovsky L, Taylor AE, et al. Single-cell RNA sequencing reveals T helper cells synthesizing steroids de novo to contribute to immune homeostasis. Cell Rep. 2014;7(4):1130–42.
Maier B, Leader AM, Chen ST, Tung N, Chang C, LeBerichel J, et al. A conserved dendritic-cell regulatory program limits antitumour immunity. Nature. 2020;580(7802):257–62.
Liu Z, Zhou Q, Wang Z, Zhang H, Zeng H, Huang Q, et al. Intratumoral TIGIT(+) CD8(+) T-cell infiltration determines poor prognosis and immune evasion in patients with muscle-invasive bladder cancer. J Immunother Cancer. 2020;8(2):e000978.
Zhao X, Liu J, Ge S, Chen C, Li S, Wu X, et al. Saikosaponin A inhibits breast cancer by regulating Th1/Th2 balance. Front Pharmacol. 2019;10:624.
De Monte L, Reni M, Tassi E, Clavenna D, Papa I, Recalde H, et al. Intratumor T helper type 2 cell infiltrate correlates with cancer-associated fibroblast thymic stromal lymphopoietin production and reduced survival in pancreatic cancer. J Exp Med. 2011;208(3):469–78.
Dai Y, Ma Y, Li S, Banerjee S, Liang S, Liu Q, et al. An LMNB1 duplication caused adult-onset autosomal dominant leukodystrophy in chinese family: clinical manifestations, neuroradiology and genetic diagnosis. Front Mol Neurosci. 2017;10:215.
Padiath QS, Saigoh K, Schiffmann R, Asahara H, Yamada T, Koeppen A, et al. Lamin B1 duplications cause autosomal dominant leukodystrophy. Nat Genet. 2006;38(10):1114–23.
Hirahara K, Nakayama T. CD4+ T-cell subsets in inflammatory diseases: beyond the Th1/Th2 paradigm. Int Immunol. 2016;28(4):163–71.
Ubukata H, Motohashi G, Tabuchi T, Nagata H, Konishi S, Tabuchi T. Evaluations of interferon-γ/interleukin-4 ratio and neutrophil/lymphocyte ratio as prognostic indicators in gastric cancer patients. J Surg Oncol. 2010;102(7):742–7.
Hutchison CJ. B-type lamins in health and disease. Semin Cell Dev Biol. 2014;29(100):158–63.
Lord CJ, Ashworth A. PARP inhibitors: synthetic lethality in the clinic. Science (New York, NY). 2017;355(6330):1152–8.
Li H, Liu ZY, Wu N, Chen YC, Cheng Q, Wang J. PARP inhibitor resistance: the underlying mechanisms and clinical implications. Mol Cancer. 2020;19(1):107.
This work was supported by grants from the National Natural Science Foundation of China (81972388, 81772710).
Ethics approval and consent to participate
This study was conducted in accordance with the Declaration of Helsinki principles. It was approved by the Ethics Committee of Drum Tower Hospital, Medical School of Nanjing University, and was undertaken with the understanding and written consent of each subject.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
: Figure S1. LMNB1 expression in human tissues, cells and plasma. Figure S2. Correlation between LMNB1 mRNA expression and CD4+ Th2 infiltration in TCGA cancers. Figure S3. Relevance between mRNA expression of LMNB1 and its correlated genes. Figure S4. GO term enrichment analysis of LMNB1-related genes across all TCGA cancers. Figure S5. GSEA of LMNB1 expression in the PARD cohort of TCGA datasets.
Primer sequences used in qRT-PCR. Table S2. shRNA sequence for gene knockdown.
GESA with HALLMARK gene sets of PRAD cohort in TCGA.
GESA with KEGG gene sets of PRAD cohort in TCGA.
GESA with GO_BP gene sets of PRAD cohort in TCGA.
About this article
Cite this article
Qin, H., Lu, Y., Du, L. et al. Pan-cancer analysis identifies LMNB1 as a target to redress Th1/Th2 imbalance and enhance PARP inhibitor response in human cancers. Cancer Cell Int 22, 101 (2022). https://doi.org/10.1186/s12935-022-02467-4
- Immune infiltration
- Homologous recombination repair