- Open Access
A comprehensive prognostic and immune analysis of enhancer RNA identifies IGFBP7-AS1 as a novel prognostic biomarker in Uterine Corpus Endometrial Carcinoma
Biological Procedures Online volume 24, Article number: 9 (2022)
Long non-coding RNAs (lncRNA) have been implicated in a hand of studies that supported an involvement and co-operation in Uterine Corpus Endometrial Carcinoma (UCEC). Enhancer RNAs (eRNA), a functional subtype of lncRNA, have a key role throughout the genome to guide protein production, thus potentially associated with diseases.
In this study, we mainly applied the Cancer Genome Atlas (TCGA) dataset to systematically discover crucial eRNAs involving UCEC. For the key eRNAs in UCEC, we employed RT-qPCR to compare eRNA expression levels in tumor tissues and paired normal adjacent tissues from UCEC patients for validation. Furthermore, the relationships between the key eRNAs and immune activities were measured from several aspects, including the analysis for tumor microenvironment, immune infiltration cells, immune check point genes, tumor mutation burden, and microsatellite instability, as well as m6A related genes. Finally, the key eRNAs were verified by a comprehensive pan-cancer analysis.
IGFBP7 Antisense RNA 1 (IGFBP7-AS1) was identified as the key eRNA for its expression patterns of low levels in tumor tissues and favorable prognostic value in UCEC correlated with its target gene IGFBP7. In RT-qPCR analysis, IGFBP7-AS1 and IGFBP7 had down-regulated expression in tumor tissues, which was consistent with previous analysis. Moreover, IGFBP7-AS1 was found closely related with immune response in relevant immune analyses. Besides, IGFBP7-AS1 and its target gene IGFBP7 correlated with a multi-omics pan-cancer analysis.
Finally, we suggested that IGFBP7-AS1 played a key role in impacting on clinical outcomes of UCEC patients for its possible influence on immune activity.
Uterine Corpus Endometrial Carcinoma (UCEC) is the second most prevalent type of malignancy among women worldwide and ranks as the third leading cause of gynecological malignancy death in females. It was estimated there would be 65,620 new cases (following breast cancer, 276,480) and 12,590 deaths (following breast cancer, 42,170; Ovary cancer, 13,940) in 2020 . The death rate for endometrial cancer has increased exceeding 100% during the past 20 years, rising by 8% since 2008 . In China, a total of 9,213 patients died from corpus uteri cancer, accounting for 2.20% of cancer-related deaths from malignancies in females . In recent years, the incidence of and mortality resulted from endometrial cancer has been growing quickly due to the increased cases of obesity and the subsequence hyperinsulinemia . The patients with prognosis for stage I disease who received surgery treatment attended good results, while the treatments for patients with advanced disease (stage III or IV) don’t work well, with overall survival (OS) rates in 5 years ranging from 47 to 69 and 15% to 17%, respectively. Moreover, therapeutic options are still limited except for first-line chemotherapy (Y. C. . Therefore, improvement in UCEC diagnosis is urgently needed.
Long non-coding RNAs (lncRNA) are defined as an important novel class of non-coding RNAs with lengths more than 200 nucleotides that are not translated into protein. lncRNAs are crucial players in regulating chromatin dynamics, gene expresson, growth, differentiation and development . The discoveries of huge numbers of lncRNA, including their wide-range of expression patterns in different types of cancer, their tumor-specificity, and stability in circulating body fluids (plasma and urine) provide a new foundation for developing cancer prognosis and therapies . Also, lncRNAs are also useful in survival outcome prediction of patients. While a lot of lncRNAs have been reported related with UCEC prognosis in many association studies [31, 51], there is little understanding of their function in causing diseases.
Enhancer RNA (eRNA) is a functional subtype of lncRNA transcribed by enhancers, which interacts with target promoter to enhance the transcription of target gene (Z. . eRNAs are widely recognized as bi-direction non-coding RNAs generated by enhancers prevalently bound by RNA polymerase II (Pol II) . There are many eRNAs that show remarkable expression difference between tumor and adjacent normal samples, which raise a potential to target eRNAs for cancer therapies (H. ,J. H. . The important function of eRNA in oncogene deregulation and cancer initiation has been established in many types of cancer (W. . However, current knowledge about the function and mechanistic roles of eRNA in UECE is limited.
In this study, we aimed to find out the key eRNA and its target gene in UCEC and we identified eRNA IGFBP7-AS1, with its target gene IGFBP7, was significantly correlated with survival outcomes of UCEC patients. Up-regulated levels of IGFBP7-AS1 correlated with increased level of IGFBP7, elevated expression of most immune-related genes and better immune response, which indicated that IGFBP7-AS1 was an immune-related eRNA in favor of clinical outcomes.
Data collection and processing
The RNA expression profile data in pan-cancers (Workflow Type: HTSeq-FPKM) were acquired from The Cancer Genome Atlas (TCGA) database, as well as relevant clinical and survival information . Relevant eRNA information was searched and summarized from relative eRNA references . Then the human gene annotation file was applied to convert the Ensemble transcript IDs of RNA expression profile data into gene symbols. The 1580 eRNAs expression profiles were thus extracted from the RNA expression profile data of UCEC by matching the eRNA gene symbols. The survival-associated eRNAs were identified by Kaplan–Meier tests, taking the P < 0.05 as standard cut-off values. The co-expression data was used to calculate the correlations between prognostic eRNAs and corresponding target genes on the basis of Spearman rank correlation coefficients (Spearman’s rank correlation coefficient R > 0.5, p < 0.05). The most significant ones which closely correlated with both OS survival and predicted target genes were considered as the key eRNAs in UCEC.
Gene Expression Profiling Interactive Analysis (GEPIA) is a well-known online tool for visual analysis from TCGA dataset (http://gepia.cancer-pku.cn/). Besides, it also includes normal and tumor samples of RNA sequencing data from Genotype-Tissue Expression projects . In this study, the comparisons of IGFBP7-AS1 and IGFBP7 expression in normal and tumor samples were visualized by GEPIA.
Reverse transcription-quantitative (RT-q) PCR for IGFBP7-AS1
Reverse transcription-quantitative (RT-q) PCR was applied to confirm the RNA-sequencing data of IGFBP7-AS1 and IGFBP7. The RNA was extracted from paired tumor and adjacent samples from 12 UCEC patients by TRIzol reagent (Thermo Fisher Scientific). Before reverse transcription to cDNA, 4 × gDNA wiper Mix (Vazyme R323-01), DEPC and total RNA (1 μg) were resuspended and reacted at 42 °C for 2 min to remove the residual genomic DNA from total RNA. Agilent Bioanalyzer 2100 (Agilent Technologies) with RNA 6000 Nano kit was applied to check the integrity of extracted RNA. PrimeScript RT reagent kit was used to synthesize the complementary RNA and SYBR Green PCR Kit (Thermo Fisher Scientific) was employed to conduct real-time quantitative analysis. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), an enzyme in glycolysis, which is widely distributed in cells of various tissues. GAPDH gene expression changes little in different states, so it was severed as the reference gene, and relative gene expression was estimated by the 2-ΔΔCt (ΔCt = Cttarget gene − Ctinternal control) method. The PCR primer sequences are presented in Table 1.
We performed K-M test to investigate the relationship between the expression of IGFBP7-AS1 and the overall survival time of UCEC patients, as well as pan-cancer patients. Also, K-M curves as the results of survival analysis were generated using “survival” package in R, showing the different survival condition in high and low IGFBP7-AS1 expression groups. Specifically, univariate and multivariate Cox proportional hazard ratio analysis for IGFBP7-AS1 levels and other clinical variables was also applied here with “survival” package, in order to identify the real effects of IGFBP7-AS1 on the prognostic survival of UCEC patients.
Gene enrichment analysis
To learn about the biological function and related signal pathways of IGFBP7-AS1, GO terms for biology process, cell component, and molecular function and KEGG enrichment pathways were all analyzed by “clusterprofiler” package of R and the results were then visualized into heatmaps through R package “enrichplot” .
Immune microenvironment and immune infiltration analysis
Stromal and immune scores were calculated by ESTIMATE algorithm for each pan-cancer sample and cancer samples were then divided into high and low scores groups by medium stromal and immune scores . In addition, Cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT), a deconvolution algorithm based on RNA-seq data, was applied to estimate the proportion of 22 immune-infiltration cell types so as to reveal the correlation between IGFBP7-AS1 and immune response . The different infiltration of degrees of immune infiltration lymphocyte types, immune-related pathway and functions in different IGFBP7-AS1 expression groups of UCEC patients were explored by ssGSEA using R package “GSVA”, in order to reflect IGFBP7-AS1’s effect on immune infiltration .
TMB and MSI calculation
Tumor mutation burden (TMB) and microsatellite instability (MSI) information was downloaded from TCGA database. TMB was calculated as somatic mutation incidences per million base pair in the code region of tumor cell genome , and MSI was counted by the times of insertion or deletion events that occurred in repeating sequences of genes (Cortes-Ciriano, Lee, Park, Kim, & Park, ).
Drug sensitivity analysis
CellMiner database (https://discover.nci.nih.gov/cellminer/home.do) was used to obtain the drug sensitivity data of different cancer cell lines, which was further filtered by the standard of clinical laboratory verification and FDA certification . Then we used Pearson correlation to figure the connection between IGFBP7-AS1 expression with drug sensitivity.
Statistical analysis was completed with R version 3.6.3 software and related R packages. The gene expression comparisons between different groups were tested by Wilcoxon rank sum test. Spearman rank correlation were used to measure the correlation between IGFBP7-AS1 expression and some targets of research, including target gene expression, TMB, MSI immune cell score and so on. A two-sided P-value lower than 0.05 was considered statistically significant. Most of the graphs were plotted through R packages “ggplot2”, and the gene expression comparison plots of paired samples in RT-qPCR were drawn by GraphPad Prism 8.
Screening prognosis-associated key eRNAs and target genes in UCEC
One thousand five hundred eighty eRNA expression data was extracted from the samples of UCEC in TCGA database, as well as the survival and clinical phenotypes dataset, and 120 eRNA gene expressions were found associated with the UCEC patients’ overall survival according to the results from the Kaplan–Meier log-rank test (p < 0.05). Spearman rank correlation coefficients were then used here to identify the significant correlations between the prognostic eRNAs and their predicted target genes thus helping with the further screening (Spearman’s rank correlation coefficient R > 0.5, p < 0.05). Finally, the relevant information of the 48 key eRNAs and their target genes was summarized in Table 2. Especially, the UCEC patients with high expression of eRNA IGFBP7-AS1 showed better overall survival performance compared to those low IGFBP7-AS1 expression ones (Fig. 1A, p < 0.05). In addition, as indicated in Fig. 1B, IGFBP7-AS1 had a closely positive correlation with IGFBP7 (R = 0.51, p < 0.001).
Data validation of IGFBP7-AS1 and IGFBP7 levels by reverse transcription quantitative PCR (RT-qPCR)
qRT-PCR was used to measure the expression levels of IGFBP7-AS1 and IGFBP7 in UCEC samples and paired adjacent samples. Down-regulation of IGFBP7-AS1 and IGFBP7 was seen in tumors when comparing it with adjacent samples, in accordance with its expression trend in the TCGA dataset (Fig. 2A, B). A significantly positive correlation was seen between these two genes in tumor tissues (R = 0.776, p < 0.001, Fig. 2C).
Correlation of IGFBP7-AS1 expression with clinical features and prognosis survival of UCEC
Using data from GEPIA database, we compared the IGFBP7-AS1 and IGFBP7 expression levels across UCEC normal and cancer tissues, and it was clear that they all had decreased expression in tumor tissue (Fig. 3A, B). Subsequently, we explored the relationship of IGFBP7-AS1 and IGFBP7 expression with clinicopathologic variables of UCEC and found that the expression level of IGFBP7-AS1 was significantly different in different age (Fig. 3C), tumor grade (Fig. 3D) and histological type (Fig. 3E) group; however, clinical stage showed no significant correlation with IGFBP7-AS1 levels (Fig. 3F). Besides, IGFBP7 was found related to patient age (Fig. 3G), tumor grade (Fig. 3H), but had no evidence for its correlation with histological type (Fig. 3I) and clinical stage (Fig. 3J). Using Cox regression, we further investigated the relationship between IGFBP7-AS1 and OS, along with the clinicopathologic characters before. The results were listed in Table 3, and some parameters, including clinical stage, tumor grade and IGFBP7-AS1 levels, turned out to be linked with OS both in univariate and multivariate regression. Figure 3K, showing a forest boxplot from Multivariate analysis, suggested that mRNA IGFBP7-AS1 was an independent favorable prognostic factor to UCEC. Since the variable age and histological type were insignificant in the multivariate analysis, another multivariate Cox model with only the significant ones of clinical stage, tumor grade and IGFBP7-AS1 levels was conducted to verify the results, in which all the three variables remained to be significant, supporting the finding of IGFBP7-AS1’s independent prognosis for UCEC (Supplementary Table 1).
GO enrichment analysis and KEGG pathway analysis were performed here so as to provide more information about the biological function of IGFBP7-AS1. Figure 4A, B were the bubble plots about the GO terms and KEGG signaling pathways of IGFBP7-AS1. Top 10 terms for biological process (BP), cellular component (CC), molecular function (MF) was extracted as the most important ones associated with IGFBP7-AS1——leukocyte migration in BP and extracellular matrix in CC were demonstrated to have strong correlation with IGFBP7-AS1. The highest-ranking signaling pathways in KEGG were PI3K − Akt signaling pathway, Neuroactive ligand − receptor interaction, Cytokine − cytokine receptor interaction and MAPK signaling pathway, which means IGFBP7-AS1 may regulate leukocyte migration through these pathways to influence on the occurrence and development of tumors. Moreover, Neuroactive ligand − receptor interaction and Cytokine − cytokine receptor interaction were all immune related pathways, which suggested the connection between IGFBP7-AS1 and immune activities.
IGFBP7-AS1 expression correlated with immune status
The tumor immune microenvironment (TME) plays an essential role in activating heterogeneity among cancer cells, thus stimulating multidrug resistance and resulting in the occurrence, development, and metastasis of tumors. The tumor immune microenvironment has been a critical part in tumor research so we calculated the stromal score, immune score, ESTIMATE score and tumor purity of each tumor sample by ESTIMATE and compared them in high and low IGFBP7-AS1 and IGFBP7 expression groups. A high stromal and immune score or low tumor purity usually contributed to better prognosis and response to immunotherapy. The results indicated that increased IGFBP7-AS1 and IGFBP7 expression was closely related with higher stromal score, immune score, estimate score and lower tumor purity (Fig. 5A-H). So high levels of IGFBP7-AS1 and IGFBP7 may be beneficial for immunotherapy effect of UCEC.
To have a better knowledge of their association with immune related cells, we employed CIBERSORT to further process the biological role of IGFBP7-AS1 in TME. Results in boxplot of Fig. 6A indicated that T cells CD4 memory activated, T cells follicular helper, T cells regulatory (Tregs), T cells gamma delta, Macrophages M1, Macrophages M2 and Mast cell resting were the immune cells affected by IGFBP7-AS1 expression. Among them, T cells CD4 memory activated, T cells follicular helper, T cells gamma delta, Macrophages M1 and Macrophages M2 showed low proportion in the high IGFBP7-AS1 expression groups compared to the low. In contrast, T cells regulatory (Tregs) proportion and Mast cell resting were clearly up-regulated in high-level IGFBP7-AS1 group.
We further investigated the detailed correlation between immune cells and expression of IGFBP7-AS1 and IGFBP7. Levels of IGFBP7-AS1 expression were inversely correlated with memory B cells, Macrophages M1, Macrophages M2, T cells follicular helper and T cells gamma delta, whereas native B cells, Mast cells resting and T cells regulatory (Tregs) were positively associated with IGFBP7-AS1 expression (Fig. 6B-I). At the same time, IGFBP7 expression levels were positively correlated with levels of Mast cells resting, Plasma cells, T cells CD4 memory resting and T cells gamma delta, and diversely correlated with memory B cells, Dendritic cells activated, Macrophages M0 and T cells follicular helper (Fig. 6J-Q).
Tumor-infiltrating lymphocytes (TILs) play a key role in sentinel lymph node status and overall survival rate prediction . Immune infiltration scores of TIL types and immune-related substances were calculated in order to determine if the IGFBP7-AS1 and IGFBP7 connected with the immune infiltration. We assessed differences between high- and low- levels in immune infiltration cell and marked differences were shown in aDCs, CCR, CD8 + T cells, Check point, Cytolytic activity, HLA, iDCs, Mast cells, Neutrophils, T helper cells, TIL, Treg and Type II IFN Reponses between IGFBP7-AS1 subgroups. Except DCs, MHC class I, NK cells, Th2 cells and Type I IFN Response, IGFBP7 expressed significantly different in other immune cells (Fig. 6R, S).
IGFBP7-AS1 expression correlated with immune checkpoint genes expression
Then, we analyzed the immune checkpoint genes expression with IGFBP7-AS1 and IGFBP7. We could find in the pictures plotted in Fig. 7A, B, that these immune checkpoints related genes (CD200, NRP1, LAIR1, CD244, CD40LG, CD48, CD28, HAVCR2, TNFSF14, HHLA2, CD70, CD27, TNFRSF4, TNFSF15, TNFRSF9) expressed differently between high- and low-expression groups of both IGFBP7-AS1 and IGFBP7. PD-1 is a protein on the surface of T and B cells that acts as an essential role in regulating the immune system's response to the cells in the human body by down-regulating the immune system and developing self-tolerance by suppression of T cell inflammatory activity . CTLA4 is a protein receptor that works as a member of immune checkpoints and has effect on immune responses downregulation . They were essential immune checkpoints in carcinoma proliferation research and immune checkpoint blockade therapy by targeting PD-1and CTLA4 revealed promising clinical effects. These key modulators were detected to express more accompanied with higher IGFBP7-AS1 levels (Fig. 7C, D).
IGFBP7-AS1 Had a Significant Negative Correlation with TMB
TMB (tumor mutational burden) is a genetic characteristic of tumorous tissue that can provide valuable information for our cancer study and relevant treatment. It is defined as the number of non-inherited mutations per Mb of investigated genomic sequence and plays a key role in cancer prognosis and response to tumor immunotherapy treatment . We calculated the TMB and discovered its connection with IGFBP7-AS1 and IGFBP7. The results suggested that both of IGFBP7-AS1 and IGFBP7 had a significant negative correlation with TMB. The more expression of IGFBP7-AS1 or IGFBP7, the less expression of TMB (Fig. 7E, F).
IGFBP7-AS1 expression correlated with m6A genes expression
m6A(6-methyladenine), a common type of mRNA methylation, impacts on the cancer development for its versatile functions in various physiological processes . The connections between m6A and numerous cancer types have been indicated in many reports involving breast cancer, prostate cancer, stomach cancer and so on. There are some genes (HNRNPC, RBM15, METTL14, YTHDC2, WTAP, YTHDF1, YTHDC1, FTO, YTHDF2) proved to regulate the modification levels of m6A and then affect the cancer cell proliferation. These genes all expressed variously in different IGFBP7-AS1 and IGFBP7 levels groups (Fig. 7G, H). Except for FTO with IGFBP7-AS1, most of the m6A genes were in down-expressed condition with high levels of IGFBP7-AS1 and IGFBP7.
Pan-cancer verification of IGFBP7-AS1 and IGFBP7 levels
Comprehensive omics analysis of IGFBP7-AS1 across 33 cancers was conducted to verify IGFBP7-AS1’s role in diverse cancers and except to provide robust evidence for potential tumor research. We obtained the expression of IGFBP7-AS1 and IGFBP7 across 33 cancer samples and respective normal tissues from TCGA project. IGFBP7-AS1 and IGFBP7 expression differences could be found in BLCA, BRCA, CESC, CHOL, COAD, GBM, HNSC, KICH, KIRC, KIRP, LUAD, LUSC, PRAD, THCA and UCEC. Besides, IGFBP7 levels differed in SARC, STAD and IGFBP7-AS1 differed in LIHC as well (Supplementary Fig. 1A, B).
Next, we identified the prognostic value of IGFBP7 − AS1 for pan-cancer. Patients in the high-expression group survived longer than those in the low-expression group in LAML, LUAD and UCEC and the outcomes of LGG, MESO and STAD were opposite in Kaplan–Meier survival curves (Supplementary Fig. 1C-H). Furthermore, univariate Cox proportional hazard regressions were modeled to understand the altered expression of IGFBP7 − AS1 and IGFBP7 with patient overall survival, and the directions of prognostic effect were manifested varied depending on cancer types (Supplementary Fig. 1I).
We first investigated the association of IGFBP7 − AS1 and IGFBP7 expression with six subtypes of immune infiltration: C1 (wound healing), C2 (INF-r dominant), C3 (inflammation), C4 (lymphopenia dominant), C5 (immunologically quiet), and C6 (TGFβ dominant) . Both IGFBP7 − AS1 and IGFBP7 were strongly correlated with immune subtypes in pan-cancer (Fig. 8A). Next, we evaluated whether a link between IGFBP7 − AS1 and IGFBP7 and expression of genes recognized as checkpoint components existed. Co-expression with immune checkpoint related genes in pan-cancers was exhibited in heatmaps of Fig. 8B, C. High positive correlations could be seen in tumors like CHOL, ESCA, LGG, LIHC, LUAD, LUSC and so on while negative ones were mainly appeared in MESO, LAML and THYM.
Also, the relationship with immune infiltration cells were analyzed here, and results were presented in Fig. 8D, E, which suggested that IGFBP7 − AS1 and IGFBP7 were strongly associated with immune cells and most of the correlations were positive.
We have done TME analysis of IGFBP7 − AS1 and IGFBP7 in UCEC, so in this section we intended to confirm IGFBP7 − AS1 and IGFBP7’s role in tumor progression and immune response in pan-cancer. Stromal, immune, ESTIMATE score and tumor purity of 33 pan-cancer were shown in Fig. 9A-D. IGFBP7 − AS1 and IGFBP7 expression were positively relevant to stromal score in majority pan-cancers with exceptions of SARC and STAD. In immune score analysis, IGFBP7 − AS1 and IGFBP7 still showed positive correlation with immune score across most of cancers, with the obvious exception of THYM, LAML, and IGFBP7 − AS1 in MESO, OV, SARC, SKCM. Up-level of ESTIMATE score and down-level tumor purity usually came with high levels of IGFBP7 − AS1 and IGFBP7(exceptions were TGCT, LAML, SARC in ESTIMATE score or tumor purity).
After that, we focused on the association with TMB and MSI across 33 pan-cancer and the correlation was exhibited in radar plot Fig. 9E-H. Most of the cancer including BLCA, BRCA, CESC, COAD, DLBC, HNSC, KIRC, LIHC, LUAD, PAAD, PRAD, SARC, STAD, UCEC, UVM had a significantly negative correlation with IGFBP7 − AS1 in terms of TMB, apart from THYM, which showed reverse result with significance. As for IGFBP7, its expression was significant associated with TMB in 21 out of 33 cancers, and most of the trends were similar with IGFBP7 − AS1. The results were consistent with the TMB analysis in UCEC.
MSI (microsatellite instability) is a good marker for determining a prognosis for cancer treatments, and the elevated MSI may be an indicator of higher tumor risk for the reason that the tumor’s disrupted function increases the gene instability . The MSI radar plots reflected that in COAD, LIHC, STAD, as well as UCEC, IGFBP7 − AS1 was acted as negative correlation with MSI, contrary to the condition in HNSC and MESO. While for IGFBP7, the significant associations with MSI turned out to be all negative in DLBC, GBM, LIHC, LUAD, LUAC, PAAD, STAD and UCEC.
Expression-based RNA stemness score (RNAss) and methylation-based DNA stemness score (DNAss) were used to reflect the features of tumor stem cells . Higher stemness index values are associated with stimulate biological processes in tumor stem cells and promoted tumor dedifferentiation . At the same time, we found IGFBP7 − AS1 and IGFBP7 were negatively correlated with RNAss and DNAss across most of the cancers (Fig. 9I, J), and a decrease in tumor stemness was usually linked with better survival. All these evidences can help to explain the IGFBP7 − AS1 and IGFBP7’s mechanisms in modulating the overall survival.
Using CellMiner database, we assessed the influence of IGFBP7 − AS1 and IGFBP7 on drug sensitivity, and such research could help us to provide better precision treatment for patients. We chose the significantly relevant ones according to P values here. With the increase in the expression of IGFBP7 − AS1, the less sensitive the cells were to chemotherapeutic drugs involving Veliparib, Mocetinostat, CYC-116 and Alisertib (Supplementary Fig. 2).
Insulin-like growth factor-binding protein 7(IGFBP7) is a new member of the IGFBPs gene family and IGFBP7 Antisense RNA 1(IGFBP7-AS1) is a related RNA gene affiliated with lncRNA class. There were already several researches providing resources for exploring the mechanism of IGFBP7-AS1 in tumorigenesis. One of them revealed IGFBP7-AS1 might have influence on glioma cell survival by regulating tumor growth and migration. In addition, further functional supported experiments demonstrated the knockdown of IGFBP7-AS1 inhibited the invasion and viability of U87 and U251 glioma cells (D. . There were also results indicating that Epstein Barr-virus (EBV) may dysregulate host lncRNAs like IGFBP7-AS1 to regulate its own replication, thus facilitating tumorigenesis . Although the IGFBP7-AS1’s intervention in the progression of tumor proliferation and migration of few specific cancer types has been confirmed, we still lack of the mechanistic explanations for the prognostic ability of IGFBP7-AS1 in majority cancers, not to mention in UCEC. Also, in-depth researches from a clinical perspective are still waiting to be conducted. It may therefore become a novel target for future molecular therapy development.
In our study, IGFBP7-AS1 was determined as the key eRNA of UCEC with its target gene IGFBP7 according to the criteria we set. High expression group of IGFBP7-AS1 was related to better clinical traits, and IGFBP7-AS1 had decreased expression in tumor samples. IGFBP7-AS1 was also identified as an independent protective factor for the clinical outcome of UCEC patients by both univariate and multivariate cox regression, in accordance with the results from Kaplan–Meier survival curves. All this evidence supported that reduced IGFBP7-AS1 expression was associated with unfavorable prognosis in UCEC patients. UCEC patients with low IGFBP7-AS1 expression were more likely to present a more advanced tumor grade and age group, which was the same with IGFBP7. Thus, we suggested that IGFBP7-AS1 had a potential influence in UCEC progression.
Insulin Growth Factor Binding Protein 7 (IGFBP7) is a protein belonging to IGFBP superfamily. There is evidence indicating that IGFBPs function as transporters of the insulin-like growth factors (IGF), lengthen their half-time and regulate their access to their receptors, thereby affecting cell growth, differentiation and metabolic processes in human body . IGFBP7 is a crucial player responsible for tumorigenesis and tumor progression. Accumulating studies clarified that IGFBP7 regulated varied in diverse types of malignancies. Its importance in many tumors like breast cancer, gastric cancer and lung tumors has been revealed in a lot of previous studies (Y. [13, 21, 39].
GO terms analysis revealed upregulated IGFBP7-AS1 was to be primarily connected with leukocyte migration, skeletal system development and extracellular matrix. Meanwhile, KEGG pathways analysis showed five pathways including PI3K − Akt signaling pathway, Neuroactive ligand − receptor interaction, Cytokine − cytokine receptor interaction, MAPK signaling pathway and Ras signaling pathways were showed significantly correlated with IGFBP7-AS1. PI3K − Akt signaling pathway has been demonstrated dispensable to cell proliferation and apoptosis and it has been confirmed that inhibition of PI3K/Akt can help to suppress the proliferation of tumor cell . Based on the results before, we can infer that IGFBP7-AS1 may regulate PI3K − Akt signaling pathway to alter the occurrence and development of tumors.
Notably, another essential finding in this study was the crucial role of IGFBP7-AS1 in immune response. TME plays a key role in tumor cell sustained growth, invasion and metastasis, and stromal cells and immune cells within TME represent attractive therapeutic targets with decreased risk of resistance and tumor recurrence . Our study demonstrated that higher the expression of IGFBP7-AS1 and IGFBP7, the content of stromal and immune cells was higher, and the tumor purity was lower, which indicated that patient with higher expression level of IGFBP7-AS1 and IGFBP7 may have better therapeutic effect from immunotherapy.
When it comes to the correlation with immune related cells, IGFBP7-AS1 expression was negatively related with score of Macrophages M1 and M2, in accordance with the results from fraction variance of Macrophages M1 and M2 in different IGFBP7-AS1 groups. Many studies have shown on the one hand macrophage abundance in a tumor had an inverse correlation with patient survival, on the other hand it had a positive linkage with the resistance to chemotherapy . Then, the increase in IGFBP7-AS1 and IGFBP7 expression indicated a higher presence of resting Mast cells. Mast cells could induce a privileged immunological environment with proliferation and the invasive ability of the cancer cells, as well as favoring the growth of a tumoral tissue with the development of its related vascular stroma (Aller, Arias, Arias, & Arias, ). Given that the significant infiltration of immune cells like macrophages, microglia, neutrophils, and monocytes may result in unfavorable prognosis in a number of cancer types , high IGFBP7-AS1 expression may indicate better UCEC overall survival. Also, high expressed IGFBP7-AS1 was accompanied with high score of TILs and immune-related substances like CCR, iDCs and HLA. These findings verified the prognostic ability of IGFBP7-AS1 in UCEC.
In recent years, significant advance has been made in studies focusing on T cells and related immune checkpoint, for their essential roles in immune infiltration and therapy. Blockade of inhibitory immune checkpoints can promote the recovery of T cell activation and prevent immune escape of cancer cells within the tumor microenvironment, while activation of stimulatory immune checkpoints can augment the effect of immune response (Y. . To data, drugs that block immune checkpoints like PD1-PDL1 or CTLA4 have gradually become the mainstream of oncology. In our study, IGFBP7-AS1 and IGFBP7 expression were observed significant connected with various types of T cell, involving CD8 + T cells, T helper cells, T cell follicular helper, T cells gamma delta, Tregs, Th1 cells, and multiple immune checkpoints including PD1 and CTLA4 we mentioned before, suggesting that targeting IGFBP7-AS1 according to its traits in immune regulation of UCEC deserved research focus.
As previous studies demonstrated, several cancer types with high TMB may be attractive aims for immuno-oncology treatment development . So, the results that high levels of IGFBP7-AS1 and IGFBP7 were both significantly associated with low TMB expression may provide implications for guidance on related immune therapy.
Finally, we also conducted pan-cancer analysis for verification, and drew the similar conclusion in most cancers. IGFBP7-AS1 expression was significantly associated with TMB, MSI, TME and Immune infiltration cells in pan-cancers, and impacted on the survival outcomes. Specifically, correlation between IGFBP7-AS1’s high expression and the upregulation of immune subtype C5 in pan-cancers was found. In previous study, IDH mutations has been reported enriched in C5, which associated with TME composition and decreased leukocyte chemotaxis, contributing to fewer tumor-associated immune cells and better clinical effect . In drug sensitivity analysis, chemotherapeutic drug sensitivity of Veliparib, Mocetinostat, CYC-116 and Alisertib was negatively correlated with IGFBP7-AS1 expression, which may illustrate IGFBP7-AS1’s potential connection with their targets like PARP, HDAC and so on.
Besides, as the target gene of IGFBP7-AS1, IGFBP7 has been confirmed in previous study that inhibits IGF signaling by the way of binding the IGF-1 receptor (IGF-1R). IGF signaling activation is one of the major oncogenic events in diverse cancers, so IGFBP7 has been reported as a candidate tumor suppressor in several cancers ,Y. [26, 35]. Moreover, IGFBP7 is manifested able to alter cell sensitivity to chemotherapeutic drugs, suggesting its beneficial value in anticancer therapies . Considering the results from our study, we suggested IGFBP7-AS1 may be likely to influence on IGFBP7 expression, thereby affecting the UCEC patient’s survival.
The disadvantage of this study is that the molecule function of IGFBP7-AS1 in UCEC has still been uncovered and the results need further clinical verification. Therefore, a specific clinical cohort is needed to further assess the diagnostic and prognostic potential of IGFBP7-AS1 in UCEC. Moreover, we demonstrated that IGFBP7-AS1 has roles in multiple kinds of cancer, including LUAD, STAD and UCEC, and the emerging roles of IGFBP7-AS1 in these cancers are also waiting to be characterized.
Taking together, increased IGFBP7-AS1 expression was related to positive prognosis of UCEC. Our comprehensive omics analysis revealed IGFBP7-AS1 may improve the tumor immune microenvironment, and impact on immune lymphocytes like T cells, Mast cells, Macrophages M1 and M2, including key immune checkpoints with its target gene IGFBP7, so that help to induce protective immunity in UCEC and provide a new reference for potential UCEC immunotherapy.
Availability of data and materials
The data and materials used to support the findings of this study are available from the corresponding author on reasonable request.
Akiel M, Guo C, Li X, Rajasekaran D, Mendoza RG, Robertson CL, Sarkar D. IGFBP7 Deletion Promotes Hepatocellular Carcinoma. Cancer Res. 2017;77(15):4014–25. https://doi.org/10.1158/0008-5472.CAN-16-2885.
Aller MA, Arias A, Arias JI, Arias J. Carcinogenesis: the cancer cell-mast cell connection. Inflamm Res. 2019;68(2):103–16. https://doi.org/10.1007/s00011-018-1201-4.
Amankulor NM, Kim Y, Arora S, Kargl J, Szulzewsky F, Hanke M, Holland EC. Mutant IDH1 regulates the tumor-associated immune system in gliomas. Genes Dev. 2017;31(8):774–86. https://doi.org/10.1101/gad.294991.116.
Barbie DA, Tamayo P, Boehm JS, Kim SY, Moody SE, Dunn IF, Hahn WC. Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1. Nature. 2009;462(7269):108–12. https://doi.org/10.1038/nature08460.
Ben-Porath I, Thomson MW, Carey VJ, Ge R, Bell GW, Regev A, Weinberg RA. An embryonic stem cell-like gene expression signature in poorly differentiated aggressive human tumors. Nat Genet. 2008;40(5):499–507. https://doi.org/10.1038/ng.127.
Bhan A, Mandal SS. LncRNA HOTAIR: A master regulator of chromatin dynamics and cancer. Biochim Biophys Acta. 2015;1856(1):151–64. https://doi.org/10.1016/j.bbcan.2015.07.001.
Bhan A, Soleimani M, Mandal SS. Long Noncoding RNA and Cancer: A New Paradigm. Cancer Res. 2017;77(15):3965–81. https://doi.org/10.1158/0008-5472.CAN-16-2634.
Braun MM, Overbeek-Wager EA, Grumbo RJ. Diagnosis and Management of Endometrial Cancer. Am Fam Physician. 2016;93(6):468–74.
Campbell BB, Light N, Fabrizio D, Zatzman M, Fuligni F, de Borja R, Shlien A. Comprehensive Analysis of Hypermutation in Human Cancer. Cell. 2017;171(5):1042–56. https://doi.org/10.1016/j.cell.2017.09.048 (e1010).
Chalmers ZR, Connelly CF, Fabrizio D, Gay L, Ali SM, Ennis R, Frampton GM. Analysis of 100,000 human cancer genomes reveals the landscape of tumor mutational burden. Genome Med. 2017;9(1):34. https://doi.org/10.1186/s13073-017-0424-2.
Chan TA, Yarchoan M, Jaffee E, Swanton C, Quezada SA, Stenzinger A, Peters S. Development of tumor mutation burden as an immunotherapy biomarker: utility for the oncology clinic. Ann Oncol. 2019;30(1):44–56. https://doi.org/10.1093/annonc/mdy495.
Chen H, Li C, Peng X, Zhou Z, Weinstein JN, Cancer Genome Atlas Research, N, Liang H. A Pan-Cancer Analysis of Enhancer Expression in Nearly 9000 Patient Samples. Cell. 2018;173(2):386–99. https://doi.org/10.1016/j.cell.2018.03.027 (e312).
Chen Y, Pacyna-Gengelbach M, Ye F, Knosel T, Lund P, Deutschmann N, Petersen I. Insulin-like growth factor binding protein-related protein 1 (IGFBP-rP1) has potential tumour-suppressive activity in human lung cancer. J Pathol. 2007;211(4):431–8. https://doi.org/10.1002/path.2132.
Cheng L, Han T, Chen B, Nie K, Peng W. TBX5-AS1, an enhancer RNA, is a potential novel prognostic biomarker for lung adenocarcinoma. BMC Cancer. 2021;21(1):794. https://doi.org/10.1186/s12885-021-08517-w.
Cortes-Ciriano I, Lee S, Park WY, Kim TM, Park PJ. A molecular portrait of microsatellite instability across multiple cancers. Nat Commun. 2017;8:15180. https://doi.org/10.1038/ncomms15180.
De’Angelis GL, Bottarelli L, Azzoni C, De’Angelis N, Leandro G, Di Mario F, Negri F. Microsatellite instability in colorectal cancer. Acta Biomed. 2018;89(9-S):97–101. https://doi.org/10.23750/abm.v89i9-S.7960.
De Santa F, Barozzi I, Mietton F, Ghisletti S, Polletti S, Tusi BK, Natoli G. A large fraction of extragenic RNA pol II transcription sites overlap enhancers. PLoS Biol. 2010;8(5):e1000384. https://doi.org/10.1371/journal.pbio.1000384.
Genin M, Clement F, Fattaccioli A, Raes M, Michiels C. M1 and M2 macrophages derived from THP-1 cells differentially modulate the response of cancer cells to etoposide. BMC Cancer. 2015;15:577. https://doi.org/10.1186/s12885-015-1546-9.
Hall M, Liu H, Malafa M, Centeno B, Hodul PJ, Pimiento J, Sarnaik AA. Expansion of tumor-infiltrating lymphocytes (TIL) from human pancreatic tumors. J Immunother Cancer. 2016;4:61. https://doi.org/10.1186/s40425-016-0164-7.
Han Y, Liu D, Li L. PD-1/PD-L1 pathway: current researches in cancer. Am J Cancer Res. 2020;10(3):727–42.
Jin L, Shen F, Weinfeld M, Sergi C. Insulin Growth Factor Binding Protein 7 (IGFBP7)-Related Cancer and IGFBP3 and IGFBP7 Crosstalk. Front Oncol. 2020;10:727. https://doi.org/10.3389/fonc.2020.00727.
Lee JH, Xiong F, Li W. Enhancer RNAs in cancer: regulation, mechanisms and therapeutic potential. RNA Biol. 2020;17(11):1550–9. https://doi.org/10.1080/15476286.2020.1712895.
Lee YC, Lheureux S, Oza AM. Treatment strategies for endometrial cancer: current practice and perspective. Curr Opin Obstet Gynecol. 2017;29(1):47–58. https://doi.org/10.1097/GCO.0000000000000338.
Li D, Lu J, Li H, Qi S, Yu L. Identification of a long noncoding RNA signature to predict outcomes of glioblastoma. Mol Med Rep. 2019;19(6):5406–16. https://doi.org/10.3892/mmr.2019.10184.
Li W, Notani D, Rosenfeld MGJNRG. Enhancers as non-coding RNA transcription units: recent insights and future perspectives. Nat Rev Genet. 2016;17(4):207–207.
Li Y, Xi Y, Zhu G, Jia J, Huang H, Liu Y, Liu L. Downregulated IGFBP7 facilitates liver metastasis by modulating epithelialmesenchymal transition in colon cancer. Oncol Rep. 2019;42(5):1935–45. https://doi.org/10.3892/or.2019.7303.
Liu Y, Hu Z, Zhang Y, Wang C. Long non-coding RNAs in Epstein-Barr virus-related cancer. Cancer Cell Int. 2021;21(1):278. https://doi.org/10.1186/s12935-021-01986-w.
Malta TM, Sokolov A, Gentles AJ, Burzykowski T, Poisson L, Weinstein JN, Wiznerowicz M. Machine Learning Identifies Stemness Features Associated with Oncogenic Dedifferentiation. Cell. 2018;173(2):338–54. https://doi.org/10.1016/j.cell.2018.03.034.
Moore K, Brewer MA. Endometrial Cancer: Is This a New Disease? Am Soc Clin Oncol Educ Book. 2017;37:435–42. https://doi.org/10.14694/EDBK_175666 (10.1200/EDBK_175666).
Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, Alizadeh AA. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7. https://doi.org/10.1038/nmeth.3337.
Ouyang D, Li R, Li Y, Zhu X. A 7-lncRNA signature predict prognosis of Uterine corpus endometrial carcinoma. J Cell Biochem. 2019;120(10):18465–77. https://doi.org/10.1002/jcb.29164.
Pai CS, Simons DM, Lu X, Evans M, Wei J, Wang YH, Fong L. Tumor-conditional anti-CTLA4 uncouples antitumor efficacy from immunotherapy-related toxicity. J Clin Invest. 2019;129(1):349–63. https://doi.org/10.1172/JCI123391.
Quail DF, Joyce JA. Microenvironmental regulation of tumor progression and metastasis. Nat Med. 2013;19(11):1423–37. https://doi.org/10.1038/nm.3394.
Rajagopal T, Talluri S, Akshaya RL, Dunna NR. HOTAIR LncRNA: A novel oncogenic propellant in human cancer. Clin Chim Acta. 2020;503:1–18. https://doi.org/10.1016/j.cca.2019.12.028.
Ruan W, Xu E, Xu F, Ma Y, Deng H, Huang Q, Lai M. IGFBP7 plays a potential tumor suppressor role in colorectal carcinogenesis. Cancer Biol Ther. 2007;6(3):354–9. https://doi.org/10.4161/cbt.6.3.3702.
Shankavaram UT, Varma S, Kane D, Sunshine M, Chary KK, Reinhold WC, Weinstein JN. Cell Miner: a relational database and query tool for the NCI-60 cancer cell lines. BMC Genomics. 2009;10:277. https://doi.org/10.1186/1471-2164-10-277.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70(1):7–30. https://doi.org/10.3322/caac.21590.
Sun T, Wu R, Ming L. The role of m6A RNA methylation in cancer. Biomed Pharmacother. 2019;112:108613. https://doi.org/10.1016/j.biopha.2019.108613.
Swisshelm K, Ryan K, Tsuchiya K, Sager R. Enhanced expression of an insulin growth factor-like binding protein (mac25) in senescent human mammary epithelial cells and induced expression with retinoic acid. Proc Natl Acad Sci U S A. 1995;92(10):4472–6. https://doi.org/10.1073/pnas.92.10.4472.
Tamborero D, Rubio-Perez C, Muinos F, Sabarinathan R, Piulats JM, Muntasell A, Gonzalez-Perez A. A Pan-cancer Landscape of Interactions between Solid Tumors and Infiltrating Immune Cell Populations. Clin Cancer Res. 2018;24(15):3717–28. https://doi.org/10.1158/1078-0432.CCR-17-3509.
Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–102. https://doi.org/10.1093/nar/gkx247.
Thorsson V, Gibbs DL, Brown SD, Wolf D, Bortone DS, Ou Yang TH, Shmulevich I. The Immune Landscape of Cancer. Immunity. 2019;51(2):411–2. https://doi.org/10.1016/j.immuni.2019.08.004.
Tomimaru Y, Eguchi H, Wada H, Noda T, Murakami M, Kobayashi S, Nagano H. Insulin-like growth factor-binding protein 7 alters the sensitivity to interferon-based anticancer therapy in hepatocellular carcinoma cells. Br J Cancer. 2010;102(10):1483–90. https://doi.org/10.1038/sj.bjc.6605669.
Wang Z, Jensen MA, Zenklusen JC. A Practical Guide to The Cancer Genome Atlas (TCGA). Methods Mol Biol. 2016;1418:111–41. https://doi.org/10.1007/978-1-4939-3578-9_6.
Xi Y, Wang W, Chen W, Han K, Qiao L, Chen W. Incidence and mortality of corpus uteri cancer in China, 2008–2012. Chin J Cancer Res. 2019;31(3):435–42. https://doi.org/10.21147/j.issn.1000-9604.2019.03.05.
Xie Y, Shi X, Sheng K, Han G, Li W, Zhao Q, Gu Y. PI3K/Akt signaling transduction pathway, erythropoiesis and glycolysis in hypoxia (Review). Mol Med Rep. 2019;19(2):783–91. https://doi.org/10.3892/mmr.2018.9713.
Yoshihara K, Shahmoradgoli M, Martinez E, Vegesna R, Kim H, Torres-Garcia W, Verhaak RG. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612. https://doi.org/10.1038/ncomms3612.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7. https://doi.org/10.1089/omi.2011.0118.
Zhang Y, Zheng J. Functions of Immune Checkpoint Molecules Beyond Immune Evasion. Adv Exp Med Biol. 2020;1248:201–26. https://doi.org/10.1007/978-981-15-3266-5_9.
Zhang Z, Lee JH, Ruan H, Ye Y, Krakowiak J, Hu Q, Han L. Transcriptional landscape and clinical utility of enhancer RNAs for eRNA-targeted therapy in cancer. Nat Commun. 2019;10(1):4562. https://doi.org/10.1038/s41467-019-12543-5.
Zhou M, Zhang Z, Zhao H, Bao S, Sun J. A novel lncRNA-focus expression signature for survival prediction in endometrial carcinoma. BMC Cancer. 2018;18(1):39. https://doi.org/10.1186/s12885-017-3983-0.
Ethics approval and consent to participate
This study was approved by Ethics Committee of Nanjing Medical University.
Consent for publication
The authors have no conflict of interest to declare.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Multivariate analysis result for the stage, grade and IGFBP7-AS1 expression only with OS of UCEC patients.
Expression levels of IGFBP7-AS1 and IGFBP7 and association with overall survival in pan-cancer. (A, B) The expression of IGFBP7-AS1 and IGFBP7 in tumor and normal tissues across 33 pan-cancers. (C) The K-M survival curves of IGFBP7-AS1 expression groups in LAML, (D) LUAD, (E) MESO, (F) LGG, (G) STAD and (H) UCEC. (I) The forest plot for overall survival with 95% confidence intervals for 33 different cancer types.
The relationship between gene expression and drug sensitivity.
About this article
Cite this article
Liu, J., Yin, J., Wang, Y. et al. A comprehensive prognostic and immune analysis of enhancer RNA identifies IGFBP7-AS1 as a novel prognostic biomarker in Uterine Corpus Endometrial Carcinoma. Biol Proced Online 24, 9 (2022). https://doi.org/10.1186/s12575-022-00172-0
- Enhancer RNA
- Tumor immunity