MIR22HG acts as a tumor suppressor via TGFβ/SMAD signaling and facilitates immunotherapy in colorectal cancer
Molecular Cancer volume 19, Article number: 51 (2020)
Long noncoding RNAs (lncRNAs) are emerging as critical regulatory elements and play fundamental roles in the biology of various cancers. However, we are still lack of knowledge about their expression patterns and functions in human colorectal cancer (CRC).
Differentially expressed lncRNAs in CRC were identified by bioinformatics screen and the level of MIR22HG in CRC and control tissues were determined by qRT-PCR. Cell viability and migration capacities were examined by MTT and transwell assay. Mouse model was used to examine the function and rational immunotherapy of MIR22HG in vivo.
We systematically investigated the expression pattern of lncRNAs and revealed MIR22HG acts as a tumor suppressor in CRC. The expression of MIR22HG was significantly decreased in CRC, which was mainly driven by copy number deletion. Reduced expression of MIR22HG was significantly associated with poor overall survival. Silencing of MIR22HG promoted cell survival, proliferation and tumor metastasis in vitro and in vivo. Mechanistically, MIR22HG exerts its tumor suppressive activity by competitively interacting with SMAD2 and modulating the activity of TGFβ pathway. Decreased MIR22HG promoted the epithelial-mesenchymal transition in CRC. Importantly, we found that MIR22HG expression is significantly correlated with CD8A and overexpression of MIR22HG triggers T cell infiltration, enhancing the clinical benefits of immunotherapy.
MIR22HG acts as a tumor suppressor in CRC. Our data provide mechanistic insights into the regulation of MIR22HG in TGFβ pathway and facilitates immunotherapy in cancer.
Colorectal cancer (CRC) is one of the most common malignant tumors with increasing incidence worldwide. Despite therapeutic advances over decades, the majority of CRC patients are diagnosed at an advanced stage with poor prognosis . Thus, it is still a challenge to find novel biomarkers for CRC as well as identifying effective therapeutic targets.
With the development of next-generation sequencing technology, it is realized that the most part of human genome is transcribed to noncoding RNAs . Of these newly discovered RNA regulatory elements, long noncoding RNAs (lncRNAs) have been demonstrated to play critical roles in diverse cellular processes [3,4,5]. LncRNAs can serve as signal mediators, scaffold or molecular decoys to function in epigenetic, transcriptional and post-transcriptional regulation. For example, HOTAIR is found to be associated with poor prognosis in CRC by regulating polycomb-dependent chromatin modification . CCAT1, CCAT2 and CCAT1-L have been shown to function as oncogenes by up-regulating c-Myc, promoting tumor cell proliferation and migration [7,8,9]. In addition, several lncRNAs play their roles through competing binding sites of microRNAs, such as HULC  and CRNDE . LncRNA SNHG1 can regulate CRC cell growth through interactions with EZH2 . Linc00659 acts as an oncogene in regulating cell growth in CRC by impairing cell cycle progression . Despite the growing list of annotated lncRNAs in CRC, the experimentally verified is small overall. The genome-wide expression pattern and function of lncRNAs in CRC still need to be analyzed.
Moreover, immunotherapy is rapidly becoming a novel therapy for many cancers, including CRC . However, some patients had been generally resistant to immunotherapy. To improve the response rates, it is critical to identify more specific biomarkers and immune checkpoint inhibitors. LncRNAs are emerging as critical regulators of gene expression in the immune system and play critical roles in the development and progression of cancer [15, 16]. However, only a few immune-related lncRNAs have been conducted in cancer so far [17,18,19]. Therefore, further studies on lncRNAs in immune regulation will be essential to identify the immunotherapy targets in cancer.
Here, we have discovered a novel lncRNA-MIR22HG, by analysis of genome-wide transcriptome across > 1500 CRC samples. We illustrated that MIR22HG is significantly downregulated in CRC, which is mainly driven by copy number deletion. Silencing of MIR22HG is associated with poor prognosis in CRC patients. Furthermore, we found that knockdown with shRNAs of MIR22HG significantly promote cell survival, proliferation and tumor metastasis in vitro and in vivo. MIR22HG functions through regulating the TGFβ pathway via competing interacting with SMAD2. Decrease of MIR22HG promotes the epithelial-mesenchymal transition in CRC. Finally, we found that overexpression of MIR22HG can enhance the response of immunotherapy by increasing the CD8 T cell infiltration, suggesting a rational combinational therapy strategy in CRC.
We collected 163 CRC patients and their clinical information in China. In total, 79 paired human CRC tissues and normal controls were collected from Affiliated Hospital of Jiangnan University. Another 84 paired cancer and normal tissues were collected at Fudan University Shanghai Cancer Center. There were 21 CRC patients with liver metastasis among these patients. All patient materials were obtained with informed consent and this study was carried out under the permission of the Clinical Research Ethics Committees.
The CRC cell lines LoVo and HCT116 were purchased from American Type Culture Collection (ATCC) and cultured following the recommended instructions. These cells were characterized by Genewiz Inc. (China) using short tandem repeat (STR) markers and were confirmed to be Mycoplasma-free (last tested in year 2017).
Genome-wide expression of lincRNAs in CRC
Genome-wide gene expression across 478 colon adenocarcinoma (COAD) patients and 41 normal controls, and 166 rectum adenocarcinoma (READ) patients and 10 normal controls were downloaded from TCGA project . Gene expression was measured by Fragments Per Kilobase of transcript per Million mapped reads (FPKM). Next, we download the gene and lncRNA annotation from GENCODE project . Here, we focused on the long intergenic noncoding RNAs (lincRNAs) and obtained 7520 lincRNAs for further analysis.
In addition, we searched Gene Expression Omnibus (GEO) by using ‘colorectal cancer’ as keywords and obtained nine datasets across 1311 samples (Table 1). The expression of MIR22HG was extracted in our analysis. The difference of MIR22HG expression in cancer and normal was compared by Wilcox rank sum test.
Identification of differentially expressed lincRNAs
The expression of lincRNAs was first log transformed. To identify the lincRNAs that are perturbed in colorectal cancer, we used Wilcoxon rank sum test to identify the differentially expressed lincRNAs in TCGA COAD and READ datasets. The p-values were adjusted by Benjamini-Hochberg procedure. The lincRNAs with adjusted p-values (false discovery rate, FDR) less than 0.005 and fold change (FC) greater than 3-folds were considered as differentially expressed in colorectal cancer.
Literature curation of lincRNAs
To investigate whether the differentially expressed lincRNAs are involved in cancer, we used lincRNA symbols and ‘cancer’ as keywords to query the PubMed. This process was performed by the R package ‘RISmed’ (https://rdrr.io/cran/RISmed/). The up-regulated and down-regulated lincRNAs are ranked based on the number of literature, separately.
Gene Set Enrichment Analysis
To identify the pathways potentially regulated by MIR22HG, we first calculated the expression correlation coefficient for each gene with MIR22HG. Here, we analyzed the gene expression from COAD and READ of TCGA, separately . We next ranked all the genes based on the correlation coefficient (R) and these genes are subjected to Gene Set Enrichment Analysis (GSEA) analysis [22, 23]. The pathways with adjusted p-values< 0.01 were identified. We further calculated the proportion of leading edges genes in each enriched pathway. The cancer hallmark-related signaling pathways from MsigDB were considered .
Immune-related scores estimation and immune cell deconvolution
We calculated three immune response-related scores estimated from gene expression of CRC patients. The immune scores were estimated based on ESTIMATE  and MHC scores were estimated from expression of “core” MHC-I set . A simple and quantitative measure of immune cytolytic activity (‘CYT’) were calculated based on expression levels of granzyme A (GZMA) and perforin (PRF1) . In addition, TIMER was used to estimate the abundances of member cell types in a mixed cell population . The mutation burden of each CRC patient was calculated as non-silent mutations per Mb.
All results are presented as the mean values+SEM. Wilcox rank sum test, student t test, and the Mann-Whitney U test were used to compared the differences among different groups. The Kaplan-Meier method and log-rank test were used to determine the differences in survival rates between MIR22HG low and high expression groups. All statistical analyses were carried out using R 3.5.1. P values < 0.05 were considered statistically significant.
Supporting materials and methods
For details regarding the quantitative RT-PCR, vector construction and siRNA, overexpression and knockdown, cell proliferation, colony formation, cell migration and invasion assays, tumour formation and metastasis assays, RNA pull-down and mass spectrometry, RNA immunoprecipitation, western blot, immunohistochemistry and other related procedures, TGFβ signaling pathway inhibition, refer to the Supporting materials (Additional file 1: Supporting Materials and Methods).
Integrative analysis reveals MIR22HG as a tumor suppressor in CRC
To identify lncRNAs that paly critical roles in CRC, we first analyzed the genome-wide expression data in COAD and READ of TCGA project (Fig. 1a and Additional file 2: Figure 1a). Particularly, we focused on the 7520 long intergenic noncoding RNAs (lincRNAs). Based on the fold changes > 3-fold and false discovery rate (FDR) < 0.005, we identified 232/26 up-regulated and 274/184 down-regulated lincRNAs in COAD and READ, separately (Fig. 1a). There are 23 up-regulated and 126 down-regulated lincRNAs were identified in both COAD and READ. We found that the expression of these common lincRNAs can effectively distinguish the cancer patients from the normal controls in COAD and READ (Fig. 1b).
Next, we investigated to what extent these lincRNAs were studied in cancer. We found that there are higher number of literature for up-regulated PVT1, CCAT1 and CRNDE (Fig. 1c), suggesting that they have been mostly investigated in cancer. Evidence has demonstrated that PVT1 and CCAT1-encoded lncRNAs have oncogenic functions [29, 30], which is consistent with their higher expression in cancer (Additional file 2: Figure 1b). In the contrast, there is less number of literatures for the down-regulated lincRNAs (Fig. 1d). We found that MIR22HG is the mostly investigated one in cancer and has been show to function as a tumor suppressor in hepatocellular carcinoma [31, 32], lung cancer  and thyroid cancer . However, we are still lack of knowledge of its roles in CRC. We found that this lincRNA is not only low expressed in CRC (Fig. 1e, FDR = 1.23E-19 and 0.003), but also shows lower expression in liver and lung cancer (Additional file 2: Figure 1c).
Although these results suggest that MIR22HG might play a tumor-suppressor role in cancer, little is known about the regulation of MIR22HG. We thus next investigated the copy number variation (CNV) of MIR22HG across 33 cancer types . We found that the genomic region of MIR22HG shows prevalent CNV loss across cancer types (Fig. 1f and Additional file 2: Figure 1d). Together, these observations suggest that CNV loss might be the main driver of its lower expression in cancer.
Validation of MIR22HG in independent cohorts of CRC
To further verify the MIR22HG expression pattern in CRC, we next surveyed the public CRC datasets across 1311 samples. We found that MIR22HG expression levels were significantly lower in CRC as compared with colon tissues in four independent datasets (Fig. 2a and Additional file 2: Figure 2a-d, p-values< 0.01). Interestingly, we found that metastatic cancer patients have even lower levels of MIR22HG as compared with primary cancer (Fig. 2b-c). These results indicated that MIR22HG expression is generally lower in CRC and may also be an indicator of metastasis of CRC. Moreover, we investigated the correlation of MIR22HG expression with the patients’ survival. Kaplan-Meier survival curve and log-rank test suggest that lower expression of MIR22HG was significantly associated with poor overall survival and disease-free survival of patients (Fig. 2d-e and Additional file 2: Figure 2e, p = 0.0008 and 0.0009).
In addition, we collected 163 CRC patients and used qRT-PCR to investigate the expression of MIR22HG. We found that MIR22HG was significantly low expressed in cancer as compared with normal controls (Fig. 2f, p < 0.001). Patients with larger tumor size were significantly enriched in the MIR2HG low expression group (Additional file 3: Table S1). Particularly, liver metastasis patients were also with lower expression of MIR22HG as compared with their own primary cancer (Fig. 2g, p < 0.001). Finally, we explored the effects of MIR22HG on prognosis of CRC patients and revealed that lower expression of MIR22HG is significantly associated with poor survival of CRC (Fig. 2h, p = 0.007). All these data demonstrated that down-regulation of MIR22HG in CRC may play a tumor suppressor role and could serve as a novel prognostic marker.
Silencing MIR22HG promotes CRC growth and metastasis in vitro
To investigate the possible function of MIR22HG in CRC pathogenesis, we firstly explored the cellular localization and coding potential of MIR22HG. The MIR22HG expression in both nuclear and cytoplasmic fractions from LoVo and HCT116 were measured by qRT-PCR. We observed that the expression of MIR22HG were in both cytoplasm and nuclear (Additional file 2: Figure 3a-b), suggesting that it may play a regulatory function in both localizations. In addition, we found that MIR22HG didn’t express to protein but play its roles in RNA (Additional file 2: Figure 3c).
We next used siRNA technology to knock down the expression of MIR22HG in HCT116 cell line and measured the colony formation, cell proliferation and migration. We found that the colony formation was increased by two-folds as compared with control siRNA (Fig. 3a, p < 0.001). Moreover, we further confirmed that cell proliferation was affected by MIR22HG knockdown and there are approximate 70% increase at 96 to 120 h (Fig. 3b, p < 0.05). Transwell assays were employed to explore whether cell migration were affected after silencing of MIR22HG. We found that MIR22HG knockdown significantly increased cell migration as compared with control (Fig. 3c, p < 0.001).
Moreover, we over-expressed MIR22HG in LoVo and further measured the colony formation, cell proliferation and migration. We found that overexpression of MIR22HG will significantly decrease the number of colony formation (Fig. 3d, p < 0.001). Overexpression of MIR22HG also significantly decreased the cell proliferation capability (Fig. 3e) and cell migration (Fig. 3f, p < 0.001). Taken together, these results suggested that MIR22HG plays tumor suppressor role by affecting colony formation, proliferation and cell migration.
Overexpression of MIR22HG suppresses cell proliferation and metastasis in vivo
Furthermore, we examined the effect of MIR22HG overexpression on tumor growth in vivo using mouse subcutaneous xenograft models. We found that tumor xenografts derived from MIR22HG-overexpressing cells exhibited smaller volumes and lower weights than those from empty vector-transduced cells (Fig. 4a, p < 0.001). Additionally, we assessed the effects of MIR22HG on colorectal cancer metastasis. Three metastasis models were considered, including an intestine metastasis, a lung metastasis mouse model and an orthotopic hepatic metastasis mouse model. We found that the number of metastatic nodules are significantly reduced in MIR22HG-overexpressing groups (Fig. 4b-e, p-values< 0.001). In addition, multiple metastases of different sizes were found in the control mouse, whereas metastatic foci were sparse and small in MIR22HG-overexpressing group. Collectively, these data indicated that overexpression of MIR22HG repressed CRC cell invasion and metastasis.
MIR22HG perturbs the TGFβ/SMADs signaling pathway in CRC
To elucidate the potential molecular mechanisms of MIR22HG in CRC, we performed the Gene Set Enrichment Analysis (GSEA). We found that the co-expressed genes of MIR22HG are significantly enriched in number of cancer-related signaling pathways (Fig. 5a), including IL6-JAK-STAT3, TNFA signaling and TGFβ signaling pathway. TGFβ signaling plays a pivotal role in diverse cellular processes, such as cell proliferation and migration . We thus investigated the function of MIR22HG in the context of TGFβ signaling pathway. First, we calculated the expression correlation coefficient between MIR22HG and genes in TGFβ pathway. We found that the expression of MIR22HG is significantly correlated with SMAD genes, such as SMAD2 (Fig. 5b), SMAD3 and SMAD4 in both COAD and READ (Additional file 2: Figure 4a). Moreover, the expression of SMAD2 is correlated with patient survival (Additional file 2: Figure 4b). Defects in SMAD signaling can result in TGFβ resistance, causing dysregulation of cell growth . These results suggest that MIR22HG might perform its function via TGFβ/SMAD signaling pathway.
Next, we performed biotin-labeled RNA pull down accompanied by mass spectrometric assays to identify MIR22HG-interacting proteins. The results showed that SMAD2 protein bands at 58.8 kDa in the MIR22HG pull-downed sample (Additional file 2: Figure 4c). We also observed that sense but not antisense of MIR22HG, specifically interact with SMAD2 (Additional file 2: Figure 4c). The interactions were further confirmed through RNA immunoprecipitation (RIP) assays (Fig. 5c). Furthermore, we found that the fragment (674–805 nucleotides) of MIR22HG was responsible for the interaction with SMAD2 (Additional file 2: Figure 4d). MIR22HG knockdown decreased the half-life of SMAD2 protein, whereas ectopic MIR22HG expression increased the half-life of SMAD2 protein in CRC cells (Fig. 5d). Together, these results suggest that MIR22HG can bind to SMAD2 in CRC.
We next characterized the molecular consequences of the associations among MIR22HG and SMAD2. Several lines of evidence have shown that disrupt TGFβ signaling by preventing heteromeric interactions between Smads . Moreover, we found that MIR22HG exerts its biological function by inhibiting the growth and migration of colorectal cancer cells through SMAD2 and TGFβ signaling pathway. Interfering the expression of SMAD2 or the activity of TGFβ signaling pathway were likely to affect the MIR22HG effects on cell growth (Additional file 2: Figure 5a) and cell migration (Additional file 2: Figure 5b). Thus, we hypothesized that MIR22HG might play its function by perturbing the interaction between SMAD2 and SMAD4 by competitively interacting with SMAD2. To test this hypothesis, the interaction between SMAD2 and SMAD4 was examined by co-immunoprecipitation assays, in the presence and absence of MIR22HG. MIR22HG knockdown could enhance the interaction between SMAD2 and SMAD4 (Fig. 5e, left). When MIR22HG was overexpressed, there was less SMAD2 that interacted with SMAD4, as well as less SMAD4 that bound to SMAD2 (Fig. 5e, right), indicating that MIR22HG exhibits a negative regulatory effect on the interaction between SMAD2 and SMAD4. Collectively, these results demonstrated that MIR22HG exerted its tumor suppressive function through competitively binding to SMAD2, thereby preventing the interaction between SMAD2 and SMAD4 in CRC cells.
Abrogation of MIR22HG promotes epithelial–mesenchymal transition in CRC
TGF-beta signaling has been shown to play an important role in epithelial to mesenchymal transition (EMT) . We thus investigated whether MIR22HG plays roles in EMT. Based on GSEA analysis, we found that the co-expressed genes with MIR22HG were significantly enriched in EMT pathway in both COAD and READ (Fig. 6a, FDRs < 0.001). Consistent with these results, we found that the overexpression of MIR22HG led to an epithelial cells property while knockdown MIR22HG acquired mesenchymal properties (Fig. 6b-c). Meanwhile, qRT-PCR and western blot assays showed that the mRNA and protein levels of the epithelial markers (E-cadherin, ZO-1 and Occludin) were markedly decreased and accompanied by an increase in the mesenchymal marker N-cadherin, vimentin and fibronrctin when knockdown MIR22HG (Fig. 6d-e). Moreover, we found that the MIR22HG effects were potentially dependent on the expression of SMAD2 and TGFβ signaling pathway activity. Interfering the expression of SMAD2 and TGFβ signaling pathway inhibition can to a certain extent reverse the EMT gene and protein expression (Additional file 2: Figure 6). These data suggest that the silencing of MIR22HG promotes tumorigenicity in CRC cells by promoting EMT process.
Given that SNAI1 is a transcription factor important for the early step of EMT and SMAD pathway [40, 41], we investigated whether expression of SNAI1 is correlated with SMAD2/SMAD4. We found that SMAD2/SMAD4 can bind to the promoter of SNAI1 (Additional file 2: Figure 7a) and the expressions of SMAD4 are significantly correlated with SNAI1 in CRC (Additional file 2: Figure 7b). Moreover, the expression of SNAI1 is negatively correlated with CDH1 (E-cadherin) and positively correlated with CDH2 (N-cadherin) (Additional file 2: Figure 7c-d). Together, these results suggest that MIR22HG-SMAD2/4-SNAI1 axis play critical roles in CRC by promoting the EMT.
Rational combination immunotherapy of MIR22HG in CRC
Patients with immunogenicity seems to be a promising strategy to develop new drugs that target the antitumor immune response in CRC . However, some patients exhibit primary or acquired resistance. It is critical to identify the additional targets for improving the sensitivity of immunotherapy. Increasing studies have demonstrated that TGFβ pathway affect the expression of PD-L1 and CD8A [43, 44]. We found that the CRC patients with MIR22HG high expression exhibit significantly higher mutation burden, MHC scores, immune scores and CYT scores (Fig. 7a). Moreover, we curated the drug treatment information from TCGA and found that patients with higher MIR22HG expression were likely to response for therapy (Additional file 2: Figure 8). These observations suggest that patients with high MIR22HG expression may be more sensitive to immunotherapy. To further investigate the function of MIR22HG in immune, syngeneic immunocompetent mouse model C57BL/6 mice was used (Fig. 7b). Ten days after CRC implantation, the mice were treated MIR22HG and PD-L1 blockade or IgG. We found that administration of MIR22HG and PD-L1 blockade can inhibit tumor growth (Fig. 7c) and reduce the tumor weight (Fig. 7d, p < 0.05). In addition, the mice treated with the combination of MIR22HG and aPD-L1 was with a much smaller tumor volume than the other mice (Fig. 7e). Bioluminescent imaging revealed that overexpression of MIR22HG efficiently enhanced the sensitivity of aPD-L1 treatment (Fig. 7f). The mice receiving a combined treatment demonstrated a slower flux (Fig. 7f) and had a significantly prolonged lifespan (Fig. 7g). Overall, these data demonstrated that MIR22HG serves as a potential therapeutic target to overcome aPD-L1 resistance, enhancing the benefits of aPD-L1 therapy.
Overexpression of MIR22HG increases the infiltration of CD8 T cells
We next investigate the potential mechanism for the combination therapy of MIR22HG and aPD-L1. Evidence has indicated that the immune cell infiltration level is correlated with immune therapy response. We thus investigated the expression of MIR22HG with the CD8 T cell infiltration. We found that the expression of MIR22HG is significantly correlated with T cell infiltration in COAD and READ (Additional file 2: Figure 9, p-values < 0.05). In addition, the expression of MIR22HG is significantly correlated with the T cell marker gene-CD8A (Fig. 8a). Flow cytometry analysis also shows that overexpressing MIR22HG significantly increases cell-surface expression of PD-L1 and CD8A in vitro and in vivo (Fig. 8b). Moreover, we found that the effects of MIR22HG on CD8A and PD-L1 can be affected by perturbation of SMAD2 expression and TGFβ signaling pathway activity (Additional file 2: Figure 10). These observations suggest that overexpression of MIR22HG is likely to increase the infiltration of CD8 T cells.
To further validate these results, we explored the expression of Cd8a in aPD-L1 treated mice and CRC patient. We found that aPD-L1 treated mice showed increased level of Cd8a (Fig. 8c) and the adjacent normal tissue of CRC patient also showed higher expression of CD8A and SMAD2 (Fig. 8c). These data demonstrated that increase of CD8A is likely to enhance aPD-L1 therapy. We next investigated the expression of Cd8a in mice treated with MIR22HG and aPD-L1. We found that the mice with combined treatment showed increased level of Cd8a than the other mice (Fig. 8d, p < 0.001). Taken together, these results demonstrate that overexpression of MIR22HG will increase the infiltration of CD8 T cells in CRC, further enhancing the sensitivity to immune therapy.
With the advancement of high-throughput sequencing technology, increasing lncRNAs have been identified in cancer. LncRNAs have been identified to play critical regulatory function in diverse processes including cancer initiation and progression. In this study, we analyzed the genome-wide lncRNAs expression across > 1500 samples and found the consistent results showing low expression of MIR22HG in CRC. Functional analyses of MIR22HG suggest that silencing of MIR22HG promotes the interactions between SMAD2 and SMAD4 of TGFβ pathway. The SMAD complex binds to the promoter of SNAI1 and further promotes the EMT process (Fig. 8e). Increasing studies have demonstrated that TGFβ pathway can remodel tumor immune environment. Importantly, the expression of MIR22HG is significantly correlated with CD8 T cell infiltration and combination of MIR22HG with aPD-L1 will enhance the sensitivity to immune therapy (Fig. 8e). These results suggest that MIR22HG acts as a tumor suppressor and facilitates immunotherapy in CRC.
By analysis of RNA sequencing data, we found that MIR22HG is low expressed in CRC. To explore the potential mechanism of its low expression, we investigated the CNV alterations of MIR22HG and found that silence of MIR2HG might be driven by CNV deletion. However, some other regulatory issues (such as transcription factors or miRNAs) still need to be addressed for investigate the detail mechanism. In addition, the mechanisms of action of lncRNAs are diverse and include RNA-RNA, RNA-DNA, and RNA-protein interactions . Evidence has shown that MIR22HG is a target of miR-141-3p in endometrial carcinoma . MIR22HG was shown to competing with miRNAs with H19 and play important regulatory roles in breast cancer progression . MIR22HG can also competitively bind to human antigen R (HuR), resulting in weakened expression of HuR-stabilized oncogenes, such as β-catenin in hepatocellular carcinoma cells . MIR22HG was found to inhibit growth, migration and invasion through regulating the miR-10a-5p/NCOR2 axis in hepatocellular carcinoma cells . Silencing of MIR22HG can also triggers cell survival/death signaling via oncogenes YBX1, MET, and p21 in lung cancer . These results suggest the same lncRNA can play diverse roles by regulating different signaling pathways in different tumor context.
It is increasingly clear that there are widespread changes in the expression of lncRNAs during the activation of immune response . A number of lncRNAs have been found to play important roles in immunotherapy. NEAT1 knockdown can induce immune tolerance in vivo and may be a promising therapeutic target in the treatment of immune-related diseases . UCA1 overexpression had been shown to protect PDL1 expression from the repression of miRNAs and contributed to the gastric cancer cells immune escape . Oncogenic LINK-A can downregulate cancer cell antigen presentation and intrinsic tumor suppression in human triple-negative breast cancer . However, limited lncRNAs are found to play immunology function in CRC. We found MIR22HG overexpression triggers CD8 T cell infiltration, suggesting that MIR22HG might be a potential target for immunotherapy in CRC.
In conclusion, our findings broaden understanding of the function of lncRNAs in CRC. Moreover, our results demonstrate that tumor suppressor MIR22HG may act as a clinical biomarker and overexpression of MIR22HG may be a novel synergistic strategy for enhancing immunotherapy sensitivity, thereby enhancing its clinical benefits for CRC patients.
Availability of data and materials
The gene expression profiles and clinical data can be found at the GDC portal (https://portal.gdc.cancer.gov/) and GEO (https://www.ncbi.nlm.nih.gov/geo/). The copy number data were downloaded from Broad GDAC Firehose (https://gdac.broadinstitute.org/). Software and resources used for the analyses are described in each method section.
Colon Cancer Associated Transcript 1
Copy number variation
Colorectal Neoplasia Differentially Expressed
False discovery rate
Gene Set Enrichment Analysis
Long non-coding RNA
MIR22 Host Gene
Programmed death-ligand 1
Short hairpin RNA or small hairpin RNA
SMAD Family Member 2
SMAD Family Member 4
Snail Family Transcriptional Repressor 1
Yonemura K, Kajiwara Y, Ao T, Mochizuki S, Shinto E, Okamoto K, Hase K, Ueno H. Prognostic value of poorly differentiated clusters in liver metastatic lesions of colorectal carcinoma. Am J Surg Pathol. 2019.
Carninci P, Kasukawa T, Katayama S, Gough J, Frith MC, Maeda N, Oyama R, Ravasi T, Lenhard B, Wells C, et al. The transcriptional landscape of the mammalian genome. Science. 2005;309:1559–63.
Mercer TR, Dinger ME, Mattick JS. Long non-coding RNAs: insights into functions. Nat Rev Genet. 2009;10:155–9.
Li Y, Li L, Wang Z, Pan T, Sahni N, Jin X, Wang G, Li J, Zheng X, Zhang Y, et al. LncMAP: Pan-cancer atlas of long noncoding RNA-mediated transcriptional network perturbations. Nucleic Acids Res. 2018;46:1113–23.
Li Y, Zhang J, Huo C, Ding N, Li J, Xiao J, Lin X, Cai B, Zhang Y, Xu J. Dynamic organization of lncRNA and circular RNA regulators collectively controlled cardiac differentiation in humans. EBioMedicine. 2017;24:137–46.
Kogo R, Shimamura T, Mimori K, Kawahara K, Imoto S, Sudo T, Tanaka F, Shibata K, Suzuki A, Komune S, et al. Long noncoding RNA HOTAIR regulates polycomb-dependent chromatin modification and is associated with poor prognosis in colorectal cancers. Cancer Res. 2011;71:6320–6.
Nissan A, Stojadinovic A, Mitrani-Rosenbaum S, Halle D, Grinbaum R, Roistacher M, Bochem A, Dayanc BE, Ritter G, Gomceli I, et al. Colon cancer associated transcript-1: a novel RNA expressed in malignant and pre-malignant human tissues. Int J Cancer. 2012;130:1598–606.
Xiang JF, Yin QF, Chen T, Zhang Y, Zhang XO, Wu Z, Zhang S, Wang HB, Ge J, Lu X, et al. Human colorectal cancer-specific CCAT1-L lncRNA regulates long-range chromatin interactions at the MYC locus. Cell Res. 2014;24:513–31.
Ling H, Spizzo R, Atlasi Y, Nicoloso M, Shimizu M, Redis RS, Nishida N, Gafa R, Song J, Guo Z, et al. CCAT2, a novel noncoding RNA mapping to 8q24, underlies metastatic progression and chromosomal instability in colon cancer. Genome Res. 2013;23:1446–61.
Panzitt K, Tschernatsch MM, Guelly C, Moustafa T, Stradner M, Strohmaier HM, Buck CR, Denk H, Schroeder R, Trauner M, Zatloukal K. Characterization of HULC, a novel gene with striking up-regulation in hepatocellular carcinoma, as noncoding RNA. Gastroenterology. 2007;132:330–42.
Ellis BC, Molloy PL, Graham LD. CRNDE: a long non-coding RNA involved in CanceR, neurobiology, and DEvelopment. Front Genet. 2012;3:270.
Xu M, Chen X, Lin K, Zeng K, Liu X, Pan B, Xu X, Xu T, Hu X, Sun L, et al. The long noncoding RNA SNHG1 regulates colorectal cancer cell growth through interactions with EZH2 and miR-154-5p. Mol Cancer. 2018;17:141.
Tsai KW, Lo YH, Liu H, Yeh CY, Chen YZ, Hsu CW, Chen WS, Wang JH. Linc00659, a long noncoding RNA, acts as novel oncogene in regulating cancer cell growth in colorectal cancer. Mol Cancer. 2018;17:72.
Le DT, Hubbard-Lucey VM, Morse MA, Heery CR, Dwyer A, Marsilje TH, Brodsky AN, Chan E, Deming DA, Diaz LA Jr, et al. A blueprint to advance colorectal Cancer immunotherapies. Cancer Immunol Res. 2017;5:942–9.
Chen YG, Satpathy AT, Chang HY. Gene regulation in the immune system by long noncoding RNAs. Nat Immunol. 2017;18:962–72.
Lin X, Jiang T, Bai J, Li J, Wang T, Xiao J, Tian Y, Jin X, Shao T, Xu J, et al. Characterization of Transcriptome transition associates long noncoding RNAs with Glioma progression. Mol Ther Nucleic Acids. 2018;13:620–32.
Huang D, Chen J, Yang L, Ouyang Q, Li J, Lao L, Zhao J, Liu J, Lu Y, Xing Y, et al. NKILA lncRNA promotes tumor immune evasion by sensitizing T cells to activation-induced cell death. Nat Immunol. 2018;19:1112–25.
Hu Q, Ye Y, Chan LC, Li Y, Liang K, Lin A, Egranov SD, Zhang Y, Xia W, Gong J, et al. Oncogenic lncRNA downregulates cancer cell antigen presentation and intrinsic tumor suppression. Nat Immunol. 2019;20:835–51.
Lin H, Jiang M, Liu L, Yang Z, Ma Z, Liu S, Ma Y, Zhang L, Cao X. The long noncoding RNA Lnczc3h7a promotes a TRIM25-mediated RIG-I antiviral innate immune response. Nat Immunol. 2019;20:812–23.
Cancer Genome Atlas N: Comprehensive molecular characterization of human colon and rectal cancer. Nature 2012, 487:330–337.
Derrien T, Johnson R, Bussotti G, Tanzer A, Djebali S, Tilgner H, Guernec G, Martin D, Merkel A, Knowles DG, et al. The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression. Genome Res. 2012;22:1775–89.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. 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.
Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, Puigserver P, Carlsson E, Ridderstrale M, Laurila E, et al. PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003;34:267–73.
Liberzon A, Subramanian A, Pinchback R, Thorvaldsdottir H, Tamayo P, Mesirov JP. Molecular signatures database (MSigDB) 3.0. Bioinformatics. 2011;27:1739–40.
Yoshihara K, Shahmoradgoli M, Martinez E, Vegesna R, Kim H, Torres-Garcia W, Trevino V, Shen H, Laird PW, Levine DA, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612.
Lauss M, Donia M, Harbst K, Andersen R, Mitra S, Rosengren F, Salim M, Vallon-Christersson J, Torngren T, Kvist A, et al. Mutational and putative neoantigen load predict clinical benefit of adoptive T cell therapy in melanoma. Nat Commun. 2017;8:1738.
Rooney MS, Shukla SA, Wu CJ, Getz G, Hacohen N. Molecular and genetic properties of tumors associated with local immune cytolytic activity. Cell. 2015;160:48–61.
McGranahan N, Furness AJ, Rosenthal R, Ramskov S, Lyngaa R, Saini SK, Jamal-Hanjani M, Wilson GA, Birkbak NJ, Hiley CT, et al. Clonal neoantigens elicit T cell immunoreactivity and sensitivity to immune checkpoint blockade. Science. 2016;351:1463–9.
Curtis C, Shah SP, Chin SF, Turashvili G, Rueda OM, Dunning MJ, Speed D, Lynch AG, Samarajiwa S, Yuan Y, et al. The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups. Nature. 2012;486:346–52.
Ozawa T, Matsuyama T, Toiyama Y, Takahashi N, Ishikawa T, Uetake H, Yamada Y, Kusunoki M, Calin G, Goel A. CCAT1 and CCAT2 long noncoding RNAs, located within the 8q.24.21 'gene desert', serve as important prognostic biomarkers in colorectal cancer. Ann Oncol. 2017;28:1882–8.
Wu Y, Zhou Y, Huan L, Xu L, Shen M, Huang S, Liang L. LncRNA MIR22HG inhibits growth, migration and invasion through regulating the miR-10a-5p/NCOR2 axis in hepatocellular carcinoma cells. Cancer Sci. 2019;110:973–84.
Zhang DY, Zou XJ, Cao CH, Zhang T, Lei L, Qi XL, Liu L, Wu DH. Identification and functional characterization of long non-coding RNA MIR22HG as a tumor suppressor for hepatocellular carcinoma. Theranostics. 2018;8:3751–65.
Su W, Feng S, Chen X, Yang X, Mao R, Guo C, Wang Z, Thomas DG, Lin J, Reddy RM, et al. Silencing of long noncoding RNA MIR22HG triggers cell survival/death signaling via oncogenes YBX1, MET, and p21 in lung Cancer. Cancer Res. 2018;78:3207–19.
Qin L, Luo JZ, Tang XL, Han CG. Identification of long noncoding RNA MIR22HG as a novel biomarker in thyroid Cancer. Pathol Oncol Res. 2019;25:703–10.
Hoadley KA, Yau C, Hinoue T, Wolf DM, Lazar AJ, Drill E, Shen R, Taylor AM, Cherniack AD, Thorsson V, et al. Cell-of-origin patterns dominate the molecular classification of 10,000 tumors from 33 types of Cancer. Cell. 2018;173:291–304 e296.
David CJ, Massague J. Contextual determinants of TGFbeta action in development, immunity and cancer. Nat Rev Mol Cell Biol. 2018;19:419–35.
Colak S, Ten Dijke P. Targeting TGF-beta signaling in Cancer. Trends Cancer. 2017;3:56–71.
Xu J, Attisano L. Mutations in the tumor suppressors Smad2 and Smad4 inactivate transforming growth factor beta signaling by targeting Smads to the ubiquitin-proteasome pathway. Proc Natl Acad Sci U S A. 2000;97:4820–5.
Xu J, Lamouille S, Derynck R. TGF-beta-induced epithelial to mesenchymal transition. Cell Res. 2009;19:156–72.
Miyazono K. Transforming growth factor-beta signaling in epithelial-mesenchymal transition and progression of cancer. Proc Jpn Acad Ser B Phys Biol Sci. 2009;85:314–23.
Kaufhold S, Bonavida B. Central role of Snail1 in the regulation of EMT and resistance in cancer: a target for therapeutic intervention. J Exp Clin Cancer Res. 2014;33:62.
Gang W, Wang JJ, Guan R, Yan S, Shi F, Zhang JY, Li ZM, Gao J, Fu XL. Strategy to targeting the immune resistance and novel therapy in colorectal cancer. Cancer Med. 2018;7:1578–603.
Park BV, Freeman ZT, Ghasemzadeh A, Chattergoon MA, Rutebemberwa A, Steigner J, Winter ME, Huynh TV, Sebald SM, Lee SJ, et al. TGFbeta1-mediated SMAD3 enhances PD-1 expression on antigen-specific T cells in Cancer. Cancer Discov. 2016;6:1366–81.
Zhao J, Feng M, Liu D, Liu H, Shi M, Zhang J, Qu J. Antagonism between HTRA3 and TGFbeta1 contributes to metastasis in non-small cell lung Cancer. Cancer Res. 2019;79:2853–64.
Li Y, Wang Z, Wang Y, Zhao Z, Zhang J, Lu J, Xu J, Li X. Identification and characterization of lncRNA mediated transcriptional dysregulation dictates lncRNA roles in glioblastoma. Oncotarget. 2016;7:45027–41.
Cui Z, An X, Li J, Liu Q, Liu W. LncRNA MIR22HG negatively regulates miR-141-3p to enhance DAPK1 expression and inhibits endometrial carcinoma cells proliferation. Biomed Pharmacother. 2018;104:223–8.
Wang P, Ning S, Zhang Y, Li R, Ye J, Zhao Z, Zhi H, Wang T, Guo Z, Li X. Identification of lncRNA-associated competing triplets reveals global patterns and prognostic markers for cancer. Nucleic Acids Res. 2015;43:3478–89.
Heward JA, Lindsay MA. Long non-coding RNAs in the regulation of the immune response. Trends Immunol. 2014;35:408–19.
Zhang M, Zheng Y, Sun Y, Li S, Chen L, Jin X, Hou X, Liu X, Chen Q, Li J, et al. Knockdown of NEAT1 induces tolerogenic phenotype in dendritic cells by inhibiting activation of NLRP3 inflammasome. Theranostics. 2019;9:3425–42.
Wang CJ, Zhu CC, Xu J, Wang M, Zhao WY, Liu Q, Zhao G, Zhang ZZ. The lncRNA UCA1 promotes proliferation, migration, immune escape and inhibits apoptosis in gastric cancer by sponging anti-tumor miRNAs. Mol Cancer. 2019;18:115.
The authors gratefully thank the TCGA Research Network for providing data for this work.
This work was supported by the National Natural Science Foundation of China (grants 31970646, 31571331, 31871338, 81602033 and 31601065), Natural Science Foundation for Distinguished Young Scholars of Heilongjiang Province (JQ2019C004); and University Nursing Program for Young Scholars with Creative Talents in Heilongjiang Province (UNPYSCT-2016189 and UNPYSCT-2017060). The funders played no roles in the study design, decision to publish, or preparation of the manuscript.
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.
About this article
Cite this article
Xu, J., Shao, T., Song, M. et al. MIR22HG acts as a tumor suppressor via TGFβ/SMAD signaling and facilitates immunotherapy in colorectal cancer. Mol Cancer 19, 51 (2020). https://doi.org/10.1186/s12943-020-01174-w