Role of ferroptosis in neuroimmunity and neurodegeneration in multiple sclerosis revealed by multi-omics data.

Previous studies have found that ferroptosis plays an important role in a variety of neurological diseases. However, the precise role of ferroptosis in the multiple sclerosis patients remains uncertain. We defined and validated a computational metric of ferroptosis levels. The ferroptosis scores were computed using the AUCell method, which reflects the enrichment scores of ferroptosis-related genes through gene ranking. The reliability of the ferroptosis score was assessed using various methods, involving cells induced to undergo ferroptosis by six different ferroptosis inducers. Through a comprehensive approach integrating snRNA-seq, spatial transcriptomics, and spatial proteomics data, we explored the role of ferroptosis in multiple sclerosis. Our findings revealed that among seven sampling regions of different white matter lesions, the edges of active lesions exhibited the highest ferroptosis score, which was associated with activation of the phagocyte system. Remyelination lesions exhibit the lowest ferroptosis score. In the cortex, ferroptosis score were elevated in neurons, relevant to a variety of neurodegenerative disease-related pathways. Spatial transcriptomics demonstrated a significant co-localization among ferroptosis score, neurodegeneration and microglia, which was verified by spatial proteomics. Furthermore, we established a diagnostic model of multiple sclerosis based on 24 ferroptosis-related genes in the peripheral blood. Ferroptosis might exhibits a dual role in the context of multiple sclerosis, relevant to both neuroimmunity and neurodegeneration, thereby presenting a promising and novel therapeutic target. Ferroptosis-related genes in the blood that could potentially serve as diagnostic and prognostic markers for multiple sclerosis.

MSGD: a manually curated database of genomic, transcriptomic, proteomic and drug information for multiple sclerosis.

Multiple sclerosis (MS) is the most common inflammatory demyelinating disease of the central nervous system. 'Omics' technologies (genomics, transcriptomics, proteomics) and associated drug information have begun reshaping our understanding of multiple sclerosis. However, these data are scattered across numerous references, making them challenging to fully utilize. We manually mined and compiled these data within the Multiple Sclerosis Gene Database (MSGD) database, intending to continue updating it in the future. We screened 5485 publications and constructed the current version of MSGD. MSGD comprises 6255 entries, including 3274 variant entries, 1175 RNA entries, 418 protein entries, 313 knockout entries, 612 drug entries and 463 high-throughput entries. Each entry contains detailed information, such as species, disease type, detailed gene descriptions (such as official gene symbols), and original references. MSGD is freely accessible and provides a user-friendly web interface. Users can easily search for genes of interest, view their expression patterns and detailed information, manage gene sets and submit new MS-gene associations through the platform. The primary principle behind MSGD's design is to provide an exploratory platform, aiming to minimize filtration and interpretation barriers while ensuring highly accessible presentation of data. This initiative is expected to significantly assist researchers in deciphering gene mechanisms and improving the prevention, diagnosis and treatment of MS. Database URL: http://bio-bigdata.hrbmu.edu.cn/MSGD.

Comprehensive quantitative radiogenomic evaluation reveals novel radiomic subtypes with distinct immune pattern in glioma.

BACKGROUND: Accurate classification of gliomas is critical to the selection of immunotherapy, and MRI contains a large number of radiomic features that may suggest some prognostic relevant signals. We aim to predict new subtypes of gliomas using radiomic features and characterize their survival, immune, genomic profiles and drug response. METHODS: We initially obtained 341 images of 36 patients from the CPTAC dataset for the development of deep learning models. Further 1812 images of 111 patients from TCGA_GBM and 152 images of 53 patients from TCGA_LGG were collected for testing and validation. A deep learning method based on Mask R-CNN was developed to identify new subtypes of glioma patients and compared the survival status, immune infiltration patterns, genomic signatures, specific drugs, and predictive models of different subtypes. RESULTS: 200 glioma patients (mean age, 33 years ± 19 [standard deviation]) were enrolled. The accuracy of the deep learning model for identifying tumor regions achieved 88.3 % (98/111) in the test set and 83 % (44/53) in the validation set. The sample was divided into two subtypes based on radiomic features showed different prognostic outcomes (hazard ratio, 2.70). According to the results of the immune infiltration analysis, the subtype with a poorer prognosis was defined as the immunosilencing radiomic (ISR) subtype (n = 43), and the other subtype was the immunoactivated radiomic (IAR) subtype (n = 53). Subtype-specific genomic signatures distinguished celllines into ISR celllines (n = 9) and control celllines (n = 13), and identified eight ISR-specific drugs, four of which were validated by the OCTAD database. Three machine learning-based classifiers showed that radiomic and genomic co-features better predicted the radiomic subtypes of gliomas. CONCLUSIONS: These findings provide insights into how radiogenomic could identify specific subtypes that predict prognosis, immune and drug sensitivity in a non-invasive manner.

The landscape of immune checkpoint-related long non-coding RNAs core regulatory circuitry reveals implications for immunoregulation and immunotherapy responses.

Long non-coding RNAs (lncRNAs) could modulate expression of immune checkpoints (ICPs) by cooperating with immunity genes in tumor immunization. However, precise functions in immunity and potential for predicting ICP inhibitors (ICI) response have been described for only a few lncRNAs. Here we present an integrated framework that leverages network-based analyses and Bayesian network inference to identify the regulated relationships including lncRNA, ICP and immunity genes as ICP-related LncRNAs mediated Core Regulatory Circuitry Triplets (ICP-LncCRCTs) that can make robust predictions. Hub ICP-related lncRNAs such as MIR155HG and ADAMTS9-AS2 were highlighted to play central roles in immune regulation. Specific ICP-related lncRNAs could distinguish cancer subtypes. Moreover, the ICP-related lncRNAs are likely to significantly correlated with immune cell infiltration, MHC, CYT. Some ICP-LncCRCTs such as CXCL10-MIR155HG-ICOS could better predict one-, three- and five-year prognosis compared to single molecule in melanoma. We also validated that some ICP-LncCRCTs could effectively predict ICI-response using three kinds of machine learning algorithms follow five independent datasets. Specially, combining ICP-LncCRCTs with the tumor mutation burden (TMB) improves the prediction of ICI-treated melanoma patients. Altogether, this study will improve our grasp of lncRNA functions and accelerating discovery of lncRNA-based biomarkers in ICI treatment.

DNA methylation heterogeneity attributable to a complex tumor immune microenvironment prompts prognostic risk in glioma.

Gliomas are malignant tumours of the human nervous system with different World Health Organization (WHO) classifications, glioblastoma (GBM) with higher grade and are more malignant than lower-grade glioma (LGG). To dissect how the DNA methylation heterogeneity in gliomas is influenced by the complex cellular composition of the tumour immune microenvironment, we first compared the DNA methylation profiles of purified human immune cells and bulk glioma tissue, stratifying three tumour immune microenvironmental subtypes for GBM and LGG samples from The Cancer Genome Atlas (TCGA). We found that more intermediate methylation sites were enriched in glioma tumour tissues, and used the Proportion of sites with Intermediate Methylation (PIM) to compare intertumoral DNA methylation heterogeneity. A larger PIM score reflected stronger DNA methylation heterogeneity. Enhanced DNA methylation heterogeneity was associated with stronger immune cell infiltration, better survival rates, and slower tumour progression in glioma patients. We then created a Cell-type-associated DNA Methylation Heterogeneity Contribution (CMHC) score to explore the impact of different immune cell types on heterogeneous CpG site (CpGct) in glioma tissues. We identified eight prognosis-related CpGct to construct a risk score: the Cell-type-associated DNA Methylation Heterogeneity Risk (CMHR) score. CMHR was positively correlated with cytotoxic T-lymphocyte infiltration (CTL), and showed better predictive performance for IDH status (AUC = 0.96) and glioma histological phenotype (AUC = 0.81). Furthermore, DNA methylation alterations of eight CpGct might be related to drug treatments of gliomas. In conclusion, we indicated that DNA methylation heterogeneity is associated with a complex tumour immune microenvironment, glioma phenotype, and patient's prognosis.

A comprehensive database of exosome molecular biomarkers and disease-gene associations.

Exosomes play a crucial role in intercellular communication and can be used as biomarkers for diagnostic and therapeutic clinical applications. However, systematic studies in cancer-associated exosomal nucleic acids remain a big challenge. Here, we developed ExMdb, a comprehensive database of exosomal nucleic acid biomarkers and disease-gene associations curated from published literature and high-throughput datasets. We performed a comprehensive curation of exosome properties including 4,586 experimentally supported gene-disease associations, 13,768 diagnostic and therapeutic biomarkers, and 312,049 nucleic acid subcellular locations. To characterize expression variation of exosomal molecules and identify causal factors of complex diseases, we have also collected 164 high-throughput datasets, including bulk and single-cell RNA sequencing (scRNA-seq) data. Based on these datasets, we performed various bioinformatics and statistical analyses to support our conclusions and advance our knowledge of exosome biology. Collectively, our dataset will serve as an essential resource for investigating the regulatory mechanisms of complex diseases and improving the development of diagnostic and therapeutic biomarkers.

Identification of somatic mutation-driven enhancers and their clinical utility in breast cancer.

Somatic mutations contribute to cancer development by altering the activity of enhancers. In the study, a total of 135 mutation-driven enhancers, which displayed significant chromatin accessibility changes, were identified as candidate risk factors for breast cancer (BRCA). Furthermore, we identified four mutation-driven enhancers as independent prognostic factors for BRCA subtypes. In Her2 subtype, enhancer G > C mutation was associated with poorer prognosis through influencing its potential target genes FBXW9, TRIR, and WDR83. We identified aminoglutethimide and quinpirole as candidate drugs targeting the mutated enhancer. In normal subtype, enhancer G > A mutation was associated with poorer prognosis through influencing its target genes ALOX15B, LINC00324, and MPDU1. We identified eight candidate drugs such as erastin, colforsin, and STOCK1N-35874 targeting the mutated enhancer. Our findings suggest that somatic mutations contribute to breast cancer subtype progression by altering enhancer activity, which could be potential candidates for cancer therapy.

Identification of cuproptosis-related diagnostic biomarkers in idiopathic pulmonary fibrosis.

Idiopathic pulmonary fibrosis (IPF) is a progressive and fatal lung disease with clinical and pathological heterogeneity. Recent studies have identified cuproptosis as a novel cell death mechanism. However, the role of cuproptosis-related genes in the pathogenesis of IPF is still unclear. Two IPF datasets of the Gene Expression Omnibus database were studied. Mann-Whitney U test, correlation analysis, functional enrichment analyses, single-sample gene set enrichment analysis, CIBERSORT, unsupervised clustering, weighted gene co-expression network analysis, and receiver operating characteristic curve analysis were used to conduct our research. The dysregulated cuproptosis-related genes and immune responses were identified between IPF patients and controls. Two cuproptosis-related molecular clusters were established in IPF, the high immune score group (C1) and the low immune score group (C2). Significant heterogeneity in immunity between clusters was revealed by functional analyses results. The module genes with the strongest correlation to the 2 clusters were identified by weighted gene co-expression network analysis results. Seven hub genes were found using the Cytoscape software. Ultimately, 2 validated diagnostic biomarkers of IPF, CDKN2A and NEDD4, were obtained. Subsequently, the results were validated in GSE47460. Our investigation illustrates that CDKN2A and NEDD4 may be valid biomarkers that were useful for IPF diagnosis and copper-related clustering.

Metabolic-Related Gene Prognostic Index for Predicting Prognosis, Immunotherapy Response, and Candidate Drugs in Ovarian Cancer.

Ovarian cancer (OC) is a highly heterogeneous disease, with patients at different tumor staging having different survival times. Metabolic reprogramming is one of the key hallmarks of cancer; however, the significance of metabolism-related genes in the prognosis and therapy outcomes of OC is unclear. In this study, we used weighted gene coexpression network analysis and differential expression analysis to screen for metabolism-related genes associated with tumor staging. We constructed the metabolism-related gene prognostic index (MRGPI), which demonstrated a stable prognostic value across multiple clinical trial end points and multiple validation cohorts. The MRGPI population had its distinct molecular features, mutational characteristics, and immune phenotypes. In addition, we investigated the response to immunotherapy in MRGPI subgroups and found that patients with low MRGPI were prone to benefit from anti-PD-1 checkpoint blockade therapy and exhibited a delayed treatment effect. Meanwhile, we identified four candidate therapeutic drugs (ABT-737, crizotinib, panobinostat, and regorafenib) for patients with high MRGPI, and we evaluated the pharmacokinetics and safety of the candidate drugs. In summary, the MRGPI was a robust clinical feature that could predict patient prognosis, immunotherapy response, and candidate drugs, facilitating clinical decision making and therapeutic strategy of OC.

LncPCD: a manually curated database of experimentally supported associations between lncRNA-mediated programmed cell death and diseases.

Programmed cell death (PCD) refers to controlled cell death that is conducted to keep the internal environment stable. Long noncoding RNAs (lncRNAs) participate in the progression of PCD in a variety of diseases. However, no specialized online repository is available to collect and store the associations between lncRNA-mediated PCD and diseases. Here, we developed LncPCD, a comprehensive database that provides information on experimentally supported associations of lncRNA-mediated PCD with diseases. The current version of LncPCD documents 6666 associations between five common types of PCD (apoptosis, autophagy, ferroptosis, necroptosis and pyroptosis) and 1222 lncRNAs in 331 diseases. We also manually curated a wealth of information: (1) 7 important lncRNA regulatory mechanisms, (2) 310 PCD-associated cell types in three species, (3) detailed information on lncRNA subcellular locations and (4) clinical applications for lncRNA-mediated PCD in diseases. Additionally, 10 single-cell sequencing datasets were integrated into LncPCD to characterize the dynamics of lncRNAs in diseases. Overall, LncPCD is an extremely useful resource for understanding the functions and mechanisms of lncRNA-mediated PCD in diseases. Database URL:  http://spare4.hospital.studio:9000/lncPCD/Home.jsp.

Proteomic characterization of epithelial ovarian cancer delineates molecular signatures and therapeutic targets in distinct histological subtypes.

Clear cell carcinoma (CCC), endometrioid carcinoma (EC), and serous carcinoma (SC) are the major histological subtypes of epithelial ovarian cancer (EOC), whose differences in carcinogenesis are still unclear. Here, we undertake comprehensive proteomic profiling of 80 CCC, 79 EC, 80 SC, and 30 control samples. Our analysis reveals the prognostic or diagnostic value of dysregulated proteins and phosphorylation sites in important pathways. Moreover, protein co-expression network not only provides comprehensive view of biological features of each histological subtype, but also indicates potential prognostic biomarkers and progression landmarks. Notably, EOC have strong inter-tumor heterogeneity, with significantly different clinical characteristics, proteomic patterns and signaling pathway disorders in CCC, EC, and SC. Finally, we infer MPP7 protein as potential therapeutic target for SC, whose biological functions are confirmed in SC cells. Our proteomic cohort provides valuable resources for understanding molecular mechanisms and developing treatment strategies of distinct histological subtypes.

Tumor biology, immune infiltration and liver function define seven hepatocellular carcinoma subtypes linked to distinct drivers, survival and drug response.

BACKGROUND & AIMS: Tumor heterogeneity is jointly determined by the components of the tumor ecosystem (TES) including tumor cells, immune cells, stromal cells, and non-cellular components. We aimed to identify subtypes using TES-related genes and determine subtype specific drivers and treatments for hepatocellular carcinoma (HCC). METHODS: We collected 68 genesets depicting tumor biology, immune infiltration, and liver function, totaling 2831 genes, and collected mRNA profiles and clinical data for over 6000 tumors from 65 datasets in the GEO, TCGA, ICGC, and several other databases. We designed a three-step clustering pipeline to identify subtypes. The microenvironment, genomic alteration, and drug response differences were systematically compared among subtypes. RESULTS: Seven subtypes (TES-1/2/3/4/5/6/7) were revealed in 159 tumors from the CHCC-HBV cohort. We constructed a single sample classifier using paired genes (sscpgsTES). TES subtypes were significantly associated with multiple clinical variables including etiology, and survival in 14 of 17 cohorts and the meta-cohort. TES-1 had the poorest prognosis and highest proliferation level. Both TES-2 and TES-7 were immune-enriched, however, TES-2 had a significantly worse prognosis, and hypoxic and immunosuppressive microenvironment. TES-4 had activated Wnt pathway, driven by CTNNB1 mutation. Good prognosis TES-6 exhibited the best differentiation. TES-5 and TES-3 were considered as novel subclasses by comparing with ten previous subtyping systems. TES-5 tumors had high AFP but good overall survival, and ∼45% of them harbored AXIN1 mutation. TES-3 was immune and stromal desert, may be driven by high copy number alteration burden, and had the poorest response to immune checkpoint inhibitor. TES-1 and TES-2 had significantly lower response to transarterial chemoembolization, but they showed significantly higher sensitivity to compound YM-155. CONCLUSIONS: Tumor ecosystem subtypes expand existing HCC subtyping systems, have distinct drivers, prognosis, and treatment vulnerabilities.

Identifying immune checkpoint-related lncRNA biomarkers for immunotherapy response and prognosis in cancers.

Long non-coding RNAs (lncRNAs) could modulate expression of immune checkpoints (ICPs) in tumor-immune. However, precise functions in immunity and potential for predicting ICP inhibitors (ICI) response have been described for only a few lncRNAs. Here, a multiple-step pipeline was developed to identify cancer- and immune-context ICP and lncRNA cooperative regulation pairs (ICPaLncCRPs) across cancers. Immune-related ICPs and lncRNAs were extracted follow immune cell lines and immunologic constant of rejection groups. ICPaLncCRP networks were constructed, which likely to modulate tumor-immune by specific patterns. Common and specific hub ICPaLncs such as MIR155HG, TRG-AS1 and PCED1B-AS1 maybe play central roles in prognosis and circulating. Moreover, these hub ICPaLncs were significantly correlated with immune cell infiltration based on bulk and single-cell RNA sequencing data. Some ICPaLncCRPs such as IDO1-MIR155HG could predict three- and five-year prognosis of melanoma in two independent datasets. We also validated that some ICPaLncCRPs could effectively predict ICI-response follow six independent datasets. Collectively, this study will enhance our understanding of lncRNA functions and accelerate discovery of lncRNA-based biomarkers in ICI treatment.

Reconstructing the immunosenescence core pathway reveals global characteristics in pan-cancer.

Immunosenescence has been demonstrated to play an important role in tumor progression. However, there is lacking comprehensive analyses of immunosenescence-related pathways. Meanwhile, the sex disparities of immunosenescence in cancer are still poorly understood. In this study, we analyzed the multi-omics data of 12,836 tumor samples, including genomics, transcriptomics, epigenomics, proteomics, and metabolomics. We systematically identified immunosenescence pathways that were disordered across cancer types. The mutations and copy number variations of immunosenescence pathways were found to be more active in pan-cancer. We reconstructed the immunosenescence core pathways (ISC-pathways) to improve the ability of prognostic stratification in 33 cancer types. We also found the head and neck squamous carcinoma (HNSC) contained abundant sex-specific immunosenescence features and showed sex differences in survival. We found that OSI-027 was a potential sex-specific drug in HNSC tumors, which tended to be more effective in male HNSC by targeting the MTOR gene in the PI3K-Akt signaling pathway. In conclusion, our study provided a systematic understanding of immunosenescence pathways and revealed the global characteristics of immunosenescence in pan-cancer. We highlighted MTOR gene could be a powerful immunosenescence biomarker of HNSC that helps to develop sex-specific immunosenescence drugs.

CRS: a circadian rhythm score model for predicting prognosis and treatment response in cancer patients.

BACKGROUND: Circadian rhythm regulates complex physiological activities in organisms. A strong link between circadian dysfunction and cancer has been identified. However, the factors of dysregulation and functional significance of circadian rhythm genes in cancer have received little attention. METHODS: In 18 cancer types from The Cancer Genome Atlas (TCGA), the differential expression and genetic variation of 48 circadian rhythm genes (CRGs) were examined. The circadian rhythm score (CRS) model was created using the ssGSEA method, and patients were divided into high and low groups based on the CRS. The Kaplan-Meier curve was created to assess the patient survival rate. Cibersort and estimate methods were used to identify the infiltration characteristics of immune cells between different CRS subgroups. Gene Expression Omnibus (GEO) dataset is used as verification queue and model stability evaluation queue. The CRS model's ability to predict chemotherapy and immunotherapy was assessed. Wilcoxon rank-sum test was used to compare the differences of CRS among different patients. We use CRS to identify potential "clock-drugs" by the connective map method. RESULTS: Transcriptomic and genomic analyses of 48 CRGs revealed that most core clock genes are up-regulated, while clock control genes are down-regulated. Furthermore, we show that copy number variation may affect CRGs aberrations. Based on CRS, patients can be classified into two groups with significant differences in survival and immune cell infiltration. Further studies showed that patients with low CRS were more sensitive to chemotherapy and immunotherapy. Additionally, we identified 10 compounds (e.g. flubendazole, MLN-4924, ingenol) that are positively associated with CRS, and have the potential to modulate circadian rhythms. CONCLUSIONS: CRS can be utilized as a clinical indicator to predict patient prognosis and responsiveness to therapy, and identify potential "clock-drugs".

Pan-cancer transcriptomic analysis identified six classes of immunosenescence genes revealed molecular links between aging, immune system and cancer.

Aging is a complex process that significantly impacts the immune system. The aging-related decline of the immune system, termed immunosenescence, can lead to disease development, including cancer. The perturbation of immunosenescence genes may characterize the associations between cancer and aging. However, the systematical characterization of immunosenescence genes in pan-cancer remains largely unexplored. In this study, we comprehensively investigated the expression of immunosenescence genes and their roles in 26 types of cancer. We developed an integrated computational pipeline to identify and characterize immunosenescence genes in cancer based on the expression profiles of immune genes and clinical information of patients. We identified 2218 immunosenescence genes that were significantly dysregulated in a wide variety of cancers. These immunosenescence genes were divided into six categories based on their relationships with aging. Besides, we assessed the importance of immunosenescence genes in clinical prognosis and identified 1327 genes serving as prognostic markers in cancers. BTN3A1, BTN3A2, CTSD, CYTIP, HIF1AN, and RASGRP1 were associated with ICB immunotherapy response and served as prognostic factors after ICB immunotherapy in melanoma. Collectively, our results furthered the understanding of the relationship between immunosenescence and cancer and provided insights into immunotherapy for patients.

Radiogenomic analysis of prediction HER2 status in breast cancer by linking ultrasound radiomic feature module with biological functions.

BACKGROUND: Human epidermal growth factor receptor 2 (HER2) overexpressed associated with poor prognosis in breast cancer and HER2 has been defined as a therapeutic target for breast cancer treatment. We aimed to explore the molecular biological information in ultrasound radiomic features (URFs) of HER2-positive breast cancer using radiogenomic analysis. Moreover, a radiomics model was developed to predict the status of HER2 in breast cancer. METHODS: This retrospective study included 489 patients who were diagnosed with breast cancer. URFs were extracted from a radiomics analysis set using PyRadiomics. The correlations between differential URFs and HER2-related genes were calculated using Pearson correlation analysis. Functional enrichment of the identified URFs-correlated HER2 positive-specific genes was performed. Lastly, the radiomics model was developed based on the URF-module mined from auxiliary differential URFs to assess the HER2 status of breast cancer. RESULTS: Eight differential URFs (p < 0.05) were identified among the 86 URFs extracted by Pyradiomics. 25 genes that were found to be the most closely associated with URFs. Then, the relevant biological functions of each differential URF were obtained through functional enrichment analysis. Among them, Zone Entropy is related to immune cell activity, which regulate the generation of calcification in breast cancer. The radiomics model based on the Logistic classifier and URF-module showed good discriminative ability (AUC = 0.80, 95% CI). CONCLUSION: We searched for the URFs of HER2-positive breast cancer, and explored the underlying genes and biological functions of these URFs. Furthermore, the radiomics model based on the Logistic classifier and URF-module relatively accurately predicted the HER2 status in breast cancer.

Immune Pathways with Aging Characteristics Improve Immunotherapy Benefits and Drug Prediction in Human Cancer.

(1) Background: Perturbation of immune-related pathways can make substantial contributions to cancer. However, whether and how the aging process affects immune-related pathways during tumorigenesis remains largely unexplored. (2) Methods: Here, we comprehensively investigated the immune-related genes and pathways among 25 cancer types using genomic and transcriptomic data. (3) Results: We identified several pathways that showed aging-related characteristics in various cancers, further validated by conventional aging-related gene sets. Genomic analysis revealed high mutation burdens in cytokines and cytokines receptors pathways, which were strongly correlated with aging in diverse cancers. Moreover, immune-related pathways were found to be favorable prognostic factors in melanoma. Furthermore, the expression level of these pathways had close associations with patient response to immune checkpoint blockade therapy in melanoma and non-small cell lung cancer. Applying a net-work-based method, we predicted immune- and aging-related genes in pan-cancer and utilized these genes for potential immunotherapy drug discovery. Mapping drug target data to our top-ranked genes identified potential drug targets, FYN, JUN, and SRC. (4) Conclusions: Taken together, our systematic study helped interpret the associations among immune-related pathways, aging, and cancer and could serve as a resource for promoting clinical treatment.

RNA2Immune: A Database of Experimentally Supported Data Linking Non-coding RNA Regulation to The Immune System.

Non-coding RNAs (ncRNAs), such as microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs), have emerged as important regulators of the immune system and are involved in the control of immune cell biology, disease pathogenesis, as well as vaccine responses. A repository of ncRNA-immune associations will facilitate our understanding of ncRNA-dependent mechanisms in the immune system and advance the development of therapeutics and prevention for immune disorders. Here, we describe a comprehensive database, RNA2Immune, which aims to provide a high-quality resource of experimentally supported database linking ncRNA regulatory mechanisms to immune cell function, immune disease, cancer immunology, and vaccines. The current version of RNA2Immune documents 50,433 immune-ncRNA associations in 42 host species, including (1) 6690 ncRNA associations with immune functions involving 31 immune cell types; (2) 38,672 ncRNA associations with 348 immune diseases; (3) 4833 ncRNA associations with cancer immunology; and (4) 238 ncRNA associations with vaccine responses involving 26 vaccine types targeting 22 diseases. RNA2Immune provides a user-friendly interface for browsing, searching, and downloading ncRNA-immune system associations. Collectively, RNA2Immune provides important information about how ncRNAs influence immune cell function, how dysregulation of these ncRNAs leads to pathological consequences (immune diseases and cancers), and how ncRNAs affect immune responses to vaccines. RNA2Immune is available at http://bio-bigdata.hrbmu.edu.cn/rna2immune/home.jsp.