TBX5-AS1, an enhancer RNA, is a potential novel prognostic biomarker for lung adenocarcinoma

Characterization of an endoplasmic reticulum stress-associated lncRNA prognostic signature and the tumor-suppressive role of RP11-295G20.2 knockdown in lung adenocarcinoma

Endoplasmic reticulum stress (ERS) is commonly induced by accumulating misfolded or unfolded proteins in tumor microenvironment. Long non-coding RNAs (lncRNAs) play important roles in ERS response and lung adenocarcinoma (LUAD) progression. However, the role of ERS-related lncRNAs in LUAD remains unknown. In this study, we aimed to identify ERS-associated lncRNAs with prognostic value in LUAD and characterize their clinical implications. Cox and least absolute shrinkage and selection operator regression analyses identified nine ERS-related lncRNAs with independent prognostic abilities, including five protective factors (CROCCP2, KIAA0125, LINC0996, RPARP-AS1 and TBX5-AS1) and four risk factors (LINC0857, LINC116, RP11-21L23.2 and RP11-295G20.2). We developed an ERS-related lncRNA risk prediction model in predicting overall survival of LUAD patients, which classified TCGA cohorts into high-risk (HS) and low-risk (LS) groups. Comprehensive bioinformatic analyses revealed HS patients featured with late-stage tumors, greater mutation burdens, weaker anti-tumor immunity/responses, and lower sensitivity to targeted drugs compared to LS patients, contributing to tumor progression and a poor prognosis. Functional enrichment analysis implicated these ERS-related lncRNAs in cell migration, cell death, and immunity. Furthermore, expression of the most significantly upregulated risk lncRNA, RP11-295G20.2, was validated at the mRNA level using clinical LUAD samples. Knockdown of RP11-295G20.2 obviously reduced ERS and suppressed proliferation, invasion, and migration of LUAD cells. This novel ERS-related lncRNA signature provides a new biomarker for prognostic prediction, and ERS-associated RP11-295G20.2 serves as a potential therapeutic target in LUAD.

Prognostic and onco-immunological value of immune-related eRNAs-driven genes in lung adenocarcinoma

BACKGROUND

We aimed to comprehensively analyze the clinical value of immune-related eRNAs-driven genes in lung adenocarcinoma (LUAD) and find the potential biomarkers for prognosis and therapeutic response to improve the survival of this malignant disease.

MATERIALS AND METHODS

Pearson's correlation analysis was performed to identify the immune-related eRNAs-driven genes. Cox regression and least absolute shrinkage and selection operator (LASSO) analyses were used to construct this prognostic risk signature. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were used to investigate the underlying molecular mechanism. The single sample gene set enrichment analysis (ssGSEA) algorithm was conducted to evaluate the immune status based on the signature. The quantitative real-time PCR (qRT-PCR) analysis was performed to evaluate the expression value of the signature genes between LUAD tissues and adjacent lung tissues.

RESULTS

Five immune-related eRNAs-driven genes (SHC1, GDF10, CCL14, FYN, and NOD1) were identified to construct a prognostic risk signature with favorable predictive capacity. The patients with high-risk scores based on the signature were significantly associated with the malignant clinical features compared with those with low-risk scores. Kaplan-Meier analysis demonstrated that the sample in the low-risk group had a prolonged survival compared with those in the high-risk group. This risk signature was validated to have a promising predictive capacity and reliability in diverse clinical situations and independent cohorts. The functional enrichment analysis demonstrated that humoral immune response and intestinal immune network for IgA production pathway might be the underlying molecular mechanism related to the signature. The proportion of the vast majority of immune infiltrating cells in the high-risk group was significantly lower than that in the low-risk group, and the immunotherapy response rate in the low-risk group was significantly higher than that in the high-risk group. Moreover, BI-2536, sepantronium bromide, and ULK1 were the potential drugs for the treatment of patients with higher risk scores. Finally, the experiment in vivo and database analysis indicated that CCL14, FYN, NOD1, and GDF10 are the potential LUAD suppressor and SHC1 is a potential treatment target for LUAD.

CONCLUSION

Above all, we constructed a prognostic risk signature with favorable predictive capacity in LUAD, which was significantly associated with malignant features, immunosuppressive tumor microenvironment, and immunotherapy response and may provide clinical benefit in clinical decisions.

The benign nature and rare occurrence of cardiac myxoma as a possible consequence of the limited cardiac proliferative/ regenerative potential: a systematic review

BACKGROUND

Cardiac Myxoma is a primary tumor of heart. Its origins, rarity of the occurrence of primary cardiac tumors and how it may be related to limited cardiac regenerative potential, are not yet entirely known. This study investigates the key cardiac genes/ transcription factors (TFs) and signaling pathways to understand these important questions.

METHODS

Databases including PubMed, MEDLINE, and Google Scholar were searched for published articles without any date restrictions, involving cardiac myxoma, cardiac genes/TFs/signaling pathways and their roles in cardiogenesis, proliferation, differentiation, key interactions and tumorigenesis, with focus on cardiomyocytes.

RESULTS

The cardiac genetic landscape is governed by a very tight control between proliferation and differentiation-related genes/TFs/pathways. Cardiac myxoma originates possibly as a consequence of dysregulations in the gene expression of differentiation regulators including Tbx5, GATA4, HAND1/2, MYOCD, HOPX, BMPs. Such dysregulations switch the expression of cardiomyocytes into progenitor-like state in cardiac myxoma development by dysregulating Isl1, Baf60 complex, Wnt, FGF, Notch, Mef2c and others. The Nkx2-5 and MSX2 contribute predominantly to both proliferation and differentiation of Cardiac Progenitor Cells (CPCs), may possibly serve roles based on the microenvironment and the direction of cell circuitry in cardiac tumorigenesis. The Nkx2-5 in cardiac myxoma may serve to limit progression of tumorigenesis as it has massive control over the proliferation of CPCs. The cardiac cell type-specific genetic programming plays governing role in controlling the tumorigenesis and regenerative potential.

CONCLUSION

The cardiomyocytes have very limited proliferative and regenerative potential. They survive for long periods of time and tightly maintain the gene expression of differentiation genes such as Tbx5, GATA4 that interact with tumor suppressors (TS) and exert TS like effect. The total effect such gene expression exerts is responsible for the rare occurrence and benign nature of primary cardiac tumors. This prevents the progression of tumorigenesis. But this also limits the regenerative and proliferative potential of cardiomyocytes. Cardiac Myxoma develops as a consequence of dysregulations in these key genes which revert the cells towards progenitor-like state, hallmark of CM. The CM development in carney complex also signifies the role of TS in cardiac cells.

Prognostic signature of lipid metabolism associated LncRNAs predict prognosis and treatment of lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD) is the most predominant histological subtype of lung cancer. Abnormal lipid metabolism is closely related to the development of LUAD. LncRNAs are involved in the regulation of various lipid metabolism-related genes in various cancer cells including LUAD. Here, we aimed to identify lipid metabolism-related lncRNAs associated with LUAD prognosis and to propose a new prognostic signature.

METHODS

First, differentially expressed lncRNAs (DE-lncRNAs) from the TCGA-LUAD and the GSE31210 dataset were identified. Then the correlation analysis between DE-lncRNAs and lipid metabolism genes was performed to screen lipid metabolism-related lncRNAs. Cox regression analyses were performed in the training set to establish a prognostic model and the model was validated in the testing set and the validation set. Moreover, The role of this model in the underlying molecular mechanisms, immunotherapy, and chemotherapeutic drug sensitivity analysis was predicted by methods such as Gene Set Enrichment Analysis, immune infiltration, tumor mutational burden (TMB), neoantigen, Tumor Immune Dysfunction and Exclusion, chemosensitivity analysis between the high- and low-risk groups. The diagnostic ability of prognostic lncRNAs has also been validated. Finally, we validated the expression levels of selected prognostic lncRNAs by quantitative real-time polymerase chain reaction (qRT-PCR).

RESULTS

The prognostic model was constructed based on four prognostic lncRNAs (LINC00857, EP300-AS1, TBX5-AS1, SNHG3) related to lipid metabolism. The receiver operating characteristic curve (ROC) and Kaplan Meier (KM) curves of the risk model showed their validity. The results of Gene Set Enrichment Analysis suggested that differentially expressed genes in high- and low-risk groups were mainly enriched in immune response and cell cycle. There statistical differences in TMB and neoantigen between high- and low-risk groups. Drug sensitivity analysis suggested that patients with low risk scores may have better chemotherapy outcomes. The results of qRT-PCR were suggesting that compared with the normal group, the expressions of EP300-AS1 and TBX5-AS1 were down-regulated in the tumor group, while the expressions of LINC00857 and SNHG3 were up-regulated. The four prognostic lncRNAs had good diagnostic capabilities, and the overall diagnostic model of the four prognostic lncRNAs was more effective.

CONCLUSION

A total of 4 prognostic lncRNAs related to lipid metabolism were obtained and an effective risk model was constructed.

Emerging Role of Enhancer RNAs as Potential Diagnostic and Prognostic Biomarkers in Cancer

Enhancers are distal cis-acting elements that are commonly recognized to regulate gene expression via cooperation with promoters. Along with regulating gene expression, enhancers can be transcribed and generate a class of non-coding RNAs called enhancer RNAs (eRNAs). The current discovery of abundant tissue-specific transcription of enhancers in various diseases such as cancers raises questions about the potential role of eRNAs in disease diagnosis and therapy. This review aimed to demonstrate the current understanding of eRNAs in cancer research with a focus on the potential roles of eRNAs as prognostic and diagnostic biomarkers in cancers.

Enhancer RNA (eRNA) in Human Diseases

Enhancer RNAs (eRNAs), a class of non-coding RNAs (ncRNAs) transcribed from enhancer regions, serve as a type of critical regulatory element in gene expression. There is increasing evidence demonstrating that the aberrant expression of eRNAs can be broadly detected in various human diseases. Some studies also revealed the potential clinical utility of eRNAs in these diseases. In this review, we summarized the recent studies regarding the pathological mechanisms of eRNAs as well as their potential utility across human diseases, including cancers, neurodegenerative disorders, cardiovascular diseases and metabolic diseases. It could help us to understand how eRNAs are engaged in the processes of diseases and to obtain better insight of eRNAs in diagnosis, prognosis or therapy. The studies we reviewed here indicate the enormous therapeutic potency of eRNAs across human diseases.

Comprehensive transcriptomic analysis of immune-related eRNAs associated with prognosis and immune microenvironment in melanoma

Background

Recent evidence suggests that enhancer RNAs (eRNAs) play key roles in cancers. Identification of immune-related eRNAs (ireRNAs) in melanoma can provide novel insights into the mechanisms underlying its genesis and progression, along with potential therapeutic targets.

Aim

To establish an ireRNA-related prognostic signature for melanoma and identify potential drug candidates.

Methods

The ireRNAs associated with the overall survival (OS-ireRNAs) of melanoma patients were screened using data from The Cancer Genome Atlas (TCGA) via WGCNA and univariate Cox analysis. A prognostic signature based on these OS-ireRNAs was then constructed by performing the least absolute shrinkage and selection operator (LASSO) Cox regression analysis. The immune landscape associated with the prognostic model was evaluated by the ESTIMATE algorithm and CIBERSORT method. Finally, the potential drug candidates for melanoma were screened through the cMap database.

Results

A total of 24 OS-ireRNAs were obtained, of which 7 ireRNAs were used to construct a prognostic signature. The ireRNAs-related signature performed well in predicting the overall survival (OS) of melanoma patients. The risk score of the established signature was further verified as an independent risk factor, and was associated with the unique tumor microenvironment in melanoma. We also identified several potential anti-cancer drugs for melanoma, of which corticosterone ranked first.

Conclusions

The ireRNA-related signature is an effective prognostic predictor and provides reliable information to better understand the mechanism of ireRNAs in the progression of melanoma.

A comprehensive prognostic and immune analysis of enhancer RNA identifies IGFBP7-AS1 as a novel prognostic biomarker in Uterine Corpus Endometrial Carcinoma

BACKGROUND

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.

METHODS

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.

RESULTS

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.

CONCLUSIONS

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.

A Two-eRNA-Based Signature Can Impact the Immune Status and Predict the Prognosis and Drug Sensitivity of Lung Adenocarcinoma

Enhancer RNAs (eRNAs) are intergenic long noncoding RNAs (lncRNAs) participating in the development of malignant cancers via targeting cancer-associated genes and immune checkpoints. Immune infiltration of the tumor microenvironment was positively associated with overall survival (OS) in lung adenocarcinoma (LUAD). In this study, we aimed to explore the clinical significance of PCBP1-AS1 in LUAD and developed a novel prognostic signature based on two eRNAs. Our team discovered that the expression of PCBP1-AS1 was distinctly downregulated in LUAD specimens compared with nontumor specimens. Lower PCBP1-AS1 expression was related to advanced clinical stages and poor prognosis. KEGG analysis unveiled that the coexpression genes of PCBP1-AS1 were involved in the regulation of several tumor-related pathways. In addition, remarkable associations were observed between the expression of PCBP1-AS1 and the levels of several immune cells. Then, we used PCBP1-AS1 and TBX5-AS1 to develop a prognostic model. Survival assays unveiled that patients with higher risk scores exhibited a shorter OS in contrast to patients with lower risk scores. In addition, multivariable Cox regressive analysis indicated that the risk score was an independent prediction factor in LUAD sufferers. The anticancer drug sensitivity analysis indicated that risk score had a positive relationship with several anticancer drugs. Taken together, our findings indicated PCBP1-AS1 as a function modulator in LUAD development. In addition, we constructed a robust immune-related eRNA signature which might be a clinical prognosis factor for LUAD patients.

Identification of Enhancer RNA CDK6-AS1 as a Potential Novel Prognostic Biomarker in Gastric Cancer

Background: This study aimed to confirm the role of enhancer RNAs (eRNAs) in gastric cancer and their clinical utility.

Methods: We used Cox survival and relevance analysis to identify the candidate eRNAs in gastric cancer and performed Gene Ontology and Reactome pathway enrichment to determine the potential functions of eRNAs. Correlation between eRNA, tumor-infiltrating immune cells, and drug sensitivity was then analyzed.

Results:CDK6-AS1, a long non-coding RNA cyclin-dependent kinase 6, may serve as a poor potential prognostic biomarker candidate in gastric cancer with a positive correlation with its target gene CDK6. The low CDK6-AS1 expression group showed more frequent mutated driver genes than the high expression group. Moreover, CDK6-AS1 is involved in a key oncogenic pathway of the cell cycle and RNA transcription. CDK6-AS1 also shows dysregulations and associations with prognosis at the pan-cancer level. This eRNA may also be associated with immune cell infiltration and drug sensitivity.

Conclusion:CDK6-AS1 may be a potential prognostic biomarker and chemotherapeutic drug sensitivity predictor in gastric cancer.

Tobacco Smoke and Electronic Cigarette Vapor Alter Enhancer RNA Expression That Can Regulate the Pathogenesis of Lung Squamous Cell Carcinoma

Simple Summary

It is well established that tobacco smoke is the key player in lung squamous cell carcinoma (LUSC) pathogenesis, and there is growing evidence that electronic cigarette (e-cigarette) vapor may also cause LUSC. Recently, several studies have associated tobacco smoke with differential enhancer RNA (eRNA) expression. However, the effects of tobacco smoke and e-cigarette vapor on eRNA expression in correlation to LUSC outcomes have not been fully elucidated. This study demonstrates that tobacco smoke and e-cigarette vapor may decrease DNA methylation and increase chromosomal alterations at key sites, which ultimately upregulate the expression of oncogenic eRNAs and downregulate the expression of tumor-suppressing eRNAs. Subsequently, we demonstrate that these eRNAs may have altered interactions with immune cells to promote LUSC pathogenesis and reduced patient survival. We hope our results can be validated in future studies, and the key eRNAs we identified may be used as effective targets for more specialized treatments for smoking-mediated LUSC.

Abstract

Tobacco is the primary etiologic agent in worsened lung squamous cell carcinoma (LUSC) outcomes. Meanwhile, it has been shown that etiologic agents alter enhancer RNAs (eRNAs) expression. Therefore, we aimed to identify the effects of tobacco and electronic cigarette (e-cigarette) use on eRNA expression in relation to LUSC outcomes. We extracted eRNA counts from RNA-sequencing data of tumor/adjacent normal tissue and before/after e-cigarette tissue from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO), respectively. Tobacco-mediated LUSC eRNAs were correlated to patient survival, clinical variables, and immune-associated elements. eRNA expression was also correlated to mutation rates through the Repeated Evaluation of Variables Conditional Entropy and Redundance (REVEALER) algorithm and methylated sites through methylationArrayAnalysis. Differential expression analysis was then completed for the e-cigarette data to compare with key tobacco-mediated eRNAs. We identified 684 downregulated eRNAs and 819 upregulated eRNAs associated with tobacco-mediated LUSC, specifically, with the cancer pathological stage. We also observed a decrease in immune cell abundance in tobacco-mediated LUSC. Yet, we found an increased association of eRNA expression with immune cell abundance in tobacco-mediated LUSC. We identified 16 key eRNAs with significant correlations to 8 clinical variables, implicating these eRNAs in LUSC malignancy. Furthermore, we observed that these 16 eRNAs were highly associated with chromosomal alterations and reduced CpG site methylation. Finally, we observed large eRNA expression upregulation with e-cigarette use, which corresponded to the upregulation of the 16 key eRNAs. Our findings provide a novel mechanism by which tobacco and e-cigarette smoke influences eRNA interactions to promote LUSC pathogenesis and provide insight regarding disease progression at a molecular level.