Long noncoding RNA LINC01132 enhances immunosuppression and therapy resistance via NRF1/DPP4 axis in hepatocellular carcinoma
Journal of Experimental & Clinical Cancer Research volume 41, Article number: 270 (2022)
Long noncoding RNAs (lncRNAs) are emerging as critical regulators of gene expression and play fundamental roles in various types of cancer. Current developments in transcriptome analyses unveiled the existence of lncRNAs; however, their functional characterization remains a challenge.
A bioinformatics screen was performed by integration of multiple omics data in hepatocellular carcinoma (HCC) prioritizing a novel oncogenic lncRNA, LINC01132. Expression of LINC01132 in HCC and control tissues was validated by qRT-PCR. Cell viability and migration activity was examined by MTT and transwell assays. Finally, our results were confirmed in vivo mouse model and ex vivo patient derived tumor xenograft experiments to determine the mechanism of action and explore LINC01132-targeted immunotherapy.
Systematic investigation of lncRNAs genome-wide expression patterns revealed LINC01132 as an oncogene in HCC. LINC01132 is significantly overexpressed in tumor and associated with poor overall survival of HCC patients, which is mainly driven by copy number amplification. Functionally, LINC01132 overexpression promoted cell growth, proliferation, invasion and metastasis in vitro and in vivo. Mechanistically, LINC01132 acts as an oncogenic driver by physically interacting with NRF and enhancing the expression of DPP4. Notably, LINC01132 silencing triggers CD8+ T cells infiltration, and LINC01132 knockdown combined with anti-PDL1 treatment improves antitumor immunity, which may prove a new combination therapy in HCC.
LINC01132 functions as an oncogenic driver that induces HCC development via the NRF1/DPP4 axis. Silencing LINC01132 may enhance the efficacy of anti-PDL1 immunotherapy in HCC patients.
Hepatocellular carcinoma (HCC) is a highly aggressive hepatic malignancy with poor survival rate . The most common therapeutic is surgical resection or liver transplantation . However, patients are usually diagnosed at an advanced stage and are not suitable for surgical treatment. Therefore, systematic investigation of the underlying mechanism associated with HCC development and progression is of high clinical significance and may lead the development of novel clinical options .
With the development of high-throughput sequencing technologies, numerous molecular markers have been linked to the development of HCC. Genes from various signaling pathways are frequently mutated in HCC, such as Wnt, P53, AKT and MAP kinase pathways . Copy number alterations have been found to be involved in the development of HCC, such as MET and PEG10 amplification [5, 6]. Cancer-related mutations have been found to perturb the RNA regulatory network and helped the identification of potential biomarkers . However, the molecular pathogenesis of HCC is still not fully understood, and novel cancer-promoting genes must be identified and characterized.
Current progress in transcriptome analysis exposed a large portion of the human genome that does not encode for proteins . Long noncoding RNAs (lncRNAs) have been discovered as a major type of regulatory RNA with important roles in cancer development . Accumulating evidence has shown that lncRNAs are involved in a wide range of biological processes, acting as scaffolds or miRNA sponges [10, 11]. For example, LINC01138 has been found to drive malignancies via activating arginine methyltransferase 5 in HCC ; to physically interact with the MYC protein and increase its stability in cancer ; and inflammation induced LINC00665 is involved in the NF-kB signaling activation in HCC . Many lncRNAs were identified as novel regulatory RNAs in HCC but their biological functions and underlying mechanism in pathogenesis remain largely unclear.
The crosstalk between the immune system and tumor cells is critical in cancer development and progression . Following the success story of immune checkpoint blockers (ICBs) therapy in various cancer types, increasing efforts were devoted to investigate novel ICBs approaches in HCC patients . Despite this breakthrough, a subset of patients remains as non- ideal candidates for immunotherapy due to the lack of efficacy [17, 18]. Therefore, numerous studies focused on modifying the expression of noncoding RNAs in combination with immunotherapy to improve the response and overall survival. LIMIT is an immunogenic lncRNA in cancer immunity targetable for cancer immunotherapy , while NKILA, another lncRNA, can promote tumor immune evasion by sensitizing T cells and activate induced cell death . Furthermore, pan-cancer analysis of immune-related lncRNAs has prioritized cancer-related lncRNAs and the identification of immune subtypes . Therefore, exploring the roles of lncRNAs in immune regulation is essential to identify additional immunotherapy targets in cancer.
In this study, we identified a new oncogenic long intergenic noncoding RNA (lincRNA) in HCC (LINC01132). LINC01132 is overexpressed in HCC, and significantly associated with malignant clinical features and poor outcomes in the clinic. Mechanistically, LINC01132 promotes cell growth, proliferation, invasion and metastasis through the LINC0132/NRF1/DPP4 axis. Finally, combinatory therapy targeting LINC01132 inhibition and anti-PDL1 blockade synergistically improves antitumor immunity in HCC in vivo and ex vivo models. In summary, our data suggests that LINC01132 is a potential biomarker and therapeutic target for HCC.
Materials and methods
Patients and ethical statement
In total, 121 HCC patients’ tumor samples and corresponding adjacent normal liver tissues were obtained from the surgical specimen archives of the Zhongshan Hospital , Shanghai, China. The patients were informed, and signed consent forms acknowledging the use of their resected tissues for research purposes, which has been previously approved .
The HCC cell lines HepG2 and Hep1–6 were purchased from American Type Culture Collection (ATCC, Manassas, Virginia, USA) and cultured following the recommended guidelines. These cells were characterized by Genewiz Inc. cultured in Dulbecco’s Modified Eagle’s Medium (DMEM) (Thermo Fisher Scientific, California, USA) with 10% fetal bovine serum (FBS) and antibiotics. Huh-7 cells were was purchased from JCRB cell bank (Tokoyo, Japan) and cultured in Roswell Park Memorial Institute (RPMI) 1640 Medium (Thermo Fisher Scientific) with 10% new born calf serum. Hep3B and SUN-449 were purchased from the Cell Bank of Chinese Academy of Sciences (Shanghai, China). Cells were maintained in Dulbecco’s modified eagle medium (DMEM, Invitrogen, Carlsbad, CA, USA) supplemented with 10% fetal bovine serum (FBS, HyClone, Logan, UT, USA), 1% penicillin/streptomycin (pen/strep, Invitrogen), and 8 mg/L of the antibiotic, tylosin tartrate, for mycoplasma (Sigma-Aldrich, St. Louis, Missouri, USA), at 37 °C in 5% CO2 (v/v). All cell lines were authenticated by autosomal STR profiling and thawed afresh every 2 months, to test for mycoplasma. None of the cell lines used was found in the database of commonly misidentified cell lines, maintained by the International Cell Line Authentication Committee.
In vivo mouse models
All studies were supervised and approved by the Shanghai University of Traditional Chinese Medicine Institutional Animal Care and Use Committee (IACUC). Female mice were used as models to study liver cancer. Power analysis indicated an n value of 5 mice per group to identify the expected effects with 90% confidence.
Total RNA was extracted from liver samples or HepG2 or Huh-7 cells with TRIzol Reagent (Thermo Fisher Scientific, California, USA). Quantitative real-time PCR (qPCR) was performed with an iQ5 machine and SYBR Premix Ex Taq™ II (TaKaRa Bio, Tokyo, Japan). Data were normalized to β-actin or IgG control (RNA pull-down assay and RIP assay). Relative genomic level of tumor tissues was compared with normal liver tissue. Primers used for RT-qPCR and RT-PCR are described in additional file 1: Table S1.
Cytoplasmic and nuclear RNA isolation
Cytoplasmic and nuclear RNA was extracted using Thermo Fisher BioReagents (Thermo Fisher Scientific) according to the manufacturer’s instructions. QRT-PCR analysis was performed using SYBR® Green Master Mix (Invitrogen, New York, USA) to amplify the localization of LINC01132 assay and β-actin U6 were used as cytoplasmic and nucleus controls. Primers are listed in additional file 1: Table S1.
Lentivirus construction and infection
The LINC01132 sequence (listed in additional file 1: Table S1) was amplified from normal genomic RNA and cloned into the pWPXL lentiviral vector to generate pWPXL-LINC01132. Virus particles were harvested 48 hrs after HEK 293 T cells were transfected with pWPXL-LINC01132, with the packaging plasmid psPAX2 and the VSV-G envelope plasmid pMD2.G using Lipofectamine 2000 reagent (Invitrogen). Huh-7, HepG2, Hep3B, SNU-449 and Hep1–6, cells were infected with recombinant lentivirus-transducing units in the presence of 1 μg/ml polybrene (Sigma-Aldrich, Missouri, USA).
LINC01132 adenovirus interference vector (AD-LINC01132-SHRNA)
Three pairs of LINC01132-shRNA oligonucleotide sequences were designed and three adv4-U6-CMV-LINC01132-shRNA adenovirus shuttle plasmids were constructed. The constructed ADV4-U6-CMV-LINC01132-shRNA and pHBAd-BHG plasmids were co-transfected into 293A cells and packaged with AD-LINC01132-shRNA. The virus titer was detected by microcytic whole-cell assay. Ad-linc01132-shrna was injected subcutaneously and the mRNA expression of LINC01132 was detected by real-time fluorescence quantitative PCR.
In vitro cell proliferation and colony formation assays
The cell proliferation assay was measured using the Cell Counting Kit-8 (CCK-8) (Dojindo, Kyushu, Japan) at 1, 3, and 5 days after LINC01132 or mock infection. Cells were seeded into the 96-well plate at a density of 103 cells/well, and 10 μl of CCK-8 was added to 90 μl of the cell culture medium per well. Cells were subsequently incubated at 37 °C for 2 hours and the optical density measured at 450 nm. For the colony formation assay, 1–1.5 *103 cells/well were plated in a 6-well plate and incubated at 37 °C for 2 weeks. The colonies were fixed and stained with 0.1% crystal violet dye in 20% methanol, and the number of colonies macroscopically counted. All assays were performed in triplicate.
In vitro migration and invasion assays
Migration assays were performed in a Transwell chemotaxis 24-well chamber (BD Biosciences, Franklin Lakes, NJ). Briefly, 2 × 104 cells were plated in the upper chamber with a non-coated membrane. For the invasion assays, 5 × 104 cells were placed into the upper chamber with a Matrigel-coated membrane. After 16 hours of incubation at 37 °C, migrating or invading cells were fixed and stained with 0.1% crystal violet dye in 20% methanol. Migrated or invaded cells were counted and imaged with an inverted microscope (Olympus, Tokyo, Japan).
In vivo assays for metastasis
For the in vivo metastasis assays, 2 × 106 Hep1–6 cells infected with the pWPXL-LINC01132 or pWPXL-GFP were resuspended in 0.2 mL of serum-free DMEM and subcutaneously injected into C57BL/6 mice liver. After 40 days mice humanly euthanized and liver, lungs and intestines were collected, fixed with phosphate buffered neutral formalin and prepared for standard histological examination. The numbers of metastatic foci in liver tissue sections were counted by H&E staining under a binocular microscope (Leica, USA).
Construction of PDX mouse model of HCC
PDX model of liver cancer was established from patient’s liver cancer samples. These were first sectioned into small tissue for subcutaneous tumor formation in nude mice. Two weeks later, the subcutaneous tumor bearing tissue was excised, sectioned and and transplanted again into the subcutaneous skin of nude mice. After subcutaneous tumor formation, shNC and shLINC01132 were injected subcutaneously and tumor diameter and weight were measured every 2 days. Animals were humanly euthanized after 30 days and tumor tissue removed for terminal endpoint analysis.
RNA pull-down assays
LINC01132 was transcribed in vitro with biotin RNA labelling mix and T7 RNA polymerase according to the manufacturer’s instructions (Invitrogen). In total, 40 μl streptavidin-linked magnetic beads (ThermoFisher Scientific) were used to pull down the biotinylated RNA at room temperature for 2 hrs. The beads-RNA-proteins were then washed with 1× binding washing buffer (5 mM Tris-HCl, 1 M NaCl, 0.5 mM EDTA, and 0.005% Tween 20) four times. The proteins were precipitated and diluted in 60 μl protein lysis buffer. Finally, the retrieved proteins were measured on SDS-PAGE gels for mass spectrometry or Western blot. Western blot in RNA pull-down assay was performed with mouse anti-NRF1 and anti-KDM5B antibodies (Cell Signaling Technology, CST, Danvers, Massachusetts, USA, 1:500) and mouse anti-β-actin antibody (CST, 1:1000). Antibodies information is available on additional file 1: Table S2.
RNA-seq and computational analysis
RNA-seq was performed at the Sequencing and Non-Coding RNA Program at the RiboBio (Guangzhou, China) using Hiseq3000(Illumina, USA). The hisat2 , StringTie  and Ballgown were used to align the reads to the genome, generate raw counts corresponding to each known gene, and calculate the RPKM (reads per kilobase per million) values .
LINC01132 levels were measured by northern blot using an Ambion Northern Max-Gly Kit (Austin, TX, USA). Total RNA was electrophoresed and siphoned to a positively charged nylon membrane (NC). RNA was then fixed to the NC membrane using UV cross-linking. In brief, the cross-linked membrane was then prehybridized with ULTRAhyb, and RNA was detected with an LINC01132-specific oligonucleotide probe (primers provided in additional file 1: Table S1) labeled with digoxigenin-ddUTP using a DIG Oligonucleotide 3′-End Labeling Kit (Roche Diagnostics, Indianapolis, IN, USA) in roller bottles.
Proteins were separated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and transferred to nitrocellulose membrane (Bio-Rad, Hercules, CA). Non-specific binding was blocked with 5% nonfat milk and subsequently incubated with indicated primary antibodies, followed by horseradish peroxidase-conjugated secondary antibodies. Immunoreactivity was visualized with chemiluminescence ECL reagents (Pierce, Rockford, IL) and imaged with ChemiDoc imaging system (12003153-s). Densitometry analysis was performed with Image-Pro Plus 6.0 (Media Cybernetics).
Mass spectrometry analysis
Specific bands were excised for proteomics screening by mass spectrometry analysis (Shanghai Applied Protein Technology, Shanghai, China). Protein identification was retrieved from the human RefSeq protein database (National Center for Biotechnology Information), using Mascot version 2.4.01 (Matrix Science, London, UK). The retrieved protein was detected by western blot.
RNA immunoprecipitation (RIP)
RIP experiments were performed using the Magna RIP™ RNA-Binding Protein Immunoprecipitation Kit (Millipore, Massachusetts, USA) according to the manufacturer’s instructions. The co-precipitated RNAs were detected by reverse transcription PCR. Total RNA (input controls) and normal mouse IgG controls were assayed simultaneously to validate RNA specificity to NRF1 and KDM5B (n = 3 for each experiment). Gene-specific primers for LINC01132 are provided in additional file 1: Table S1).
Huh-7 and HepG2 cells infected with the lentivirus expressing LINC01132 or Huh-7 cells transfected with LINC01132 siRNA were lysed with RIPA buffer (Beyotime Biotechnology) with protease (Thermo Fisher Scientific Inc.) and RNase (Thermo Fisher Scientific Inc.) inhibitors, and then centrifuged at 16,400 g for 15 min. Supernatants were collected and the amount of NRF1 and DPP4 or KDM5B and DPP4 protein examined by immunoblotting to normalized for DPP4 and NRF1 loading. Supernatants were then incubated with the indicated antibody-coated Protein G Dynabeads (Life Technologies) overnight at 4 °C with gentle rotation. The beads were washed five times with NT2 buffer (50 mM Tris-HCL [pH 7.4], 150 mM NaCl, 1 mM MgCl2, 0.05% Nonidet P-40) containing protease and RNase Inhibitor, and then three times with PBS again containing protease and RNase Inhibitor. After washing, proteins were eluted and immuno-complexes were analyzed by immunoblotting (antibodies provided in additional file 1: Table S2).
Genome-wide expression of lncRNAs in HCC
Genome-wide expression of lncRNAs in HCC was collected from Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA). Gene expression profile of HCC was collected from GEO under the accession number GSE104310, which provided the expression of tumors and paired non-tumor tissues collected from Sun Yat-sen University Cancer Center. In total, eight normal and 12 tumor samples were included in our analysis. Moreover, genome-wide expression of liver hepatocellular carcinoma (LIHC) were obtained from TCGA project , which included 374 tumor and 50 normal samples. RNA sequencing of 21 Hepatitis B virus-HCC patients with non-neoplastic liver and tumor tissues were obtained from GSE94660 .
Identification of differentially expressed lncRNAs in HCC
Wilcoxon’s rank sum test was used to identify differentially expressed lncRNAs in HCC tumor compared to normal tissues. Only long intergenic non-coding RNAs (lincRNAs) were considered in our analyses. We identified the lincRNAs with fold changes > 2 or < 0.5 and adjusted p-values < 0.05 as differentially expressed lincRNAs.
Predict the functions of lncRNA
Spearman correlation coefficients (SCC) were used to predict the function of lncRNA and protein-coding gene expression. All genes were ranked based on SCC and subjected into Gene Set Enrichment Analysis (GSEA) . The Biocarta pathway dataset from MsigDB was used for this analysis . Pathways with normalized enrichment score > 1.96 or < − 1.96 and false discovery rate (FDR) < 0.01 were identified.
Transcription regulation analysis
To identify the transcription factors (TFs) that bind to the promoter region of gene (e.g., DPP4), we queried the ChIPBase database . The TFs than can bind to the 5 kb upstream to 1 kb downstream of transcript start sites were identified.
Integrative analyses identify onco-lncRNA LINC01132 in HCC
To identify the potential oncogenic lncRNAs in HCC, we first analyzed the genome-wide expression profile of 12 tumors and 8 adjacent non-tumor tissues. In total, we identified 22 up-regulated long intergenic noncoding RNAs (lincRNAs) and 35 down-regulated lincRNAs in HCC (Fig. 1A and additional file 1: Table S3). Numerous tumor associated lincRNAs were identified, such as PVT1 , HIF1A-AS1  and MIAT . Moreover, the expression of these lincRNAs can effectively distinguish the tumor tissue from the normal controls in HCC (Fig. 1B).
Next, we investigated the genetic alterations of differentially expressed lincRNAs. The majority of highly expressed lincRNAs were associated with copy number amplification in HCC (Fig. 1C). For example, approximate 11% of HCC patients had PVT1 CNV amplification, which might induce the expression of PVT1 in HCC. In the contrast, lower number of patients had copy number deletion for the down-regulated lincRNAs (Additional file 2: Fig. S1A). We then focused on the lincRNA-LINC01132, which exhibited the second higher copy number amplification (7%) in HCC (Fig. 1C) and was highly expressed in HCC tissue. Moreover, LINC01132 exhibited higher copy number amplification across various cancer types (Additional file 2: Fig. S1B).
To further validate the oncogenic roles of LINC01132, we analyzed the expression in another two cohorts of HCC. We found that LINC01132 exhibited significantly higher expression in tumor patients (Additional file 2: Fig. S1C-D). We also evaluated the expression of LINC01132 in 121 paired tumor tissues and corresponding non-cancerous tissues (NCT), showing that LINC01132 was highly expressed in cancer (Fig. 1D, p < 0.0001). Patients with more than one cancer embolus (Fig. 1E, p = 0.0032) or larger tumor size (Fig. 1F, p = 0.0321) had even higher LINC01132 expression levels. Finally, we explored the effects of LINC01132 on prognosis in HCC patients and observed that high LINC01132 expression was associated with poor survival of HCC (Fig. 1G, p = 0.0143). Altogether, these results suggested that LINC01132 plays an oncogenic role in HCC.
LINC01132 promotes HCC cell growth and metastasis in vitro
To investigate the possible roles of LINC01132 in HCC pathogenesis, we first explored the expression in HCC cell lines. We found that LINC01132 was with relatively high expressed in HepG2 and but not in Huh-7 (Additional file 2: Fig. S2A). We also observed that the expression of LINC01132 was in both cytoplasm and nuclear (Additional file 2: Fig. S2B), suggesting that it may play a regulatory function in both cellular localizations. Moreover, we found that LINC01132 didn’t express to protein but functioned as RNA (Additional file 2: Fig. S2C-D).
We next knocked down or overexpressed LINC01132 in HCC cell lines and performed colony and cell proliferation assays. We found that the knocked down or overexpressed LINC01132 can effectively alter its expression in cell lines (Additional file 2: Fig. S2E). LINC01132 knocked down significantly decreased the number of colony formation (Fig. 2A, p < 0.001) while overexpression of LINC01132 significantly increased the number of colony formation (Fig. 2B, p < 0.001). Moreover, knockdown of LINC01132 significantly decreased cell proliferation (Fig. 2C, p < 0.001), while LINC01132 overexpression significantly promoted the cancer cells proliferation (Fig. 2D, p < 0.001).
We next explored the roles of LINC01132 in cell invasion and migration. LINC01132 knockdown significantly decreased cell migration and invasion (Fig. 2E-F, p-values < 0.001). In the contrast, overexpression of LINC01132 significantly increased cell migration and invasion (Fig. 2G, p-values < 0.001). In addition, we validated the biological roles of LINC01132 in another two cell lines (Hep3B and SNU-449). We found that the results were consistent across different cell lines (Additional file 2: Fig. S3). Taken together, these results demonstrated that LINC01132 significantly promotes in vitro cell growth, proliferation, migration and invasion in HCC.
LINC01132 promotes HCC cell growth and metastasis in vivo
To confirm the functions of LINC01132 on the tumorigenicity of HCC, LINC01132 knockdown cells and control cells derived from cell lines were subcutaneously injected into nude immunodeficient mice. Tumor xenografts derived from LINC01132-knockdown cells exhibited smaller volumes and lower weights than those from empty vector-transduced cells (Fig. 3A-C). In addition, we assessed the effects of LINC01132 on metastasis, observing that the numbers of metastatic nodules were significantly increased in LINC01132-overexpressing groups (Fig. 3D, p = 0.0092).
The patient-derived xenograft (PDX) mouse model has been shown to recapitulate multiple characteristics of human cancer biological context. Therefore, we next silenced LINC01132 on PDX mouse model 7 days post inoculation (Fig. 3E). LINC01132 knockdown significantly inhibited cancer growth as shown by decreased tumor volumes and weights (Fig. 3E, p = 0.0028 and 0.039). As the transfection efficiency of plasmid vector was low, and adenovirus vector was constructed to validate the function of LINC01132. We found that silence of LINC01132 induced significantly decreased tumor volumes and weights (Additional file 2: Fig. S4). Together, these results indicated that elevated LINC01132 expression may contribute to the development and progression of HCC.
LINC01132 potentially regulates DPP4 in HCC
To elucidate the potential molecular mechanisms of LINC01132 in HCC, we calculated the correlation between protein expression and expression of LINC01132 in TCGA HCC cohort . LINC01132 expression was significantly positively correlated with SETD2, PRDX1 and CD26, while negatively correlated with AKT and BRD4 (Fig. 4A). Deficiency of histone methyltransferase SET Domain-Containing 2 (SETD2) in liver has been demonstrated to lead to abnormal lipid metabolism and HCC . Similarly, PRDX1 can act as a pro-cancer protein in HCC HepG2 cells . The expression of LINC01132 was consistently correlated with CD26 protein and RNA expression in HCC (Fig. 4B). CD26/dipeptidyl peptidase (DPP) 4 is a membrane-bound protein found in many cell types and has been suggested as a potential biomarker and target for cancer therapy .
Next, we performed GSEA analysis based on the expression correlation between LINC01132 and all protein coding genes and found that LINC01132 expression was positively correlated with EIF2 and ETC pathways, while negatively correlated with TCRA and NKT pathways (Fig. 4C). EIF2 pathway has been associated with various types of cancer, including HCC [37, 38]. The observed negative correlation with immune-related pathways suggests that LINC01132 might suppress immune pathways to promote the cancer development and progression.
Moreover, we in two cell lines and performed RNA-Seq for genome-wide expression profiling. Potential correlated genes such as PRDX1, DPP4 (CD26) and AKT2 exhibited significant expression changes in, LINC01132 knockdown cells (Fig. 4D). GSEA analysis revealed that highly expressed genes LINC01132 knockdown cells were significantly involved in COMP and STEM pathways (Fig. 4E-F, false discovery rates = 0.030 and < 0.001). In particular, we found that LINC01132 knockdown increased the expression of numerous immune-related genes, such as CD8A, IL6 and IL11 (Fig. 4F). We further investigated the expression levels of DPP4 in LINC01132-knockdown or LINC01132-overexpressing HCC cells. DPP4 transcriptional and protein levels were increased LINC01132-overexpressing and decreased in LINC01132-knockdown HCC cells, respectively (Fig. 4G-H). These results suggest that LINC01132 might exert its oncogene functions by modulating the DPP4 signaling pathway.
LINC01132 interferes with NRF1 binding to DPP4 in HCC
To further explore the molecular mechanism underlying the oncogenic activity of LINC01132 in HCC, we performed RNA pull-down assay to identify the proteins associated with LINC01132. LINC01132 pull-down identified 598 potential interacting proteins (Fig. 5A). Of these proteins, NRF1 and KDM5B, and one transcription cofactor (CDK8) have binding sites around the transcription start site (TSS) of DPP4 gene (Fig. 5A). Moreover, we confirmed that sense but not antisense LINC01132 was specifically associated with NRF1 and KDM5B (Fig. 5B-C). RIP assays further showed that NRF1 and KDM5B antibodies could significantly enrich LINC01132 whereas the GAPDH antibody and IgG control could not (Fig. 5D-E). The expression of NRF1 and KDM5B was then investigated showing that higher grade HCC patients exhibited significantly higher expression of NRF1 (Additional file 2: Fig. S5A-B). Survival analysis revealed that poor survival of patients with higher expression of NRF1, KDM5B, LINC01132 and DPP4 (log-rank p = 0.038, Additional file 2: Fig. S5C).
Furthermore, we retrieved the public chromatin immunoprecipitation (ChIP) sequencing data from ChIPBase  and found that the NRF1 and KDM5B can bind to the promoter region of DPP4 across cell lines (Fig. 5F and Additional file 2: Fig. S6). Additionally, LINC01132 knockdown significantly blocked the interaction between NRF1 and DPP4 (Fig. 5G). However, the effects were not observed for KDM5B. These results indicated that LINC01132 acts as a scaffold in the interaction between NRF1 and DPP4.
LINC01132 improves the response to anti-PD1 immunotherapy in HCC
Previous studies have identified DPP4 decreased chemokines and other immune molecules . Our data revealed that LINC01132 overexpression was significantly associated with decreased immune pathway activity. We also identified a negative relationship between DPP4 and CD8+ T cell infiltration levels in HCC tissues (Fig. 6A). Combination therapy with anti-PDL1 blockade immunotherapy and other therapies have been shown to improve the efficacy of the tumor-specific T-cell response [40, 41]. Based on the above results of LINC01132, we predicted that LINC01132 knockdown could enhance lymphocyte trafficking and improve tumor responses to PDL1 blockage in HCC. Thus, we investigated the combination therapy of LINC01132 knockdown and PDL1 inhibitor in the Hep1–6-shLINC01132 tumor model (Fig. 6B).
We found that shLINC01132 therapy resulted in delayed tumor growth, smaller tumor volume and weight (Fig. 6B-D). Importantly, much clearer tumor regression was observed in the group treated with shLINC01132 plus anti-PDL1 (Fig. 6B-D). The number of CD8+ T cells was significantly increased in the tumors of mice treated with shLINC01132 compared with control (Fig. 6E, p = 0.012). Moreover, the increase was even higher in the tumor of mice treated with the shLINC01132 and PDL1 blockage (Fig. 6E, p = 0.0098). Thus, the above results further indicated that LINC01132/NRF1/DPP4 axis is involved in the immunosuppression of HCC (Fig. 6F) and suggested that knockdown LINC01132 could improve the efficacy of PDL1 blockage immunotherapy.
Rapid progresses in high throughput sequencing technologies have successfully identified a large number of lncRNAs. Moreover, emerging evidence has indicated that lncRNAs are expressed in a cell type-specific or tissue-specific patterns [42, 43], suggesting important roles in diverse biological processes. Expression perturbation of lncRNAs has been observed in various cancer types. By integration analysis of genome-wide lncRNA expression and genetic alterations, we revealed that LINC01132 is up-regulated in HCC tissues and that its high expression might be driven by copy number amplification. Moreover, overexpression of LINC01132 was associated with poor prognosis of HCC patients.
LINC01132 promoted cell growth, proliferation, invasion and metastasis in vitro and in vivo. Functional analysis revealed that LINC01132 physically interacted with NRF1 and KDM5B and promoted the expression of DPP4. To date, expression perturbation of LINC01132 has only been associated with oncogenic activities in ovarian cancer by regulating miR-431-5p/SOX9 axis  and involved in hypoxia regulation in glioblastoma . Therefore, the roles of LINC01132 in the regulation of NRF1/DPP4 axis in HCC are here first described. These results suggested that the same lncRNA can play diverse functions in regulating various pathways in different tumor context.
CD26/DPP4 is a membrane-bound protein, and its higher expression has been found in a wide variety of tumor pathologies . CD26 expression was significantly increased in tumor HCC specimens and was associated with larger tumor size . CD26 has been proven as a pro-oncogenic gene in HCC and a potential therapeutic target. We also found that DPP4 interacts with number of cancer-related genes (Additional file 2: Fig. S7). Moreover, DPP4 inhibition improves antitumor effect of PD1 in HCC by enhancing CD8+ T cell infiltration . It is increasingly clear that there are widespread changes in lncRNA expressions during the immune response. Numerous lncRNAs, such as NEAT1, UCA1, MIR22HG, and LINK-A, have involved in immune regulation in cancer [48,49,50]. Here, we demonstrated the inhibition of LINC01132 can achieve the similar antitumor effects. Altogether we expose the critical role of LINC01132/NRF1/DPP4 promoting the development of HCC.
In conclusion, our results demonstrated that LINC01132 may act as an onco-lncRNA and overexpression of LINC01132 promoted HCC development via the NRF1/DPP4 signaling axis. Contrastingly, LINC01132 silencing may be a novel synergistic strategy to improve the efficacy of PDL1 inhibitor therapy in a subgroup of ICB resistant HCC patients.
Availability of data and materials
All data generated or analysed during this study are included in this published article. The gene expression profiles and clinical data can be found at the GDC portal and GEO (https://www.ncbi.nlm.nih.gov/geo/). Software and resources used for the analyses are described in each method section.
Dipeptidyl Peptidase 4
Gene Set Enrichment Analysis
Long non-coding RNAs
Lysine Demethylase 5B
Long Intergenic Non-Protein Coding RNA 1132
Nuclear Respiratory Factor 1
Patient-derived tumor xenograft mouse model
Short hairepin RNA or small hairpin RNA
SRY-Box Transcription Factor 9
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71:209–49.
Maluccio M, Covey A. Recent progress in understanding, diagnosing, and treating hepatocellular carcinoma. CA Cancer J Clin. 2012;62:394–9.
Ozer Etik D, Suna N, Boyacioglu AS. Management of Hepatocellular Carcinoma: prevention, surveillance, diagnosis, and staging. Exp Clin Transplant. 2017;15:31–5.
Rebouissou S, Nault JC. Advances in molecular classification and precision oncology in hepatocellular carcinoma. J Hepatol. 2020;72:215–29.
Bang H, Ha SY, Hwang SH, Park CK. Expression of PEG10 is associated with poor survival and tumor recurrence in hepatocellular carcinoma. Cancer Res Treat. 2015;47:844–52.
Ke AW, Shi GM, Zhou J, Wu FZ, Ding ZB, Hu MY, et al. Role of overexpression of CD151 and/or c-met in predicting prognosis of hepatocellular carcinoma. Hepatology. 2009;49:491–503.
Li Y, McGrail DJ, Xu J, Li J, Liu NN, Sun M, et al. MERIT: systematic analysis and characterization of mutational effect on RNA Interactome topology. Hepatology. 2019;70:532–46.
Frankish A, Diekhans M, Jungreis I, Lagarde J, Loveland JE, Mudge JM, et al. Gencode 2021. Nucleic Acids Res. 2021;49:D916–23.
Goodall GJ, Wickramasinghe VO. RNA in cancer. Nat Rev Cancer. 2021;21:22–36.
Peng WX, Koirala P, Mo YY. LncRNA-mediated regulation of cell signaling in cancer. Oncogene. 2017;36:5661–7.
Huang Z, Zhou JK, Peng Y, He W, Huang C. The role of long noncoding RNAs in hepatocellular carcinoma. Mol Cancer. 2020;19:77.
Li Z, Zhang J, Liu X, Li S, Wang Q, Di C, et al. The LINC01138 drives malignancies via activating arginine methyltransferase 5 in hepatocellular carcinoma. Nat Commun. 2018;9:1572.
Zhang J, Li S, Zhang L, Xu J, Song M, Shao T, et al. RBP EIF2S2 promotes tumorigenesis and progression by regulating MYC-mediated inhibition via FHIT-related enhancers. Mol Ther. 2020;28:1105–18.
Ding J, Zhao J, Huan L, Liu Y, Qiao Y, Wang Z, et al. Inflammation-induced long intergenic noncoding RNA (LINC00665) increases malignancy through activating the double-stranded RNA-activated protein kinase/nuclear factor kappa B pathway in hepatocellular carcinoma. Hepatology. 2020;72:1666–81.
Schreiber RD, Old LJ, Smyth MJ. Cancer immunoediting: integrating immunity's roles in cancer suppression and promotion. Science. 2011;331:1565–70.
Pinter M, Scheiner B, Peck-Radosavljevic M. Immunotherapy for advanced hepatocellular carcinoma: a focus on special subgroups. Gut. 2021;70:204–14.
Zongyi Y, Xiaowu L. Immunotherapy for hepatocellular carcinoma. Cancer Lett. 2020;470:8–17.
Llovet JM, Castet F, Heikenwalder M, Maini MK, Mazzaferro V, Pinato DJ, et al. Immunotherapies for hepatocellular carcinoma. Nat Rev Clin Oncol. 2022;19:151–72.
Li G, Kryczek I, Nam J, Li X, Li S, Li J, et al. LIMIT is an immunogenic lncRNA in cancer immunity and immunotherapy. Nat Cell Biol. 2021;23:526–37.
Huang D, Chen J, Yang L, Ouyang Q, Li J, Lao L, et al. NKILA lncRNA promotes tumor immune evasion by sensitizing T cells to activation-induced cell death. Nat Immunol. 2018;19:1112–25.
Li Y, Jiang T, Zhou W, Li J, Li X, Wang Q, et al. Pan-cancer characterization of immune-related lncRNAs identifies potential oncogenic biomarkers. Nat Commun. 2020;11:1000.
Kim D, Paggi JM, Park C, Bennett C, Salzberg SL. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. 2019;37:907–15.
Pertea M, Pertea GM, Antonescu CM, Chang TC, Mendell JT, Salzberg SL. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol. 2015;33:290–5.
Pertea M, Kim D, Pertea GM, Leek JT, Salzberg SL. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat Protoc. 2016;11:1650–67.
Cancer Genome Atlas Research Network. Electronic address wbe, Cancer genome atlas research N: comprehensive and integrative genomic characterization of hepatocellular carcinoma. Cell. 2017;169:1327–1341 e1323.
Yoo S, Wang W, Wang Q, Fiel MI, Lee E, Hiotis SP, et al. A pilot systematic genomic comparison of recurrence risks of hepatitis B virus-associated hepatocellular carcinoma with low- and high-degree liver fibrosis. BMC Med. 2017;15:214.
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 U S A. 2005;102:15545–50.
Liberzon A, Birger C, Thorvaldsdottir H, Ghandi M, Mesirov JP, Tamayo P. The molecular signatures database (MSigDB) hallmark gene set collection. Cell Syst. 2015;1:417–25.
Zhou KR, Liu S, Sun WJ, Zheng LL, Zhou H, Yang JH, et al. ChIPBase v2.0: decoding transcriptional regulatory networks of non-coding RNAs and protein-coding genes from ChIP-seq data. Nucleic Acids Res. 2017;45:D43–50.
Luo Z, Cao P. Long noncoding RNA PVT1 promotes hepatoblastoma cell proliferation through activating STAT3. Cancer Manag Res. 2019;11:8517–27.
Hong F, Gao Y, Li Y, Zheng L, Xu F, Li X. Inhibition of HIF1A-AS1 promoted starvation-induced hepatocellular carcinoma cell apoptosis by reducing HIF-1alpha/mTOR-mediated autophagy. World J Surg Oncol. 2020;18:113.
Huang X, Gao Y, Qin J, Lu S. lncRNA MIAT promotes proliferation and invasion of HCC cells via sponging miR-214. Am J Physiol Gastrointest Liver Physiol. 2018;314:G559–65.
Chen MM, Li J, Wang Y, Akbani R, Lu Y, Mills GB, et al. TCPA v3.0: an integrative platform to explore the Pan-Cancer analysis of functional proteomic data. Mol Cell Proteomics. 2019;18:S15–25.
Steen CB, Luca BA, Esfahani MS, Azizi A, Sworder BJ, Nabet BY, et al. The landscape of tumor cell states and ecosystems in diffuse large B cell lymphoma. Cancer Cell. 2021;39:1422–1437 e1410.
Aguilar-Melero P, Prieto-Alamo MJ, Jurado J, Holmgren A, Pueyo C. Proteomics in HepG2 hepatocarcinoma cells with stably silenced expression of PRDX1. J Proteome. 2013;79:161–71.
Enz N, Vliegen G, De Meester I, Jungraithmayr W. CD26/DPP4 - a potential biomarker and target for cancer therapy. Pharmacol Ther. 2019;198:135–59.
Ji P, Wang H, Cheng Y, Liang S. Prognostic prediction and gene regulation network of EIF2S2 in hepatocellular carcinoma based on data mining. J Gastrointest Oncol. 2021;12:3061–78.
Zhang J, Li S, Zhang L, Xu J, Song M, Shao T, et al. RBP EIF2S2 promotes tumorigenesis and progression by regulating MYC-mediated inhibition via FHIT-related enhancers. Mol Ther. 2021;29:886.
Hollande C, Boussier J, Ziai J, Nozawa T, Bondet V, Phung W, et al. Inhibition of the dipeptidyl peptidase DPP4 (CD26) reveals IL-33-dependent eosinophil-mediated control of tumor growth. Nat Immunol. 2019;20:257–64.
Hayashi H, Nakagawa K. Combination therapy with PD-1 or PD-L1 inhibitors for cancer. Int J Clin Oncol. 2020;25:818–30.
Xu J, Shao T, Song M, Xie Y, Zhou J, Yin J, et al. MIR22HG acts as a tumor suppressor via TGFbeta/SMAD signaling and facilitates immunotherapy in colorectal cancer. Mol Cancer. 2020;19:51.
Lv D, Xu K, Jin X, Li J, Shi Y, Zhang M, et al. LncSpA: LncRNA spatial atlas of expression across normal and cancer tissues. Cancer Res. 2020;80:2067–71.
Xu K, Cai Y, Zhang M, Zou H, Chang Z, Li D, et al. Pan-cancer characterization of expression and clinical relevance of m (6) A-related tissue-elevated long non-coding RNAs. Mol Cancer. 2021;20:31.
Zhu W, Xiao X, Chen J. Silencing of the long noncoding RNA LINC01132 alleviates the oncogenicity of epithelial ovarian cancer by regulating the microRNA4315p/SOX9 axis. Int J Mol Med. 2021;48.
Bao L, Chen Y, Lai HT, Wu SY, Wang JE, Hatanpaa KJ, et al. Methylation of hypoxia-inducible factor (HIF)-1alpha by G9a/GLP inhibits HIF-1 transcriptional activity and cell migration. Nucleic Acids Res. 2018;46:6576–91.
Kawaguchi T, Kodama T, Hikita H, Makino Y, Saito Y, Tanaka S, et al. Synthetic lethal interaction of combined CD26 and Bcl-xL inhibition is a powerful anticancer therapy against hepatocellular carcinoma. Hepatol Res. 2015;45:1023–33.
Huang XY, Zhang PF, Wei CY, Peng R, Lu JC, Gao C, et al. Circular RNA circMET drives immunosuppression and anti-PD1 therapy resistance in hepatocellular carcinoma via the miR-30-5p/snail/DPP4 axis. Mol Cancer. 2020;19:92.
Heward JA, Lindsay MA. Long non-coding RNAs in the regulation of the immune response. Trends Immunol. 2014;35:408–19.
Wang CJ, Zhu CC, Xu J, Wang M, Zhao WY, Liu Q, et al. The lncRNA UCA1 promotes proliferation, migration, immune escape and inhibits apoptosis in gastric cancer by sponging anti-tumor miRNAs. Mol Cancer. 2019;18:115.
Zhang M, Zheng Y, Sun Y, Li S, Chen L, Jin X, et al. Knockdown of NEAT1 induces tolerogenic phenotype in dendritic cells by inhibiting activation of NLRP3 inflammasome. Theranostics. 2019;9:3425–42.
This work was supported by the Hainan Province Science and Technology Special Fund [ZDYF2021SHFZ051]; Hainan Provincial Natural Science Foundation of China [820MS053]; Bioinformatics for Mjaor Diseases Science Innovation Group of Hainan Medical University, National Natural Science Foundation of China [82173106, 82130115,31871338, 31970646, 61873075, 32060152, 32070673, 32170676]; Shanghai Pujiang talent program [20PJ1413000], Major Science and Technology Program of Hainan Province [ZDKJ202003]; HMU MarshalI Initiative Funding [HMUMIF-21024]; Marshal Initiative Funding of Hainan Medical University [JBGS202103]; project supported by Hainan Province Clinical Medical Center [QWYH202175]; National Key R&D Program of China [2018YFC2000100]; Natural Science Foundation for Distinguished Young Scholars of Heilongjiang Province [JQ2019C004]; Heilongjiang Touyan Innovation Team Program; Hainan Provincial Key Laboratory of Carcinogenesis and Intervention [JCKF2021003].
Ethics approval and consent to participate
Patient data we used were acquired by publicly available datasets that were collected with patients’ informed consent.
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.
Primers, probes and siRNA/shRNA used in this study. Table S2. Antibodies for western blot and IHC and in vivo used in this study. Table S3. Differentially expressed lncRNAs in HCC.
Genomic and transcriptome alterations of lncRNAs in cancer. Fig. S2. Expression of LINC01132 in HCC cell lines. Fig. S3. LINC01132 increases cancer cell growth, proliferation, invasion and metastasis in vitro. Fig. S4. PDX models shLINC01132. Fig. S5. Expression and clinical association of NRF1 and KDM5B. Fig. S6. Integrative Genomics Viewer of KDM5B binding around the TSS of DPP4 in HCC cell lines. Fig. S7. Interactions of DPP4 in protein-protein interaction networks.
About this article
Cite this article
Zhang, J., Pan, T., Zhou, W. et al. Long noncoding RNA LINC01132 enhances immunosuppression and therapy resistance via NRF1/DPP4 axis in hepatocellular carcinoma. J Exp Clin Cancer Res 41, 270 (2022). https://doi.org/10.1186/s13046-022-02478-z