ImmReg: the regulon atlas of immune-related pathways across cancer types

The Oncogenic Role of VWA8-AS1, a Long Non-Coding RNA, in Epstein-Barr Virus-Associated Oral Squamous Cell Carcinoma: An Integrative Transcriptome and Functional Analysis

Dysregulated long non-coding RNA (lncRNA) expression is linked to various cancers and may be influenced by oncogenic Epstein-Barr virus (EBV) infection, a known and detectable risk factor in oral squamous cell carcinoma (OSCC) patients. However, research on the oncogenic role of EBV-induced lncRNAs in OSCC is limited. To identify lncRNA-associated EBV infection and OSCC carcinogenesis, the differential expression of RNA-seq datasets from paired normal adjacent and OSCC tissues, and microarray data from EBV-negative and EBV-positive SCC25 cells, were identified and selected, respectively, for interaction, functional analysis, and CCK-8 cell proliferation, wound healing, and invasion Transwell assays. In OSCC tissues, 6731 differentially expressed lncRNAs were identified when compared to normal tissues from RNA-seq datasets, with 295 linked to EBV-induced OSCC carcinogenesis from microarray datasets. The EBV-induced lncRNA VWA8-AS1 showed significant upregulation in EBV-positive SCC25 cells and EBV-infected adjacent and OSCC tissue samples. VWA8-AS1 potentially promotes OSCC via the lncRNA-miRNA-mRNA axis or direct protein interactions, affecting various cellular processes. Studies in OSCC cell lines revealed that elevated VWA8-AS1 levels enhanced cell migration and invasion. This study demonstrates VWA8-AS1's contribution to tumor progression and possible interactions with its targets in OSCC, offering insights for future research on functional mechanisms and therapeutic targets in EBV-associated OSCC.

Pan-cancer analyses identify MIR210HG overexpression, epigenetic regulation and oncogenic role in human tumors and its interaction with the tumor microenvironment

BACKGROUND

Molecular entities showing dysregulation in multiple cancers may hold great biomarker or therapeutic potential. There is accumulating evidence that highlights the dysregulation of a long non-coding RNA, MIR210HG, in various cancers and its oncogenic role. However, a comprehensive analysis of MIR210HG expression pattern, molecular mechanisms, diagnostic or prognostic significance or evaluation of its interaction with tumor microenvironment across various cancers remains unstudied.

METHODS

A systematic pan-cancer analysis was done using multiple public databases and bioinformatic tools to study the molecular role and clinical significance of MIR210HG. We have analyzed expression patterns, genome alteration, transcriptional and epigenetic regulation, correlation with patient survival, immune infiltrates, co-expressed genes, interacting proteins, and pathways associated with MIR210HG.

RESULTS

The Pan cancer expression analysis of MIR210HG through various tumor datasets demonstrated that MIR210HG is significantly upregulated in most cancers and increased with the tumor stage in a subset of them. Furthermore, prognostic analysis revealed high MIR210HG expression is associated with poor overall and disease-free survival in specific cancer types. Genetic alteration analysis showed minimal alterations in the MIR210HG locus, indicating that overexpression in cancers is not due to gene amplification. The exploration of SNPs on MIR210HG suggested possible structural changes that may affect its interactions with the miRNAs. The correlation of MIR210HG with promoter methylation was found to be significantly negative in nature in majority of cancers depicting the possible epigenetic regulation of expression of MIR210HG. Additionally, MIR210HG showed negative correlations with immune cells and thus may have strong impact on the tumor microenvironment. Functional analysis indicates its association with hypoxia, angiogenesis, metastasis, and DNA damage repair processes. MIR210HG was found to interact with several proteins and potentially regulate chromatin modifications and transcriptional regulation.

CONCLUSIONS

A first pan-can cancer analysis of MIR210HG highlights its transcriptional and epigenetic deregulation and oncogenic role in the majority of cancers, its correlation with tumor microenvironment factors such as hypoxia and immune infiltration, and its potential as a prognostic biomarker and therapeutic target in several cancers.

SORC: an integrated spatial omics resource in cancer

The interactions between tumor cells and the microenvironment play pivotal roles in the initiation, progression and metastasis of cancer. The advent of spatial transcriptomics data offers an opportunity to unravel the intricate dynamics of cellular states and cell-cell interactions in cancer. Herein, we have developed an integrated spatial omics resource in cancer (SORC, http://bio-bigdata.hrbmu.edu.cn/SORC), which interactively visualizes and analyzes the spatial transcriptomics data in cancer. We manually curated currently available spatial transcriptomics datasets for 17 types of cancer, comprising 722 899 spots across 269 slices. Furthermore, we matched reference single-cell RNA sequencing data in the majority of spatial transcriptomics datasets, involving 334 379 cells and 46 distinct cell types. SORC offers five major analytical modules that address the primary requirements of spatial transcriptomics analysis, including slice annotation, identification of spatially variable genes, co-occurrence of immune cells and tumor cells, functional analysis and cell-cell communications. All these spatial transcriptomics data and in-depth analyses have been integrated into easy-to-browse and explore pages, visualized through intuitive tables and various image formats. In summary, SORC serves as a valuable resource for providing an unprecedented spatially resolved cellular map of cancer and identifying specific genes and functional pathways to enhance our understanding of the tumor microenvironment.

LncSEA 2.0: an updated platform for long non-coding RNA related sets and enrichment analysis

Long non-coding RNAs (lncRNAs) possess a wide range of biological functions, and research has demonstrated their significance in regulating major biological processes such as development, differentiation, and immune response. The accelerating accumulation of lncRNA research has greatly expanded our understanding of lncRNA functions. Here, we introduce LncSEA 2.0 (http://bio.liclab.net/LncSEA/index.php), aiming to provide a more comprehensive set of functional lncRNAs and enhanced enrichment analysis capabilities. Compared with LncSEA 1.0, we have made the following improvements: (i) We updated the lncRNA sets for 11 categories and extremely expanded the lncRNA scopes for each set. (ii) We newly introduced 15 functional lncRNA categories from multiple resources. This update not only included a significant amount of downstream regulatory data for lncRNAs, but also covered numerous epigenetic regulatory data sets, including lncRNA-related transcription co-factor binding, chromatin regulator binding, and chromatin interaction data. (iii) We incorporated two new lncRNA set enrichment analysis functions based on GSEA and GSVA. (iv) We adopted the snakemake analysis pipeline to track data processing and analysis. In summary, LncSEA 2.0 offers a more comprehensive collection of lncRNA sets and a greater variety of enrichment analysis modules, assisting researchers in a more comprehensive study of the functional mechanisms of lncRNAs.

Dissecting the tumour immune microenvironment in merkel cell carcinoma based on a machine learning framework

Merkel cell carcinoma (MCC) is a primary cutaneous neoplasm of neuroendocrine carcinoma of the skin, which is characterized by molecular heterogeneity with diverse tumour microenvironment (TME). However, we are still lack knowledge of the cellular states and ecosystems in MCC. Here, we systematically identified and characterized the landscape of cellular states and ecotypes in MCC based on a machine learning framework. We obtained 30 distinct cellular states from 9 immune cell types and investigated the B cell, CD8 T cell, fibroblast, and monocytes/macrophage cellular states in detail. The functional profiling of cellular states were investigated and found the genes highly expressed in cellular states were significantly enriched in immune- and cancer hallmark-related pathways. In addition, four ecotypes were further identified which were with different patient compositions. Transcriptional regulation analysis revealed the critical transcription factors (i.e. E2F1, E2F3 and E2F7), which play important roles in regulating the TME of MCC. In summary, the findings of this study may provide rich knowledge to understand the intrinsic subtypes of MCCs and the pathways involved in distinct subtype oncogenesis, and will further advance the knowledge in developing a specific therapeutic strategy for these MCC subtypes.

Pan-cancer landscape of immunology PIWI-interacting RNAs

PIWI-interacting RNAs (piRNAs), an emergent type of non-coding RNAs during oncogenesis, play critical roles in regulating tumor microenvironment. Systematic analysis of piRNAs' roles in modulating immune pathways is important for tumor immunotherapy. In this study, in-depth analysis of piRNAs was performed to develop an integrated computational algorithm, the immunology piRNA (ImmPI) pipeline, for uncovering the global expression landscape of piRNAs and identifying their regulatory roles in immune pathways. The immunology piRNAs show a tendency towards overexpression patterns in immune cells, causing perturbations in tumors, being significantly associated with infiltration of immune cells, and having prognostic value. The ImmPI score can contribute to prioritizing tumor-related piRNAs and distinguish two subtypes of SKCM (immune-cold and hot phenotypes), as characterized by different prognoses, immunogenicity and antitumor immunity. Finally, we developed an interactive web resource (ImmPI portal: http://www.hbpding.com/ImmPi) for the biomedical research community, with several useful modules to browse, visualize, and download the results of immunology piRNAs analysis. Overall, our work provides a comprehensive landscape of piRNAs across multiple cancer types and sheds light on their regulatory and functional roles in tumor immunity. These findings pave the way for future research and development of piRNA-based immunotherapies for cancer treatment.

Pan-cancer characterization of ncRNA synergistic competition uncovers potential carcinogenic biomarkers

Non-coding RNAs (ncRNAs) act as important modulators of gene expression and they have been confirmed to play critical roles in the physiology and development of malignant tumors. Understanding the synergism of multiple ncRNAs in competing endogenous RNA (ceRNA) regulation can provide important insights into the mechanisms of malignant tumors caused by ncRNA regulation. In this work, we present a framework, SCOM, for identifying ncRNA synergistic competition. We systematically construct the landscape of ncRNA synergistic competition across 31 malignant tumors, and reveal that malignant tumors tend to share hub ncRNAs rather than the ncRNA interactions involved in the synergistic competition. In addition, the synergistic competition ncRNAs (i.e. ncRNAs involved in the synergistic competition) are likely to be involved in drug resistance, contribute to distinguishing molecular subtypes of malignant tumors, and participate in immune regulation. Furthermore, SCOM can help to infer ncRNA synergistic competition across malignant tumors and uncover potential diagnostic and prognostic biomarkers of malignant tumors. Altogether, the SCOM framework (https://github.com/zhangjunpeng411/SCOM/) and the resulting web-based database SCOMdb (https://comblab.cn/SCOMdb/) serve as a useful resource for exploring ncRNA regulation and to accelerate the identification of carcinogenic biomarkers.

Roles of non-coding RNAs in eye development and diseases

The prevalence of ocular disorders is dramatically increasing worldwide, especially those that cause visual impairment and permanent loss of vision, including cataract, glaucoma, age-related macular degeneration, and diabetic retinopathy. Extensive evidence has shown that ncRNAs are key regulators in various biogenesis and biological functions, controlling gene expression related to histogenesis and cell differentiation in ocular tissues. Aberrant expression and function of ncRNA can lead to dysfunction of visual system and mediate progression of eye disorders. Here, we mainly offer an overview of the role of precise modulation of ncRNAs in eye development and function in patients with eye diseases. We also highlight the challenges and future perspectives in conducting ncRNA studies, focusing specifically on the role of ncRNAs that may hold expanded promise for their diagnostic and therapeutic applications in various eye diseases. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA in Disease and Development > RNA in Disease RNA in Disease and Development > RNA in Development.

WITHDRAWN: JUN and ATF3 are deficient in prostate cancer patients and their delivery in vivo via lipid nanoparticles has therapeutic efficacy by enhancing immune surveillance

This article has been withdrawn at the request of the author(s) and/or editor. The Publisher apologizes for any inconvenience this may cause. The full Elsevier Policy on Article Withdrawal can be found at https://www.elsevier.com/about/policies/article-withdrawal.

Deep immunophenotyping reveals clinically distinct cellular states and ecosystems in large-scale colorectal cancer

Determining the diverse cell types in the tumor microenvironment (TME) and their organization into cellular communities, is critical for understanding the biological heterogeneity and therapy of cancer. Here, we deeply immunophenotype the colorectal cancer (CRC) by integrative analysis of large-scale bulk and single cell transcriptome of 2350 patients and 53,137 cells. A rich landscape of 42 cellular states and 7 ecosystems in TMEs is uncovered and extend the previous immune classifications of CRC. Functional pathways and potential transcriptional regulators analysis of cellular states and ecosystems reveal cancer hallmark-related pathways and several critical transcription factors in CRC. High-resolution characterization of the TMEs, we discover the potential utility of cellular states (i.e., Monocytes/Macrophages and CD8 T cell) and ecosystems for prognosis and clinical therapy selection of CRC. Together, our results expand our understanding of cellular organization in TMEs of CRC, with potential implications for the development of biomarkers and precision therapies.

Identification of immunity-related lncRNAs and construction of a ceRNA network of potential prognostic biomarkers in acute myeloid leukemia

Objective: Using bioinformatics analyses, this study aimed to identify lncRNAs related to the immune status of acute myeloid leukemia (AML) patients and ascertain the potential impact in immunity-related competing endogenous RNA (ceRNA) networks on AML prognosis. Methods: AML-related RNA-seq FPKM data, AML-related miRNA expression microarray data, and gene sets associated with immunity-related pathways were, respectively, obtained from the TCGA, GEO, and ImmReg databases. An immunity-related ceRNA network was then constructed according to the predicted interactions between AML-related mRNAs, lncRNAs, and miRNAs. After performing LASSO and multivariate Cox regression analyses, lncRNAs in the ceRNA network were used to establish an AML prognostic model. According to mutual regulatory relationships and consistent trends of expression among candidate ceRNAs, two ceRNA subnetworks related to the AML prognostic model were determined. Finally, the correlation between the expression levels of mRNAs, lncRNAs, and miRNAs in each ceRNA subnetwork and immune cell infiltration (assessed by combining the ESTIMATE and CIBERSORT methods and ssGSEA) was analyzed. Results: A total of 424 immunity-related differentially expressed (IR-DE) mRNAs (IR-DEmRNAs), 191 IR-DElncRNAs, and 69 IR-DEmiRNAs were obtained, and a ceRNA network of 20 IR-DElncRNAs, 6 IR-DEmRNAs, and 3 IR-DEmiRNAs was established. Univariate Cox regression analysis was conducted on 20 IR-DElncRNAs, and 7 of these were identified to be significantly correlated with the overall survival (OS) time in AML patients. Then, two IR-DElncRNAs (MEG3 and HCP5) were screened as independent OS-related factors by LASSO and multivariable Cox regression analyses, and a prognostic model was constructed to evaluate the survival risk in AML patients. Survival analyses indicated that the OS of patients was often poor in the high-risk group. Additionally, from this model, two ceRNA regulatory pathways, namely, MEG3/miR-125a-5p/SEMA4C and HCP5/miR-125b-5p/IL6R, which were potentially involved in the immune regulation of AML prognosis were identified. Conclusion: lncRNAs HCP5 and MEG3 may act as key ceRNAs in the pathogenesis in AML by regulating immune cell representation as part of the regulatory lncRNA-miRNA-mRNA axes. The candidate mRNAs, lncRNAs, and miRNAs included in the ceRNA network identified here may serve as useful prognostic biomarkers and immunotherapeutic targets for AML.

RNAseq profiling of blood from patients with coronary artery disease: Signature of a T cell imbalance

Background

Cardiovascular disease had a global prevalence of 523 million cases and 18.6 million deaths in 2019. The current standard for diagnosing coronary artery disease (CAD) is coronary angiography either by invasive catheterization (ICA) or computed tomography (CTA). Prior studies employed single-molecule, amplification-independent RNA sequencing of whole blood to identify an RNA signature in patients with angiographically confirmed CAD. The present studies employed Illumina RNAseq and network co-expression analysis to identify systematic changes underlying CAD.

Methods

Whole blood RNA was depleted of ribosomal RNA (rRNA) and analyzed by Illumina total RNA sequencing (RNAseq) to identify transcripts associated with CAD in 177 patients presenting for elective invasive coronary catheterization. The resulting transcript counts were compared between groups to identify differentially expressed genes (DEGs) and to identify patterns of changes through whole genome co-expression network analysis (WGCNA).

Results

The correlation between Illumina amplified RNAseq and the prior SeqLL unamplified RNAseq was quite strong (r = 0.87), but there was only 9 % overlap in the DEGs identified. Consistent with the prior RNAseq, the majority (93 %) of DEGs were down-regulated ~1.7-fold in patients with moderate to severe CAD (>20 % stenosis). DEGs were predominantly related to T cells, consistent with known reductions in Tregs in CAD. Network analysis did not identify pre-existing modules with a strong association with CAD, but patterns of T cell dysregulation were evident. DEGs were enriched for transcripts associated with ciliary and synaptic transcripts, consistent with changes in the immune synapse of developing T cells.

Conclusions

These studies confirm and extend a novel mRNA signature of a Treg-like defect in CAD. The pattern of changes is consistent with stress-related changes in the maturation of T and Treg cells, possibly due to changes in the immune synapse.

Data analysis guidelines for single-cell RNA-seq in biomedical studies and clinical applications

The application of single-cell RNA sequencing (scRNA-seq) in biomedical research has advanced our understanding of the pathogenesis of disease and provided valuable insights into new diagnostic and therapeutic strategies. With the expansion of capacity for high-throughput scRNA-seq, including clinical samples, the analysis of these huge volumes of data has become a daunting prospect for researchers entering this field. Here, we review the workflow for typical scRNA-seq data analysis, covering raw data processing and quality control, basic data analysis applicable for almost all scRNA-seq data sets, and advanced data analysis that should be tailored to specific scientific questions. While summarizing the current methods for each analysis step, we also provide an online repository of software and wrapped-up scripts to support the implementation. Recommendations and caveats are pointed out for some specific analysis tasks and approaches. We hope this resource will be helpful to researchers engaging with scRNA-seq, in particular for emerging clinical applications.

Development and Validation of Lactate Metabolism-Related lncRNA Signature as a Prognostic Model for Lung Adenocarcinoma

BACKGROUND

Lung cancer has been a prominent research focus in recent years due to its role in cancer-related fatalities globally, with lung adenocarcinoma (LUAD) being the most prevalent histological form. Nonetheless, no signature of lactate metabolism-related long non-coding RNAs (LMR-lncRNAs) has been developed for patients with LUAD. Accordingly, we aimed to develop a unique LMR-lncRNA signature to determine the prognosis of patients with LUAD.

METHOD

The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases were utilized to derive the lncRNA expression patterns. Identification of LMR-lncRNAs was accomplished by analyzing the co-expression patterns between lncRNAs and LMR genes. Subsequently, the association between lncRNA levels and survival outcomes was determined to develop an effective signature. In the TCGA cohort, Cox regression was enlisted to build an innovative signature consisting of three LMR-lncRNAs, which was validated in the GEO validation cohort. GSEA and immune infiltration analysis were conducted to investigate the functional annotation of the signature and the function of each type of immune cell.

RESULTS

Fourteen differentially expressed LMR-lncRNAs were strongly correlated with the prognosis of patients with LUAD and collectively formed a new LMR-lncRNA signature. The patients could be categorized into two cohorts based on their LMR-lncRNA signatures: a low-risk and high-risk group. The overall survival of patients with LUAD in the high-risk group was considerably lower than those in the low-risk group. Using Cox regression, this signature was shown to have substantial potential as an independent prognostic factor, which was further confirmed in the GEO cohort. Moreover, the signature could anticipate survival across different groups based on stage, age, and gender, among other variables. This signature also correlated with immune cell infiltration (including B cells, neutrophils, CD4⁺ T cells, CD8⁺ T cells, etc.) as well as the immune checkpoint blockade target CTLA-4.

CONCLUSION

We developed and verified a new LMR-lncRNA signature useful for anticipating the survival of patients with LUAD. This signature could give potentially critical insight for immunotherapy interventions in patients with LUAD.