- 1Deparment of Nutrition and Food Hygiene, Medical College of Soochow University, Suzhou, China
- 2Department of Obstetrics and Gynecology, The First Affiliated Hospital of Soochow University, Suzhou, China
- 3Department of Occupational and Environmental Health, School of Public Health, Medical College of Soochow University, Suzhou, China
- 4Department of Obstetrics and Gynecology, The Second Affiliated Hospital of Soochow University, Suzhou, China
- 5State Key Laboratory of Radiation Medicine and Protection, School of Radiation Medicine and Protection, Soochow University, Suzhou, China
- 6Department of Pathology, The First Affiliated Hospital of Soochow University, Suzhou, China
- 7Department of Epidemiology and Biostatistics, School of Public Health, Soochow University, Suzhou, China
Endometrial carcinoma (EC) is one of the common gynecological cancers with increasing incidence and revived mortality recently. Given the heterogeneity of tumors and the complexity of lncRNAs, a panel of lncRNA biomarkers might be more precise and stable for prognosis. In the present study, we developed a new lncRNA model to predict the prognosis of patients with EC. EC-associated differentially expressed long noncoding RNAs (lncRNAs) were identified from The Cancer Genome Atlas (TCGA). Univariate COX regression and least absolute shrinkage and selection operator (LASSO) model were selected to find the 8-independent prognostic lncRNAs of EC patient. Furthermore, the risk score of the 3-lncRNA signature for overall survival (OS) was identified as CTD-2377D24.6 expression × 0.206 + RP4-616B8.5 × 0.341 + RP11-389G6.3 × 0.343 by multivariate Cox regression analysis. According to the median cutoff value of this prognostic signature, the EC samples were divided into two groups, high-risk set (3-lncRNAs at high levels) and low-risk set (3-lncRNAs at low levels), and the Kaplan–Meier survival curves demonstrated that the low-risk set had a higher survival rate than the high-risk set. In addition, the 3-lncRNA signature was closely linked with histological subtype (p = 0.0001), advanced clinical stage (p = 0.011), and clinical grade (p < 0.0001) in EC patients. Our clinical samples also confirmed that RP4-616B8.5, RP11-389G6.3, and CTD-2377D24.6 levels were increased in tumor tissues by qRT-PCR and in situ hybridization. Intriguingly, the p-value of combined 3-lncRNAs was lower than that of each lncRNA, indicating that the 3-lncRNA signature also showed higher performance in EC tissue than paracancerous. Functional analysis revealed that cortactin might be involved in the mechanism of 3-lncRNA signatures. These findings provide the first hint that a panel of lncRNAs may play a critical role in the initiation and metastasis of EC, indicating a new signature for early diagnosis and therapeutic strategy of uterine corpus endometrial carcinoma.
Introduction
Globally, the incidence for uterine corpus endometrial carcinoma (UCEC) persistently increased with 1.3% per year from 2007–2016, in part due to continued declines in the fertility rate as well as increased obesity (Siegel et al., 2020). In China, the incidence of EC was also increasing from 2014, which ranked second in female reproductive malignancies on account of the increased risk factors such as diabetes and obesity (Chen et al., 2019). Although EC has a good prognosis with 5-year overall survival (OS) of 74–91%, the advanced or metastatic EC patients still have a poor prognosis due to tumor metastasis and poor differentiation (Piulats et al., 2017). Histological classification and the International Federation of Gynecology and Obstetrics (FIGO) staging system are the traditional treatment guideline and prognostic indicators (Pecorelli, 2009; Morice et al., 2016). However, distinct molecular characteristics have been demonstrated in the same stage and histology of cancers (Murali et al., 2014; Yang et al., 2016). With the development of precision medicine, a new therapeutic approach according to molecular profiling has been provided. In 2021, to improve outcomes of EC patients, molecular classification was recommended to select appropriate treatment regimens by the National Comprehensive Cancer Network (NCCN). Four molecular subgroups have been classified in 2013 based on the integrated genomic data of 373 endometrial carcinomas (Levine et al., 2013). Nevertheless, the integrated classification had limited application due to high expense and complex procedures. Therefore, identifying an efficient prognostic and diagnostic signature to guide clinical practice for EC is urgent.
Noncoding RNA was initially recognized as simply leaky transcription noise because they are not translated into proteins. However, numerous noncoding RNAs showed specific functions in cellular processes, as well as the dysregulation in human pathologies. Long noncoding RNA (lncRNA) is a class of noncoding transcripts with more than 200 nucleotides in length. Compelling studies reported that lncRNAs were associated with various human diseases including cancer by participating in biological processes widely (Schmitt and Chang, 2016; Peng et al., 2017; Yang et al., 2019). Meanwhile, accumulating evidence supported the potential ability of lncRNAs as cancer biomarkers (Lim et al., 2019; Li et al., 2020; Xie et al., 2020) and the prognostic value of lncRNAs (He et al., 2014; Sun et al., 2021). For example, Liu et al. systematically discussed the EC-related lncRNAs and their roles in different cancer hallmarks, including tumor growth, metastasis, maintenance of cancer stem cells, and chemoresistance (Liu H et al., 2019). Until now, some biomarkers for EC have been identified using gene expression profile data. However, these models are limited to a specific stage or grade of EC. For example, one study identified a prognostic model for patients with early-stage EC using reverse-phase protein arrays (Yang et al., 2016). Others found a prognostic value of immune, metabolic, or autophagy-related coding and noncoding lncRNAs for EC (Ouyang et al., 2019; Gao et al., 2020; Li and Wan, 2020; Wang et al., 2021). However, given the heterogeneity of EC and the complexity of lncRNAs, a panel of lncRNA biomarkers might be more precise and stable for predicting prognosis rather other one lncRNA. Therefore, it is timely to investigate the new lncRNA biomarkers by combining The Cancer Genome Atlas (TCGA) data with UCEC-specific data.
In the present study, we obtained the lncRNA expression profile and clinical information of UCEC patients from the datasets of TCGA project. By bioinformatic approaches, a potential 3-lncRNA signature was identified in EC, and the association between the signature and clinical characteristics was confirmed. Furthermore, clinical samples were used to demonstrate that 3-lncRNA signature has a much better performance than independent 3 lncRNAs, providing a new signature for early diagnosis and therapeutic strategy of EC.
Results
Identification of Differentially Expressed Long Noncoding RNAs Associated with Uterine Corpus Endometrial Carcinoma from The Cancer Genome Atlas
We obtained lncRNA expression profiles in 548 UCEC tissues and 35 normal tissues from TCGA datasets to screen DElncRs. To obtain reliable and stable results, lncRNA expression data were downloaded and performed using “DEseq,” “edgeR,” and “limma” R package separately in the R software (Figures 1A–C), and the intersections were acquired. Among the acquired lncRNAs, a set of 233 lncRNAs, including 93 upregulated and 140 downregulated, was abundantly expressed in all the uterine corpus endometrial carcinoma (Figures 1D–F, Supplementary Table S1). These results indicated the role of differentially expressed lncRNAs in the initiation and progression of uterine corpus endometrial carcinoma.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/814456/fcell-09-814456-HTML/image_m/fcell-09-814456-g001.jpg)
FIGURE 1. Identification of differentially expressed lncRNAs associated with UCEC from TCGA. (A-C) Volcano plot of DElncRs (different expressions of lncRNAs) between tumors and normal tissue by using “DEseq” (A), “edgeR” (B), and “limma” (C) R package. An absolute log2 fold change (FC) > mean (abs (logFC)) + 2*sd (abs (logFC) and an adjusted p-value of < 0.05 cutoff was used to defined DElncRs. (D) Venn diagram of upregulated lncRNAs. (E) Venn diagram of downregulated lncRNAs. (F) The heat map of DElncRs in TCGA–UCEC.
Validation of Prognostic Long Noncoding RNA Signature
Subsequently, univariate Cox regression analysis was conducted to estimate the prognostic relationship between DElncRs and EC patient OS, and 31 prognostic lncRNAs were obtained with a p < 0.05 (Figure 2A). Furthermore, to minimize prediction errors, 9 lncRNAs were screened out using the LASSO regression method. Kaplan–Meier survival curves were used to further analyze the relationship between the 9 lncRNAs and the OS of EC patients. Ultimately, 8 lncRNAs were identified to be related with OS (Figures 2B–I). Multivariable Cox regression analysis revealed the hazard ratios of 8 lncRNAs for OS of endometrium carcinoma (Figure 2J). The area under the ROC curve (AUC) for OS was 0.71 (Figure 2K). These results implied that the 8-lncRNA model could efficiently identify the risk of EC prognosis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/814456/fcell-09-814456-HTML/image_m/fcell-09-814456-g002.jpg)
FIGURE 2. Validation of prognostic lncRNA signature. (A) LASSO coefficient profiles of 31 lncRNAs. (B-I) Overall survivals of the 8 lncRNAs in TCGA–UCEC. Kaplan–Meier curves were plotted for overall survival in patients with either high or low levels of lncRNAs, classified according to median of each lncRNA. p values were calculated using log-rank test. Y-axis, overall survival probability; X-axis, time in months. (J) Hazard ratios of the signature for overall survival of 8 lncRNAs. (K) ROC curve analysis of 8-lncRNA model for survival prediction of UCEC.
Assessment of Prognostic Risk in Uterine Corpus Endometrial Carcinoma Patients Using a 3-Long Noncoding RNA Model
To better identify the prognostic signature for UCEC, 3 lncRNAs (hazard ratios for CTD-2377D24.6 = 1.229, RP4-616B8.5 = 1.407, and RP11-389G6.3 = 1.409) with the lowest p-value (p < 0.1) were picked out for further investigation (Figure 3A). Based on the coefficients of 3 prognostic lncRNAs from multivariate Cox regression analysis (Liu Y et al., 2019; Jiang et al., 2021), the risk score of the 3-lncRNA signature for OS was identified as CTD-2377D24.6 expression × 0.206 + RP4-616B8.5 × 0.341 + RP11-389G6.3 × 0.343. According to the median cutoff value of this prognostic signature, patients were divided into low-risk and high-risk sets. The survival results demonstrated that the low-risk set had a higher survival rate than that of the high-risk set (p < 0.0001, Figure 3B). To assess the potential prediction of 3-lncRNAs for overall survival of UCEC patients, the AUC analysis was performed to test the 3-lncRNA signature compared with each lncRNA. The results showed that the 3-lncRNA signature insignificantly showed an excellent performance than that of each lncRNA and two lncRNAs (Supplementary Figure S1, S2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/814456/fcell-09-814456-HTML/image_m/fcell-09-814456-g003.jpg)
FIGURE 3. Assessment of prognostic risk in UCEC patients using a 3-lncRNA model. (A) Hazard ratios of the signature for overall survival of 3 lncRNAs. (B) Kaplan–Meier plot of OS according to risk score of the 3-lncRNA signature in TCGA–UCEC. TCGA patients were divided into high risk (3-lncRNAs expressions higher than median) and low risk (3-lncRNA expressions lower than median), classified according to formula of signature, CTD-2377D24.6 expression × 0.206 + RP4-616B8.5 × 0.341 + RP11-389G6.3 × 0.343.
Correlation Between the 3-Long Noncoding RNA Signature and Clinical Characteristics of The Cancer Genome Atlas–Uterine Corpus Endometrial Carcinoma
To better understand the prognostic value of the 3-lncRNA signature, we further evaluated the relationships between the 3-lncRNA signature and traditional clinical characteristics. According to the median expressions of CTD-2377D24.6, RP4-616B8.5, RP11-389G6.3, and the 3-lncRNA signature risk score, the UCEC samples were divided into two sets. Pearson chi-square or Fisher’s exact tests revealed that the 3-lncRNA signature was closely linked with histological subtype (p < 0.0001), advanced clinical stage (p = 0.011), and clinical-grade (p < 0.0001) (Table 1). Compared to low-risk sets, the high-risk set tended to be serous adenocarcinoma (SAC), a histopathological type with worse differentiation and distant metastasis. These results demonstrated that 3-lncRNA signature was closely related to the conventional prognostic indicators.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/814456/fcell-09-814456-HTML/image_m/fcell-09-814456-t001.jpg)
TABLE 1. Relationships between 3-lncRNA set with histological subtype, clinical stage, and grade of EC patients.
Expressions of 3 Long Noncoding RNAs in Paracancerous and Tumor Tissues of Uterine Corpus Endometrial Carcinoma Patients
In addition, we validated the expressions of 3 lncRNAs in 30 paired paracancerous and tumor tissues of UCEC patients. First, the transcript abundances of RP4-616B8.5, RP11-389G6.3, and CTD-2377D24.6 were evaluated by qRT-PCR, and the results indicated that the expressions of RP11-389G6.3 and CTD-2377D24.6 were significantly higher in tumor tissues with p-values of 0.023 and 0.002, respectively, while the expression of RP4-616B8.5 did not show significant difference between tumor and paracancerous tissues with a p-value of 0.087 (Figure 4A), and the p-value of combined 3-lncRNAs was 0.027 using Hotelling T2 test (F = 3.56) (Figure 4B). Furthermore, in situ hybridization assay was also utilized to confirm the expressions of lncRNAs in paracancerous and tumor tissues of UCEC patients (Figure 4C). The staining scores of RP4-616B8.5, RP11-389G6.3, and CTD-2377D24.6 in EC tissues were significantly higher than those in paracancerous tissues with p-values of 0.042, 0.005, and 0.011, respectively (Figure 4D), and the p-value of combined 3-lncRNAs was 0.0002 using the Fisher’s methods (χ2 = 25.84) (Figure 4E). These results revealed that 3-lncRNA signature exhibited a better performance than the independent 3 lncRNAs for EC diagnosis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/814456/fcell-09-814456-HTML/image_m/fcell-09-814456-g004.jpg)
FIGURE 4. Validation of 3 lncRNAs expressions in clinical UCEC samples. (A) lncRNA RP4-616B8.5, RP11-389G6.3, and CTD-2377D24.6 mRNA expression in 30 paired paracanerous and tumor tissues of UCEC patients from the First Affiliated Hospital of Soochow University. (B) The combined effect of 3-lncRNAs using the Hotelling T2 test. (C) Representative images of in situ hybridization. (D) Staining scores of lncRNA RP4-616B8.5, RP11-389G6.3, and CTD-2377D24.6 levels in UCEC tissues and adjacent tissues by in situ hybridization. All data were presented as the mean ± SD. (E) The combined effect of 3-lncRNAs using Fisher’s method.
Functional Analysis of 3-Long Noncoding RNA Signature in Uterine Corpus Endometrial Carcinoma
To explore the potential roles of 3-lncRNA signature in UCEC, differentially expressed mRNAs (DeRNAs) between the high-risk (3-lncRNAs at high levels) and low-risk (3-lncRNAs at low levels) groups were identified in endometrium carcinoma patients (Supplementary Table S2), and KEGG and GO analysis were conducted. Functional enrichment analysis revealed that these DeRNAs were significantly enriched in 5 KEGG pathways, including carcinogenesis, drug metabolism and resistance, fluid shear stress, and steroid hormone biosynthesis (Figure 5A), 20 GO terms in biological processes, 10 GO terms in cellular components, and 10 GO terms in molecular functions (Figure 5B), indicating that drug metabolism, chemical carcinogenesis, and cell motility–related pathways might be involved. Given the importance of cortactin for invadopodia formation, cancer cell migration, and metastasis (Ji et al., 2020), we examined the mRNA expression levels and location of cortactin by qRT-PCR and immunohistochemical staining. The results showed that two cortactin-encoding genes, CTTN and HCLS1, were markedly increased in tumor tissues (Figures 5C,D). Interestingly, immunohistochemical staining revealed that cortactin exhibited in gland duct cells, but not in supporting cells (Figure 5I). In addition, the greatest differentially expressed mRNAs between the high-risk (3-lncRNA signature at high levels) and low-risk (3-lncRNA signature at low levels) groups were also determined. DNAH5, LTF, and Ezrin were significantly increased in tumor tissues (p < 0.05, Figures 5E–G), and WNT7A displayed a slight increase (p = 0.0669, Figure 5H). These results indicated that cortactin might be associated with the function of 3-lncRNA signature.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/814456/fcell-09-814456-HTML/image_m/fcell-09-814456-g005.jpg)
FIGURE 5. Functional analysis of 3-lncRNA signature. KEGG pathway (A) and Gene ontology (B) analysis (p < 0.05) of the 3-lncRNA signature related functional mRNAs. (C-H) qRT-PCR of indicated DeRNA expressions in 16 paired paracanerous and tumor tissues of UCEC patients from the First Affiliated Hospital of Soochow University. (I) Representative images of cortactin expression and location using immunohistochemical staining in UCEC tissues.
Discussion
Currently, a growing number of literatures demonstrated that dysregulated lncRNAs were involved in various diseases, as well as cancers (Evans et al., 2016). lncRNA might be a promising biomarker for cancer diagnosis, treatment, and prognosis prediction. Due to the heterogeneity of the tumor, a panel of lncRNA signature was more precise than a single lncRNA. In the present study, using Cox and LASSO regression, a 3-lncRNA signature was identified for predicting OS of EC patients. According to the median cutoff value of 3-lncRNA model, we demonstrated that the high-risk set displayed a poor survival, a higher clinical stage, and clinical grade and tended to be serous adenocarcinoma, a histopathological type with worse differentiation and distant metastasis. Our clinical samples also confirmed that 3-lncRNA, RP4-616B8.5, RP11-389G6.3, and CTD-2377D24.6 levels were increased in EC tissues than in paracancerous tissues by qRT-PCR and in situ hybridization. These findings provide an important hint that the 3-lncRNA signature has the potential performance for EC diagnosis and prognosis.
Mechanically, lncRNAs play a crucial role in EC progression by multiple patterns such as signaling, decoying, scaffolding, and guidance (Liu H et al., 2019). For example, Dong et al. reported that DLEU2 potentially interacted with miR-455 and miR-181a to promote epithelial-to-mesenchymal transition (EMT) and aerobic glycolysis (Dong et al., 2021). NEAT1 initiates a miR-361–mediated network to drive aggressive EC progression (Dong et al., 2019). lncRNA SOCS2-AS1 was found to suppress EC progression by promoting Aurora kinase A (AURKA) degradation via the ubiquitin–proteasome pathway (Jian et al., 2021). Additionally, Zhao et al. revealed that DLX6-AS1 promoted EC progression by recruiting p300/E2F1 in the DLX6 promoter region (Zhao and Xu, 2020). However, we found that little was known about the roles of RP4-616B8.5, RP11-389G6.3, and CTD-2377D24.6 in tumorigenesis and progression. One study reported that CTD-2377D24.6 was significantly induced by heat shock (Kim et al., 2020). Another study showed that CTD-2377D24.6 was a predictive factor in HCC patients with cirrhosis (Ye et al., 2020). No public study reported other two lncRNAs (RP4-616B8.5 and RP11-389G6.3). In the present study, we verified that the expressions of RP4-616B8.5, RP11-389G6.3, and CTD-2377D24.6 were higher in EC tissues than paracancerous tissues by qRT-PCR and in situ hybridization assays. According to the median cutoff value of 3-lncRNA signature, low-risk and high-risk sets were divided, and DeRNAs were identified. KEGG and GO analysis found that drug metabolism, chemical carcinogenesis, and cell motility–related pathway were enriched, indicating the potential roles of a panel of lncRNAs in initiation, metastasis, and chemoresistance of EC. The extensive quantity of published reports suggested that cell motility at an early stage in cancer correlated with metastasis (Lambert et al., 2017). In particular, the importance of cortactin for invadopodia formation, cancer cell migration, and metastasis has been proven (Schnoor et al., 2018; Ji et al., 2020). However, the link between lncRNAs and cortactin in endometrial carcinoma remains unclear. Here, we demonstrated that cortactin was markedly increased in UCEC tumor tissues, and especially exhibited in gland duct cells. In addition, the greatest differentially expressed mRNAs between the high-risk and low-risk groups, such as DNAH5, LTF, and Ezrin, were significantly increased in tumors. These hinted that more comprehensive studies about the molecular mechanism of 3-lncRNAs will remain to be lucubrated.
Traditional therapeutic strategies and risk stratification for EC patients are based on clinical and histological characteristics. However, the conventional classification does not adequately depict tumor biology owing to the high heterogeneity of EC. Recently, molecular or genomic classification has drawn attention as a promising approach to predict cancer prognosis. Levine et al. assessed the genome, transcriptome, and proteome of 373 endometrial carcinomas. Based on integrated genomic data, they were classified into four subgroups: POLE ultramutated, microsatellite instability hypermutated, copy-number low, and copy-number high. Subsequently, another molecular classification for EC termed “ProMisE” was identified (Talhouk et al., 2015). A similar integrated risk profile was established by the TransPORTEC international consortium (Stelloo et al., 2015; Stelloo et al., 2016). However, compared with genome sequencing, our 3-lncRNA signature was more suitable for clinical diagnosis and classification due to its higher stability and lower cost.
By bioinformatic approaches and verification of clinical samples, we demonstrated that the 3-lncRNA signature might be a reliable prognostic biomarker. However, there are several limitations in our study. First, the 3-lncRNA signature was constructed by the TCGA–UCEC datasets, in which the Caucasian race was the main patient. So, the prognostic value in other races is needed to be validated. Second, we detected the independent difference of 3 lncRNAs between the paracancerous tissues and UCEC tissues in our clinical samples, while the prognostic value of the signature was not analyzed due to insufficient prognostic data. Third, the 3 lncRNAs were rarely reported, and their potential function was unclear. Although functional enrichment analysis based on the DeRNAs in high- and low-risk signatures was performed, the potential mechanisms should be further experimentally investigated.
In conclusion, we revealed a potential 3-lncRNA signature that could accurately predict outcomes for UCEC patients. Meanwhile, we found that the 3-lncRNA signature was closely associated with clinical characteristics. Furthermore, we validated the different expressions of the 3 lncRNAs in our clinical samples, indicating that a panel of 3-lncRNAs exhibited better performance for EC diagnosis. These findings provide the first hint that the set of lncRNAs may play a critical role in the initiation and metastasis of EC, indicating a new signature for early diagnosis and therapeutic strategy of uterine corpus endometrial carcinoma.
Materials and Methods
Data Source and Differentially Expressed Long Noncoding RNA Screening
The lncRNA expression data and corresponding clinical information of UCEC patients were open-accessed from TCGA data portal (https://www.cancer.gov/tcga), including 548 UCEC tissues and 35 normal tissues. The differentially expressed lncRNAs (DElncRs) were identified by using the R package “DEseq,” “edgeR,” and “limma” with |Log2 fold-change (FC)| > mean|Log2FC| ± 2 × sd |Log2FC| and adjusted p < 0.05. A volcano plot of DElncRs was obtained by using the R package (Anders and Huber, 2010; Robinson et al., 2010; Ritchie et al., 2015).
Construction and Assessment of Long Noncoding RNA-Based Prognostic Signature
First, univariate COX regression and least absolute shrinkage and selection operator (LASSO) model were selected to find the independent prognostic lncRNAs. The LASSO method was performed by the package “glmnet” in the R software. Subsequently, multivariate Cox regression analysis was used to construct the prognostic signature. A risk score formula was constructed as follows: gene 1 × b1 + gene 2 × b2 + gene 3 × b3 +···gene n × bn, in which b represented the respective coefficient of genes. Gene represented the expression level of each gene. Subsequently, the risk score of prognostic signature formula was calculated as follows: CTD-2377D24.6 expression × 0.206 + RP4-616B8.5 × 0.341 + RP11-389G6.3 × 0.343. According to the median of risk score, the TCGA–UCEC patients were divided into the high-risk set and the low-risk set. To evaluate the prognostic signature of lncRNAs, the Kaplan–Meier and time-dependent receiver operating characteristic (ROC) curve analysis were performed.
Functional Enrichment Analysis
Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology analysis were performed with the “clusterProfiler” R package to identify the function of lncRNA-based signature (Yu et al., 2012). Significant functional categories were filtered into the meaning of p-value and false discovery rate (FDR) values <0.05.
RNA Extraction and Quantitative Real-Time Polymerase Chain Reaction
Total RNA was extracted using TRIzol (Life Technologies, NY, United States). The concentration and integrity of RNA were verified by a NanoDrop spectrophotometer (Thermo Fisher Scientific, Waltham, MA, United States). Afterward, the total RNA was reverse-transcribed into cDNA using the PrimeScript RT reagent kit (Takara, Dalian, China). The expressions of 3 lncRNAs were measured by qRT-PCR using the Hiff qPCR SYBR Green Master Mix (Yeasen Biotech Co., Shanghai, China) in the QuantStudio 6 system (Applied Biosystems, Waltham, MA). The primers synthesized are listed in Table 2.
Ethics Statement
Publicly available TCGA datasets were analyzed in this study, and approval from a local Ethics Committee was not necessary. For human subjects, all procedures were carried out according to Helsinki Declaration and institutional guidelines and were approved by the Ethics Committee at the First Affiliated Hospital of Soochow University.
In situ Hybridization Assay
The paraffin embedded UCEC and adjacent normal tissues were stained to detect the lncRNA expression. The lncRNA probes were designed and produced by SimaifuBio (Suzhou, Jiangsu, China). The probe sequences are presented in Table 3. In brief, sections were deparaffinized, digested, and blocked with 3% methanol-H2O2; then, the sections were dropped with prehybridization solution and incubated for 1 h in the incubator at 37°C. With the absorption of the excess liquid, the hybridization solution containing indicated lncRNA probes was added and then incubated in the incubator at 42°C overnight. Next day, after washing, the samples were dropped with block solution and incubated for 30 min at room temperature. After that, digoxigenin-labeled peroxidase antibody was added to incubate for 40 min in the incubator at 37°C. Afterward, the sections were added with DAB coloration, and the positive signal appeared brown–yellow. Hematoxylin staining solution was used to stain the nucleus. CaseViewer 2.2.1 (3DHISTECH Ltd.) and Image Pro Plus 6 were used for image capture and analysis, respectively.
Statistical Analysis
All of the expression profiles and clinical information were obtained from TCGA by R software. All statistical analyses were carried out using SPSS23.0 (SPSS, Chicago, IL, United States) or R software. For continuous variables, Student’s t-test was used to compare the difference between the two groups. For categorical variables, χ2 test was used to compare the differences among groups. Fisher’s method and Hotelling T2 test were used to combine p value. p < 0.05 was considered to be statistically significant.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics Statement
The studies involving human participants were reviewed and approved by The Ethics Committee at the First Affiliated Hospital of Soochow University. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
Conceptualization: BL and YC; methodology: HD, FJ, LD, WS, and JW; software: FJ and PW; validation: HD, LD, and JW; formal analysis: HD and LD; investigation: HD, FJ, and JL; resources: LD and JW; data curation: LD and JW; writing—original draft preparation: HD, FJ, LD, and JW; writing—review and editing: LC, JL, YP, YZ, and ZZ; visualization: PW and MJ; supervision: BL and YC; project administration: BL and YC; funding acquisition: BL, JL, YP, and FJ.
Funding
This study was partially supported by the National Natural Science Foundation of China (81872622, 81703209, and 31771417), the Natural Science Foundation of Jiangsu Province (BK20210737), and a sub-project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD) and Science and Technology Program of Suzhou (SKY2021035).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.814456/full#supplementary-material
References
Anders, S., and Huber, W. (2010). Differential Expression Analysis for Sequence Count Data. Genome Biol. 11 (10), R106. doi:10.1186/gb-2010-11-10-r106
Chen, W., Xia, C., Zheng, R., Zhou, M., Lin, C., Zeng, H., et al. (2019). Disparities by Province, Age, and Sex in Site-specific Cancer burden Attributable to 23 Potentially Modifiable Risk Factors in China: a Comparative Risk Assessment. Lancet Glob. Health. 7 (2), e257–e269. doi:10.1016/S2214-109X(18)30488-1
Dong, P., Xiong, Y., Konno, Y., Ihira, K., Kobayashi, N., Yue, J., et al. (2021). Long Non-Coding RNA DLEU2 Drives EMT and Glycolysis in Endometrial Cancer through HK2 by Competitively Binding with miR-455 and by Modulating the EZH2/miR-181a Pathway. J. Exp. Clin. Cancer Res. 40 (1), 216. doi:10.1186/s13046-021-02018-1
Dong, P., Xiong, Y., Yue, J., Xu, D., Ihira, K., Konno, Y., et al. (2019). Long Noncoding RNA NEAT1 Drives Aggressive Endometrial Cancer Progression via miR-361-Regulated Networks Involving STAT3 and Tumor Microenvironment-Related Genes. J. Exp. Clin. Cancer Res. 38 (1), 295. doi:10.1186/s13046-019-1306-9
Evans, J. R., Feng, F. Y., and Chinnaiyan, A. M. (2016). The Bright Side of Dark Matter: lncRNAs in Cancer. J. Clin. Invest. 126 (8), 2775–2782. doi:10.1172/JCI84421
Gao, L., Nie, X., Zhang, W., Gou, R., Hu, Y., Qi, Y., et al. (2020). Identification of Long Noncoding RNA RP11-89K21.1 and RP11-357H14.17 as Prognostic Signature of Endometrial Carcinoma via Integrated Bioinformatics Analysis. Cancer Cel Int. 20, 268. doi:10.1186/s12935-020-01359-9
He, X., Bao, W., Li, X., Chen, Z., Che, Q., Wang, H., et al. (2014). The Long Non-Coding RNA HOTAIR Is Upregulated in Endometrial Carcinoma and Correlates with Poor Prognosis. Int. J. Mol. Med. 33 (2), 325–332. doi:10.3892/ijmm.2013.1570
Ji, R., Zhu, X.-J., Wang, Z.-R., and Huang, L.-Q. (2020). Cortactin in Epithelial-Mesenchymal Transition. Front. Cel Dev. Biol. 8, 585619. doi:10.3389/fcell.2020.585619
Jian, F., Che, X., Zhang, J., Liu, C., Liu, G., Tang, Y., et al. (2021). The Long-Noncoding RNA SOCS2-AS1 Suppresses Endometrial Cancer Progression by Regulating AURKA Degradation. Cell Death Dis. 12 (4), 351. doi:10.1038/s41419-021-03595-x
Jiang, X., Xia, Y., Meng, H., Liu, Y., Cui, J., Huang, H., et al. (2021). Identification of a Nuclear Mitochondrial-Related Multi-Genes Signature to Predict the Prognosis of Bladder Cancer. Front. Oncol. 11, 746029. doi:10.3389/fonc.2021.746029
Kim, W. S., Ghassemi Nejad, J., Peng, D. Q., Jung, U. S., Kim, M. J., Jo, Y. H., et al. (2020). Identification of Heat Shock Protein Gene Expression in Hair Follicles as a Novel Indicator of Heat Stress in Beef Calves. Animal. 14 (7), 1502–1509. doi:10.1017/S1751731120000075
Lambert, A. W., Pattabiraman, D. R., and Weinberg, R. A. (2017). Emerging Biological Principles of Metastasis. Cell. 168 (4), 670–691. doi:10.1016/j.cell.2016.11.037
Levine, D. A., Kandoth, C., Schultz, N., Cherniack, A. D., Akbani, R., Liu, Y., et al. (2013). Integrated Genomic Characterization of Endometrial Carcinoma. Nature. 497 (7447), 67–73. doi:10.1038/nature12113
Li, B. L., and Wan, X. P. (2020). Prognostic Significance of Immune Landscape in Tumour Microenvironment of Endometrial Cancer. J. Cel Mol Med. 24 (14), 7767–7777. doi:10.1111/jcmm.15408
Li, W., Liu, J., and Zhao, H. (2020). Identification of a Nomogram Based on Long Non-Coding RNA to Improve Prognosis Prediction of Esophageal Squamous Cell Carcinoma. Aging. 12 (2), 1512–1526. doi:10.18632/aging.102697
Lim, L. J., Wong, S. Y. S., Huang, F., Lim, S., Chong, S. S., Ooi, L. L., et al. (2019). Roles and Regulation of Long Noncoding RNAs in Hepatocellular Carcinoma. Cancer Res. 79 (20), 5131–5139. doi:10.1158/0008-5472.CAN-19-0255
Liu, H., Wan, J., and Chu, J. (2019). Long Non-Coding RNAs and Endometrial Cancer. Biomed. Pharmacother. 119, 109396. doi:10.1016/j.biopha.2019.109396
Liu, Y., Liu, B., Jin, G., Zhang, J., Wang, X., Feng, Y., et al. (2019). An Integrated Three-Long Non-Coding RNA Signature Predicts Prognosis in Colorectal Cancer Patients. Front. Oncol. 9, 1269. doi:10.3389/fonc.2019.01269
Morice, P., Leary, A., Creutzberg, C., Abu-Rustum, N., and Darai, E. (2016). Endometrial Cancer. The Lancet. 387 (10023), 1094–1108. doi:10.1016/S0140-6736(15)00130-0
Murali, R., Soslow, R. A., and Weigelt, B. (2014). Classification of Endometrial Carcinoma: More Than Two Types. Lancet Oncol. 15 (7), e268–e278. doi:10.1016/S1470-2045(13)70591-6
Ouyang, D., Li, R., Li, Y., and Zhu, X. (2019). A 7‐lncRNA Signature Predict Prognosis of Uterine Corpus Endometrial Carcinoma. J. Cel Biochem. 120 (10), 18465–18477. doi:10.1002/jcb.29164
Pecorelli, S. (2009). Revised FIGO Staging for Carcinoma of the Vulva, Cervix, and Endometrium. Int. J. Gynecol. Obstet. 105 (2), 103–104. doi:10.1016/j.ijgo.2009.02.012
Peng, W.-X., Koirala, P., and Mo, Y.-Y. (2017). LncRNA-Mediated Regulation of Cell Signaling in Cancer. Oncogene. 36 (41), 5661–5667. doi:10.1038/onc.2017.184
Piulats, J. M., Guerra, E., Gil-Martín, M., Roman-Canal, B., Gatius, S., Sanz-Pamplona, R., et al. (2017). Molecular Approaches for Classifying Endometrial Carcinoma. Gynecol. Oncol. 145 (1), 200–207. doi:10.1016/j.ygyno.2016.12.015
Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). Limma Powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies. Nucleic Acids Res. 43 (7), e47. doi:10.1093/nar/gkv007
Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor Package for Differential Expression Analysis of Digital Gene Expression Data. Bioinformatics. 26 (1), 139–140. doi:10.1093/bioinformatics/btp616
Schmitt, A. M., and Chang, H. Y. (2016). Long Noncoding RNAs in Cancer Pathways. Cancer Cell. 29 (4), 452–463. doi:10.1016/j.ccell.2016.03.010
Schnoor, M., Stradal, T. E., and Rottner, K. (2018). Cortactin: Cell Functions of A Multifaceted Actin-Binding Protein. Trends Cel Biol. 28 (2), 79–98. doi:10.1016/j.tcb.2017.10.009
Siegel, R. L., Miller, K. D., and Jemal, A. (2020). Cancer Statistics, 2020. CA A. Cancer J. Clin. 70 (1), 7–30. doi:10.3322/caac.21590
Stelloo, E., Bosse, T., Nout, R. A., MacKay, H. J., Church, D. N., Nijman, H. W., et al. (2015). Refining Prognosis and Identifying Targetable Pathways for High-Risk Endometrial Cancer; a TransPORTEC Initiative. Mod. Pathol. 28 (6), 836–844. doi:10.1038/modpathol.2015.43
Stelloo, E., Nout, R. A., Osse, E. M., Jürgenliemk-Schulz, I. J., Jobsen, J. J., Lutgens, L. C., et al. (2016). Improved Risk Assessment by Integrating Molecular and Clinicopathological Factors in Early-Stage Endometrial Cancer-Combined Analysis of the PORTEC Cohorts. Clin. Cancer Res. 22 (16), 4215–4224. doi:10.1158/1078-0432.CCR-15-2878
Sun, Y., Gao, X., Li, P., Song, L., and Shi, L. (2021). LncRNA ZFAS1, as a Poor Prognostic Indicator, Promotes Cell Proliferation and Epithelial-Mesenchymal Transition in Endometrial Carcinoma. Personalized Med. 18 (1), 43–53. doi:10.2217/pme-2020-0014
Talhouk, A., McConechy, M. K., Leung, S., Li-Chang, H. H., Kwon, J. S., Melnyk, N., et al. (2015). A Clinically Applicable Molecular-Based Classification for Endometrial Cancers. Br. J. Cancer. 113 (2), 299–310. doi:10.1038/bjc.2015.190
Wang, X., Dai, C., Ye, M., Wang, J., Lin, W., and Li, R. (2021). Prognostic Value of an Autophagy-Related Long-Noncoding-RNA Signature for Endometrial Cancer. Aging. 13 (4), 5104–5119. doi:10.18632/aging.202431
Xie, S., Chang, Y., Jin, H., Yang, F., Xu, Y., Yan, X., et al. (2020). Non-Coding RNAs in Gastric Cancer. Cancer Lett. 493, 55–70. doi:10.1016/j.canlet.2020.06.022
Yang, J.-Y., Werner, H. M. J., Li, J., Westin, S. N., Lu, Y., Halle, M. K., et al. (2016). Integrative Protein-Based Prognostic Model for Early-Stage Endometrioid Endometrial Cancer. Clin. Cancer Res. 22 (2), 513–523. doi:10.1158/1078-0432.CCR-15-0104
Yang, Z., Jiang, S., Shang, J., Jiang, Y., Dai, Y., Xu, B., et al. (2019). LncRNA: Shedding Light on Mechanisms and Opportunities in Fibrosis and Aging. Ageing Res. Rev. 52, 17–31. doi:10.1016/j.arr.2019.04.001
Ye, J., Li, H., Wei, J., Luo, Y., Liu, H., Zhang, J., et al. (2020). Risk Scoring System Based on lncRNA Expression for Predicting Survival in Hepatocellular Carcinoma with Cirrhosis. Asian Pac. J. Cancer Prev. 21 (6), 1787–1795. doi:10.31557/APJCP.2020.21.6.1787
Yu, G., Wang, L.-G., Han, Y., and He, Q.-Y. (2012). clusterProfiler: an R Package for Comparing Biological Themes Among Gene Clusters. OMICS: A J. Integr. Biol. 16 (5), 284–287. doi:10.1089/omi.2011.0118
Keywords: endometrial carcinoma, lncRNA signature, prognosis, diagnostic, biomarker
Citation: Ding H, Jiang F, Deng L, Wang J, Wang P, Ji M, Li J, Shi W, Pei Y, Li J, Zhang Y, Zhang Z, Chen Y and Li B (2022) Prediction of Clinical Outcome in Endometrial Carcinoma Based on a 3-lncRNA Signature. Front. Cell Dev. Biol. 9:814456. doi: 10.3389/fcell.2021.814456
Received: 13 November 2021; Accepted: 31 December 2021;
Published: 01 February 2022.
Edited by:
Feng Wang, Affiliated Hospital of Nantong University, ChinaReviewed by:
Yongye Huang, Northeastern University, ChinaJinpeng He, Institute of Modern Physics (CAS), China
Dong Pan, Duke University, United States
Chunlin Shao, Fudan University, China
Copyright © 2022 Ding, Jiang, Deng, Wang, Wang, Ji, Li, Shi, Pei, Li, Zhang, Zhang, Chen and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Youguo Chen, chenyouguo@suda.edu.cn; Bingyan Li, bingyanli@suda.edu.cn
†These authors have contributed equally to this work