scholarly journals Alternative Splicing Events in Tumor Immune Infiltration in Colorectal Cancer

2021 ◽  
Vol 11 ◽  
Author(s):  
Jian-yu Shi ◽  
Yan-yan Bi ◽  
Bian-fang Yu ◽  
Qing-feng Wang ◽  
Dan Teng ◽  
...  

Despite extensive research, the exact mechanisms involved in colorectal cancer (CRC) etiology and pathogenesis remain unclear. This study aimed to examine the correlation between tumor-associated alternative splicing (AS) events and tumor immune infiltration (TII) in CRC. We analyzed transcriptome profiling and clinical CRC data from The Cancer Genome Atlas (TCGA) database and lists of AS-related and immune-related signatures from the SpliceSeq and Innate databases, respectively to develop and validate a risk model of differential AS events and subsequently a TII risk model. We then conducted a two-factor survival analysis to study the association between TII and AS risk and evaluated the associations between immune signatures and six types of immune cells based on the TIMER database. Subsequently, we studied the distribution of six types of TII cells in high- and low-risk groups for seven AS events and in total. We obtained the profiles of AS events/genes for 484 patients, which included 473 CRC tumor samples and 41 corresponding normal samples, and detected 22581 AS events in 8122 genes. Exon Skip (ES) (8446) and Mutually Exclusive Exons (ME) (74) exhibited the most and fewest AS events, respectively. We then classified the 433 patients with CRC into low-risk (n = 217) and high-risk (n = 216) groups based on the median risk score in different AS events. Compared with patients with low-risk scores (mortality = 11.8%), patients with high-risk scores were associated with poor overall survival (mortality = 27.6%). The risk score, cancer stage, and pathological stage (T, M, and N) were closely correlated with prognosis in patients with CRC (P < 0.001). We identified 6479 differentially expressed genes from the transcriptome profiles of CRC and intersected 468 differential immune-related signatures. High-AS-risk and high-TII-risk predicted a poor prognosis in CRC. Different AS types were associated with different TII risk characteristics. Alternate Acceptor site (AA) and Alternate Promoter (AP) events directly affected the concentration of CD4T cells, and the level of CD8T cells was closely correlated with Alternate Terminator (AT) and Exon Skip (ES) events. Thus, the concentration of CD4T and CD8T cells in the CRC immune microenvironment was not specifically modulated by AS. However, B cell, dendritic cell, macrophage, and neutrophilic cell levels were strongly correlated with AS events. These results indicate adverse associations between AS event risk levels and immune cell infiltration density. Taken together, our findings show a clear association between tumor-associated alternative splicing and immune cell infiltration events and patient outcome and could form a basis for the identification of novel markers and therapeutic targets for CRC and other cancers in the future.

2021 ◽  
Vol 19 (1) ◽  
Author(s):  
Carly A. Conran ◽  
Zhuqing Shi ◽  
William Kyle Resurreccion ◽  
Rong Na ◽  
Brian T. Helfand ◽  
...  

Abstract Background Genome-wide association studies have identified thousands of disease-associated single nucleotide polymorphisms (SNPs). A subset of these SNPs may be additively combined to generate genetic risk scores (GRSs) that confer risk for a specific disease. Although the clinical validity of GRSs to predict risk of specific diseases has been well established, there is still a great need to determine their clinical utility by applying GRSs in primary care for cancer risk assessment and targeted intervention. Methods This clinical study involved 281 primary care patients without a personal history of breast, prostate or colorectal cancer who were 40–70 years old. DNA was obtained from a pre-existing biobank at NorthShore University HealthSystem. GRSs for colorectal cancer and breast or prostate cancer were calculated and shared with participants through their primary care provider. Additional data was gathered using questionnaires as well as electronic medical record information. A t-test or Chi-square test was applied for comparison of demographic and key clinical variables among different groups. Results The median age of the 281 participants was 58 years and the majority were female (66.6%). One hundred one (36.9%) participants received 2 low risk scores, 99 (35.2%) received 1 low risk and 1 average risk score, 37 (13.2%) received 1 low risk and 1 high risk score, 23 (8.2%) received 2 average risk scores, 21 (7.5%) received 1 average risk and 1 high risk score, and no one received 2 high risk scores. Before receiving GRSs, younger patients and women reported significantly more worry about risk of developing cancer. After receiving GRSs, those who received at least one high GRS reported significantly more worry about developing cancer. There were no significant differences found between gender, age, or GRS with regards to participants’ reported optimism about their future health neither before nor after receiving GRS results. Conclusions Genetic risk scores that quantify an individual’s risk of developing breast, prostate and colorectal cancers as compared with a race-defined population average risk have potential clinical utility as a tool for risk stratification and to guide cancer screening in a primary care setting.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Tuo Liang ◽  
Jiarui Chen ◽  
GuoYong Xu ◽  
Zide Zhang ◽  
Jiang Xue ◽  
...  

AbstractWe established a relationship among the immune-related genes, tumor-infiltrating immune cells (TIICs), and immune checkpoints in patients with osteosarcoma. The gene expression data for osteosarcoma were downloaded from UCSC Xena and GEO database. Immune-related differentially expressed genes (DEGs) were detected to calculate the risk score. “Estimate” was used for immune infiltrating estimation and “xCell” was used to obtain 64 immune cell subtypes. Furthermore, the relationship among the risk scores, immune cell subtypes, and immune checkpoints was evaluated. The three immune-related genes (TYROBP, TLR4, and ITGAM) were selected to establish a risk scoring system based on their integrated prognostic relevance. The GSEA results for the Hallmark and KEGG pathways revealed that the low-risk score group exhibited the most gene sets that were related to immune-related pathways. The risk score significantly correlated with the xCell score of macrophages, M1 macrophages, and M2 macrophages, which significantly affected the prognosis of osteosarcoma. Thus, patients with low-risk scores showed better results with the immune checkpoints inhibitor therapy. A three immune-related, gene-based risk model can regulate macrophage activation and predict the treatment outcomes the survival rate in osteosarcoma.


2021 ◽  
Vol 8 ◽  
Author(s):  
Mingqin Ge ◽  
Jie Niu ◽  
Ping Hu ◽  
Aihua Tong ◽  
Yan Dai ◽  
...  

Objective: This study aimed to construct a prognostic ferroptosis-related signature for thyroid cancer and probe into the association with tumor immune microenvironment.Methods: Based on the expression profiles of ferroptosis-related genes, a LASSO cox regression model was established for thyroid cancer. Kaplan-Meier survival analysis was presented between high and low risk groups. The predictive performance was assessed by ROC. The predictive independency was validated via multivariate cox regression analysis and stratified analysis. A nomogram was established and verified by calibration curves. The enriched signaling pathways were predicted via GSEA. The association between the signature and immune cell infiltration was analyzed by CIBERSORT. The ferroptosis-related genes were validated in thyroid cancer tissues by immunohistochemistry and RT-qPCR.Results: A ferroptosis-related eight gene model was established for predicting the prognosis of thyroid cancer. Patients with high risk score indicated a poorer prognosis than those with low risk score (p = 1.186e-03). The AUCs for 1-, 2-, and 3-year survival were 0.887, 0.890, and 0.840, respectively. Following adjusting other prognostic factors, the model could independently predict the prognosis (p = 0.015, HR: 1.870, 95%CI: 1.132–3.090). A nomogram combining the signature and age was constructed. The nomogram-predicted probability of 1-, 3-, and 5-year survival approached the actual survival time. Several ferroptosis-related pathways were enriched in the high-risk group. The signature was distinctly associated with the immune cell infiltration. After validation, the eight genes were abnormally expressed between thyroid cancer and control tissues.Conclusion: Our findings established a prognostic ferroptosis-related signature that was associated with the immune microenvironment for thyroid cancer.


2021 ◽  
Author(s):  
Chen-jie Qiu ◽  
Xue-bing Wang ◽  
Zi-ruo Zheng ◽  
Chao-zhi Yang ◽  
Kai Lin ◽  
...  

Abstract Background: The purpose of this study was to identify ferroptosis-related genes (FRGs) associated with the prognosis of pancreatic cancer and to construct a prognostic model based on FRGs. Methods: Based on pancreatic cancer data obtained from The Cancer Genome Atlas database, we established the prognostic model from 232 FRGs. A nomogram was constructed by combining the prognostic model and clinicopathological features. Gene Expression Omnibus datasets and tissue samples obtained from our center were utilized to validate the model. Relationship between risk score and immune cell infiltration was explored by CIBERSORT and TIMER.Results: The prognostic model was established based on four FRGs (ENPP2, ATG4D, SLC2A1 and MAP3K5) and can be an independent risk factor in pancreatic cancer (HR 1.648, 95% CI 1.335-2.035, p < 0.001). Based on the median risk score, patients were divided into a high-risk group and a low-risk group. The prognosis of the low-risk group was significantly better than that of the high-risk group. In the high-risk group, patients treated with chemotherapy had a better prognosis. The nomogram showed that the model was the most important element. Gene set enrichment analysis identified three key pathways, namely, TGFβ signaling, HIF signaling pathway and adherens junction. The prognostic model can also affect the immune cell infiltration, such as macrophages M0, M1, CD4+T cell and CD8+T cell. Conclusion: A ferroptosis-related prognostic model can be employed to predict the prognosis of pancreatic cancer. Ferroptosis can be an important marker and immunotherapy can be a potential therapeutic target for pancreatic cancer.


Author(s):  
Liuxing Wu ◽  
Xin Hu ◽  
Hongji Dai ◽  
Kexin Chen ◽  
Ben Liu

Despite robust evidence for the role of m6A in cancer development and progression, its association with immune infiltration and survival outcomes in melanoma remains obscure. Here, we aimed to develop an m6A-related risk signature to improve prognostic and immunotherapy responder prediction performance in the context of melanoma. We comprehensively analyzed the m6A cluster and immune infiltration phenotypes of public datasets. The TCGA (n = 457) and eleven independent melanoma cohorts (n = 758) were used as the training and validation datasets, respectively. We identified two m6A clusters (m6A-clusterA and m6A-clusterB) based on the expression pattern of m6A regulators via unsupervised consensus clustering. IGF2BP1 (7.49%), KIAA1429 (7.06%), and YTHDC1 (4.28%) were the three most frequently mutated genes. There was a correlation between driver genes mutation statuses and the expression of m6A regulators. A significant difference in tumor-associated immune infiltration between two m6A clusters was detected. Compared with m6A-clusterA, the m6A-clusterB was characterized by a lower immune score and immune cell infiltration but higher mRNA expression-based stemness index (mRNAsi). An m6A-related risk signature consisting of 12 genes was determined via Cox regression analysis and divided the patients into low- and high-risk groups (IL6ST, MBNL1, NXT2, EIF2A, CSGALNACT1, C11orf58, CD14, SPI1, NCCRP1, BOK, CD74, PAEP). A nomogram was developed for the prediction of the survival rate. Compared with the high-risk group, the low-risk group was characterized by high expression of immune checkpoints and immunophenoscore (IPS), activation of immune-related pathways, and more enriched in immune cell infiltrations. The low-risk group had a favorable prognosis and contained the potential beneficiaries of the immune checkpoint blockade therapy and verified by the IMvigor210 cohort (n = 298). The m6A-related signature we have determined in melanoma highlights the relationships between m6A regulators and immune cell infiltration. The established risk signature was identified as a promising clinical biomarker of melanoma.


2021 ◽  
Vol 12 ◽  
Author(s):  
Luping Zhang ◽  
Shaokun Wang ◽  
Yachen Wang ◽  
Weidan Zhao ◽  
Yingli Zhang ◽  
...  

BackgroundImbalanced nutritional supply and demand in the tumor microenvironment often leads to hypoxia. The subtle interaction between hypoxia and immune cell behavior plays an important role in tumor occurrence and development. However, the functional relationship between hypoxia and the tumor microenvironment remains unclear. Therefore, we aimed to investigate the effect of hypoxia on the intestinal tumor microenvironment.MethodWe extracted the names of hypoxia-related genes from the Gene Set Enrichment Analysis (GSEA) database and screened them for those associated with colorectal cancer prognosis, with the final list including ALDOB, GPC1, ALDOC, and SLC2A3. Using the sum of the expression levels of these four genes, provided by The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, and the expression coefficients, we developed a hypoxia risk score model. Using the median risk score value, we divided the patients in the two databases into high- and low-risk groups. GSEA was used to compare the enrichment differences between the two groups. We used the CIBERSORT computational method to analyze immune cell infiltration. Finally, the correlation between these five genes and hypoxia was analyzed.ResultThe prognosis of the two groups differed significantly, with a higher survival rate in the low-risk group than in the high-risk group. We found that the different risk groups were enriched by immune-related and inflammatory pathways. We identified activated M0 macrophages in TCGA and GEO databases and found that CCL2/4/5, and CSF1 contributed toward the increased infiltration rate of this immune cell type. Finally, we observed a positive correlation between the five candidate genes’ expression and the risk of hypoxia, with significant differences in the level of expression of each of these genes between patient risk groups.ConclusionOverall, our data suggest that hypoxia is associated with the prognosis and rate of immune cell infiltration in patients with colorectal cancer. This finding may improve immunotherapy for colorectal cancer.


2021 ◽  
Author(s):  
Yangyang Yin ◽  
Ying Liu ◽  
Yun Pan ◽  
Tianyu Yu ◽  
Bin Liu ◽  
...  

Abstract Background: Colorectal cancer (CRC) is the third most commonly diagnosed cancer worldwide. Epithelial-mesenchymal transition (EMT) has been viewed to play a vital role in immune regulation and treatment response shaping. EMT is associated with an invasive or metastatic phenotype in CRC and is also related to patient prognosis and treatment responsiveness. We aimed to explore an EMT gene panel with potential application in patient classification and precise treatment and assist clinicians to generate individualized immunotherapeutic strategies for CRC patients.Methods: TCGA, GEO, STRING, TRRUST, TragetScan, miRTarBase, miRDB, cBioPortal, StarBase databases were utilized in this study.Results: In this study, EMT factors were screened in three different ways(EMT factors, EMT-related pathways, EMT genotyping genes).And the prognosis-related modules were screened by the WGCNA method. Then, Cox single factor regression analysis was performed on the module hub gene, combined with multivariate Cox regression, the prognosis model was established, and the risk score of each sample was calculated. Then the samples were divided into high-risk samples and low-risk samples according to the risk score, and the differences in immune cell infiltration, mutation, CNV, clinical characteristics between high-risk and low-risk samples were compared. The risk model can effectively predict the prognosis of samples, which is verified by two external data, and it can also effectively predict the prognosis of CRC samples in the other two digestive tract cancers (liver cancer and gastric cancer) and has a good indication for the effect of chemotherapy treatment response and immunotherapy. Conclusion: The prognosis model can effectively predict the prognosis of samples and may be an effective tool for treatment guidance in CRC patients.


2021 ◽  
Vol 2021 ◽  
pp. 1-14
Author(s):  
Ruixiang Luo ◽  
Mengjun Huang ◽  
Yinhuai Wang

Background. Prostate cancer (PC) is one of the most critical cancers affecting men’s health worldwide. The development of many cancers involves dysregulation or mutations in key transcription factors. This study established a transcription factor-based risk model to predict the prognosis of PC and potential therapeutic drugs. Materials and Methods. In this study, RNA-sequencing data were downloaded and analyzed using The Cancer Genome Atlas dataset. A total of 145 genes related to the overall survival rate of PC patients were screened using the univariate Cox analysis. The Kdmist clustering method was used to classify prostate adenocarcinoma (PRAD), thereby determining the cluster related to the transcription factors. The support vector machine-recursive feature elimination method was used to identify genes related to the types of transcription factors and the key genes specifically upregulated or downregulated were screened. These genes were further analyzed using Lasso to establish a model. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were used for the functional analysis. The TIMER algorithm was used to quantify the abundance of immune cells in PRAD samples. The chemotherapy response of each GBM patient was predicted based on the public pharmacogenomic database, Genomics of Drug Sensitivity in Cancer (GDSC, http://www.cancerrxgene.org). The R package “pRRophetic” was applied to drug sensitivity (IC50) value prediction. Results. We screened 10 genes related to prognosis, including eight low-risk genes and two high-risk genes. The receiver operating characteristic (ROC) curve was 0.946. Patients in the high-risk score group had a poorer prognosis than those in the low-risk score group. The average area under the curve value of the model at different times was higher than 0.8. The risk score was an independent prognostic factor. Compared with the low-risk score group, early growth response-1 (EGR1), CACNA2D1, AC005831.1, SLC52A3, TMEM79, IL20RA, CRACR2A, and FAM189A2 expressions in the high-risk score group were decreased, while AC012181.1 and TRAPPC8 expressions were increased. GO and KEGG analyses showed that prognosis was related to various cancer signaling pathways. The proportion of B_cell, T_cell_CD4, and macrophages in the high-risk score group was significantly higher than that in the low-risk score group. A total of 25 classic immune checkpoint genes were screened out to express abnormally high-risk scores, and there were significant differences. Thirty mutant genes were identified; in the high- and low-risk score groups, SPOP, TP53, and TTN had the highest mutation frequency, and their mutations were mainly missense mutations. A total of 36 potential drug candidates for the treatment of PC were screened and identified. Conclusions. Ten genes of both high-and low-risk scores were associated with the prognosis of PC. PC prognosis may be related to immune disorders. SPOP, TP53, and TTN may be potential targets for the prognosis of PC.


2020 ◽  
Author(s):  
Carly Ann Conran ◽  
Zhuqing Shi ◽  
William Kyle Resurreccion ◽  
Rong Na ◽  
Brian T. Helfand ◽  
...  

Abstract Background: Genome-wide association studies have identified thousands of disease-associated single nucleotide polymorphisms (SNPs). A subset of these SNPs may be additively combined to generate genetic risk scores (GRSs) that confer risk for a specific disease. Although the clinical validity of GRSs to predict risk of specific diseases has been well established, there is still a great need to determine their clinical utility by applying GRSs in primary care for cancer risk assessment and targeted intervention.Methods: This clinical study involved 281 primary care patients without a personal history of breast, prostate or colorectal cancer who were 40-70 years old. DNA was obtained from a pre-existing biobank at NorthShore University HealthSystem. GRSs for colorectal cancer and breast or prostate cancer were calculated and shared with participants through their primary care provider. Additional data was gathered using questionnaires as well as electronic medical record information. A t-test or Chi-square test was applied for comparison of demographic and key clinical variables among different groups.Results: The median age of the 281 participants was 58 years and the majority were female (66.6%). One hundred one (36.9%) participants received 2 low risk scores, 99 (35.2%) received 1 low risk and 1 average risk score, 37 (13.2%) received 1 low risk and 1 high risk score, 23 (8.2%) received 2 average risk scores, 21 (7.5%) received 1 average risk and 1 high risk score, and no one received 2 high risk scores. Before receiving GRSs, younger patients and women reported significantly more worry about risk of developing cancer. After receiving GRSs, those who received at least one high GRS reported significantly more worry about developing cancer. There were no significant differences found between gender, age, or GRS with regards to participants’ reported optimism about their future health neither before nor after receiving GRS results.Conclusions: Genetic risk scores that quantify an individual’s risk of developing breast, prostate and colorectal cancers as compared with a race-defined population average risk have potential clinical utility as a tool for risk stratification and to guide cancer screening in a primary care setting.


Author(s):  
Shuo Hong ◽  
Yueming Zhang ◽  
Manming Cao ◽  
Anqi Lin ◽  
Qi Yang ◽  
...  

Objective: Resistance to immune checkpoint inhibitors (ICIs) has been a massive obstacle to ICI treatment in metastatic urothelial carcinoma (MUC). Recently, increasing evidence indicates the clinical importance of the association between hypoxia and immune status in tumor patients. Therefore, it is necessary to investigate the relationship between hypoxia and prognosis in metastatic urothelial carcinoma.Methods: Transcriptomic and clinical data from 348 MUC patients who underwent ICI treatment from a large phase 2 trial (IMvigor210) were investigated in this study. The cohort was randomly divided into two datasets, a training set (n = 213) and a testing set (n = 135). Data of hypoxia-related genes were downloaded from the molecular signatures database (MSigDB), and screened by univariate and multivariate Cox regression analysis to construct a prognosis-predictive model. The robustness of the model was evaluated in two melanoma cohorts. Furthermore, an external validation cohort, the bladder cancer cohort, from the Cancer Genome Atlas (TCGA) database, was t used to explore the mechanism of gene mutation, immune cell infiltration, signaling pathway enrichment, and drug sensitivity.Results: We categorized patients as the high- or low- risk group using a four-gene hypoxia risk model which we constructed. It was found that patients with high-risk scores had significantly worse overall survival (OS) compared with those with low-risk scores. The prognostic model covers 0.71 of the area under the ROC curve in the training set and 0.59 in the testing set, which is better than the survival prediction of MUC patients using the clinical characteristics. Mutation analysis results showed that deletion mutations in RB1, TP53, TSC1 and KDM6A were correlated with hypoxic status. Immune cell infiltration analysis illustrated that the infiltration T cells, B cells, Treg cells, and macrophages was correlated with hypoxia. Functional enrichment analysis revealed that a hypoxic microenvironment activated inflammatory pathways, glucose metabolism pathways, and immune-related pathways.Conclusion: In this investigation, a four-gene hypoxia risk model was developed to evaluate the degree of hypoxia and prognosis of ICI treatment, which showed a promising clinical prediction value in MUC. Furthermore, the hypoxia risk model revealed a close relationship between hypoxia and the tumor immune microenvironment.


Sign in / Sign up

Export Citation Format

Share Document