- 1Department of Urology, The Second People’s Hospital of Wuhu, Wuhu, China
- 2School of Pediatrics, Nanjing Medical University, Nanjing, China
- 3The First Clinical Medical College of Nanjing Medical University, Nanjing, China
Purpose: Epigenetic modifications play a crucial role in cancer development, and our study utilized public data to analyze which leads to the discovery of significant epigenetic abnormalities in lncRNAs, offering valuable insights into prognosis and treatment strategies for renal carcinoma.
Methods: Public data were obtained from the Cancer Genome Atlas (TCGA), International Cancer Genome Consortium (ICGC) and Gene Expression Omnibus (GEO) database. The analysis of the online public data was all completed in R software.
Results: We discovered a great number of epigenetic abnormalities of lncRNA in renal cancer, which is achieved by comparing the following modification and methylation of histone region changes on the promoter and enhancer of lncRNA: H3K27ac, H3K4me1, H3K4me3. As a result, 12 specific epigenetic disorders of lncRNA genes in renal cancer were identified. Finally, based on this lncRNA, we investigated the prognosis of renal cancer samples, among which 8 lncRNA can be seen as markers of prognosis in renal cancer, which had great prediction ability for ccRCC prognosis. Meanwhile, high risk score may pose response better to axitinib and nilotinib, but not sorafenib or sunitinib. Beyond, we observed an elevated level of risk score in immunotherapy non-responders. Further, biological enrichment and immuno-infiltration analysis was conducted to investigate the fundamental differences between patients categorized as high or low risk.
Conclusion: Our research improves the understanding in the function of epigenetic dysregulated long non-coding RNAs in renal carcinoma.
1 Introduction
Renal cell carcinoma is predominantly represented by clear cell renal cell carcinoma (ccRCC), constituting over three-quarter of all reported cases, making it the most significant histological subtype among renal carcinomas (Capitanio and Montorsi, 2016). Early-stage patients with renal cancer who undergo surgery might benefit from long-term therapeutic outcomes (Campbell et al., 2021). Nonetheless, for those with advanced stage or distant metastasis, treatment options are often limited. Consequently, identifying novel targets that have the potential to serve as ones of guiding renal carcinoma therapy is meaningful.
Epigenetic modification is a common physiological process in organisms, like DNA methylation, genomic imprinting, gene silencing, and so on (Miranda Furtado et al., 2019). Generally, normal and orderly epigenetics help maintain normal physiological functions of cells and the body. However, there is increasing evidence that epigenetic regulation abnormalities are closely linked to cancer occurrence and development (Sun et al., 2022). Becker et al. (2020) found that epigenetic reprogramming can maintain the glycolytic phenotype of breast carcinoma cells by regulating HIF-1α and glycolytic enzymes, further enhancing cancer progression. Hamamoto et al. (2004) revealed that SMYD3 was overexpressed in colon and liver cancer, and can facilitate cancer development based on its histone methyltransferase activity. Meanwhile, Ulanovskaya et al. (2013) indicated that the nicotinamide N-methyltransferase NNMT can affect tumor metabolism by impairing the methylation potential of cancer cells. Long non-coding RNAs (lncRNAs) represent a category of noncoding RNAs that exceed 200 nucleotides in length, famous for their broad regulatory effect (Bhan et al., 2017). Moreover, lncRNAs are widely known of being involved in cancer epigenetics. For example, in colonic carcinoma, Ni et al. (2019) noticed that lncRNA GAS5 has direct interaction with WW, subsequently facilitating the phosphorylation and ubiquitin-mediated degradation of YAP. Xiu et al. (2019) demonstrated that the lncRNA LINC02273 can form a complex with hnRNPL. This complex can upregulate AGR2 transcription by influencing epigenetic modification of the AGR2 promoter region and therefore promote cancer development. Therefore, exploring the lncRNAs involved in epigenetic modification is important.
In this research, we identified the lncRNAs involved in epigenetic modification in ccRCC and comprehensively explored their role. Moreover, based on the identified epi-lncRNAs, a prognosis model was established with good prediction ability of patients’ prognosis. In addition, drug-susceptibility analysis illustrated that patients who got high-risk scores may have a more positive response to axitinib and nilotinib than to sorafenib or sunitinib. Meanwhile, we noticed a higher risk score among non-responders with immunotherapy. Further, biological enrichment and immuno-infiltration analysis was published to investigate the fundamental differences between patients categorized as high or low risk. Our result can provide novel insight into the identified epi-lncRNAs.
2 Materials and methods
2.1 Public data collection
In the research, the gene expression profile of ccRCC and its corresponding normal samples were downloaded from International Cancer Genome Consortium (ICGC) and Cancer Genome Atlas (TCGA) database, including FPKM, count number expression profile data and clinical information. Then we converted FPKM to TPM based on the author’s R code. GENCODE’s gene annotation file (version 24) was used to differentiate between long noncoding RNAs and protein-coding genes (PCGs). We also transformed the Ensembl ID of the chosen genes into symbol form. All expression profiles have been processed before analysis. Also, 450k microarray chip data were downloaded from the TCGA-KIRC project. Besides, all the missing values in renal cancer data were completed with the KNN method. Based on the Gene Expression Omnibus (GEO) database, the GSE86091 dataset with para-cancerous and tumor samples was downloaded, including h3k4me1, h3k4me3, h3k27ac histone information (Yao et al., 2017).
2.2 Identification of lncRNAs and PCGs associated with the epigenetic disorder
Then, the lncRNAs and PCGs which expressed differentially were identified by R package deseq2 with p-value determined by the Benjamin Hochberg method. LncRNAs or PCGs with FDR < 0.05 and |logFC| > 1 were considered significant. Next, we selected the unique peak in RCC in accordance with the physical location of the histone modified peak. Only the peak with a p-value < 0.05 was retained as the difference peak. Furthermore, we combined with the GTF file of GENCODE to get differentially histone-modified genes as well as identity of gene enhancers from the human enhancer database FANTOM5. The promoter of the gene is defined as 2 kb upstream and 0.5 kb downstream of the transcription initiation site (TSS), which can be identified by the R package ChIPseeker and the human lincRNAsTranscripts database. Finally, Using R-packet CHAMP, the Bumphunter method is used to identify DMR and the regions with Bumphunter DMR. p-values < 0.01 are considered significant DMRs. In the end, we made a definition that lncRNAs meet the following conditions as epiRNAs: i) there was differential expression between tumor and normal samples; ii) promoter or enhancer has more than one histone modified or methylated region.
2.3 Genomic characterization
To figure out the genomic characteristics of epigenetic dysregulation and none dysregulation lncRNA/PCG, we made comparations among the number and length of exons, transcripts, genes of epi-lncRNA, non-epi-lncRNA, epi-PCG, non-epi-PCG.
2.4 Genome mapping
To investigate the epigenetic features of lncRNA caused by histone modification, we analyzed the distribution characteristics of promoter and enhancer among epi-lncRNA modified through different histones in genome level.
2.5 Identification of LncRNA with renal carcinoma-specific candidate epigenetic characteristic
To learn more about the relation between epigenetic disorders lncRNA and carcinoma, we made an analysis of epi-lncRNA enrichment and cancer lncRNAs with lnc2cancer v3.0. Besides, epi-lncRNA related to ccRCC was screened out for the sake of exploring the regulation of epi-lncRNA. In addition, considering that most of the genes that affect the process of the diseases tend to show expression disorder, we identified the epi-lncRNA in ccRCC through the TCGA TPM expression profile.
2.6 Enrichment analysis
Pathway enrichment and quantification of specific pathways activity were completed by Gene Set Enrichment Analysis (GSEA) and Single Sample GSEA (ssGSEA) algorithm, including the Kyoto Encyclopedia of Genes and Genomes (KEGG) and epigenetically dysregulated lncRNAs score (Subramanian et al., 2005; Pang Y et al., 2022).
2.7 Prognosis analysis
The Kaplan-Meier method was applied to estimate the overall survival (OS) rate and survival time, and based on the method of log-rank test, the differences among different sample subgroups were obtained. Finally, the potential prognostic markers of epi-lncRNA were identified. A risk model was built on basis of multivariate Cox regression analysis.
2.8 Immune infiltration analysis
The CIBERSORT algorithm (version 1.03) and the LM22 signature matrix were used to analyze immune infiltration according to the expression profile (Chen et al., 2018).
2.9 Drug sensitivity and immunotherapy analysis
The Genomics of Drug Sensitivity in Cancer Database (GDSC) was used to analyze drug susceptibility (Yang et al., 2013). The evaluation of patients taken immunotherapy was made by application of the Tumor Immune Dysfunction and Exclusion (TIDE) algorithm (Fu et al., 2020).
2.10 Statistical analysis
All analysis was performed with the programs of R software (version 4.2.1). Statistical significance was determined using two-tailed p-value and less than 0.05 was considered statistically significant.
2.11 Clinical samples and tissue
This study was approved by the Institutional Review Board of the First Affiliated Hospital of Nanjing Medical University. All patients participated in the study voluntarily donated samples for clinical research and signed informed consents before surgery. We collected samples from patients with renal cell carcinoma from March 2024 to June 2024. All the samples were from patients took partial nephrectomy in the Department of Urology Department of the First Affiliated Hospital of Nanjing Medical University. Samples were collected under the guidance of one experienced pathologist, and immediately stored in a liquid nitrogen tank.
2.12 RNA preparation and quantitative real-time PCR (qRT-PCR)
TRIZOL reagent (Invitrogen) was used to extract total RNA in both tissues and cells according to the manufacturer’s instruction. The cDNA was obtained from the purified RNA using a PrimeScript RT Reagent Kit (Takara). SYBR Premix Ex Taq II Kit (Takara) was used for qRT-PCR assays. Results were normalized to β-actin expression. All experiments were performed at least three times. All specific primers used in this study have been listed in Supplementary Table S1.
3 Results
3.1 Identification of the epigenetically dysregulated lncRNAs
Genes that are differentially expressed were identified using DESeq2, including 5906 PCGs and 5370 lncRNAs. The combination of histone modifications and 450K methylation microarray data enabled us to identify 149 epiRNAs, 13426 non-epiRNAs, 2241 epi-PCGs, and 17200 non-epi-PCGs. Moreover, lncRNAs in ccRCC are significantly less frequently than in PCG (Figure 1A). A comparison of epi-lncRNA transcripts and exons with non-epi-lncRNA transcripts, epi-PCG transcripts, and non-epi-PCG transcripts was followed to represent the genomic characteristics of epigenetically dysregulated lncRNAs. Epi-lncRNAs appear with a long transcript length and large numbers of transcripts (Figures 1B, C) and epi-PCG showed similar structural characteristics as epi-lncRNAs. Simultaneously, epi-lncRNAs were presented with more exons and epi-PCG owns more exons and length (Figures 1D, E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1406150/fgene-15-1406150-HTML/image_m/fgene-15-1406150-g001.jpg)
Figure 1. Identification of the epi-lncRNAs in ccRCC. Notes: (A) The proportion of epi-lncRNAs and epi-PCGs to total lncRNAs and PCGs in the genome; (B–E) Genomic characteristics comparison between epigenetic dysregulated lncRNA/PCGs and non-epigenetic dysregulated lncRNA/PCGs.
3.2 Genomic landscape of epigenetically dysregulated lncRNAs
Systematic analysis of epigenetic disorders of epi-lncRNA in renal cell carcinoma revealed the landscape of epi-lncRNAs resulted from various histone modifications and different methylation regions (Figure 2A). Besides, in these lncRNAs, we found several abnormal histone modifications including H3K4me3, H3K4me1, and H3K27ac. Also, most of the lncRNAs with apparent disorders coexist with various histone modifications. Moreover, we also found that the abnormal histone modifications were mostly located in the promoter region (Figure 2B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1406150/fgene-15-1406150-HTML/image_m/fgene-15-1406150-g002.jpg)
Figure 2. Genomic landscape of epi-lncRNAs. Notes: (A) Differential methylation regions and histone modifications lead to a distinct epi-lncRNA genomic landscape; (B) Type distribution in epi-lncRNAs with several apparent disorders.
3.3 Biological enrichment of dysregulated epi-lncRNAs
Next, we completed an analysis of the relation between the expression of epi-lncRNA and the pathway in ccRCC. The ssGSEA was applied to calculate each sample enrichment score in all the lncRNAs affected by histone modifications. We observed that the scores in cancer samples were remarkably higher than adjacent samples, suggesting that these epi-lncRNAs may own a great possibility of being carcinogenic (Figure 3A). Further, we analyzed the relationship between these enrichment scores and OS and found that high scores were associated with poor prognosis (Figure 3B). Additionally, we observed that there are 18 pathways in the KEGG pathway related to six sorts of epi-lncRNA scores closely, suggesting that there is a certain consistency among various epi-lncRNA-related pathways. Of these 18 pathways, the tumor-related pathway of BLADDER_CANCER and the immune-related pathways of PRIMARY_IMMUNODEFICIENCY, CYTOKINE_CYTOKINE_RECEPTOR_INTERACTION, NATURAL_KILLER_CELL_MEDIATED_CYTOTOXICITY and AUTOIMMUNE_THYROID_DISEASE have Positive correlations to the epi-lncRNAs score. On the contrary, the metabolism-related pathways of SELENOAMINO_ACID_METABOLISM, CITRATE_CYCLE_TCA_CYCLE, PYRUVATE_METABOLISM and PROPANOATE_METABOLISM are negatively correlated with the epi-lncRNAs score (Figure 3C). These results reveal that epi-lncRNA plays an important role in tumor development and immunity.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1406150/fgene-15-1406150-HTML/image_m/fgene-15-1406150-g003.jpg)
Figure 3. Functional analysis of epi-lncRNA. Notes: (A) Differences between cancer and adjacent cancers in six types of apparent dysregulated lncRNAs; (B) Relationship between different histone modifications of six epi-lncRNAs and overall survival; (C) Six kinds of apparent dysregulated lncRNAs are most related to KEGG pathway.
3.4 Relationship between epigenetically dysregulated lncRNA and RNA modification
We calculated the association of the enrichment scores of six epi-lncRNAs and m6a, m5c and m1A genes, indicating that they are significantly related. The enhancer and promoter of H3K4me3, H3K4me1 and H3K27ac not only have a similar correlation with these genes but also have their unique correlation (Figure 4A). The results indicate that there might exist unknown regulatory modes of epi-lncRNA resulted from enhancer and promoter histone modification, with a tight link between epi-lncRNA and RNA modification.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1406150/fgene-15-1406150-HTML/image_m/fgene-15-1406150-g004.jpg)
Figure 4. Epi-lncRNAs in ccRCC. Notes: (A) Correlation between m6a, m5c, and m1A genes and enrichment scores of six epi-lncRNAs; (B) Intersection of disease-related lncRNA and renal cell carcinoma epi-lncRNA; (C) Specific epi-lncRNA expression associated with renal cell carcinoma is different between normal samples and tumors.
3.5 Identification of renal cancer-specific epigenetic dysregulated lncRNA
Using Lnc2Cancer V3.0, it can be discovered that 51 epi-lncRNAs were associated with normal human cancers (Figure 4B). To explore how candidate lncRNAs play the regulatory role in ccRCC, the list of ccRCC-related lncRNAs was collected from Lnc2Cancer v3.0. Of 51 epi-lncRNAs, 12 lncRNAs were identified that have direct relationships with ccRCC. We subsequently calculated the differences of 12 epi-lncRNAs between ccRCC samples and normal samples, respectively (Figure 4C). Of that, 11 epi-lncRNAs are presented with significant differences, except H19. Except that the expression of HOXA11-AS and LINC00461 in normal samples is higher than that in tumor samples, the other 9 lncRNAs are highly expressed in tumor samples. Among these epi-lncRNAs, HOTAIR, HOTTIP, XIST, H19, HCP5, and LINC00160 were only related to histone-modified promoters; LINC00461 is only related to DMR enhancers; PVT1 is associated with histone modified enhancer and promoter (Figure 5A). Then, we collected a clinical ccRCC cohrt (normal group: 8, disease group: 8) and performed PCR validation on the relevant lncRNAs, and found significant differences between the disease group and the normal group (Figure 5B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1406150/fgene-15-1406150-HTML/image_m/fgene-15-1406150-g005.jpg)
Figure 5. Epigenetic modification profiles of identified epi-lncRNAs. Notes: (A) Epigenetic modification profiles of PVT1, HOTAIR, LINC00160, HOTTIP, H19, LINC00461, HCP5 and XIST, (B) quantitative real-time PCR validation.
3.6 Epi-lncRNAs are associated with the prognosis of ccRCC
Considering the potential prognostic value of epi-lncRNA is unclear yet, an analysis of the 12 epi-lncRNAs survival mentioned above based on the sample in TCGA-KIRC was made. Results indicated that 8 epi-lncRNAs are significantly associated with survival based on the univariate Cox regression analysis (Figure 6A). Kaplan-Meier survival curves illustrated that the epi-lncRNAs HOTAIR, HOTTIP, HOXA11-AS1, LINC00461, LUCAT1, PVT1, XIST, H19, LINC00160, MIAT, GAS6-AS1 were correlated with poor prognosis, while HCP5 was correlated with better prognosis (Figures 6B–M).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1406150/fgene-15-1406150-HTML/image_m/fgene-15-1406150-g006.jpg)
Figure 6. Prognosis value of identified epi-lncRNAs. Notes: (A) Univariate analysis of 12 kidney cancer-associated epi-lncRNA genes; (B–M) Kaplan-Meier survival curve for samples with high and low expression of 12 renal cell carcinoma-related EPI lncRNA genes.
3.7 Establishment of prognosis signature consisting of eight epi-lncRNAs
Multivariate survival analysis was used to establish the 8-EpiLncRNA signature on basis of the TCGA dataset, in which patients with high-risk score had higher mortality (Figure 7A). Receiver operating characteristic curve (ROC) analysis shows that the designed signature owned a remarkable prediction ability of patient OS (Figure 7B, 1-year AUC of 0.75, 3-year AUC of 0.67, and 5-year AUC of 0.72). Similarly, Kaplan-Meier survival curves indicated that the high-risk score patients have a worse OS performance tendency (Figure 7C). Then, we selected the RECA-EU with prognostic information from the ICGC database for model validation. In the ICGC cohort, we also observed satisfactory performance (Figure 7D). The ROC result of 1-, 3- and 5- years were 0.75, 0.68, and 0.68, respectively (Figure 7E). Kaplan-Meier survival curves indicated that the high-risk score ones also incline to have poor prognosis (Figure 7F). Univariate and multivariate analysis suggested that the new model we built can be considered as an independent risk factor of other clinical features (Supplementary Figures S1A, B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1406150/fgene-15-1406150-HTML/image_m/fgene-15-1406150-g007.jpg)
Figure 7. Prognosis model based on the epi-lncRNAs. Notes: (A) Overview of our model in TCGA cohort; (B) Kaplan-Meier survival curve of high and low risk patients in TCGA cohort; (C) ROC curves of our model in TCGA cohort; (D) Overview of our model in ICGC cohort; (E) Kaplan-Meier survival curve for ICGC cohort patients at high and low risk; (F) ROC curves of our model in ICGC cohort.
3.8 Drug sensitivity and immunotherapy analysis
Next, we investigated the fundamental distinction in Supplementary Figures S1A, B drug sensitivity between patients at high and low risk. Results indicated patients with high-risk score might respond better to axitinib and nilotinib, but not sorafenib and sunitinib (Figure 8). Moreover, A positive association was found between the TIDE score and the risk score. (Supplementary Figure S2A, cor = 0.274, p < 0.001). Meanwhile, it was worth noted that immunotherapy non-responders had a much higher score (Supplementary Figure S2B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1406150/fgene-15-1406150-HTML/image_m/fgene-15-1406150-g008.jpg)
Figure 8. Drug sensitivity analysis. Notes: (A–L) Drug sensitivity differences between people at high and low risk for specific target drugs.
3.9 Pathway enrichment and immune infiltration
Based on the Hallmark gene collection, GSEA analysis revealed that the pathways of epithelial-mesenchymal transition, P53 pathway, UV response up, E2F targets, myogenesis, pancreas beta cells, and G2M checkpoint were triggered in patients at high risk (Supplementary Figure S3). Analysis of immune infiltration revealed that there’s a chance that the high-risk patients have more plasma cell infiltration, activated CD4+ memory T cells, M0 macrophages, follicular helper T cells, and regulatory T cells (Tregs), while a lower level of monocytes, resting mast cells and eosinophils (Supplementary Figure S4).
4 Discussion
Among the most dangerous tumors of the urinary system, renal cancer poses a long-term medical and health burden for patients (Turajlic et al., 2018). Radiotherapy and chemotherapy are ineffective in treating patients with advanced ccRCC patients, therefore limiting their treatment options to some extent (Turajlic et al., 2018). Nowadays, researchers had focused on exploring the target affecting the response of ccRCC patients to specific target drugs. Wei et al. (2021) suggested that MX2 could be a prognosis marker of ccRCC and can affect sunitinib sensibility.
The regulatory effect of lncRNAs has made them possible targets for intervention in cancer. As well as participating in many important physiological processes, lncRNA also contributes to cell remodeling. Tan et al. (2021) found that the posttranslational modifications medicated by lncRNAs can affect the metabolism reprogramming of cancer. Numerous studies have identified dysregulated lncRNAs in RCC tissues compared to normal kidney tissues. These lncRNAs participate in diverse biological processes such as cell proliferation, apoptosis, invasion, and metastasis by modulating chromatin remodeling, transcriptional regulation, and post-transcriptional gene expression. For example, Liu et al. (2022) found that AC010973.2 can affect tumor stemness and is associated with ccRCC patients’ prognosis. Qu et al. (2016) indicated that the exosome-transmitted lncRNA lncARSR can upregulate the level of c-MET and AXL expression in RCC cells through a competitive endogenous RNA (ceRNA) mechanism, partly responsible for the cancer progression and sunitinib resistance. In terms of epigenetic modification, lncRNAs exert their regulatory functions through interactions with epigenetic modifiers including DNA methyltransferases, histone methyltransferases, and chromatin remodeling complexes. They can act as scaffolds, guides, or decoys to target specific genomic loci and regulate gene expression in cis or trans (Gu et al., 2021; Hu et al., 2021; Yang et al., 2021). In conclusion, while significant strides have been made in elucidating the epigenetic landscape of RCC, challenges such as tumor heterogeneity and technical limitations necessitate ongoing research efforts. Future directions should focus on overcoming these challenges through innovative research methodologies and collaborative multiomics and single-cell approaches, ultimately translating scientific discoveries into clinical benefits for RCC patients.
In our research, we discovered the lncRNAs involved in the epigenetic modification in ccRCC. A total of 51 lncRNAs were identified associated with the enriched score of promoter or enhancer of H3K27ac, H3K4me1 and H3K4me3. Among these lncRNAs, 12 epi-lncRNAs were selected for further analysis for their role in ccRCC previously reported, including HOTAIR, HOTTIP, HOXA11-AS, LINC00461, LUCAT1, PVT1, XIST, H19, HCP5, LINC00160, MIAT and GAS6-AS1. In renal cancer, Bai et al. (2021) found that HOTAIR can coordinate with the androgen receptor to enhance the transcription of GLI2, further promoting tumor stemness and angiogenesis. Wang et al. (2018) uncovered that the lncRNA HOTTIP may accelerate the development of renal cancer through a ceRNA mechanism (miR-615/IGF-2 axis). According to reports, all these lncRNAs play a role in renal cancer development. Our result indicated that these lncRNAs were associated with abnormal epigenetic modifications, which can broaden the research direction of these lncRNAs.
Also, based on the identified epi-lncRNAs, a prognosis model was established with good prediction ability of patients’ prognosis. Moreover, the results of drug sensitivity analysis showed that Axitinib and nilotinib may be far more effective for patients with high-risk scores, but not sunitinib and sorafenib. When treating advanced renal cancer, sunitinib is a highly potent targeted drug (Powles et al., 2020). Researchers have focused on the factors affecting the target drug of ccRCC. Xiong et al. (2021) demonstrated how RRM2 can influence PD-1 blockade and sunitinib sensitivity by promoting the AKT pathway and stabilizing ANXA1. Based on our results, assessing the model molecules’ relative expression level might reveal the patient risk group, further guiding the therapy option for renal cancer.
Some shortcomings also need to be addressed. Initially, the patients we registered were mainly Western people. The race bias brought by this can reduce the effectiveness of our findings. Also, the sample counts in the validation cohort ICGC are less than in the TCGA cohort. In the future, advances in high-throughput sequencing technologies and computational algorithms will facilitate integrated analyses of lncRNA expression profiles, epigenetic modifications, and genomic alterations in urologic tumors (Pang et al., 2020; Pang K et al., 2022). This holistic approach promises to unravel complex regulatory networks and identify actionable targets for precision medicine. Moreover, Continued efforts in RNA biology and nanomedicine hold promise for developing innovative RNA-targeted therapies, including lncRNA-specific inhibitors or mimics, for personalized treatment strategies in RCC.
Data availability statement
All the data used for this study was open-accessed through the following link: https://portal.gdc.cancer.gov/ (TCGA databse); https://icgc.org/ (ICGC databse); https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE86091 (GSE86091). All original data can be obtained from the corresponding author based on reasonable requirements.
Ethics statement
The studies involving humans were approved by the Institutional Review Board of the First Affiliated Hospital of Nanjing Medical University. The studies were conducted in accordance with the local legislation and institutional requirements. The human samples used in this study were acquired from primarily isolated as part of your previous study for which ethical approval was obtained. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.
Author contributions
JW: Writing–original draft. PD: Writing–original draft, Writing–review and editing. YS: Writing–review and editing. GW: Investigation, Software, Writing–review and editing. HL: Supervision, Writing–review and editing. LT: Conceptualization, Supervision, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was supported by the 2022 Wuhu Science and Technology Plan Project (Project ID: 2022cg30).
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/fgene.2024.1406150/full#supplementary-material
References
Bai, J. Y., Jin, B., Ma, J. B., Liu, T. J., Yang, C., Chong, Y., et al. (2021). Hotair and androgen receptor synergistically increase Gli2 transcription to promote tumor angiogenesis and cancer stemness in renal cell carcinoma. Cancer Lett. 498, 70–79. Epub 2020/11/07. doi:10.1016/j.canlet.2020.10.031
Becker, L. M., O'Connell, J. T., Vo, A. P., Cain, M. P., Tampe, D., Bizarro, L., et al. (2020). Epigenetic reprogramming of cancer-associated fibroblasts deregulates glucose metabolism and facilitates progression of breast cancer. Cell. Rep. 31 (9), 107701. doi:10.1016/j.celrep.2020.107701
Bhan, A., Soleimani, M., and Mandal, S. S. (2017). Long noncoding RNA and cancer: a new paradigm. Cancer Res. 77 (15), 3965–3981. doi:10.1158/0008-5472.CAN-16-2634
Campbell, S. C., Clark, P. E., Chang, S. S., Karam, J. A., Souter, L., and Uzzo, R. G. (2021). Renal mass and localized renal cancer: evaluation, management, and follow-up: AUA guideline: Part I. J. Urol. 206 (2), 199–208. doi:10.1097/JU.0000000000001911
Capitanio, U., and Montorsi, F. (2016). Renal cancer. Ren. cancer. Lancet. 387 (10021), 894–906. doi:10.1016/S0140-6736(15)00046-X
Chen, B., Khodadoust, M. S., Liu, C. L., Newman, A. M., and Alizadeh, A. A. (2018). Profiling tumor infiltrating immune cells with cibersort. Methods Mol. Biol. Clift. NJ 1711, 243–259. Epub 2018/01/19. doi:10.1007/978-1-4939-7493-1_12
Fu, J., Li, K., Zhang, W., Wan, C., Zhang, J., Jiang, P., et al. (2020). Large-scale public data reuse to model immunotherapy response and resistance. Genome Med. 12 (1), 21. Epub 2020/02/28. doi:10.1186/s13073-020-0721-z
Gu, Y., Niu, S., Wang, Y., Duan, L., Pan, Y., Tong, Z., et al. (2021). Dmdrmr-mediated regulation of M(6)a-modified cdk4 by M(6)a reader Igf2bp3 drives ccrcc progression. Cancer Res. 81 (4), 923–934. Epub 2020/12/10. doi:10.1158/0008-5472.Can-20-1619
Hamamoto, R., Furukawa, Y., Morita, M., Iimura, Y., Silva, F. P., Li, M., et al. (2004). SMYD3 encodes a histone methyltransferase involved in the proliferation of cancer cells. Nat. Cell. Biol. 6 (8), 731–740. doi:10.1038/ncb1151
Hu, G., Ma, J., Zhang, J., Chen, Y., Liu, H., Huang, Y., et al. (2021). Hypoxia-induced lnchilar promotes renal cancer metastasis via cerna for the mir-613/206/1-1-3p/jagged-1/notch/cxcr4 signaling pathway. Mol. Ther. J. Am. Soc. Gene Ther. 29 (10), 2979–2994. Epub 2021/06/01. doi:10.1016/j.ymthe.2021.05.020
Liu, Y., Wang, J., Li, L., Qin, H., Wei, Y., Zhang, X., et al. (2022). Ac010973.2 promotes cell proliferation and is one of six stemness-related genes that predict overall survival of renal clear cell carcinoma. Sci. Rep. 12 (1), 4272. Epub 2022/03/13. doi:10.1038/s41598-022-07070-1
Miranda Furtado, C. L., Dos Santos Luciano, M. C., Silva Santos, R. D., Furtado, G. P., Moraes, M. O., and Pessoa, C. (2019). Epidrugs: targeting epigenetic marks in cancer treatment. Epigenetics 14 (12), 1164–1176. doi:10.1080/15592294.2019.1640546
Ni, W., Yao, S., Zhou, Y., Liu, Y., Huang, P., et al. (2019). Long noncoding RNA GAS5 inhibits progression of colorectal cancer by interacting with and triggering YAP phosphorylation and degradation and is negatively regulated by the m6A reader YTHDF3. Mol. Cancer 18 (1), 143. doi:10.1186/s12943-019-1079-y
Pang, K., Dong, Y., Hao, L., Shi, Z. D., Zhang, Z. G., Chen, B., et al. (2022a). ERH interacts with EIF2α and regulates the eif2α/ATF4/CHOP pathway in bladder cancer cells. Front. Oncol. 12, 871687. Published 2022 Jun 14. doi:10.3389/fonc.2022.871687
Pang, K., Li, M. L., Hao, L., Shi, Z. D., Feng, H., Chen, B., et al. (2022b). ERH gene and its role in cancer cells. Front. Oncol. 12, 900496. Published 2022 May 23. doi:10.3389/fonc.2022.900496
Pang, K., Wang, W., Qin, J. X., Shi, Z. D., Hao, L., Ma, Y. Y., et al. (2020). Role of protein phosphorylation in cell signaling, disease, and the intervention therapy. MedComm 3 (4), e175. Published 2022 Nov 3. doi:10.1002/mco2.175
Powles, T., Plimack, E. R., Soulières, D., Waddell, T., Stus, V., Gafanov, R., et al. (2020). Pembrolizumab plus axitinib versus sunitinib monotherapy as first-line treatment of advanced renal cell carcinoma (Keynote-426): extended follow-up from a randomised, open-label, phase 3 trial. Lancet Oncol. 21 (12), 1563–1573. Epub 2020/12/08. doi:10.1016/s1470-2045(20)30436-8
Qu, L., Ding, J., Chen, C., Wu, Z. J., Liu, B., Gao, Y., et al. (2016). Exosome-transmitted lncarsr promotes sunitinib resistance in renal cancer by acting as a competing endogenous rna. Cancer Cell. 29 (5), 653–668. Epub 2016/04/28. doi:10.1016/j.ccell.2016.03.004
Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene Set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U. S. A. 102 (43), 15545–15550. Epub 2005/10/04. doi:10.1073/pnas.0506580102
Sun, L., Zhang, H., and Gao, P. (2022). Metabolic reprogramming and epigenetic modifications on the path to cancer. Protein Cell. 13 (12), 877–919. doi:10.1007/s13238-021-00846-7
Tan, Y. T., Lin, J. F., Li, T., Li, J. J., Xu, R. H., and Ju, H. Q. (2021). Lncrna-mediated posttranslational modifications and reprogramming of energy metabolism in cancer. Cancer Commun. Lond. Engl. 41 (2), 109–120. Epub 2020/10/30. doi:10.1002/cac2.12108
Turajlic, S., Swanton, C., and Boshoff, C. (2018). Kidney cancer: the next decade. J. Exp. Med. 215 (10), 2477–2479. Epub 2018/09/16. doi:10.1084/jem.20181617
Ulanovskaya, O. A., Zuhl, A. M., and Cravatt, B. F. (2013). NNMT promotes epigenetic remodeling in cancer by creating a metabolic methylation sink. Nat. Chem. Biol. 9 (5), 300–306. doi:10.1038/nchembio.1204
Wang, Q., Wu, G., Zhang, Z., Tang, Q., Zheng, W., Chen, X., et al. (2018). Long non-coding rna hottip promotes renal cell carcinoma progression through the regulation of the mir-615/igf-2 pathway. Int. J. Oncol. 53 (5), 2278–2288. Epub 2018/09/19. doi:10.3892/ijo.2018.4539
Wei, Y., Chen, X., Ren, X., Wang, B., Zhang, Q., Bu, H., et al. (2021). Identification of Mx2 as a novel prognostic biomarker for sunitinib resistance in clear cell renal cell carcinoma. Front. Genet. 12, 680369. Epub 2021/07/27. doi:10.3389/fgene.2021.680369
Xiong, W., Zhang, B., Yu, H., Zhu, L., Yi, L., and Jin, X. (2021). Rrm2 regulates sensitivity to sunitinib and Pd-1 blockade in renal cancer by stabilizing Anxa1 and activating the akt pathway. Adv. Sci. Weinheim, Baden-Wurttemberg, Ger. 8 (18), e2100881. Epub 2021/07/29. doi:10.1002/advs.202100881
Xiu, B., Chi, Y., Liu, L., Chi, W., Zhang, Q., Chen, J., et al. (2019). Linc02273 drives breast cancer metastasis by epigenetically increasing Agr2 transcription. Mol. cancer 18 (1), 187. Epub 2019/12/21. doi:10.1186/s12943-019-1115-y
Yang, L., Chen, Y., Liu, N., Shi, Q., Han, X., Gan, W., et al. (2021). Low expression of traf3ip2-as1 promotes progression of nono-tfe3 translocation renal cell carcinoma by stimulating N(6)-methyladenosine of Parp1 mrna and downregulating pten. J. Hematol. Oncol. 14 (1), 46. Epub 2021/03/21. doi:10.1186/s13045-021-01059-5
Yang, W., Soares, J., Greninger, P., Edelman, E. J., Lightfoot, H., Forbes, S., et al. (2013). Genomics of drug sensitivity in cancer (gdsc): a resource for therapeutic biomarker discovery in cancer cells. Nucleic acids Res. 41 (Database issue), D955–D961. Epub 2012/11/28. doi:10.1093/nar/gks1111
Keywords: lncRNAs, ccRCC, sensitivity, epigenetic, prognosis
Citation: Wang J, Dou P, Sun Y, Zheng J, Wu G, Liu H and Tao L (2024) Epigenetic dysregulated long non-coding RNAs in renal cell carcinoma based on multi-omics data and their influence on target drugs sensibility. Front. Genet. 15:1406150. doi: 10.3389/fgene.2024.1406150
Received: 24 March 2024; Accepted: 16 July 2024;
Published: 02 August 2024.
Edited by:
Giovanni Nassa, University of Salerno, ItalyCopyright © 2024 Wang, Dou, Sun, Zheng, Wu, Liu and Tao. 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: Heqian Liu, urologylhq1220@163.com; Lingsong Tao, taolingsong@sina.com
†These authors have contributed equally to this work and share first authorship